text:
# Magnetohydrostatic (MHS) atmospheres

## 1 Introduction

In our observation of solar ‘threads’ (November & Koutchmy 1996, hereafter NK96) we puzzled over whether coronal flux tubes should expand with height to maintain pressure balance with the ambient gas-pressure decrease or follow the lines of the external magnetic field. After all, if the corona is independently heated, the gas and magnetic height variations should be generally different. The many observations of EUV and X-ray loops as well as the purely density-sensitive observations of white-light $`>15`$ Mm ‘voids’ (MacQueen et al. 1974), $`>2`$ Mm threads, and scintillation $`>1`$ kilometer filamentary microstructure (Coles & Harmon 1978; Woo et al. 1995) all suggest that the large-scale solar coronal magnetic field is interspersed with flux tubes with a great range of sizes. Figure 1 shows dark and bright threads in a square region about a solar radius on a side over the west limb taken with the Canada-France-Hawaii Telescope (CFHT) at the unique total-eclipse opportunity on Mauna Kea on July 11, 1991. The dark and bright threads, which were essentially unchanged over the 4-minute eclipse duration, appear to be aligned and organized in arched and radial surfaces, probably current sheets, overlapping in their projection onto the plane of the sky in the line-of-sight view.

We know that isolated thin flux tubes must expand with height in balance with the ambient gas-pressure decrease. The flux-tube gas pressure difference, external minus internal, is related to the magnetic pressure difference by the static equilibrium relation

$$p_ep_i=\frac{1}{8\pi }\left(B_i^2B_e^2\right).$$
(1)
The Wilson depression in sunspots shows that at least large-scale photospheric flux tubes are nearly evacuated. Evacuated flux tubes exhibit the maximum magnetic flux concentration and field strength possible. In evacuated flux tubes $`p_i=0`$ and the photospheric external field $`B_e=0`$, so the internal field strength must follow the ambient gas-pressure decrease $`B_i(r)^2=8\pi p_e(r)`$, giving flux tubes that expand with height to conserve their total internal magnetic flux. Since the scale height for $`p_e(r)`$ is small in the relatively cool photosphere and chromosphere, a dramatic flux-tube expansion occurs, which is the accepted explanation for the solar ‘magnetic canopy’ (Giovanelli 1982). Both voids and a range of sizes of dark threads were also found to have contrasts consistent with what is expected for fully evacuated cylindrical flux tubes (MacQueen et al. 1983; NK96). In the corona, where the external magnetic field is not zero, the magnetic flux inside evacuated flux tubes must decrease with height to follow the total external ambient gas plus magnetic pressures, or $`B_i(r)^2=8\pi p_e(r)+B_e(r)^2`$. This gives a small flux-tube expansion with height in the large gas-pressure scale height of the hot corona. A slight flux-tube expansion with height seems evident in some of the threads in Fig. 1.

We can argue too that the radial expansion of isolated thin flux tubes must conform to the external magnetic field. Parker (1972) shows that

$$\left(B_e\right)B_i=0,$$
(2)
in MHS equilibria without gravity, where $`B_e`$ represents a uniform background magnetic field, and the field in flux tubes $`B_i`$ contains all of the spatial variations. The theorem is locally applicable everywhere in perturbed MHS equilibria with gravity, i.e. Eq. (2) with $`B_e`$ allowed to vary slowly spatially, as we discuss in Sect. 2. Thus thin flux tubes in a gravitational field must spread to follow the lines of the relatively uniform background magnetic field; the radial functional variations must be directly related, that is $`|B_i(r)||B_e(r)|`$ along flux tubes.

Therefore both ways of looking at the problem are essentially correct. The coronal flux-tube field must comply with the external field $`B_i(r)^2B_e(r)^2`$, and the radial variation of the external gas pressure outside isolated thin evacuated coronal flux tubes is proportional to the magnetic pressure difference $`p_e(r)B_i(r)^2B_e(r)^2`$. Hence, the external gas and magnetic pressures must be coupled with $`p_e(r)B_e(r)^2`$ outside and along isolated thin evacuated flux tubes. The background magnetic field is dominant in lower stellar coronae and must determine the expansion of embedded flux tubes. The expansion of evacuated flux tubes then defines the radial variation of the ambient gas pressure for equilibrium, the gas pressure being horizontally stratified in the external force-free or potential field. More generally, flux tubes of arbitrary base pressure but in thermal balance with the local ambient atmosphere satisfy $`p_i(r)p_e(r)`$, giving the same external atmosphere $`p_e(r)B_e(r)^2`$. Thus an atmosphere containing isolated thin non-force-free flux tubes in different forms of local balance must exhibit a temperature profile defined by the dominant background magnetic field for a perfect gas in hydrostatic equilibrium.

The properties of flux tubes reflect a natural degeneracy in the form of non-force-free magnetic fields, and imply the existence of an underlying mathematical relationship between an equilibrium atmosphere and large-scale magnetic field. In Sect. 2, we show that the height variations of the gas pressure and magnetic field amplitude are directly coupled in the most-significant-order terms of the gravitational perturbation in general non-force-free MHS equilibria, with smaller field divergence and curl effects. In low-$`\beta `$ systems, with strong magnetic fields, the external gas pressure corresponds to a small difference $`p_e(r)=(B_i(r)^2B_e(r)^2)/(8\pi )`$, which allows large pressure deviations for a given relative flexibility of the large-scale magnetic field around individual flux tubes. Magnetic fields can exhibit and persist in significant deviations from equilibrium, however the coupling $`p_e(r)B_e(r)^2`$ describes the general equilibrium condition for arbitrary $`\beta `$.

Historical MHS equilibrium solutions with gravity have also suggested that a general relationship exists between the large magnetic and gas-pressure scales. Dungey (1953) derived plane-parallel MHS equilibria for a magnetic field with its degenerate axis directed perpendicular to gravity (see also: Low 1975; Cheng & Choe 1998). His solutions are represented in a modified Grad-Shafranov (GS) Equation that retains the gas-pressure radial scale-height variation as a separable factor in the pressure term. Dungey showed that when the magnetic field is also zero in its degenerate axis, the separable factor can be absorbed into the coordinates, giving the classical GS Equation in a remapped coordinate system. The new GS coordinates are defined by the gas-pressure scale height, meaning that magnetic fields in a prespecified atmosphere must be spatially distorted. Horizontal flux tubes exhibit substantial radial deformation in a solar corona of given constant temperature (Zweibel & Hundhausen 1982; Low 1992). Depending upon their size compared to the gas-pressure scale height, flux tubes are either compressed or expanded in the gravitational direction, and approximately undistorted only in a specific atmosphere.

In Sect. 3, we examine equilibrium radial profiles of temperature and flux-tube diameter using an approximate solution to the equilibrium equations. The profiles reproduce the salient features of the solar atmosphere: the rapid expansion of flux tubes in the photosphere where the atmospheric gas pressure is dominant, a transition region of reasonable thickness located just above the height where the gas and magnetic energy densities are equal, and a corona of correct base temperature for the sun’s main magnetic spherical harmonic. Steady flows do not alter the coupling, but may produce a Bernoulli pressure, which modifies the gas pressure, and leads to possible accelerated steady wind solutions.

Although the dynamics of the system can be quite complicated, the equilibrium state characterizes the average in many types of MHD systems, e.g. ‘quasi-steady’ (Low 1980). Hydrodynamics in such systems tend to maintain the MHS atmosphere, as we discuss in Sect. 4. The possible importance of MHS atmospheres for other problems in physics and astrophysics beyond the formation of solar and stellar coronae and in more marginal non-quasi-steady physical conditions needs to be considered, but discussion lies beyond the scope of this paper.

## 2 MHS equilibria with gravity

The equations for static equilibrium in a gravitational field are the MHS equation and Gauss’ Law

$`p+{\displaystyle \frac{p}{h}}r={\displaystyle \frac{1}{4\pi }}(\times B)\times B,`$ (3)
$`B=0.`$ (4)
The static contribution to the pressure scale height $`h`$ is directly related to the temperature for a perfect gas, $`hp/(\rho g(r))=kT/(\overline{m}g(r))`$. The variables have their usual meanings: $`B`$ is the magnetic field vector in Gaussian units, $`p`$ the gas pressure, $`\rho `$ the density, $`T`$ the temperature, $`\overline{m}`$ the mean particle mass, $`k`$ Boltzman’s constant, $`c`$ the speed of light (used below); $`g(r)=GM_s/r^2`$ is the gravitational acceleration, with $`G`$ the universal gravitational constant and $`M_s`$ the stellar mass. The coordinate $`r`$ measures the distance from the system or stellar center, and the formal notation $`r`$ is adopted for algebraic convenience to denote the radially directed unit vector.

The set of four vector-element Eqs. (3) and (4) contain five unknowns, $`p`$, $`h`$, and the three vector elements of $`B`$. However, the solution space for $`B`$ is more restrictive than would be obtained with linear relations, as we discuss in this Section and elaborate in Appendix A.

The gravitational term in Eq. (3) compared to the pressure gradient is of order $`d/h`$, where $`d`$ is a characteristic flux-tube thickness; for the main power in threads in the solar corona $`d/h<10^3`$. In the vicinity of non-force-free fields gas-pressure changes must be mainly magnetically determined, and only far away can the gravitational gradient be significant. Where the Lorentz term is negligible, the gas pressure is hydrostatic and defined by a single radial temperature profile. Since the gravitational term is small in non-force-free conditions, the classical equilibrium solutions without gravity described in Appendices A and B are applicable, and Eqs. (1) and (2) are valid in every small volume, but with the allowance that quantities may vary slowly spatially. Thus we obtain coupled atmospheric and magnetic height variations, at least if certain types of flux tubes are present as described in the Introduction.

Conditions along field lines are represented by a parallel-field equation

$$B\left(p+\frac{p}{h}r\right)=0,$$
(5)
which is the projected component of the MHS Eq. (3) along $`B`$. The relation says that gas pressure changes along field lines must be hydrostatic, which can be seen by introducing the integrating factor $`\varphi `$ into the gas pressure

$$p=\varphi \widehat{p},$$
(6)
for an unrestricted function $`\widehat{p}`$. Upon substitution

$$B\left(\varphi \widehat{p}+\widehat{p}\left(\varphi +\frac{\varphi }{h}r\right)\right)=0,$$
(7)
so taking for the integrating factor

$$\mathrm{ln}\varphi =\frac{1}{h}r,$$
(8)
we obtain $`B\widehat{p}=0`$, or $`\widehat{p}`$ constant along field lines. The only solutions to Eq. (8) are the 1D $`\varphi =\varphi (r)`$ and $`h=h(r)`$ or $`T=T(r)`$, represented in the hydrostatic relation

$$\frac{\mathrm{ln}\varphi (r)}{r}=\frac{1}{h(r)},$$
(9)
which can be integrated to give

$$\varphi (r)=\mathrm{exp}\left(_{r_s}^r\frac{\mathrm{d}r}{h(r)}\right),$$
(10)
where $`r_s`$ denotes an arbitrary base height for the radial variations. Thus pressure variations along arbitrarily directed field lines follow the 1D hydrostatic scale-height function $`\varphi (r)`$ from a constant base pressure $`\widehat{p}`$; 3D effects are possible because the base pressure $`\widehat{p}`$ and the scale-height function $`\varphi (r)`$ or temperature profile $`T(r)`$ can vary from field line to field line. One scale-height function $`\varphi (r)`$ with temperature profile $`T(r)`$ applies throughout force-free regions and in isolated non-force-free flux tubes in thermal balance with their local surroundings.

