in euv images of the solar corona , coronal loops appear as long graceful arcs of bright plasma that trace magnetic field lines through the atmosphere .
these loops can be preferentially illuminated because localized heating and inefficient cross - field diffusion lead to hot plasma spreading along individual field lines . despite the obvious magnetic nature of these structures
, it has proven challenging to measure through spectroscopic means the magnetic - field strength within the corona .
of course measuring the magnetic - field strength is equivalent to measuring the energy density , and is therefore a key constraint in the modeling of energetic and eruptive phenomena such as flares and cmes .
coronal loops are sometimes observed to vacillate back and forth with a regular frequency .
the identification of these oscillations as resonant mhd kink waves , trapped between the loop s footpoints ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) , launched the field of coronal loop seismology .
coronal seismology promises the opportunity to measure the magnetic - field strength along the loop through inversion of the kink - mode eigenfrequencies ( e.g. , * ? ? ?
however , before such inversions can be performed the ability to construct stable , curved coronal loop models with realistic density and magnetic profiles is needed .
the state of the art in the modeling of static coronal loops is nonlinear force - free field ( nlfff ) models .
typically such models have used vector magnetograph measurements in the photosphere as an observational constraint . a substantial weakness
to such an approach is that the force - free assumption is rather inappropriate within the chromosphere and the low corona @xcite ; thus , the region in which the force - free assumption is valid and the region in which the observational constraint is applied are disjoint .
this leads to a substantial mismatch between the model field and the actual magnetic field @xcite .
more recent work has partially overcome this difficulty by applying additional constraints higher in the corona .
these constraints have taken the form of euv images of bright coronal loops from instruments such as stereo and the _ atmospheric imaging assembly _ ( aia ) .
initially , multiple such euv images were used , each taken from a different vantage point either using simultaneous images from the two stereo spacecraft or using solar rotation to view the presumably static magnetic structure from different angles .
solar stereoscopy was then used to reconstruct the 3-d curve traced by a loop ( see the reviews of * ? ? ?
* ; * ? ? ?
* ) . however , recent attempts to deduce the 3-d shape of a loop from just a single high - resolution euv image ( such as those from aia ) and a coeval photospheric magnetogram have had intriguing success @xcite .    despite the achievements that the nlfff models have made in reconstructing the geometry of the magnetic field , the force - free assumption decouples the thermodynamic variables from the magnetic field . certainly , if the dynamic timescales are significantly shorter than the cooling times , hydrostatic balance along field lines is still valid .
however , deducing the mass density and gas pressure requires either the specification of the temperature by fiat , or the inclusion of an energy equation that models the heating and cooling of the loop ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
here we constrain the mass density within the loop in a different manner .
in essence , we find an equilibrium solution by considering deviations from a force - free state .
this deviation only occurs within the loop and appears as a change in the field strength without a concomitant change in the field s direction .
such a field strength perturbation produces both magnetic buoyancy and a small lorentz force , both of which depend intimately on the shape or geometry of the loop and which oppose each other in equilibrium .
the establishment of this equilibrium requires that the mass along the loop redistributes itself with a timescale shorter than the cooling time .
thus , accounting for the buoyancy of the loop allows one to directly connect the mass density and other thermodynamic variables ( such as the temperature ) to the shape of the loop and the magnetic - field strength within the loop .
we will show that the temperature profile of the loop is not a function that can be freely specified , but instead has a functional form that is a direct consequence of the geometry of the loop .
our goal here is to self - consistently include curvature and buoyancy in the equilibria of coronal loops and to develop models of the temperature , density , and field strength with the geometry of the loop as the primary input . in order to permit analytic solutions we treat coronal loops as slender magnetic fibrils and adopt the thin - flux - tube approximation when deriving the force balance .
the balance of forces is characterized by a magnetic bond number which is the dimensionless ratio of the buoyancy force to the lorentz force . the shape and curvature of the loop is succinctly expressed in terms of the magnetic bond number , which may be a function of position along the loop . in
2 we will derive the equation that describes the balance of forces and from this equation we identify the magnetic bond number . then , assuming that the shape of the loop is provided by observations , we derive the temperature and mass density profiles that are consistent with this shape . in
3 we present a simple equilibrium solution for an embedded fibril which has a uniform magnetic bond number .
we discuss the atmospheric properties that are consistent with a constant magnetic bond number and derive the resulting coronal magnetic field . in
4 we demonstrate another simple solution corresponding to a semi - circular fibril with a uniform radius of curvature .
finally , in  5 we discuss the implications of our findings and summarize the results .
we will model a coronal loop as a curved , magnetic fibril embedded in a larger coronal magnetic structure .
we will further assume that the fibril is compact and thin . by _
thin _ we mean that the radius of the fibril is small compared to any other relevant length scale . even though the thin - tube approximation may not apply to all coronal loops
, it is an appropriate approximation for many .
recent high - spatial - resolution observations of the solar corona in the fe xiii 19.5 nm line by the hi - c instrument have enabled a resolution of about 150 km , sufficient to resolve the cross - section of most coronal loops . using these observations
, @xcite examined brightness cross - sections for 91 loops in the solar corona and found a distribution of radii sharply peaked at 270 km . by examining the pixel - to - pixel brightness fluctuations across loop cross - sections ,
@xcite have argued that coronal loops are unlikely to be structured on a finer unresolved spatial scale .
since the corona s pressure and density scale heights are generally a hundred times larger than this spatial scale and the lengths of loops a thousand times larger , the thin - flux - tube approximation appears to be quite relevant for a substantial fraction of coronal loops .
we also assume that the external corona is magnetically dominated and its magnetic field is a force - free field .
we adopt the notation that quantities evaluated in the external corona have a subscript ` e ' , while those within the fibril lack a subscript .
thus , the external magnetic field is @xmath0 ; whereas , the internal magnetic field of the fibril is just @xmath1 .
one way to envision the fibril is to select a bundle of field lines within the coronal magnetic field and uniformly increase or decrease the field strength within that bundle by a constant factor , @xmath2 .
we call the constant @xmath3 the field - strength deviation and it can be positive or negative depending on whether the fibril is strongly or weakly magnetized compared to its surroundings .
we consider only loops that are small deviations from the force - free state .
thus , the field - strength deviation will be a small quantity , @xmath4 .
the constant @xmath3 was chosen such that it represents the constant of proportionality between the exterior magnetic pressure and the magnetic - pressure contrast ( the difference in the magnetic pressure between the inside and outside of the fibril ) ,    @xmath5    the field - strength deviation @xmath3 must be constant along the fibril .
otherwise the fibril and surrounding corona would not have a common flux surface where they join .
since the internal magnetic field is proportional to the external magnetic field , the internal field is also force free . however , because of the discontinuity in the field strength at the edge of the fibril , there exists a current sheath that surrounds the fibril .
we neglect the spherical geometry of the solar atmosphere and assume that the corona can be treated as a plane - parallel atmosphere with constant gravity @xmath6 .
we employ a cartesian coordinate system , with the @xmath7@xmath8 plane corresponding to the photosphere and the @xmath9 coordinate increasing upwards ( i.e. , @xmath10 ) .
we restrict our attention to coronal loops that are symmetric about the origin and that are confined to the @xmath7@xmath9 plane .
such loops lack torsion .
in addition to the cartesian coordinate system , both within the fibril and within the external corona we will employ the local frenet coordinates for a field line ( illustrated in figure 1 ) . the direction tangent to the magnetic field will be denoted with the unit vector @xmath11 ( thus @xmath12 ) , and the longitudinal coordinate @xmath13 is the pathlength along a field line measured from the photosphere ( @xmath14 corresponds to the footpoint intersecting the photosphere in the region @xmath15 ) .
the curvature vector for the field line is indicated by @xmath16 , and points in the direction of the principle normal @xmath17 with a modulus equal to the reciprocal of the local radius of curvature @xmath18 of the field line .
the direction of the unit vector in the binormal direction will be indicated with @xmath19 and the torsion of the field line with the variable @xmath20 . in the equilibrium considered here ,
the loop itself lacks torsion ,
@xmath21 , and its binormal uniformly points in the @xmath8-direction , @xmath22 .
however , the exterior coronal magnetic field may be 3-d ( with spatial symmetries near the loop ) .
therefore , for completeness and to aid follow - up work we consider the more general case for the moment and specialize only as necessary .
the standard geometrical relations between these coordinate vectors , i.e. , the frenet - serret formulae , are given below for reference , assuming that the position vector of a field line is given by @xmath23 ,    @xmath24    notice the lack of ` e ' subscripts on the frenet vectors .
for a sufficiently thin tube the frenet vectors will be nearly constant across the fibril with the same value as the surrounding corona .
therefore , we avoid appending subscripted labels to the frenet vectors of the external field to simplify notation and to emphasize that the field lines have the same direction and shape inside and immediately outside the fibril .
the forces acting on an isolated , thin , magnetic fibril embedded in a _ field - free _ atmosphere have been previously derived by a variety of authors @xcite through averaging of the mhd momentum equation over the cross - sectional area of the tube .
here we rederive the equilibrium forces in the presence of a _ magnetized _ external atmosphere .
following @xcite , we begin by averaging the mhd momentum equation over the cross - section of the fibril , with the goal of deriving an equation which describes the force per unit length along the fibril . in equilibrium
, there will be a balance of these forces , and this balance will specify the shape of the fibril and its thermodynamic properties .
let @xmath25 be a cross - sectional surface of the fibril at the location @xmath13 that is everywhere perpendicular to the magnetic field ( i.e. , perpendicular to @xmath11 ) and let @xmath26 be a differential area of this surface .
we define the cross - sectional average of a general quantity @xmath27 in the natural way ,    @xmath28    we write the mhd momentum equation using a formulation of the lorentz force that directly acknowledges that the magnetic force is transverse to the field itself ,    @xmath29    where    @xmath30    in equation  ( [ eqn : mhdmomentum ] ) , the gas pressure , magnetic - field strength , and mass density within the loop are @xmath31 , @xmath32 , and @xmath33 , respectively and @xmath34 denotes the lagrangian time derivative
.    we average equation  ( [ eqn : mhdmomentum ] ) over the cross - section @xmath35 , and seek a static solution by setting the acceleration to zero ,    @xmath36 \ , da          + \frac { \overline{b^2 { { \mbox{{\boldmath$k$ } } } } } } { 4\pi } + { { \mbox{{\boldmath$g$ } } } } \bar{\rho }   = 0 \ ; .\ ] ]    the two terms involving gradients can be expressed as contour integrals around @xmath37 , the boundary of @xmath35 , by using the 2-d form of the divergence theorem appropriate for integration over an open surface ,    @xmath38    the unit vector @xmath39 is the outward normal to the tube s bounding surface . the differential @xmath40 is the differential pathlength around @xmath37 . in deriving this equation , we have decomposed the gradient of the gas pressure into longitudinal and transverse components ,    @xmath41    the total pressure must be continuous across the loop s bounding surface .
therefore , on @xmath37    @xmath42    and the pressures that appear within the integrand of the contour integral in equation  ( [ eqn : avemomentum ] ) may be replaced with those from the external fluid ,    @xmath43    buoyancy is the sum of the gravity acting on a body and the net gas - pressure force acting on the outer surface of the body . in this case , we also must consider the effects of the external magnetic pressure and magnetic tension .
we can construct an extension of archimedes principle that is relevant for our problem by performing a similar cross - sectional average of the mhd momentum equation in the exterior fluid .
this requires that we assume that the mathematical form of the external pressures and density can be analytically continued inside the loop .
the transverse component of the resulting equation provides an expression for the net external pressure force ,    @xmath44    where @xmath45 is the component of @xmath6 perpendicular to @xmath11 . since we have assumed that the tube is thin , to lowest order in the radius of the tube , the tangent vector and curvature vector are constant across a cross - section
thus , the previous equation can be rewritten ,    @xmath46    this expression can be used to eliminate the external pressure integrals in equation  ( [ eqn : forcebalance ] ) .
if we once again assume that the tube is thin and the frenet vectors do not vary significantly over the tube s cross - section , we can replace the cross - sectional averages of the densities and pressures with their axial values to obtain    @xmath47 { { { \mbox{{\boldmath$\hat{s}$ } } } } } +           \frac{b^2-b_{{\rm e}}^2}{4\pi } \ , { { \mbox{{\boldmath$k$ } } } } + \left ( \rho - \rho_{{\rm e}}\right ) { { \mbox{{\boldmath$g$}}}}_\perp = 0 \ ; .\ ] ]    where @xmath48 is the tangential component of gravity .
for simplicity we have labeled the axial values without accents or subscripts and the external magnetic field and density are to be evaluated at the location of the magnetic fibril .
note , this equation is very general ; we did not assume that either of the internal or external fields were force free , nor did we assume that the loop is torsionless .
therefore , this equation is valid for a fibril that describes a fully 3-d curve through an equilibrium atmosphere filled with a general 3-d external magnetic field .
the term in square brackets is the longitudinal component and represents hydrostatic balance along field lines .
the remaining terms are transverse and correspond to the magnetic and buoyancy forces . in some ways these two transverse terms
are analogous ; the last term , buoyancy , is a combination of gravity and the net support provided by the external gas pressure , whereas the second term , the magnetic force , is the residual lorentz force that remains once magnetic tension and the net support provided by the external magnetic pressure have been combined .
both forces are normal to the surface of the tube .
our intuition may tell us that buoyancy is aligned with gravity ; but , this is only true for bodies with closed symmetric surfaces , such as a sphere .
further , just as buoyancy can point upwards or downwards depending on whether the fibril is over- or underdense , the net magnetic force can point up or down depending on the relative strength of the magnetic field inside and outside the fibril .
we now invoke the assumption that the fibril is torsionless ( @xmath49 ) and vertically oriented , therefore lacking forces in the binormal direction . under this assumption
we can separate the mean force equation  ( [ eqn : meanforces ] ) into only two components ,    @xmath50    where we have used @xmath51 , which holds in equilibrium when the fibril is confined to the @xmath7@xmath9 plane .
equation  ( [ eqn : equil_s ] ) expresses the balance of forces in the tangential or axial direction @xmath11 , whereas equation  ( [ eqn : equil_k ] ) describes the balance in the transverse direction of the principle normal @xmath17 .
the axial equation is simply hydrostatic balance along magnetic field lines .
the transverse equation constrains the density contrast , @xmath52 , through the balance of magnetic and buoyancy forces .
this balance can be characterized by a magnetic bond number @xcite ,    @xmath53    which is a nondimensional ratio of the buoyancy to the magnetic forces .
the length scale that we have used in this definition , @xmath54 , is the footpoint separation , where @xmath55 is the location of each of the fibril s footpoints in the photosphere .
the magnetic bond number appears in calculations of the rayleigh - taylor instability when one or more of the fluid layers are filled with a horizontal field , providing the critical wavenumber below which instability ensues @xcite .
equation  ( [ eqn : equil_k ] ) expresses a condition on the radius of curvature @xmath18 of the fibril s axis in terms of the magnetic bond number @xmath56 and the geometry of the fibril ( i.e. , the direction of the principle normal relative to gravity @xmath57 ) ,    @xmath58    to proceed we need expressions for the frenet unit vectors , @xmath11 and @xmath17 , the radius of curvature , @xmath18 , and the arclength , @xmath13 .
if @xmath59 is the height of the fibril above the photosphere , given as a function of the horizontal coordinate @xmath7 , and the fibril s axis is traced by the position vector @xmath60 , these quantities can be easily derived from equations  ( [ eqn : s_hat ] ) and ( [ eqn : k_hat ] ) ,    @xmath61    where primes denote differentiation with respect to the photospheric coordinate @xmath7 .
if we insert equations  ( [ eqn : k ] ) and ( [ eqn : r ] ) into equation  ( [ eqn : shapeeqn ] ) , we obtain a nonlinear ode for the height of the fibril @xmath59 ,    @xmath62    a quick examination of this equation reveals that the fibril will be locally concave or convex depending on the sign of the magnetic bond number and that points of inflection correspond to vanishing magnetic bond number .
the magnetic bond number is proportional to the signed curvature .
a stable fibril with a single concave arch will have a negative bond number everywhere , whereas a multi - arched structure will have a magnetic bond number that changes sign .
the temperature , density and gas pressure within the loop can all differ from the surrounding corona .
we define the contrasts is these properties in the following manner :    @xmath63    technically the external variables are all functions of two spatial coordinates  e.g . , @xmath7 and @xmath9 . in all of these definitions
, however , the external variable is evaluated at the location of the fibril .
we indicate this by expressing these variables as a function of the pathlength alone . only the contrasts in the mass density @xmath64 , gas pressure @xmath65 , and magnetic pressure @xmath66 play an active role in the force balance .
the temperature contrast @xmath67 is a dependent variable that can be derived from these active variables post facto .
there are four primary equations that provide all of the inter - relations between these properties of the loop ,    @xmath68    the first of these is a direct consequence of equation  ( [ eqn : pressurecont ] ) .
the second results from transforming equation  ( [ eqn : equil_s ] ) from the pathlength variable @xmath13 to the fibril height @xmath69 . the third relation , equation  ( [ eqn:4e - transverse ] ) , arises from a combination of the definition of the magnetic bond number  ( [ eqn : bond ] ) and the equation of transverse force balance  ( [ eqn : ode ] ) , where the magnetic - pressure contrast has been replaced through use of pressure continuity , equation  ( [ eqn:4e - pcont ] ) .
the last equation is a restatement of equation  ( [ eqn : matchedflare ] ) which ensures that loop and external magnetic field share a common bounding flux surface .
the temperature contrast @xmath67 can be found by using the ideal gas law ( @xmath70 ) ,    @xmath71    where @xmath72 is the external corona s pressure scale height , @xmath73 , and @xmath74 is the scale height of the gas - pressure contrast ,    @xmath75    conveniently , these equations can be solved analytically to obtain the contrast variables as a function of the fibril s shape .
if we combine equations  ( [ eqn:4e - hydrostatics ] ) and ( [ eqn:4e - transverse ] ) , we obtain a differential equation that relates the pressure contrast to the derivative of the fibril s height ,    @xmath76    this equation can be directly integrated to obtain ,    @xmath77 \ ; , \ ] ]    where @xmath78 is a positive integration constant that can be fixed in a variety of ways .
one could apply a photospheric boundary condition or one could use spectroscopic observations to fix the temperature or density contrast at the apex of the loop . strictly speaking , the integration constant is the product @xmath79 .
we have included the field - strength deviation in this product for later convenience and in order to remind the reader that the thermodynamic contrasts arise from a small deviation from the force - free equilibrium .
if we subsequently use the equations of pressure continuity  ( [ eqn:4e - pcont ] ) and hydrostatic balance  ( [ eqn:4e - hydrostatics ] ) , we can derive equations for the mass density and magnetic - pressure contrasts ,    @xmath80 \ ; .\end{aligned}\ ] ]    thus , if we are provided the shape of the loop , we can derive all of the loop s contrast variables .
different values of the integration constant @xmath79 result in different temperature profiles along the loop .
equation  ( [ eqn : deltat ] ) can be rewritten to show how the temperature contrast depends on the geometry of the loop and the constant @xmath79 ,    @xmath81}{p_{{\rm e}}+ 2 \alpha c h_{{\rm e}}z_0^{\prime\prime } }   \ ; .\ ] ]    note that the integration constant @xmath78 was chosen such that it appears only in the combination of @xmath79 in the equations above .
thus , one needs only a single constraint to define the thermodynamic variables .
given the pressure contrast @xmath65 , we can use equations  ( [ eqn:4e - pcont ] ) , ( [ eqn:4e - flarerate ] ) , and ( [ eqn : deltap ] ) to express the external field strength at the location of the fibril in terms of the geometry of the fibril ,    @xmath82 \ ; .\ ] ]    this last equation shows that the integration constant @xmath78 is really just the value of the external magnetic pressure at the apex of the loop ( hence @xmath78 is positive ) .
further , the equation demonstrates that not all external fields are self - consistent with an embedded fibril
. if we select an arbitrary field line within a general magnetic field , equation  ( [ eqn : bext ] ) may not be satisfied for a constant @xmath78 .
this becomes readily apparent when we express the direction of the tangent vector in terms of the angle @xmath83 that it makes with the @xmath7-axis , @xmath84 . with a little basic manipulation ,
we find the following relation between the angle @xmath83 and the height function @xmath69 :    @xmath85    subsequently , one can easily show that equation  ( [ eqn : bext ] ) is equivalent to the statement that the @xmath7-component of the external magnetic field is constant along the length of the fibril ,    @xmath86    of course , not all force - free fields will have this property .
if needed , a self - consistent external magnetic field can be derived by noting that the shape of the fibril and the strength of the magnetic field at the location of the fibril provide a boundary condition for the coronal magnetic field that occupies the bulk of the domain . for simplicity
we assume that the external magnetic field is potential and 2d . therefore , since the field is solenoidal and 2d it can be described by a flux function @xmath87 ,    @xmath88    since the field is potential the flux function must be harmonic ,    @xmath89    therefore , we can find a self - consistent external field solution by solving equation  ( [ eqn : harmonic ] ) under the boundary condition    @xmath90    applied at the location of the fibril with @xmath91 fixed by equation  ( [ eqn : bext ] ) .
this ensures that the external field is both parallel to the fibril and the field strength is given by @xmath91 .
there is an infinity of such solutions , all with different photospheric boundary conditions .
for our first example solution we wish to consider one that is both simple and fully analytic .
therefore , we consider a fibril which has constant magnetic bond number @xmath56 along its length . for constant @xmath56 , the transverse force equation  ( [ eqn:4e - transverse ] ) can be integrated twice to obtain the height of the fibril @xmath59 . similarly , using the result of the first of these integrations , equation  ( [ eqn : sprime ] ) can be integrated to obtain the arclength @xmath92
we choose the three constants of integration such that the footpoints occur at @xmath93 , the fibril is symmetric about @xmath94 , and the footpoint at @xmath95 corresponds to @xmath14 ,    @xmath96    with these conditions , we find the following solution for the geometrical properties of the fibril :    @xmath97    figure  [ fig : constbond ] displays solutions for several different values of @xmath56 .
it is clear that fibril posesses a nonconstant radius of curvature , with strong curvature near the apex and weaker curvature in the legs .
the length of the fibril @xmath98 is obtained by inserting the rightmost footpoint position @xmath99 into equation  ( [ eqn : arclength ] ) ,    @xmath100    two interesting limits of this equation exist . as @xmath101 the length of the loop converges to the footpoint separation @xmath102 .
this arises because the loop becomes straight , flat , and confined to the photosphere . as @xmath103 the loop length diverges logarithmically because the height of the loop grows without bound
this can be seen by recognizing that the loop reaches its apex at its center ( @xmath94 ) , therefore achieving a maximum height of    @xmath104    one can easily see that @xmath105 corresponds to a logarithmic singularity in the height .
clearly the height @xmath59 should be a positive function for the range @xmath106 . with a little thought , from equation  ( [ eqn : z0 ] )
we can see that this is only possible if @xmath107 , once again emphasizing that a fibril with a single concave arch has a negative magnetic bond number .
the contrast variables have a common geometrical factor , @xmath108 , appearing in their functional forms . for the model with uniform magnetic bond number this factor takes on a simple form . using equations  ( [ eqn : z0 ] ) and ( [ eqn : dz0dx ] ) , we find    @xmath109    using this expression and the equations derived in subsection   [ subsec : thermo ] we can solve for all of the contrast variables and the external magnetic pressure in terms of the scale height @xmath74 that depends on the bond number , @xmath110 ,    @xmath111    the density and pressure contrasts , as well as the external magnetic pressure ,
are revealed to be exponential functions of height , all with the same constant scale height @xmath74 .
the density and temperature contrast are illustrated in figure 3 , for an isothermal exterior with a scale height of @xmath112 mm .
the constant @xmath79 has been chosen such that the fractional density contrast at the apex of the loop is @xmath113 .
the solid curves correspond to overdense loops and the dotted curves to underdense loops , with the different colors corresponding to loops with different magnetic bond numbers .
we will see in the discussion section  [ sec : discussion ] that overdense loops have been heated relative to the surrounding corona and underdense loops have been cooled
. therefore , overdense loops are probably more physically relevant .    for an isothermal exterior , where @xmath72 is a constant
, one could choose the magnetic bond number such that @xmath114 .
for such a model , the temperature contrast would be identically zero all along the length of the loop and the plasma parameter @xmath115 would be a constant both inside and outside the loop , while the density contrast remains nonzero .
the condition @xmath116 marks a transition between disparate behavior .
loops with @xmath117 have fractional density and temperature contrasts that decrease with height , whereas loops with @xmath118 have fractional density and temperature contrasts that increase with height .      in the previous subsection we discovered that a fibril with constant magnetic bond number requires that the external magnetic pressure have a constant scale height as we slide along the fibril
this suggests that we seek a solution in which the magnetic pressure is only a function of height and varies exponentially . a simple harmonic solution that satisfies this requirement is a sinusoid in the horizontal @xmath7-direction multiplied by a decaying exponential in height @xmath9 ,    @xmath119    in this solution
, @xmath120 is the horizontal repetition length of the field ( @xmath87 changes sign every @xmath121 ) . with such a flux function
, the magnetic field is as follows :    @xmath122 \ , e^{-\kappa z } \ ; .\ ] ]    since the magnetic field is periodic in the @xmath7-direction , we will only consider the arcade that is symmetric about the origin and exists in the range @xmath123 , where @xmath124 .
the magnetic pressure for such a field depends only on height and falls off exponentially with a scale height of @xmath125 ,    @xmath126    in order to satisfy equation  ( [ eqn : consteps_be ] ) , we must insist that the external field s horizontal wavenumber is proportional to the magnetic bond number and that the photospheric value of the field strength depends on @xmath78 ,    @xmath127    lines of constant @xmath87 in the @xmath7@xmath9 plane denote field lines .
the height above the photosphere , @xmath128 , of a field line with a specified value of the flux function is easily obtained by inverting equation  ( [ eqn : psi ] ) ,    @xmath129    where the footpoints intersect the photosphere at @xmath130 ,    @xmath131    we can clearly see from equation  ( [ eqn : z0 ] ) that this potential field solution for the external corona is consistent with the loop model , equation , if the magnetic bond number is inversely proportional to the repetition length of the external magnetic field , @xmath132 , which is a condition we have already imposed to ensure that the external field strength has the appropriate value along the fibril .    the field lines of this potential field are illustrated in figure  [ fig : constbondfield ] . in the beginning of this section
we discovered that stability requires that the magnetic bond number be bounded @xmath107 . in the context of this specific model
, we can see that the magnetic bond number selects the field line within the coronal field that corresponds to the axis of the fibril .
fibrils with different footpoint separations and magnetic bond numbers are shown as various colored field lines in figure  [ fig : constbondfield ] .
a small value of the magnetic bond number corresponds to a short , flat , low - lying loop confined very near the origin , while a magnetic bond number near the lower limit is a tall , long loop , with foot - points approaching the outer edge of the arcade .
in fact , in the limit @xmath133 , the loop is infinitely tall with vertical legs located at @xmath134 .
therefore , there is a direct mapping between the magnetic bond number and the value of the flux function that corresponds to the loop s axis ,    @xmath135    and one can choose to place the loop along any field line in the coronal model by dialing the magnetic bond number between @xmath136 and 0 .
another simple model to consider is a fibril with a constant radius of curvature .
this model is of course an example of a fibril with a nonuniform magnetic bond number .
however , it is still a single concave arch and therefore the bond number will be negative everywhere . with a constant radius of curvature @xmath18 , equation  ( [ eqn : shapeeqn ] ) simplifies to    @xmath137    where @xmath138 is the height above the photosphere of the fibril s center of curvature . from this expression
we can see that the magnetic bond number diverges when @xmath139 resulting in a divergence of all the contrast variables at the same location .
therefore , for a physically meaningful solution we must insist that the center of curvature lies below the photosphere , i.e. , @xmath140 .
the fibril is therefore a circular arc that is less than a full semi - circle and intersects the photosphere at oblique angles .    using the procedure outlined in subsection   [ subsec : thermo ]
we derive all of the contrast variables and the external magnetic pressure from the magnetic bond number ,    @xmath141    figure  [ fig : tandd_circ ] illustrates the density and temperature contrasts for loops with radii of curvature spanning 25125 mm .
each is embedded in an isothermal corona with a scale height of 75 mm .
the integration constant @xmath79 has been chosen such that the fractional density contrast is @xmath113 at the apex of the loop .
for all of the loops the magnitude of the density contrast decreases with height . for the parameters chosen for the figure , the overdense loops tend to be cold compared to their surroundings ,
although those that are sufficiently large ( @xmath142 ) can have a warm crest .
similarly , underdense loops tend to be hot with the tallest and widest loops possessing a cool crest .
if a change of sign occurs , it will happen at a distance of two scale heights above the center of curvature .
many external field solutions could be found that possess circular field lines ; however , most will not possess the needed functional form for the external field strength along the loop .
consider the field generated by a line current located at the center of curvature .
such a line current generates a field with concentric field lines that have constant field stength as one moves along a line .
we can see that this field fails because equation  ( [ eqn : constr_be ] ) dictates that the field strength must decrease with height along the fibril .
the functional form for @xmath143 given by that equation and the fact that the fibril is a circular arc , suggests that we seek a solution in polar coordinates with the following form ,    @xmath144    where @xmath145 and @xmath146 .
the origin of the polar coordinate system is located at the center of curvature , @xmath18 is the radius of curvature of the fibril , and @xmath147 points in the positive @xmath7-direction .
once again , @xmath148 is the value of the external magnetic field strength at the footpoints .
we can verify that this particular potential field has a circular field line at @xmath149 by noting that the flux function is constant there ( in fact , @xmath150 on the fibril ) .
direct evaluation further verifies that the magnetic pressure when evaluated at the location of the fibril has the requisite functional form , equation  ( [ eqn : constr_be ] ) , as long as    @xmath151    the field lines for this potential field are illustrated in figure  [ fig : constcurvfield ] .
the external field lines are drawn in black and the circular fibril is marked in blue .
note , unlike the constant magnetic bond number model presented in section   [ sec : constantepsilon ] where one could dial the control parameter @xmath56 to change where the fibril appears in the corona , only one field line in figure  [ fig : constcurvfield ] meets the prerequisite of constant curvature .
thus , while self - similar , this external magnetic - field model differs for fibrils with different radii of curvature .
our equilibrium model is predicated on the implicit assumption that force balance is achieved at a much faster time scale than radiative cooling and thermal conduction redistribute heat .
hydrostatic equilibrium should be established on a dynamical time scale in between the free - fall time and the acoustic crossing time .
thus , for a loop that reaches 50 mm above the photosphere and for a pressure scale height of 80 mm , the dynamical time scale lies between 1040 minutes .
transverse force balance is achieved by mass redistribution along the loop and should therefore be established on a similar time scale .
estimates of the radiative cooling and conduction times vary , but for active region loops radiative times can be on the order of 140 hr and conduction times are typically 1 hr @xcite . therefore , we expect that immediately after an impulsive heating event , a coronal loop will quickly re - establish force balance and will stay in force balance as the loop slowly cools . in the following section we discuss the implications of this balance .
the shape of a coronal loop is determined by the force - free equilibrium of the surrounding corona .
given this shape , the mass along the loop must redistribute itself such that the buoyancy forces exactly oppose the lorentz forces acting on the loop .
the resulting balance of forces can be succinctly characterized by the magnetic bond number which is fully specified by the height of the loop above the photosphere as a function of position . in equilibrium ,
the magnetic bond number is proportional to the signed curvature of the loop , and therefore concave curvature results in a negative magnetic bond number and convex curvature in a positive bond number . of course
, inflection points correspond to locations where the bond number vanishes .
therefore , a loop comprised of a single concave arch must have negative magnetic bond number everywhere if in equilibrium .
equilibrium loops that possess multiple arches will have negative bond number in the crests and positive bond number in the dips or troughs .
all loops with the same geometric shape have exactly the same profile of the magnetic bond number .
however , these loops can be achieved in a continuum of different ways , all varying in the distribution of mass and field strength .
this continuum is represented by different values of an integration constant @xmath79 which is a direct measure of the magnetic - pressure contrast at the apex of the loop .
for example , an underdense loop with enhanced magnetic - field strength can have exactly the same bond number as an overdense loop with reduced field strength . while the shape of these two loops will be the same , the
thermodynamic and magnetic properties of these two loop might be very different .      for a loop with enhanced field strength ,
the mass density is proportional to the second derivative of the height function with respect to the horizontal photospheric coordinate , equation  ( [ eqn : deltarho ] ) .
thus , we can immediately deduce that peaks or crests in a loop will have a mass deficit while dips or troughs have mass accumulation .
furthermore , if we define the total mass deviation between two points as the integrated density contrast ,          where @xmath154 is the cross - sectional area of the loop at the apex .
if we change variable from the pathlength to the photospheric @xmath7 coordinate and replace the density contrast using equation  ( [ eqn : deltarho ] ) the integral can be evaluated trivially ,      thus , the net mass deviation between a crest and a neighboring trough is zero and the crest drains mass into the dips . for a loop comprised of a single concave arch , the mass deviation between the two footpoints is negative .
mass must drain out of the loop through the photosphere into the solar interior . in the opposite case ,
a loop with a reduced field strength , the mass density is proportional to the negative of the second derivative .
the peaks in such a loop must have enhanced density so that buoyancy can oppose the upward lorentz force .
this obviously implies that the gas pressure must increase throughout the loop to increase the hydrostatic support .
the additional mass is lifted from the photosphere into the loop .
observations of the temperature along the length of a coronal loops have suggested that some loops might be nearly isothermal @xcite .
for the relatively small density contrasts that have been explored here , the temperature contrast is also relatively small and the temperature is perforce roughly isothermal as long as the external corona is isothermal .
however , we wish to point out that the temperature contrast itself is markedly constant over height for a signficant fraction of the models we have illustrated .
three of the models with constant magnetic bond number , appearing in figure 3 , have a temperature variation of less than a 30% percent from footpoint to apex .
these three correspond to the models with the squattest loops .    given a specific shape for a coronal loop
there exists a family of possible solutions , each characterized by a different value of @xmath79 , the integration constant  see equations  ( [ eqn : deltap])([eqn : deltat_alt ] ) . when illustrating our models , we chose to fix this parameter by specifying the fractional density contrast at the apex of the loop .
we wish to point out , however , that the parameter @xmath79 is really a measure of the total heat that has been input into the coronal loop . the excess heat contained by the loop per unit length
is given by      where @xmath157 is the specific heat at constant volume and @xmath158 is the adiabatic exponent .
we can obtain the total excess heat contained by the loop by integrating @xmath159 over the loop s length .
if we eliminate the gas - pressure contrast by using equation  ( [ eqn : deltap ] ) , we reduce the integral to a constant factor times a positive definite integral that depends only on the geometry of the loop ,      since the constant @xmath78 and the integral appearing on the right hand side of the previous equation are always positive , we can immediately discern two facts : ( 1 ) heated loops correspond to those with a negative field - strength deviation @xmath3 and cooled loops to those with a positive value . thus , heated loops are those that are undermagnetized and overdense and cooled loops are overmagnetized and underdense .
( 2 ) the integration constant @xmath78 is a measure of the magnitude of the heat that has been injected into or extracted from the loop . obviously , heat deposition is more likely to be physically relevant . the family of solutions defined by different values of @xmath78 form a continuum of loop models with different amounts of stored heat .
for example , after a large impulsive heating event , the heat content of the loop jumps but force balance is quickly re - established .
this requires a redistribution of mass ( and heat ) along the loop and the resulting balance is characterized by a nonzero value of @xmath78 .
then as radiative cooling and conduction slowly dissipate the loop s excess heat , the loop responds by moving through a sequence of equilibria while simultaneously maintaining force balance .
this sequence is represented by slow temporal attenuation of the parameter @xmath78 .
seismology of a coronal loop is only sensitive to the distribution of wave speed along the loop ( e.g. , * ? ? ?
if we assume that the observed waves are kink oscillations and that they propagate at the kink speed ,      then , at best , a seismic analysis can provide information about the ratio of the field strength to the density .
it is impossible from only a seismic analysis to measure the field strength profile independently of the density profile .
however , with additional constraints , provided either from observations or from basic physical principles , one might be able to disentangle field - strength and density variations along the loop .
the force - balance model constructed here provides just such a constraint .
the density contrast and magnetic - pressure contrast are related to the geometry of the loop through the magnetic bond number ,      if we assume that the interior magnetic field is nearly the same as the external field ( hence the field - strength deviation is small ) , we can use this previous equation to express the kink speed solely in terms of the density contrast , the exterior density and the loop s geometry ,      if we further assume that the exterior corona is isothermal , @xmath164 , this expression has only four free parameters : the integration constant @xmath78 , the field - strength deviation @xmath3 , the coronal scale height @xmath72 and photospheric value of the exterior density @xmath165 .
thus , even the measurement of only a few mode frequencies that provide different spatial averages of the kink speed should be sufficient to either fix these parameters or to verify a choice made using other information .
we have developed a model of curved coronal loops that self - consistently couples deviations from the force - free state to the thermodynamic properties .
the coupling is accomplished by requiring that the loop is in equilibrium and the lorentz force arising from the deviation from a force - free field is balanced by buoyancy .
this links the curvature of the loop to the loop s mass - density contrast and hence to the pressure and temperature through hydrostatic balance and the ideal gas law .
the end result is that specification of the loop s geometry is sufficient to derive the temperature and density profiles to within a single integration constant .
this integration constant can be selected in a variety of ways , either seismically through estimates of the mean kink - wave speed or spectroscopically through estimates of the density or temperature contrast at one position along the loop .
further , this integration constant is a direct assessment of the thermal history of the loop providing an estimate of the excess heat contained by the loop relative to the surrounding corona .