Substituting Eq. (6) into the MHS Eq. (3) and using Eq. (8) gives the cross-field equation

$$\varphi (r)\widehat{p}=\frac{1}{4\pi }(\times B)\times B,$$
(11)
for $`\widehat{p}`$ allowed to vary from field line to field line. We adopt the single scale-height function $`\varphi (r)`$, assuming for our argument that one temperature profile applies throughout. At a base height $`r_s`$, $`\varphi (r_s)=1`$, and the cross-field Eq. (11) resembles the classical MHS Eq. (24) without gravity. Since the base height $`r_s`$ is arbitrary, it is again evident that the classical equilibrium solutions must be applicable in every small volume.

Since the base pressure $`\widehat{p}`$ is constant along field lines, the variation in the scale-height function $`\varphi (r)`$ on the left side of Eq. (11) must be reproduced on the right side in the Lorentz term, and so in $`B^2`$ too. Thus $`B`$ must vary like $`\varphi (r)^{1/2}`$ along field lines in non-force-free fields within $`d/h`$, which is the order of the derivatives of the scale-height function. Formally we can absorb $`\varphi (r)`$ into the magnetic field on the right side by decomposing $`B`$ into the general product

$$B=\varphi _B(r)\widehat{B}=\varphi _{B2}(r)^{1/2}\widehat{B}=\varphi (r)^{1/2}\zeta (r)\widehat{B}.$$
(12)
A 1D magnetic scale-height function $`\varphi _B(r)`$ is sufficient to cancel $`\varphi (r)`$; we introduce the residual radial multiplier $`\zeta (r)`$ to allow differing scale-height functions, leaving $`\widehat{B}`$ an unrestricted vector function. Substituting back into Eqs. (11) and (4) eliminates $`\varphi (r)`$ and gives

$`\widehat{p}={\displaystyle \frac{1}{4\pi }}\left(\left({\displaystyle \frac{r}{2h(r)}}\right)\times \zeta (r)\widehat{B}\right)\times \zeta (r)\widehat{B},`$ (13)
$`\zeta (r)\widehat{B}={\displaystyle \frac{r}{2h(r)}}\zeta (r)\widehat{B}.`$ (14)
The equations can be solved by substituting perturbative expansions in powers of the small quantity $`d/h(r)`$, $`\widehat{p}=_n\widehat{p}_n(d/h(r))^n`$ and $`\widehat{B}=_n\widehat{B}_n(d/h(r))^n`$ for $`n0`$ with $`d`$ a constant characteristic flux-tube thickness. The expansions separate according to powers of $`(d/h(r))^n`$, giving the most-significant order-zero equations

$`\widehat{p}_0={\displaystyle \frac{1}{4\pi }}\left(\times \zeta (r)\widehat{B}_0\right)\times \zeta (r)\widehat{B}_0,`$ (15)
$`\zeta (r)\widehat{B}_0=0.`$ (16)
Equations (15) and (16) in $`\widehat{p}_0`$ and $`\zeta (r)\widehat{B}_0`$ are equivalent to the classical plane-parallel MHS equations without gravity, which are discussed in Appendix A. The solutions are the classical equilibria, which exhibit a degenerate direction $`z`$ in the local coordinate system $`(x,y,z)`$ taken to be arbitrarily oriented with respect to the radial $`r`$. There is no loss in generality in just taking $`\zeta (r)=1`$ and $`\varphi _{B2}(r)=\varphi (r)`$ in Eq. (12) leaving $`\widehat{B}_0`$ the unperturbed solution from Eq. (28) and constant in $`z`$. If $`z⟂̸r`$ the solutions for $`B`$ exhibit a variation along the field direction that goes like $`\varphi (r)^{1/2}`$. In the case considered by Dungey (1953) with $`zr`$, distortion effects are introduced perpendicular to the main magnetic field direction as described in the Introduction. Thus $`h_B(r)=2h(r)`$ at every $`r`$, where the magnetic scale height is defined in the usual way, $`h_B(r)=[\mathrm{ln}\varphi _B(r)/r]^1`$. The scale heights are coupled to most-significant order in the perturbing gravitational field in all non-force-free MHS equilibria containing thin flux tubes in thermal balance with their local surroundings.

Additional perturbative equations for the base functions $`\widehat{p}`$ and $`\widehat{B}`$ originate with the rescaling of the magnetic field vector in $`r`$ as discussed in Appendix C. The rescaling leads to the small term $`r/(2h(r))`$ in both Eqs. (13) and (14). The term appears on the right side of Eq. (14) and is needed since a slow radial decrease in $`|B|`$ cannot occur without a divergence of field lines; it represents a compensating creation of flux in $`\widehat{B}(r)`$ with $`r`$. The term is subtracted from the curl in the Lorentz force in Eq. (13) and corresponds to a small residual pseudo-force directed between $`\widehat{B}`$ and the outward radial $`r`$; it represents a progressive left-handed twist of field lines at the rate of 1 radian in every magnetic scale height $`2h(r)`$. In MHS equilibria with gravity, the amplitude, the divergence, and the curl of $`B`$ all vary on the large scale $`2h(r)`$.

In principle anyway, a single thin flux tube somewhere in the atmosphere, evacuated or in thermal balance with its surroundings, fixes the equilibrium gas pressure for the entire atmosphere, but all non-force-free fields must be consistent to satisfy mutually the horizontal pressure boundary condition. In reality of course, magnetic fields may exhibit local deviations, and the equilibrium state might only be reflected in the average in an atmosphere containing many flux tubes.

A flux tube in thermal balance satisfies the condition $`p_e(r)B_e(r)^2`$ whether of increased or decreased field strength compared to the local background magnetic field. A flux tube with increased field strength $`|B_i(r)|>|B_e(r)|`$, must exhibit a decreased internal gas pressure in Eq. (1), $`p_i(r)<p_e(r)`$ consistent with the extreme dark-thread evacuated case of $`p_i(r)=0`$, which requires the minimum surrounding atmosphere for a given magnetic pressure difference. Flux tubes with less field strength than the external magnetic field $`|B_i(r)|<|B_e(r)|`$, must have an increased gas pressure $`p_i(r)>p_e(r)`$, consistent with bright threads.

Threads appear to be isolated thin flux tubes aligned with the background magnetic field and embedded in current sheets in a bimodal amplitude distribution, consistent with the natural equilibrium form discussed in Appendix B. The temperature of bright threads was indeed found to be the same as the external surroundings within observational uncertainties (NK96). The solar EUV coronal temperature for unresolved quiet and active regions is relatively uniform, 1.5 – 2.1 $`\times 10^6`$ K (Withbroe 1975), and a narrow temperature range around $`2\times 10^6`$ K is found for a broad range of sizes and densities of X-ray loops for all but prominence flare loops (Davis et al. 1975; Vaiana et al. 1976; Rosner et al. 1978). Large-area coronal temperatures derived from the solar white-light scale height and forbidden line ratios give temperatures somewhat lower than the EUV and X-ray temperatures, in the range of 1 – 2 $`\times 10^6`$ K at the solar radial height $`r=1.15r_{}`$ (Guhathakurta et al. 1992).

A slightly heated (or cooled) flux tube exhibits an increased (or decreased) internal pressure $`p_i(r)`$ with height and thereby must affect a decreased (or increased) internal field strength $`|B_i(r)|`$ or an increased (or decreased) external pressure $`p_e(r)`$ to satisfy the equilibrium Eq. (1). Small changes in internal field strength might be accommodated by changes in the flux-tube size as a state of local stress in the large-scale magnetic field, and changes in the surrounding ambient pressure $`p_e(r)`$ might be modified by compensating external dynamics like we describe in Sect. 4. Hot flux tubes should tend to cool to the temperature of their immediate surroundings, bringing them back to the normal MHS equilibrium condition. A relatively cool nonevacuated flux tube with excess field strength tends to evacuation with height as a different equilibrium condition, and so might just persist out of equilibrium within about one internal scale height of its base.

Continuous 3D non-force-free equilibrium solutions are claimed (Low 1985, 1991), but these allow an arbitrary hydrostatic 1D gas pressure independent of the magnetic field (see the use of $`p_0(r)`$ and discussion around Eq. (24) of Low 1991, in Sect. III in Bogdan & Low 1986, or Sect. 4 in Neukirch 1997). In the general force balance, that is the cross-field Eq. (11) with $`\varphi (r)`$ allowed to vary from field line to field line, the pressure and magnetic variations cannot be separated except in force-free regions unconstrained by a horizontal pressure boundary condition, and so necessarily lacking isolated thin non-force-free flux tubes evacuated or in thermal balance with their surroundings.

We have confined our study here to the static solutions, but it is straightforward to broaden consideration to include certain types of steady flows. Steady uniform flows along field lines add a flow pressure to the parallel-field Eq. (5), which leads to a modified hydro-steady relation in Eqs. (9) and (10), but the same cross-field Eq. (11) results. Thus we obtain a gas pressure containing static and wind components coupled to the large-scale magnetic field outside flux tubes. For a given gas pressure variation or scale-height function $`\varphi (r)`$ as is defined by the magnetic field in a magnetically dominated environment, the static contribution to the pressure scale height $`h(r)`$ and corresponding temperature $`kT(r)=\overline{m}g(r)h(r)`$ are always larger with a steady flow than without (Parker 1960). Steady flux-tube flows can modify the coupling too.

## 3 Model atmospheres

The derivation of an explicit static equilibrium solution is a formidable task because such an analysis must consider the specific form that distorted magnetic fields take in the presence of a gas. However an approximate coupled scale-height function can be written

$$\stackrel{~}{\varphi }(r)=\varphi _a(r)+\frac{1}{\beta _s}\varphi _{B2}(r),$$
(17)
where $`\varphi _a(r)`$ and $`\varphi _{B2}(r)`$ are nominal separate gas and magnetic pressure solutions to the MHS Eqs. (3) and (4) normalized at a base height $`r_s`$; the coefficient $`1/\beta _s`$ defines the relative magnetic-field strength. The approximate coupled scale-height function $`\stackrel{~}{\varphi }(r)`$ is proportional to the total pressure $`p+B^2/(8\pi )`$ of the nominal separate solutions. The gradient of the total pressure approximates the gas-pressure gradient minus the Lorentz term in the MHS equilibrium Eq. (3) reflecting the local balance in Eq. (1), as discussed at the end of Appendix B. The coupled scale-height function $`\stackrel{~}{\varphi }(r)`$ thus gives a correct MHS solution where one pressure component is dominant, $`pB^2/(8\pi )`$ or $`B^2/(8\pi )p`$.

At a base height $`r_s`$ taken as the stellar photospheric surface, we assume that magnetic fields are relatively weak, $`\beta _s1`$ and $`\stackrel{~}{\varphi }(r)\varphi _a(r)`$. There the gas pressure is dominant and must follow its nominal atmospheric form not being much affected by the magnetic field, whereas the magnetic field lines must be highly distorted. The nominal atmospheric-pressure scale-height function $`\varphi _a(r)`$ is defined by the hydrostatic relation Eq. (10) with a temperature found using a transfer equation that includes all of the usual energy input and loss mechanisms.

The atmospheric contribution $`\varphi _a(r)`$ characteristically falls off much more rapidly than the magnetic pressure $`\varphi _{B2}(r)`$, so around some transition height $`r_t`$ where $`\varphi _a(r_t)=\varphi _{B2}(r_t)/\beta _s`$, $`\stackrel{~}{\varphi }(r)`$ changes from $`\varphi _a(r)`$ to $`\varphi _{B2}(r)/\beta _s`$ over about one gas-pressure scale height. Above $`r_t`$ the magnetic field is dominant $`\stackrel{~}{\varphi }(r)\varphi _{B2}(r)/\beta _s`$, and the field must follow its nominal form irrespective of the gas pressure. There the gas pressure profile may be very distorted from its nominal form, representing an added atmospheric heating produced by MHS restoring flows as we discuss in Sect. 4.

Both the magnetic and gas pressure variations must be distorted from their nominal separate forms around the intermediate transition height $`r_t`$. Preserving the total pressure $`p+B^2/(8\pi )`$ gives an approximate flux-conserving extrapolation below the magnetically dominated corona and an ostensible gas-pressure extrapolation above the photosphere and chromosphere. The total relative pressure around the intermediate transition height goes from $`A_{\mathrm{tot}}(\varphi _a(r)+\varphi _{B2}(r)/\beta _s)`$ for a nominal superposition of independent gas and magnetic pressures $`\varphi _a(r)`$ and $`\varphi _{B2}(r)/\beta _s`$ with $`A_{\mathrm{tot}}`$ the total surface area, to $`(A_{\mathrm{tot}}A_{\mathrm{vac}})\stackrel{~}{\varphi }(r)+A_B\stackrel{~}{\varphi }(r)`$ for the coupled pressure $`\stackrel{~}{\varphi }(r)`$ with $`A_{\mathrm{vac}}`$ the evacuated area and $`A_B`$ the magnetically filled area. The two pressure totals are equal in general only with $`\stackrel{~}{\varphi }(r)`$ from Eq. (17) and when the field-filling regions are evacuated with $`A_B=A_{\mathrm{vac}}`$. Anyway the choice of coupled scale-height function is not too critical for our demonstration, as smooth switching functions with correct asymptotic behavior in the photosphere and corona exhibit similar temperature profiles even around the transition height.

For demonstration purposes, we take $`\varphi _a(r)`$ for a polytrope atmosphere as derived in Appendix D

$$\varphi _a(r)=\left(1\frac{\mathrm{\Gamma }1}{\mathrm{\Gamma }}\frac{r_s}{h_s}\left(1\frac{r_s}{r}\right)\right)^{\frac{\mathrm{\Gamma }}{\mathrm{\Gamma }1}},$$
(18)
where $`h_sh(r_s)=kT_s/(\overline{m}g_s)`$ denotes the surface scale height and $`\mathrm{\Gamma }`$ the ratio of specific heats or polytrope adiabat.

A nominal scale-height function $`\varphi _{B2}(r)`$ is defined for a single-spherical-harmonic potential magnetic field. The potential magnetic form seems to be consistent with what is observed in the lower solar corona for $`r<1.6r_{}`$ (Altschuler & Newkirk 1969; Schatten et al. 1969). The potential magnetic field vector is the gradient of a sum of scalar spherical-harmonic component functions, each of which is the separable product of a 2D surface function and radial multiplier $`1/r^{\mathrm{}+1}`$. The resulting magnetic field in each spherical harmonic goes like $`1/r^{\mathrm{}+2}`$ in all its vector elements; a monopole $`\mathrm{}=0`$ exhibits a $`1/r^2`$ radial falloff, a dipole field $`\mathrm{}=1`$, a $`1/r^3`$ falloff, etc. For a single spherical harmonic $`\mathrm{}`$, the magnetic energy density falls off like

$$\varphi _{B2}(r)=\left(\frac{r_s}{r}\right)^{2\mathrm{}+4}.$$
(19)
With a superposition of spherical harmonics, a nonuniform radial dependence can occur, and a representative global-average spherical harmonic $`\mathrm{}`$ might be obtained by appropriately weighting flux-tube locations.

The temperature is written from the hydrostatic Eq. (9)

$$kT(r)=\overline{m}g(r)h(r)=\overline{m}g(r)\left(\frac{1}{\varphi (r)}\frac{\varphi (r)}{r}\right)^1.$$
(20)
Taking $`\varphi (r)`$ to be the approximate $`\stackrel{~}{\varphi }(r)`$ from Eq. (17) with $`\varphi _a(r)`$ from Eq. (18) and $`\varphi _{B2}(r)`$ from Eq. (19), we obtain an explicit formula for $`T(r)`$. In the corona $`\varphi (r)\varphi _{B2}(r)/\beta _s`$, and the formula exhibits the limiting scale height

$$h_{\mathrm{cor}}(r)=\frac{r}{2\mathrm{}+4},$$
(21)
and temperature

$$kT_{\mathrm{cor}}(r)=\frac{GM_s\overline{m}}{(2\mathrm{}+4)r}=\frac{\overline{m}g(r)r}{2\mathrm{}+4}.$$
(22)
The limiting base coronal temperature $`T_{\mathrm{cor}}(r_s)`$ is defined entirely by the stellar surface gravity $`g_s`$, radius $`r_s`$, mean particle mass $`\overline{m}`$, and main magnetic spherical harmonic $`\mathrm{}`$. For a toroidal field $`\mathrm{}=2`$, the scale height at the stellar surface is $`h_{\mathrm{cor}}(r_s)=r_s/8`$, which gives a solar base coronal temperature of $`T_{\mathrm{cor}}(r_s)=1.73\times 10^6`$ K (using $`\overline{m}=0.6m_p`$ for $`m_p`$ the proton mass, $`\overline{m}=1\times 10^{24}`$ g, $`r_s=r_{}=6.96\times 10^{10}`$ cm, $`g_s=2.75\times 10^4`$ cm s<sup>-2</sup>, and $`k=1.38\times 10^{16}`$ ergs K<sup>-1</sup>). With a solar base coronal temperature of $`T_{\mathrm{cor}}(r_s)=1.6\times 10^6`$ K, we obtain the spherical harmonic $`\mathrm{}=2.33`$.

Figure 2 shows some example scale-height functions $`\stackrel{~}{\varphi }(r)`$ from Eqs. (17) – (19) in the lowest part of the solar atmosphere around the transition height $`r_t`$. Scale-height functions with the polytrope adiabat $`\mathrm{\Gamma }=1.1`$ and spherical harmonic $`\mathrm{}=2.33`$ are plotted for different $`\beta _s`$, using the solar parameters and surface photospheric temperature $`T_s=`$6820 K from the VAL model atmospheres (Vernazza et al. 1981), which defines the surface scale height $`h_s=kT_s/(\overline{m}g_s)`$ for $`\varphi _a(r)`$ in Eq. (18). As all of the curves are based on the same $`\mathrm{\Gamma }`$, they coincide until the transition height $`r_t`$ for the model is reached, and then switch rapidly to the nominal magnetic $`\varphi _{B2}(r)`$, which appears to be relatively flat on the log scale.

The radial expansion of a flux tube is another way to visualize the field strength decrease in the atmosphere and the properties of the coupled scale-height function. For conservation of the total flux through the cross-sectional area of a flux tube $`\pi (d/2)^2|B|`$, the flux-tube diameter $`d(r)`$ must increase with radial distance $`r`$, $`d(r)|B|^{1/2}\varphi (r)^{1/4}`$. Figure 3 portrays the relative flux-tube diameter as a function of height $`d(r)/d(r_s)`$ for different solar atmospheric models denoted by $`\mathrm{\Gamma }`$ and $`\beta _s`$.

For a given polytrope adiabat $`\mathrm{\Gamma }`$, all flux tubes exhibit approximately the same relative shape up to a height that depends upon $`\beta _s`$, where the field lines straighten up. Stronger fields exist with a lower transition height and exhibit less overall relative expansion before straightening up. The canopy depends upon the surface distribution of the fields, but strong fields might give the appearance of a lower canopy height too consistent with observations of sunspots (Giovanelli & Jones 1982). The straightening height is somewhat below the location of the fastest temperature change due to the differing dependencies: the field strength, which goes like $`\varphi (r)^{1/2}`$, corrects from its atmospherically determined form at $`r_t`$, but the most rapid temperature change occurs significantly higher where the slope of the scale-height function flattens out, as $`T(r)(\mathrm{ln}\varphi (r)/r)^1`$.

Figure 4 illustrates the temperature variation in the lower solar atmosphere using Eq. (20) with $`T_s=`$6820 K from the VAL atmospheres. The upper panel shows the model atmospheres with $`\mathrm{}=2.33`$ and $`\beta _s=`$5000 for different polytrope adiabats $`\mathrm{\Gamma }`$, and the lower panel the model atmospheres with $`\mathrm{}=2.33`$ and $`\mathrm{\Gamma }=1.1`$ for different $`\beta _s`$. The VAL A and F model temperatures for relatively cool inner network and hot network bright points are also shown.

In the solar photosphere, the nominal atmospheric polytrope follows the VAL temperature roughly, giving the best tradeoff between photospheric and chromospheric temperatures with $`\mathrm{\Gamma }1.1`$. Of course the features of a real solar atmosphere can never be well approximated by a polytrope. The shape of the temperature function through the chromosphere and transition height, where the scale-height function $`\stackrel{~}{\varphi }(r)`$ goes from $`\varphi _a(r)`$ to $`\varphi _{B2}(r)/\beta _s`$, depends upon $`\mathrm{\Gamma }`$ and $`\beta _s`$ but not on $`\mathrm{}`$, which separately determines the base coronal temperature. Temperature profiles with $`\beta _s5000`$ or 1000 give best agreement around the transition height with the VAL A or F atmospheres, respectively, averaging approximately through the Ly$`\alpha `$ plateau.

For flux-tube evacuation in the photosphere $`B(r_s)^2=8\pi p(r_s)`$, and we obtain a maximum field strength of $`|B(r_s)|=`$ 1715 G using $`p(r_s)=1.17\times 10^5`$ dynes cm<sup>-2</sup> from the VAL solar models; $`1/\beta _s^{1/2}`$ might be taken as the fractional area covered by evacuated magnetic fields. Strictly the coefficient $`1/\beta _s`$ represents the coronal base gas or non-force-free magnetic pressure extrapolated back to $`r_s`$, which can be no larger than the average surface magnetic pressure. Thus the fractional area covered by evacuated magnetic fields at $`r_s`$ must be at least $`1/\beta _s^{1/2}`$, or $`1/71`$ for $`\beta _s=5000`$ and $`1/32`$ for $`\beta _s=1000`$, which corresponds to a minimum base coronal field strength of $`|B_{\mathrm{cor}}|=`$ 24 G for $`\beta _s=5000`$ and $`|B_{\mathrm{cor}}|=`$ 54 G for $`\beta _s=1000`$, within the range of solar observations (Lin et al. 2000). A lower transition region occurs where the photospheric flux-tube areal coverage and coronal magnetic field strength are larger near active regions, consistent with the known tendency.

The solar atmospheric model of Fontenla et al. (1990), which includes particle diffusion and conduction effects, lacks the Ly$`\alpha `$ plateau and exhibits a much more abrupt transition region than what we obtain. While radiative, diffusive, and conductive losses must be largely balanced by MHS restoring flows in the corona, as we discuss in Sect. 4, loss mechanisms can influence the detailed shape of transition-region profiles. Real radiative models coupled to more general magneto-hydro-steady equilibrium solutions need to be developed, but such work may be complicated by the intrinsic limitations and uncertainties in atmospheric modeling, especially associated with inhomogeneous magnetic structure (Ayres 1981; Carlsson & Stein 1995). It is possible that general methods based upon transition-region emissivity profiles, which have been used to determine the nonradiative atmospheric heating contribution (Craig & Brown 1976; Anderson et al. 1996), might be able to distinguish coronal heating by MHS restoring flows from other heating mechanisms.

Figure 5 illustrates the solar atmospheric temperature function $`T(r)`$ for various values of $`\mathrm{}`$ with the coronal temperatures inferred from the white-light-intensity radial gradient taken from eclipse photographs (Newkirk et al. 1970) and inferred from FeXIV 5303Å line-width measurements (Jarrett & von Klüber 1958). The limited resolution in $`r`$ in the figure hides the transition region near $`r_s`$. The coronal model for the spherical harmonic $`\mathrm{}=1`$ seems to give the best overall agreement with the measured radial profiles, but $`\mathrm{}2.3`$ better matches the base coronal temperature. However all of the model curves show a more rapid falloff than the measured temperature profiles. Models containing an outward wind should give systematically higher temperatures with height in the corona.

## 4 MHS restoring flows

Under usual conditions, e.g. quasi-steady, the large-scale and long-term system average is represented by the equilibrium state. In the quasi-steady approximation, the system evolution can be described by a sequence of nearby equilibria. Low (1977) shows that the quasi-steady approximation is applicable when magnetic adjustments are fast compared to evolutionary processes and the magnetic field is dominant over the gas, as in the lower solar corona.

A qualitative feature of MHD processes in a slowly evolving dominant magnetic field where reconnections are fast is an atmospheric restoring flow. Hydro-steady equilibria are always self-restoring. In the absence of sufficient external gas or flow pressure, a time-dependent outward gas acceleration arises outside and along rigidly held non-force-free flux tubes in the overall MHD balance, as described in Appendix E. Flows in a fixed magnetic field act generally to reestablish the equilibrium atmosphere, heating the gas compressively. As long as the timescale for magnetic evolution is relatively long, such MHD processes can have a large-scale effect since the pressure perturbation that drives the flow propagates horizontally away from flux tubes at the sound speed.

The available energy for outflows must be limited and the atmosphere has to eventually cool, distorting and geometrically stressing the coupled magnetic field, which contributes to magnetic instability. The magnetic field can be perturbed in other ways too, for example by new flux coming up from the convection zone below or by photospheric twisting and flux-tube motions. Magnetic evolution is complicated since unstable equilibria are generally possible where perturbations reduce the overall energy of the system. Also local deviations may persist within field allowances before equilibria can be resolved. However if the new magnetic equilibrium after fast reconnection is mostly force-free containing rigidly held non-force-free flux tubes, outflows again must reestablish the equilibrium gas-pressure profile of a depleted cool atmosphere, applying the argument from Appendix E.

MHS restoring flows may provide an important coronal heating mechanism, which can exist consistently along with other heating processes. By offseting the total coronal energy requirement, other heating sources should add stability to the large-scale magnetic field and thus longevity to the magnetic evolution. The total coronal energy requirement is estimated conventionally using heat-flux scaling relations for loops and homogeneous regions (Rosner et al. 1978; Hearn & Kuin 1981; Hammer 1982). If sufficient magnetic energy is not available from evolutionary processes to offset losses, then pure force-free cool equilibria may be the only possibility.

Localized heating in flux tubes, which is thought to be an important coronal heating process (Vaiana & Rosner 1978; Airapetian & Smartt 1995), may affect ambient MHS restoring flows. Heated flux tubes must exhibit field-line perturbations that stress the large-scale magnetic field and tend to produce outflows in the surroundings at least on the average, which heat the ambient atmosphere above the normal equilibrium temperature. The ambient hydrostatic imbalance can have a relatively large scale of influence since it propagates horizontally away from its source much faster than thermal diffusion effects. Such post-reconnection hydrodynamic heating acting throughout the coronal volume may help overcome the often-cited difficulty that reconnection events are so spatially localized and short-lived that only a larger event rate than what is directly inferred from observations could account for coronal heating (Heyvaerts & Priest 1984; Hudson 1991).

Thus MHD processes may lead to refilling outflows, a ‘magnetic suction’, which might be in evidence as the large-scale flows seen in the vicinity of coronal voids (Wagner et al. 1983), or be responsible for the supergranular-scale chromospheric upward velocities interpreted from post-flare solar spectral-line blue shifts (Schmieder et al. 1987; Cauzzi et al. 1996), or appear as spicule eruptions in the chromosphere. The upward mass flux produced by solar spicule events is about 100 times the total solar-wind mass flux, but represents only a small fraction of its total energy flux (Pneuman & Kopp 1978).

The remarkable feature of an MHS corona is that the temperature is essentially a geometric parameter defined by the stellar surface gravity $`g_s`$, radius $`r_s`$, mean particle mass $`\overline{m}`$, and magnetic spherical harmonic $`\mathrm{}`$

$$kT_{\mathrm{cor}}(r_s)=\frac{\overline{m}g_sr_s}{2\mathrm{}+4}.$$
(23)
The observed average solar base coronal temperature of $`1.6\times 10^6`$ K corresponds to a spherical harmonic of $`\mathrm{}=2.33`$ consistent with the large-scale solar magnetic field, which exhibits substantial power in the toroidal $`\mathrm{}=2`$ and spherical harmonics $`\mathrm{}=3`$ or 4 as evidenced by the presence of active longitudes. With a general mixture of spherical harmonics, discrepant horizontal gas-pressure boundary conditions at flux tubes can arise, and a nonuniform large-scale surface temperature distribution is possible as a deviation from the pure MHS equilibrium state.

Formulae for the average base coronal temperature like Eq. (23) have been suggested already in the hypothesis of a ‘geometric boundary condition’, relevant for different kinds of heating models where the base pressure scale height is proportional to the stellar radius (Menzel 1968; Scudder 1992). The relation appears to give reasonable coronal temperatures for a range of stellar types (Williams & Mullan 1996).

###### Acknowledgements.

The author is grateful to many for useful discussion during the course of this work, especially to Eric Priest and Ray Smartt for their careful reviews of the manuscript, and to the referee for many insightful remarks and helpful suggestions.

## Appendix A Local MHS equilibria

There are well-known solutions to the classical static problem describing a gas existing in conjunction with elongated magnetic fields without a perturbing gravitational field (Parker 1979, Chapter 6). It is helpful to revisit the classical problem here in a general way. The classical MHS equilibrium equations without gravity are the static equilibrium equation

$$p=\frac{1}{4\pi }(\times B)\times B,$$
(24)
and Gauss’ Law, Eq. (4).

Plane-parallel solutions can be developed based upon the general Cartesian vector function

$$B=(\frac{a(x)}{y},\frac{a(x)}{x}\frac{b(x)}{z},\frac{b(x)}{y}).$$
(25)
The $`x`$ and $`z`$ vector elements of $`B`$, $`a/y`$ and $`b/y`$, are taken to be arbitrary functions, and the $`y`$ element is written so that Gauss’ Law is always satisfied. It is convenient to use derivative forms in the vector elements to avoid complicating integrals in the expression.

Any solutions to Eq. (24) must satisfy the two conditions $`Bp=0`$ and $`(\times B)p=0`$, which are written with $`B`$ from Eq. (25)

$`\mathrm{C}_{xy}a\mathrm{C}_{yz}b=0,`$ (26)
$`\left(\mathrm{C}_{xz}{\displaystyle \frac{}{x}}+\mathrm{C}_{yz}{\displaystyle \frac{}{y}}\right)a+\left(\mathrm{C}_{xy}{\displaystyle \frac{}{y}}+\mathrm{C}_{xz}{\displaystyle \frac{}{z}}\right)b=0,`$ (27)
using the linear commutation operators $`\mathrm{C}_{xy}\frac{p}{x}\frac{}{y}\frac{p}{y}\frac{}{x}`$, etc., which contain $`p(x)`$ as an implicit function.

The two Eqs. (26) and (27) can be used to write two of the functions in terms of the third, for example $`a(x)`$ and $`b(x)`$ in terms of $`p(x)`$. The third relation, Eq. (24) in the vector direction perpendicular to both $`B`$ and $`\times B`$, gives a differential relation for the third quantity. The separate vanishing of every commutator term is a reduction that is satisfied with certain general functional dependencies between the variables and leads to a well-known self-consistent form for the third relation, the Grad-Shafranov Equation. Other possibilities are not considered here.

Taking every commutator term to vanish in Eqs. (26) and (27), but without other reduction, gives $`a(x)`$, $`b(x)`$, and spatial derivatives of $`a(x)`$ and $`b(x)`$ as functions of $`p(x)`$ alone, which is only possible with specific spatial forms for real functions, so other reductions are required in general. Taking $`/z=0`$ gives a consistent solution defined by the functional dependencies $`a=a(p)`$ in Eq. (26) and $`\frac{b}{y}=\frac{b}{y}(p)`$ in Eq. (27). A similar solution is obtained by taking $`/x=0`$. These two solutions or linear combinations are equivalent with rotation of the $`x`$$`z`$ axes around $`y`$. Taking $`z`$ as the degenerate direction gives the usual 2D solution, ordinarily written with $`a`$ as the implicit function

$$B=(\frac{a(x,y)}{y},\frac{a(x,y)}{x},B_z(a(x,y))).$$
(28)
The degenerate $`z`$ is a direction of field elongation. The magnetic potential function $`a(x,y)`$ defines a 2D planform for the solutions, which determines the spatial variations in the longitudinal field element $`B_z=B_z(a(x,y))`$ and gas pressure $`p=p(a(x,y))`$ as arbitrary 1D mappings. The magnetic field is everywhere perpendicular to the gradient of $`a(x,y)`$, since $`Ba(x,y)=0`$.

It is widely believed that this 2D solution is the only plane-parallel one. Using perturbative expansions, Parker (1979, Sect. 14.2) shows that no nearby 3D solutions exist for bounded quantities in an infinite spatial domain. However 3D variations do arise as large-scale deviations from plane parallel (Arendt & Schindler 1988). In the common shorthand,‘2.5D’ solutions refer to perturbed solutions in the problem with gravity, which admit large-scale variations in the degenerate direction of the magnetic field.

## Appendix B Properties of local MHS equilibria

Substituting $`B`$ from Eq. (28) into Eq. (24) gives the governing equation for the implicit magnetic potential function $`a(x,y)`$ for the static equilibrium problem without gravity

$$\left(^2a+\frac{d}{da}\left(4\pi p(a)+\frac{B_z(a)^2}{2}\right)\right)a=0.$$
(29)
Where $`a0`$, the potential function must satisfy the Grad-Shafranov (GS) Equation, written

$$^2a(x,y)=P^{}(a(x,y)),$$
(30)
defining the total pressure $`P(a)=4\pi p(a)+B_z(a)^2/2`$, where $`P^{}(a)=dP(a)/da`$ .

The total pressure $`P(a)`$ and its component functions $`p(a)`$ and $`B_z(a)`$ are all one dimensional and nonlinear, but must all be well-defined everywhere in the solution domain consistent with their forms at the boundaries. Strictly, disagreeing $`z`$ boundary conditions are inconsistent with the $`z`$ independence of the solutions; differences might produce small deviations in the solutions like twist or divergence or be a source of dynamical instability. It is popular to restrict consideration to entirely force-free solutions with $`p^{}(a)=0`$; the restriction does not change the nature of the basic GS Eq. (30) for $`a(x,y)`$ but requires the specific current density $`J=\frac{c}{4\pi }B_z^{}(a)B`$ for an arbitrary $`B_z(a)`$, as is evident by expanding $`J\frac{c}{4\pi }\times B`$ using $`B`$ from Eq. (28) with $`JB`$.

The magnetic potential function $`a`$ can be seen to be constant everywhere in the local plane of $`B`$ and $`J`$. From Eq. (24), $`Jp(a)=0`$. Then $`Jp(a)=(Ja)p^{}(a)=0`$. Thus $`Ja=0`$ at least when $`p^{}(a)0`$, and when $`p^{}(a)=0`$, $`JB`$ so $`Ja=0`$ anyway. We think of $`a`$ as constant on ribbon-like current sheets, the solution surfaces being more curved along $`J`$ as it is defined by the derivatives of $`B`$.

The GS Eq. (30) has the unusual feature that it contains the 1D filter function $`P^{}(a)`$. The Fourier transform of the GS Equation in the linear case, with $`P^{}(a)=k_0^2a`$ for $`k_0`$ a constant wavenumber, requires a transform function $`\overline{a}(k_x,k_y)`$ that is zero everywhere in its 2D wavenumber domain $`k=(k_x,k_y)`$ except on a thin annulus at $`|k|=k_0`$, where arbitrary complex values Hermitian in $`\pm k`$ are allowed. Taking the Fourier transform back gives the potential function $`a(x,y)`$ as a common Bessel-function radial kernel times azimuthal factor convolved with a spatial distribution of delta-function source points. Boundary conditions on $`a(x,y)`$ constrain the azimuthal factor, leading to possible planar solutions defined by all power at one azimuth, or axisymmetric solutions with power uniformly distributed in azimuth, consistent with our visualization of current sheets and cylindrical flux tubes. Physical arguments show that nonlinear kernels are similarly constrained by boundary conditions (Vainshtein & Parker 1986).

Random spatial distributions of a common GS kernel are nonlinear GS solutions too. We take the potential function $`a(x,y)`$ to be the convolution of a distribution of sources $`D(x,y)`$ with a common 2D kernel function $`A(x,y)`$, $`a(x,y)=D(x,y)A(x,y)`$. The distribution is written $`D(x,y)=_jc_j\delta (xx_j,yy_j)`$ counting sources $`j`$ of varying strength $`c_j`$, where $`\delta (x,y)`$ denotes the 2D delta function. For a spatially incoherent or random distribution $`D`$ of equal strength sources with all $`c_j=1`$, $`D`$ raised to a power is $`D`$ alone; for equal amplitude sources with $`c_j=\pm 1`$, $`D`$ raised to an odd power is $`D`$ alone. For an analytic nonlinear driver function $`P(a)`$ the convolution factors out $`P^{}(a)=P^{}(DA)=DP^{}(A)`$, when the distribution in $`D`$ is incoherent, disjoint, and contains suitably restricted amplitudes. Then the GS Eq. (30) reduces to the same GS Equation, but for the common kernel function $`A`$ in place of $`a`$. Dark and bright threads appear to be isolated thin flux tubes aligned with the background magnetic field organized in current sheets suggestive of a bimodal amplitude distribution of a common flux-tube kernel.

Parker (1972) considers MHS equilibria that contain a relatively strong constant background field. Expanding the Lorentz term in Eq. (24) gives an alternate form for the MHS equation

$$\left(p+\frac{B^2}{8\pi }\right)=\frac{1}{4\pi }\left(B\right)B,$$
(31)
or with the superposition $`B(x)=B_e+B_{}(x)`$ for a constant background field $`B_e`$ and a spatially varying field $`B_{}(x)`$, we obtain

$$\left(p+\frac{B^2}{8\pi }\right)=\frac{1}{4\pi }\left(B_e\right)B_{}+\frac{1}{4\pi }\left(B_{}\right)B_{}.$$
(32)
Taking the divergence, applying Gauss’s Law $`B_{}=0`$, and using the property that the total pressure $`p+B^2/(8\pi )`$ is a bounded quantity in an infinite domain, we obtain

$$\left(B_e\right)B_{}=0,$$
(33)
at least to first order in the small quantity $`|B_{}|/|B_e|`$. Thus we have

$$\left(p+\frac{B^2}{8\pi }\right)0.$$
(34)
with equality to first order in $`|B_{}|/|B_e|`$. Eq. (33) proves to be accurate to all orders in the perturbative expansions, and the equilibrium Eqs. (1) and (2) are justified at least in the presence of a relatively strong constant background field, where the internal flux-tube field is defined $`B_iB_e+B_{}(=B)`$. Additional discussion on Parker’s perturbative expansion is contained in Appendix C.

The Parker theorem Eq. (33) reflects the special features of the classical equilibria in Eq. (28) in the degenerate direction of the magnetic field $`z`$. MHS equilibria allow only certain alignments for the background field. A constant magnetic field can enter into the $`z`$ vector element $`B_z(a)`$ without producing other ramifications, whereas an added offset in the $`x`$ or $`y`$ elements of $`B`$ requires adding a uniformly inclined plane to $`a(x,y)`$. Such a plane limits the solutions and thus cannot be considered general, and even appears to be precluded by the boundedness of $`a(x,y)`$: low-wavenumber components in the Fourier domain needed to represent the added plane are at odds with the Fourier ring solutions in the linear case of the GS Eq. (30), and also contrary to the form of a superposition of spatially compact axisymmetric nonlinear kernels.

## Appendix C Perturbative expansions

The less-significant order $`n>0`$ equations in the perturbative series developed from Eqs. (13) and (14) are written for $`\zeta (r)=1`$; from the MHS Eq. (13)

$$\begin{array}{c}\widehat{p}_1\frac{1}{4\pi }\left((\times \widehat{B}_0)\times \widehat{B}_1+(\times \widehat{B}_1)\times \widehat{B}_0\right)=\hfill \\ \frac{1}{8\pi d}(r\times \widehat{B}_0)\times \widehat{B}_0,\hfill \end{array}$$
(35)
or for the general order $`n>0`$

$$\begin{array}{c}\widehat{p}_n\frac{1}{4\pi }_{j=0}^n(\times \widehat{B}_j)\times \widehat{B}_{nj}=\hfill \\ \frac{1}{8\pi d}_{j=0}^{n1}(r\times \widehat{B}_j)\times \widehat{B}_{nj1}\left(1+2j\frac{h(r)}{r}\right).\hfill \end{array}$$
(36)
From Eq. (14)

$$\widehat{B}_1=\frac{1}{2d}(r\widehat{B}_0),$$
(37)
or in general for $`n>0`$

$$\widehat{B}_n=\frac{1}{2d}(r\widehat{B}_{n1})\left(1+2(n1)\frac{h(r)}{r}\right).$$
(38)
The equations are linear in the variables $`\widehat{p}_n`$ and $`\widehat{B}_n`$ based upon order $`<n`$ variables. The order-zero base functions $`\widehat{p}_0`$ and $`\widehat{B}_0`$ are 2D and constant in the degenerate magnetic field direction $`z`$; 3D effects can be introduced with the cross-product terms or with general variation in the function $`h(r)`$ with $`r`$ not perpendicular to $`z`$. For consistency in Eqs. (13) and (14), $`\widehat{B}_0\widehat{p}_0=0`$, which is the order-zero component of $`B\widehat{p}=0`$.

A wide latitude is allowed for choosing perturbative quantities, and as long as they are relatively small almost everywhere the standard procedure remains valid. Van Ballegooijen (1985, Appendix A) questions the Parker (1972) choice of perturbative quantity and equation separation procedure, but van Ballegooijen’s assigning of approximate magnitudes to individual terms to separate his equations is not an algebraic procedure, and his resulting equations have inherent contradictions as described by Parker (1987). Parker (1972) does not ascribe a particular functional meaning to his kernel $`ϵ`$, quoting from the beginning of his Sect. II: “Expand the field $`b`$ and pressure $`p`$ in ascending powers of some parameter $`ϵ`$, which is of the order of $`ϵ=|b|/|B|`$”, where $`b`$ is the small spatially varying magnetic field on the uniform background field $`B`$. Our choice for perturbative kernel $`d/h(r)`$ with $`d`$ constant leads to a convenient separation of orders.

The perturbative procedure is a functional separation method rooted in the natural independence of the basis functions $`ϵ^n`$ irrespective of their real amplitudes. If the parameter $`ϵ`$ is small almost everywhere, the series solutions will be convergent. There could be locations where the amplitudes of components seem to be of differing orders as van Ballegooijen contends. Where a normally small expansion parameter $`ϵ`$ actually becomes relatively large series expansions may be divergent, but such behavior characteristically defines the isolated singularities of differential equations as discussed in theorems on the Frobinius method for series solutions.

## Appendix D Polytrope atmosphere

A polytrope atmosphere is defined by the pressure relation $`p(r)\rho (r)^\mathrm{\Gamma }`$ for the polytrope adiabat $`\mathrm{\Gamma }`$. Thus the 1D polytrope scale height $`h(r)=p(r)/(\rho (r)g(r))`$ can be written

$$\frac{h(r)}{h_s}=\left(\frac{p(r)}{p(r_s)}\right)^{\frac{\mathrm{\Gamma }1}{\mathrm{\Gamma }}}\frac{g_s}{g(r)}=\varphi (r)^{\frac{\mathrm{\Gamma }1}{\mathrm{\Gamma }}}\frac{g_s}{g(r)},$$
(39)
making reference to a base height $`r_s`$ taken to be the stellar surface, and introducing the scale-height function $`\varphi (r)`$ for the normalized gas pressure.

Substituting $`h(r)`$ into the hydrostatic Eq. (9) defines the scale-height function for a polytrope

$$\varphi (r)^{1/\mathrm{\Gamma }}\frac{\varphi (r)}{r}=\frac{g(r)}{g_sh_s},$$
(40)
or

$$\frac{}{r}\varphi (r)^{\frac{\mathrm{\Gamma }1}{\mathrm{\Gamma }}}=\frac{\mathrm{\Gamma }1}{\mathrm{\Gamma }}\frac{g(r)}{g_sh_s},$$
(41)
giving

$$\varphi (r)=\left[\frac{\mathrm{\Gamma }1}{\mathrm{\Gamma }}\frac{g(r)}{g_sh_s}dr\right]^{\frac{\mathrm{\Gamma }}{\mathrm{\Gamma }1}},$$
(42)
and with $`g(r)=g_s(r_s/r)^2`$, we obtain

$$\varphi (r)=\left[1\frac{\mathrm{\Gamma }1}{\mathrm{\Gamma }}\frac{r_s}{h_s}\left(1\frac{r_s}{r}\right)\right]^{\frac{\mathrm{\Gamma }}{\mathrm{\Gamma }1}},$$
(43)
choosing the constant of integration so that $`\varphi (r_s)=1`$. The part of the expression contained in brackets $`[]`$ corresponds to the radial temperature variation, the thermal lapse rate. The polytrope adiabat $`\mathrm{\Gamma }`$ ranges from 1 to $`5/3`$ for an atmosphere ranging from lossy and isothermal to adiabatically stratified. For $`\mathrm{\Gamma }=1`$, the scale-height function $`\varphi (r)`$ has a limiting pure exponential decay form. For $`\mathrm{\Gamma }>1`$, the thermal lapse rate and scale-height function have the unphysical feature that they become negative at some height, showing that a polytrope of constant $`\mathrm{\Gamma }`$ can exist only over a limited range in $`r`$. Such behavior is well known in stellar-interior polytrope models and leads to convection zones of finite extent.

## Appendix E Hydrodynamics around fixed fields

If the magnetic field is dominant over the gas, and either the field remains stationary, or the magnetic adjustment time is much shorter than the hydrodynamic adjustment time, then the Lorentz term in the MHD equation can be treated as constant in time and the hydrodynamics separate from the magnetodynamics.

Fast magnetic adjustment is a remarkable feature of solar reconnection phenomena. Sporadic crossings of coronal loops in low-energy flare events (Lin et al. 1992) suggest ‘X’-type configurations, and the instability of the ‘X’ configuration is well known to form a new current sheet at near the Alfvén velocity if a nonzero resistivity can be affected in the region (Dungey 1958; Syrovatskii 1981); such a resistivity is predicted in some models (Petschek & Thorne 1967; Priest 1972). Taking a field strength of $`|B|`$400 G as representative of chromospheric magnetic fields near neutral lines in solar active regions where flares commonly occur, we obtain an Alfvén speed $`c_a=|B|/(4\pi \rho )^{1/2}10^8`$ cm s<sup>-1</sup> (with $`\rho =N\overline{m}=10^{12}`$ g cm<sup>-3</sup>, using $`N=10^{12}`$ cm<sup>-3</sup> and $`\overline{m}=10^{24}`$ g). The Alfvén speed stays fairly constant into the corona giving a collapse time of $`\tau _B=h_B/c_a`$ 160 s for the propagation of the instability through a magnetic scale height $`h_B=2h_{\mathrm{cor}}=r_s/(\mathrm{}+2)1.6\times 10^{10}`$ cm, using a primary magnetic spherical harmonic of $`\mathrm{}=`$2.33 in Eq. (21).

The hydrodynamic response to already present or newly formed isolated thin non-force-free flux tubes can be obtained by perturbative analysis of the dynamical equations projected along the magnetic field. Taking the velocity to be aligned with the magnetic field and of lesser order than the mean quantities adds a single time-derivative term to the projected MHS Eq. (3), which is written with the continuity equation

$`\rho {\displaystyle \frac{u_1}{t}}+{\displaystyle \frac{p}{z}}+\rho g(zr)=0,`$ (44)
$`{\displaystyle \frac{\rho }{t}}+{\displaystyle \frac{}{z}}\left(\rho u_1\right)=0;`$ (45)
the equations are taken with the thermodynamic energy equation $`p\rho ^\gamma `$; $`z`$ is the direction of the magnetic field; $`u_1`$ denotes the velocity amplitude; $`\gamma `$ is the ratio of specific heats, which ranges from 1 to 5/3 for lossy isothermal refilling flows to adiabatic. We now choose to use the density $`\rho `$ as a problem variable, rather than the temperature or scale height $`h`$.

Perturbing the quantities backward in time, $`p=p_0p_1`$ and $`\rho =\rho _0\rho _1`$, from a final equilibrium state $`(p_0,\rho _0)`$ leaves the hydrostatic relation Eq. (9) with $`p_0`$ coupled to the magnetic field through the cross-field Eq. (11), and three coupled perturbative equations for $`(p_1,\rho _1,u_1)`$

$`\rho _0{\displaystyle \frac{u_1}{t}}={\displaystyle \frac{p_1}{z}}+\rho _1g\mathrm{cos}\theta ,`$ (46)
$`{\displaystyle \frac{\rho _1}{t}}={\displaystyle \frac{}{z}}\left(\rho _0u_1\right),`$ (47)
and $`p_1=\gamma (p_0/\rho _0)\rho _1`$, introducing $`\theta `$ as the angle between the magnetic field direction $`z`$ and the radial $`r`$. With small positive pressure and density deficits $`p_1`$ and $`\rho _1`$, a small positive velocity $`u_1`$ is required in Eq. (46) corresponding to an upward flow opposite the normal gravitational stratification. Refilling changes the gas-pressure profile, in effect heating the gas compressively as the dominant heating term in the thermodynamic energy equation.

The unperturbed density $`\rho _0=p_0/(gh)`$ varies on a scale similar to the gas pressure $`p_0`$, since the gravity $`g`$ and temperature contained in $`h`$ vary on much larger scales; thus the derivative is approximated $`\rho _0/z\rho _0\mathrm{cos}\theta /h`$. The perturbed quantities must vary on the scale of the magnetic field $`2h/\mathrm{cos}\theta `$, giving $`(p_1,u_1)/z(p_1,u_1)\mathrm{cos}\theta /(2h)`$. Substituting for the time derivative $`/t1/\tau _u`$ lets us solve for the hydrodynamic timescale

$$\tau _u=\frac{1}{\mathrm{cos}\theta }\left(\frac{4h}{(6+3\gamma )g}\right)^{1/2},$$
(48)
which is proportional to the gravitational free-fall time through the scale height $`h`$. For a radial field $`\theta =0`$, mixed lossy conditions $`\gamma =1.2`$, and $`h=r_s/(2\mathrm{}+4)`$, we obtain $`\tau _u=`$349 s for solar surface conditions with $`\mathrm{}=2.33`$, which is just over twice the magnetic adjustment time $`\tau _B`$. The hydrodynamic adjustment time increases with more lossy conditions and smaller $`\gamma `$, but overall depends only weakly on the thermodynamics.


format: text

text:
# Comment on the paper ’Thermodynamics of two-dimensional magneto nanoparticles’ by P. Vargas, D. Altbir, M. Knobel and D. Laroze

## Abstract

It is shown that there is no bi-stability in the system of two-dimensional nanoparticles considered by P. Vargas et al. in Ref. .

In the recent Letter  the thermodynamics of two-dimensional magnetic nanoparticles in an external magnetic field has been considered and equivalent mechanical model was presented. We want to comment on their result that the equilibrium magnetization has a maximum at finite temperatures and on their mechanical analogue result.

We start with the mechanical model which is presented in Fig. 1 of . This model is desribed by the following three dimensional Lagrangian function

$`={\displaystyle \frac{m}{2}}\left(\dot{x}^2+\dot{y}^2+\dot{z}^2\right)+m\omega \left(\dot{x}yx\dot{y}\right)+mgz`$ (1)
with a constraint

$`x^2+y^2+z^2=R^2`$ (2)
where $`m`$ is the mass of the particle, $`R`$ is the radius of the ring, $`g`$ is the gravity acceleration. The second term in Eq. (1) represents an inertia contribution due to the rotation of the ring around the $`z`$-axis. By using the Legendre transformation  we get the Hamiltonian of the system in the form

$`H={\displaystyle \frac{1}{2m}}\left(p_x^2+p_y^2+p_z^2\right)\omega \left(xp_yyp_x\right)+mgz`$ (3)
From Eq. (3) we obtain the partition function in the general form

$$Z_{mech}=N\underset{\mathrm{}}{\overset{\mathrm{}}{}}𝑑p_x𝑑p_y𝑑p_z\underset{\mathrm{}}{\overset{\mathrm{}}{}}d_xd_yd_ze^{H/kT}\delta \left(\sqrt{x^2+y^2+z^2}R\right)$$

where $`N`$ is a normalization factor, $`k`$ is the Boltzman constant and the constraint (2) was taken into account. This can be written for (3) as

$`Z_{mech}=2\pi R^2N{\displaystyle \underset{0}{\overset{\pi }{}}}𝑑\theta \mathrm{sin}\theta \mathrm{exp}\{mgR{\displaystyle \frac{𝒰}{kT}}\}`$ (4)
with

$$𝒰=\omega ^2R^2/g\mathrm{sin}^2\theta \mathrm{cos}\theta $$

where the angle $`\theta `$ is measured with respect to the z-axis. This potential function coincides with the potential function introduced in  but due to the condition that the angle $`\theta `$ varies only in the interval $`(0,\pi )`$ there is no bistability in the system. Another point which is missing in the partition function $`Z_{2D}`$ presented in  is the factor $`\mathrm{sin}\theta `$ in the Jacobian of the transformation from the cartesian- to the polar co-ordinates.

Thus we see that the interpretation of the dimensionless energy ( see Fig. 2 of ) is misleading because it does not present a double well potential and therefore there is no bi- stability of the system. This can easily be seen from the fact that the two angles for minimum energy belong to two different directions of rotation: the positive angle for the anti-clockwise, the negative angle for the clock-wise rotation. With this remarks in mind it should be also obvious that there is no bifurcation for a given rotation ( Fig. 3 in ).

When Vargas et al.  analyze the thermodynamical properties of magnetic nanoparticles they used the enrgy of the single particle in the form

$$E=mH\mathrm{cos}\theta \frac{1}{2}mH_a\mathrm{sin}^2\theta $$
(5)
where $`H`$ is an external magnetic field and $`H_a`$ is the anisotropy field. It is clearly seen that in the absence of magnetic field $`H`$ Eq. (5) gives the energy of an easy-plane magnet : the magnetic moment can freely rotate in the plane perpendicular to the z-axis. Note in passing that this is in a complete agreement with the above-considered mechanical model when in the no-gravity case the particle may be situated at any point of the circular ring of radius $`R`$. We would like to remark that the desription of an easy-axes magnet is quite another problem and does not correpond to the above mechanical model.

Exactly the same partition function as given in Eq. (4) ( with appropriate change of variables) can be obtained for a magnetic particle which is described by equation (5). The reason for the pronounced discrepancy between their partition function and that in Eq. (4) is the following. The authors of  overlooked the fact that by calculating the partition function one should integrate over the phase space $`d\mathrm{\Gamma }=dpdq`$ where $`p`$ and $`q`$ are canonically conjugated variables (see e.g. ). In magnetic systems which are characterized by the classical three dimensional magnetization vectors $`\stackrel{}{M}=(\mathrm{sin}\theta \mathrm{cos}\varphi ,\mathrm{sin}\theta \mathrm{sin}\varphi ,\mathrm{cos}\theta )`$ the $`z`$-component of the magnetization, $`m=\mathrm{cos}\theta `$, and the azimuthal angle, $`\varphi `$, are canonically conjugated variables. Thus the element of the phase space volume is given by $`d\mathrm{\Gamma }=\mathrm{sin}\theta d\theta d\varphi `$ and the partition function for magnetic particles takes the form (4). When one calculates the mean magnetization, $`m`$, and the susceptibility, $`\chi `$, along the field by using the correct partition function of the form

$`Z_{magn}=2\pi N{\displaystyle \underset{0}{\overset{\pi }{}}}𝑑\theta \mathrm{sin}\theta \mathrm{exp}\{{\displaystyle \frac{mH\mathrm{cos}\theta +mH_a\mathrm{sin}^2\theta /2}{kT}}\}`$ (6)
one obtains dependences shown in Fig.1.

Thus in contrast to the claim of Ref. the correctly defined magnetization does not manifest non-monotonic behaviour as a function of temperature. The magnetic susceptibility has a maximum for $`H_a<H`$ ( not seen in Fig.1 )in the low temperature interval but for $`H_aH`$ it is a monotonic function of $`T`$.

Acknowledgements Yu. G. is grateful for the hospitality of the University of Bayreuth where this work was performed.


text:
# Inhomogeneous Dark Radiation Dynamics on a de Sitter Brane

## 1 Introduction

In the context of the intensive search for extra dimensions the Randall and Sundrum (RS) brane world scenario is specially compelling for its simplicity and depth . In this scenario the observable Universe is a 3-brane boundary of a non-compact $`Z_2`$ symmetric 5-dimensional Anti-de Sitter (AdS) space. The matter fields are restricted to the brane but gravity exists in the whole AdS bulk and is bound to the brane by the warp of the infinite fifth dimension.

In recent years numerous studies have been conducted within the RS scenario \[see  for a recent review and notation\]. From an effective 4-dimensional point of view  the interactions existing between the brane and the bulk lead to a modification of the Einstein equations by a set of two distinct terms, namely a local high energy embedding term generated by the matter energy-momentum tensor and a non-local term induced by the bulk Weyl tensor. The resulting equations have a complex non-linear dynamics. For example the exterior vaccum of a collapsing distribution of matter on the brane is now filled with modes originated by the bulk Weyl curvature and can no longer be a static space .

Previous investigations of the RS scenario have been focused on static or homogeneous dynamical solutions. In this proceedings we analyse some aspects of the inhomogeneous dynamics of a RS brane world vaccum \[see  for more details\].

## 2 Vaccum Einstein Equations on the Brane

In the 4-dimensional effective geometric approach to the RS brane world scenario first introduced by Shiromizu, Maeda and Sasaki , the induced Einstein vaccum field equations on the brane are given by

$$G_{\mu \nu }=\mathrm{\Lambda }g_{\mu \nu }_{\mu \nu },$$
(1)
where $`\mathrm{\Lambda }`$ is the brane cosmological constant and the tensor $`_{\mu \nu }`$ is the limit on the brane of the projected 5-dimensional Weyl tensor. On one hand the Weyl symmetries ensure it is a symmetric and traceless tensor. On the other the Bianchi identities constrain it to satisfy the conservation equations

$$_\mu _\nu ^\mu =0.$$
(2)
The tensor $`_{\mu \nu }`$ can be written in the following general form

$$_{\mu \nu }=\left(\frac{\stackrel{~}{\kappa }}{\kappa }\right)^4\left[𝒰\left(u_\mu u_\nu +\frac{1}{3}h_{\mu \nu }\right)+𝒫_{\mu \nu }+𝒬_\mu u_\nu +𝒬_\nu u_\mu \right],$$
(3)
where $`u_\mu `$ such that $`u^\mu u_\mu =1`$ is the 4-velocity field and $`h_{\mu \nu }=g_{\mu \nu }+u_\mu u_\nu `$ is the tensor which projects orthogonaly to $`u_\mu `$. The forms $`𝒰`$, $`𝒫_{\mu \nu }`$ and $`𝒬_\mu `$ represent different characteristics of the effects induced on the brane by the free gravitational field in the bulk. Thus, $`𝒰`$ is interpreted as an energy density, $`𝒫_{\mu \nu }`$ as stress and $`𝒬_\mu `$ as energy flux.

Since the 5-dimensional metric is unknown, in general $`_{\mu \nu }`$ is not fully determined on the brane . As a consequence the effective 4-dimensional theory is not closed and to close it we need simplifying assumptions about the bulk degrees of freedom. For example we may consider a static and spherically symmetric brane vaccum with $`𝒬_\mu =0`$, $`𝒫_{\mu \nu }0`$ and $`𝒰0`$ to find the Reissner-Nordström black hole solution on the brane . Let us now show that it is possible to take a non-static spherically symmetric brane vaccum with $`𝒬_\mu =0`$, $`𝒰0`$, $`𝒫_{\mu \nu }0`$ and still close the system of dynamical equations.

Consider the general, spherically symmetric metric in comoving coordinates $`(t,r,\theta ,\varphi )`$,

$$ds^2=g_{\mu \nu }dx^\mu dx^\nu =e^\sigma dt^2+A^2dr^2+R^2d\mathrm{\Omega }^2,$$
(4)
where $`d\mathrm{\Omega }^2=d\theta ^2+\mathrm{sin}^2\theta d\varphi ^2`$, $`\sigma =\sigma (t,r)`$, $`A=A(t,r)`$, $`R=R(t,r)`$ and $`R`$ is interpreted as the physical spacetime radius. If the stress is isotropic then the traceless $`𝒫_{\mu \nu }`$ will have the general form

$$𝒫_{\mu \nu }=𝒫\left(r_\mu r_\nu \frac{1}{3}h_{\mu \nu }\right),$$
(5)
where $`𝒫=𝒫(t,r)`$ and $`r_\mu `$ is the unit radial vector, given in the above metric by $`r_\mu =(0,A,0,0)`$. Then

$$_\mu ^\nu =\left(\frac{\stackrel{~}{\kappa }}{\kappa }\right)^4\text{diag}(\rho ,p_r,p_T,p_T),$$
(6)
where the energy density and pressures are, respectively, $`\rho =𝒰`$, $`p_r=(1/3)\left(𝒰+2𝒫\right)`$ and $`p_T=(1/3)\left(𝒰𝒫\right)`$. Substituting in the conservation Eq. (2) we obtain the following expanded system

$$2\frac{\dot{A}}{A}\left(\rho +p_r\right)=2\dot{\rho }4\frac{\dot{R}}{R}\left(\rho +p_T\right),$$

$$\sigma ^{}\left(\rho +p_r\right)=p_{r}^{}{}_{}{}^{}+4\frac{R^{}}{R}(p_Tp_r),$$
(7)
where the dot and the prime denote, respectively, derivatives with respect to $`t`$ and $`r`$. A synchronous solution is obtained by taking the equation of state $`\rho +p_r=0`$, giving $`𝒫=2𝒰`$ with $`𝒰`$ having the dark radiation form

$$𝒰=\left(\frac{\kappa }{\stackrel{~}{\kappa }}\right)^4\frac{Q}{R^4},$$
(8)
where the dark radiation tidal charge $`Q`$ is constant. Consequently, we obtain

$$G_{\mu \nu }=\mathrm{\Lambda }g_{\mu \nu }+\frac{Q}{R^4}\left(u_\mu u_\nu 2r_\mu r_\nu +h_{\mu \nu }\right),$$
(9)
an exactly solvable closed system which depends on the two free parameters $`\mathrm{\Lambda }`$ and $`Q`$. Indeed, its solutions can be written in the LeMaître-Tolman-Bondi form

$$ds^2=dt^2+\frac{R_{}^{}{}_{}{}^{2}}{1+f}dr^2+R^2d\mathrm{\Omega }^2,$$
(10)
where $`R`$ satisfies

$$\dot{R}^2=\frac{\mathrm{\Lambda }}{3}R^2\frac{Q}{R^2}+f.$$
(11)
The function $`f=f(r)>1`$ is naturally interpreted as the energy inside a shell labelled by $`r`$ in the dark radiation vaccum and is fixed by its initial configuration.

## 3 Inhomogeneous de Sitter Solutions

Let us assume from now on that $`Q<0`$ and $`\mathrm{\Lambda }>0`$ \[note that according to the recent supernovae data $`\mathrm{\Lambda }10^{84}\text{GeV}^2`$ \- see e.g. \]. Then we have a de Sitter brane with negative dark energy density and the gravitational field is confined to the vicinity of the brane. The marginally bound models corresponding to $`f=0`$ are the static solutions

$$ds^2=\left(1+\frac{Q}{R^2}\frac{\mathrm{\Lambda }}{3}R^2\right)dT^2+\left(1+\frac{Q}{R^2}\frac{\mathrm{\Lambda }}{3}R^2\right)^1dR^2+R^2d\mathrm{\Omega }^2$$
(12)
and were used for the static exterior of a collapsing sphere of homogeneous dark radiation . When $`\mathrm{\Lambda }=0`$ this corresponds to the zero mass limit of the tidal Reissner-Nordström black hole solution on the brane .

It is for $`f0`$ that the dark radiation vaccum really becomes dynamical. The solutions are better organized by a redefinition of $`f`$, the function $`\beta =(3/\mathrm{\Lambda })[3f^2/(4\mathrm{\Lambda })+Q]`$. As an explicit example take $`\beta >0`$. Then $`|f|>2\sqrt{Q\mathrm{\Lambda }/3}`$ and

$$\left|R^2+\frac{3f}{2\mathrm{\Lambda }}\right|=\sqrt{\beta }\mathrm{cosh}\left[\pm 2\sqrt{\frac{\mathrm{\Lambda }}{3}}t+\mathrm{cosh}^1\left(\frac{\left|r^2+\frac{3f}{2\mathrm{\Lambda }}\right|}{\sqrt{\beta }}\right)\right],$$
(13)
where $`+`$ or $``$ refer respectively to expansion or collapse. Since $`R`$ is a non-factorizable function of $`t`$ and $`r`$ these are new non-static, cosmological and intrinsically inhomogeneous exact solutions for negative dark energy dynamics on a spherically symmetric de Sitter brane.

## 4 Singularities and Rebounces

The dark radiation dynamics defined by Eq. (11) can lead to the formation of shell focusing singularities at $`R=0`$ and of regular rebounce points at some $`R0`$. To analyse these issues consider

$$R^2\dot{R}^2=V(R,r)=\frac{\mathrm{\Lambda }}{3}R^4+fR^2Q.$$
(14)
If for all $`R0`$ we find $`V>0`$ then a shell focusing singularity forms at $`R=R_s=0`$. On the other hand if there exists $`R=R_{}0`$ where $`V=0`$ then a regular rebounce point forms at $`R=R_{}`$. For the dark radiation vaccum it is clear that no more than two regular rebounce epochs can be found. Since $`\mathrm{\Lambda }>0`$ there is always a phase of continuous expansion to infinity with ever increasing speed. For $`\mathrm{\Lambda }>0`$ and $`Q<0`$ other phases are possible depending on $`f(r)`$. To ilustrate consider $`\beta >0|f|>2\sqrt{Q\mathrm{\Lambda }/3}`$. If $`f>2\sqrt{Q\mathrm{\Lambda }/3}`$ then $`V`$ is positive for all $`R0`$. There are no rebounce points and the dark radiation shells may either expand continuously or collapse to a shell focusing singularity at $`R_s=0`$.

However for $`1<f<2\sqrt{Q\mathrm{\Lambda }/3}`$ (see Fig. 1) there are two rebounce epochs at $`R=R_\pm `$ with $`R_\pm ^2=3f/(2\mathrm{\Lambda })\pm \sqrt{\beta }`$. Since $`V(0,r)=Q>0`$ a singularity also forms at $`R_s=0`$. The region between the two rebounce points is forbidden because there $`V`$ is negative. The phase space of allowed dynamics is thus divided in two disconnected regions separated by the forbidden interval $`R_{}<R<R_+`$. For $`0RR_{}`$ the dark radiation shells may expand to a maximum radius $`R=R_{}`$, rebounce and then fall to the singularity. If $`RR_+`$ then there is a collapsing phase to the minimum radius $`R=R_+`$ followed by reversal and subsequent accelerated continuous expansion. The singularity at $`R_s=0`$ does not form and so the solutions are globally regular.

## 5 Conclusions

In this work we have reported some new results on the gravitational dynamics of inhomogeneous dark radiation on a RS brane. Taking an effective 4-dimensional approach we have shown that under certain simplifying but natural assumptions the Einstein field equations on the brane form a closed, solvable system. We have presented exact dynamical and inhomogeneous solutions for $`\mathrm{\Lambda }>0`$ and $`Q<0`$ showing they further depend on the energy function $`f(r)`$. We have also described the conditions under which a singularity or a regular bounce develop inside the dark radiation vaccum. Left for future research are for instance the analysis of the confinement of gravity near the brane and the introduction of an inhomogeneous collapsing distribution of matter.

Acknowledgements

We thank the Fundação para a Ciência e a Tecnologia (FCT) for financial support under the contracts SFRH/BPD/7182/2001 and POCTI/32694/FIS/2000. We also would like to thank Louis Witten and T.P. Singh for kind and helpful comments.


text:
# On the nonextensivity of the long range X-Y model

## Abstract

It will be given analytical and numerical evidence supporting that the X-Y model yields an extensive, i.e. proportional to the number of degrees of freedom $`N`$, internal energy $`U`$ for any value of the interaction range.

X-Y model. Long–range interactions. Monte Carlo simulations. Non-extensive systems.

The X-Y is a well known model of Statistical Mechanics. It is defined by a set of conjugate pairs of variables $`(\theta _i,p_i)_{i=1,\mathrm{},N}`$, $`N`$ being the number of degrees of freedom, with a Hamiltonian function

$$=\frac{1}{2}\underset{i=1}{\overset{N}{}}p_i^2+\frac{\kappa }{2}\underset{i,j=1}{\overset{N}{}}\frac{1\mathrm{cos}(\theta _i\theta _j)}{r_{ij}^\alpha }K+V,$$
(1)
where the kinetic, $`K`$, and potential, $`V`$, contributions to the Hamiltonian have been identified. In this expression $`r_{ij}`$ is the distance between the sites $`(i,j)`$ of a regular lattice of dimensionality $`d`$ and lattice spacing equal to one. For example, in the linear case that will be analysed in detail later, $`d=1`$, it is $`r_{ij}=|ij|`$. Periodic boundary conditions are assumed, and it is understood the convention of minimal distance between sites $`i`$ and $`j`$. The constant $`\kappa `$ sets the interaction strength. The definition is such that, for ferromagnetic coupling $`\kappa 0`$ (the only one considered in this paper) the minimum value of the energy is equal to zero.

The parameter $`\alpha `$ sets the interaction range and its variation allows to recover some limits of interest. For example, the case $`\alpha =\mathrm{}`$, for which $`r_{ij}^\alpha =0`$ unless $`r_{ij}=1`$, is the nearest neighbor version, whereas the opposite limit, $`\alpha =0`$, is a mean field Hamiltonian in which all sites interact equally. This mean field limit is usually studiedST71  under the condition that the interaction strength is of order $`N^1`$, i.e. by setting $`\kappa =ϵ/N`$ with $`ϵ=O(1)`$. It is known that in this case the model exhibits a phase transition from an ordered to a disordered phase at a critical value $`\beta _c=2/ϵ`$. This result followsAR95  from an analysis of the partition function

$$𝒵(\beta ,N)=_{\mathrm{}}^{\mathrm{}}𝑑p_1\mathrm{}_{\mathrm{}}^{\mathrm{}}𝑑p_N_0^{2\pi }𝑑\theta _1\mathrm{}_0^{2\pi }𝑑\theta _N\mathrm{e}^\beta $$
(2)
at inverse temperature $`\beta =1/k_BT`$, as usual.

The model has received recently some attentionTA00  for a value of the interaction strength $`\kappa =O(1)`$ as an example of a non–extensive system. It has been arguedTS95  that in this case the internal energy $`U=`$, as well as other thermodynamical potentials, should scale as $`UN^{2\alpha /d}`$ for $`\alpha <d`$ (the so–called nonextensive regime), whereas for $`\alpha >d`$, the extensive regime $`UN`$ is recovered.

In this paper, I will present arguments showing that the internal energy is an extensive function $`U(N,T)=Nu(T)`$ for any value of $`\alpha `$. The result is based upon an exact calculation of the partition function in the mean–field limit $`\alpha =0`$ as well as some numerical simulations.

Let me first sketch the proof for $`\alpha =0`$. Introducing the total “spin vector” $`\stackrel{}{M}=_i\stackrel{}{m}_i`$ with $`\stackrel{}{m}_i=(\mathrm{cos}\theta _i,\mathrm{sin}\theta _i)`$ it is possible to write the potential energy as $`V=\frac{\kappa }{2}(N^2M^2)`$. Using this relation and the Hubbard-Stratonovich transformation, one arrives at the following expression for the partition functionarfor :

$$𝒵=\left(\frac{2\pi }{\beta }\right)^{\frac{N}{2}}\mathrm{e}^{\beta \frac{\kappa }{2}N^2}\frac{(2\pi )^N}{\beta \kappa }_0^{\mathrm{}}𝑑yy\mathrm{e}^{\frac{y^2}{2\beta \kappa }}\left[I_0(y)\right]^N𝒵_K𝒵_V$$
(3)
where $`𝒵_K=\left(\frac{2\pi }{\beta }\right)^{\frac{N}{2}}`$ and $`𝒵_V`$ are, respectively, the kinetic and potential contributions and $`I_n(y)`$ is the modified Bessel function of order $`n`$. The standard calculationAR95  assumes that $`\kappa =O(N)`$ and finds the dominant saddle point contribution to the above integral in the limit $`N\mathrm{}`$. A similar calculation can be carried out assuming that $`\kappa =O(1)`$. In this case, the location of the saddle point $`y_0`$ is given as the solution of

$$\frac{1}{Ny_0}\frac{2y_0}{\beta \kappa N}+\frac{I_1(y_0)}{I_0(y_0)}=0.$$
(4)
This can be solved in the limit $`N\mathrm{}`$, with the result that $`y_0=\frac{N\beta \kappa }{2}\frac{1}{2}+O(N^1)`$. The resulting partition function is:

$$𝒵=\left(\frac{2\pi }{\beta }\right)^{\frac{N}{2}}\frac{1}{\beta \kappa }\left(\frac{2\pi }{\beta \kappa N}\right)^{\frac{N}{2}}$$
(5)
It turns out that the Helmholtz free energy $`F=k_BT\mathrm{ln}𝒵`$ is not an extensive function since it contains terms of order $`N\mathrm{ln}N`$. However, these terms disappear when computing the internal energy $`U=\frac{\mathrm{ln}𝒵}{\beta }`$. This yields $`U=(N+1)/\beta `$. In the thermodynamic limit the dominant term is:

$$U=\frac{N}{\beta }$$
(6)
showing that the internal energy is indeed an extensive function of the number of degrees of freedom. Note that the kinetic and potential terms in the Hamiltonian contribute equally, $`K=V={\displaystyle \frac{N}{2\beta }}`$, to the internal energy.

The nearest neighbor interaction, equivalent to the limit $`\alpha =\mathrm{}`$, has been studied in Ref.LPR87 ; EKL94  for the linear lattice $`d=1`$. It is found that the internal energy is again an extensive function, namely:

$$U=N\left[\frac{1}{2\beta }+\kappa \left(1\frac{I_1(\beta \kappa )}{I_0(\beta \kappa )}\right)\right]$$
(7)
Therefore, since the limits $`\alpha =0`$ and $`\alpha =\mathrm{}`$ both yield an extensive internal energy, it is natural to speculate that a similar behavior will hold for any value of the parameter $`\alpha `$. In order to check this result I have performed Monte Carlo simulations of the one dimensional system given in (1) for the values of $`\alpha =0,0.5,\mathrm{}`$ and a value of the coupling strength $`\kappa =1`$. Whereas the limiting cases $`\alpha =0,\mathrm{}`$ are only for the purposes of checking that the numerical simulations do reproduce the theoretical results, the intermediate value $`\alpha =0.5`$ can not be compared with any a priori theoretical calculation.

I have used a standard Metropolis algorithmBH88  to sample the phase space with a probability proportional to $`\mathrm{e}^\beta `$: (i) Choose randomly a lattice site $`i`$. (ii) Propose a change $`(p_i,\theta _i)(p_i^{},\theta _i^{})`$ of the variables associated to this site; the proposed values for $`p_i^{}`$ and $`\theta _i^{}`$ are selected randomly from a uniform distribution in $`(p_i\delta _p,p_i+\delta _p)`$ and $`(\theta _i\delta _\theta ,\theta _i+\delta _\theta )`$, respectively. (iii) Accept the change with a probability $`\mathrm{min}[1,\mathrm{e}^{\beta \mathrm{\Delta }}]`$, being $`\mathrm{\Delta }=^{}`$ the change in energy involved in the proposed change. The values of $`\delta _p`$ and $`\delta _\theta `$ are chosen such that acceptance probability remains close to 50%. Typically, after $`10^5`$ Monte Carlo steps per particle for equilibration, results have been averaged for $`10^6`$ Monte Carlo steps per particle for each value of the inverse temperature $`\beta `$. In all cases, the following system sizes have been taken: $`N=50,100,200`$.

Figures 1-3 plot the internal energy rescaled by the number of degrees of freedom, $`U/N`$, as a function of the inverse temperature $`\beta `$, for the differer values of $`\alpha `$ considered.. Figure 4 plots the rescaled potential energy, $`V/N`$, in the case $`\alpha =0.5`$. In all cases the data for all values of $`N`$ collapse in a single curve, showing the extensive nature of the internal energy for all the values of $`\alpha `$. Moreover, the cases $`\alpha =0`$ and $`\alpha =\mathrm{}`$, Figs. 1 and 2, follow the theoretical predictions given by Eqs. (6) and (7), respectively. Note that, as shown in Figs. 3, there is not much quantitative difference between the values of the internal energies corresponding to these two limiting cases.

As other authors have foundTA00 , it is also possible to map the results of the X-Y model with $`\kappa =O(1)`$ to the same model with $`\kappa =O(N)`$ if one allows to rescale the temperature $`T`$ by the number of degrees of freedom $`N`$. In this way it is possible to achieve a sort of non–extensive regime in which $`U(N,T)=N^2u(T/N)`$ (for $`\alpha =0`$). The results of this paper, however, concern the true thermodynamic limit, $`N\mathrm{}`$, while keeping $`T=O(1)`$. It has been shown that this limit yields an extensive internal energy and that this extensive behavior can be recovered by simple Monte Carlo simulations. It is interesting to remark that the same extensive behavior will be found in microcanonical simulations. In fact, it is possible to compute the density of states $`g(E)`$ for the potential part of the Hamiltonian as the inverse Laplace transform of the partition function $`𝒵_V`$PA96 . This gives:

$$\frac{\mathrm{log}(g(E))}{N}=\frac{1}{2}\mathrm{log}\left(\frac{E}{N^2}\right)+c$$
(8)
where $`c`$ is a constant of order one. As shown in figure 5, this result is fully consistent with a numerical calculation of the density of statesSTP02  which predicted empirically a form $`g(E)=\mathrm{e}^{N\varphi (E/N^2)}`$. Within this microcanonical ensemble, the internal (potential) energy is found by using $`\beta =\frac{\mathrm{ln}g(E)}{E}`$ or $`E=N/2\beta `$, the same result we obtain before for the average potential energy $`V`$.

The existence of this extensive regime can be considered surprising at first sight, since the interaction terms are very strong. From the mathematical point of view, the extensive regime is due to the continuous nature of the excitations of the ground state of the Hamiltonian (1). A similar analysis for the long range Ising model in which the spins can only vary by an integer amount indeed produces a non-extensive internal energy. In the X-Y model with long range interactions and $`\kappa =O(1)`$ there is no phase transition at any finite value of the temperature in the thermodynamic limit, and the system is always in an ordered situation. In fact, the order parameter is $`M=O(N)`$ for all values of the temperature.

In summary, it has been given analytical and numerical evidence showing that the internal energy derived from the X-Y Hamiltonian is an extensive function of the number of degrees of freedom, $`U(N,T)=Nu(T)`$, for any value of the parameter $`\alpha `$ setting the interaction range. Other thermodynamic potentials, such as the Helmholtz free energy, contain terms of the order $`N\mathrm{log}N`$. Non extensive regimes require to scale the temperature with $`N`$.

Acknowledgments

This work is supported by the Ministerio de Ciencia y Tecnología (Spain) and FEDER, projects BFM2001-0341-C02-01 and BFM2000-1108.