Free Access
Issue
A&A
Volume 498, Number 2, May I 2009
Page(s) 551 - 566
Section The Sun
DOI https://doi.org/10.1051/0004-6361/200810971
Published online 05 March 2009

Causes and consequences of magnetic cloud expansion

P. Démoulin1 - S. Dasso2,3

1 - Observatoire de Paris, LESIA, UMR 8109 CNRS, 92195 Meudon Principal Cedex, France
2 - Instituto de Astronomía y Física del Espacio, CONICET-UBA, CC. 67, Suc. 28, 1428 Buenos Aires, Argentina
3 - Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, 1428 Buenos Aires, Argentina

Received 17 September 2008 / Accepted 23 February 2009

Abstract
Context. A magnetic cloud (MC) is a magnetic flux rope in the solar wind (SW), which, at 1 AU, is observed $\sim$2-5 days after its expulsion from the Sun. The associated solar eruption is observed as a coronal mass ejection (CME).
Aims. Both the in situ observations of plasma velocity distribution and the increase in their size with solar distance demonstrate that MCs are strongly expanding structures. The aim of this work is to find the main causes of this expansion and to derive a model to explain the plasma velocity profiles typically observed inside MCs.
Methods. We model the flux rope evolution as a series of force-free field states with two extreme limits: (a) ideal magneto-hydrodynamics (MHD) and (b) minimization of the magnetic energy with conserved magnetic helicity. We consider cylindrical flux ropes to reduce the problem to the integration of ordinary differential equations. This allows us to explore a wide variety of magnetic fields at a broad range of distances to the Sun.
Results. We demonstrate that the rapid decrease in the total SW pressure with solar distance is the main driver of the flux-rope radial expansion. Other effects, such as the internal over-pressure, the radial distribution, and the amount of twist within the flux rope have a much weaker influence on the expansion. We demonstrate that any force-free flux rope will have a self-similar expansion if its total boundary pressure evolves as the inverse of its length to the fourth power. With the total pressure gradient observed in the SW, the radial expansion of flux ropes is close to self-similar with a nearly linear radial velocity profile across the flux rope, as observed. Moreover, we show that the expansion rate is proportional to the radius and to the global velocity away from the Sun.
Conclusions. The simple and universal law found for the radial expansion of flux ropes in the SW predicts the typical size, magnetic structure, and radial velocity of MCs at various solar distances.

Key words: Sun: coronal mass ejections (CMEs) - Sun: magnetic fields - interplanetary medium

1 Introduction

Interplanetary coronal mass ejections (ICMEs) are the manifestation of magnetized plasma structures ejected from the Sun as CMEs (e.g., Zurbuchen & Richardson 2006; Wimmer-Schweingruber et al. 2006). ICMEs are dynamic structures that move through and interact with the solar wind (SW). The size of an ICME drastically increases, typically by more than a 103 factor, during its travel from the Sun to the external heliosphere.

Magnetic clouds (MCs) are a particular subset of ICMEs that have a well-structured magnetic field. They are formed by twisted magnetic flux tubes, called flux ropes (e.g., Burlaga 1995). Though there is still no consensus, non-MC ICMEs are probably formed by magnetic structures similar to MCs. It is plausible that MCs are only observed when the spacecraft crosses the magnetic structure close to its center (Jian et al. 2006). The MC/ICME ratio increases from ${\approx} 15\%$ at solar maximum to ${\approx} 100\%$ at solar minimum (Cane & Richardson 2003). This evolution is most likely due to an observational selection effect since at solar maximum CMEs are launched from higher solar latitudes than at minimum, so a spacecraft located in the ecliptic more frequently crosses the periphery of the magnetic structure than at solar minimum (Richardson & Cane 2004). A secondary effect is that there are more interactions between ICMEs at solar maximum, so more distortion/reconnection of the magnetic structures is expected at solar maximum than minimum.

The in situ measurements are limited to the spacecraft trajectory crossing the incoming ICME. Therefore, one needs to rely on modeling to derive the global magnetic structure from local measurements. This has restricted the analysis so far to MCs. Most of the magnetic models neglect the evolution of the magnetic field during the crossing time of the observed MC, and consider a flux rope locally invariant along its axis. These different approaches are reviewed by Riley et al. (2004), Dasso et al. (2005), Forbes et al. (2006), and briefly summarized below.

The magnetic field in MCs can be relatively well modeled by the so-called Lundquist model (Lundquist 1950), which considers a static and axially-symmetric linear force-free magnetic configuration (e.g., Burlaga 1988; Lepping et al. 1990; Lynch et al. 2003). Alternately, Farrugia et al. (1999) considered a cylindrical shape for the cloud cross-section and a non-linear force-free field, while Mulligan et al. (1999), Hidalgo et al. (2002), and Cid et al. (2002) have considered a cylindrical cloud with a finite plasma pressure. Elliptical models have also been developed with a linear force-free field (Vandas & Romashets 2003) and with a non-force free field (Hidalgo 2003). They are able to describe distortions of the magnetic structures due to the interaction with its environment.

The magnetic structure of MCs has also been analyzed by solving the magneto-static force balance as a Cauchy problem (e.g., Hu et al. 2005; Hu & Sonnerup 2002). The amount of distortion from circular cross section is variable in the MCs analyzed. Some MCs are significantly distorted (within the limit that a Cauchy problem is ill-posed, since it is sensitive to modifications of the boundary conditions, so that the results can be significantly affected by the temporal resolution and the range of the data used).

Inside MCs, and more generally inside ICMEs, the measured plasma velocity typically has a linear variation along the spacecraft trajectory. A much higher velocity is present in the front than in the back, indicating expansion. Burlaga & Behannon (1982) found consistency between the expansion speed estimated from in situ observations and the increase of their typical size, obtained from measurements with different spacecraft located between 2 and 4 AUs.

Self-similar models have been used to describe the expansion of flux ropes. Some consider only a radial expansion (e.g., Farrugia et al. 1993; Osherovich et al. 1993; Farrugia et al. 1997; Nakwacki et al. 2008), while others consider an isotropic expansion (e.g., Shimazu & Vandas 2002; Berdichevsky et al. 2003; Dasso et al. 2007). The main problem of the first group of models is that, even taking a force-free configuration at some point of the evolution, the only radial expansion implies the creation of strong unbalanced magnetic forces during the evolution. Then, the force-free state is singular in these models, happening only at one single time during the evolution, while observations show that MCs typically have a low-$\beta $ plasma. A self-similar isotropic expansion does not change the force balance, so the second group of models are expected to be a better approximation. However, it is not obvious that the expansion could be isotropic since the magnetic force is anisotropic in a flux rope.

Flux rope evolution has also been investigated with MHD simulations (Manchester et al. 2004; Chané et al. 2006; Cargill et al. 2000; Wu et al. 2007; Shen et al. 2007; Cargill & Schmidt 2002; Riley et al. 2003). A flux rope model is typically inserted in a SW model close to the Sun as an initial condition. Typically, these MHD simulations found that a flux rope, moving much faster than the surrounding SW, has its cross section flattened as it moves outward. Flattening is interpreted as two effects, one geometrical (spherical expansion) and the second due mainly to the dynamic pressure of the overtaken SW (Riley et al. 2003). The magnetic tension of the flux rope is not sufficient in these simulations to maintain an approximately circular cross section.

In present work we analyze the expansion of flux ropes in the SW with a complementary approach to the self-similar expansion models and to the MHD simulations cited above. Here we model the dynamical evolution of the flux rope as a moving boundary problem (see, e.g., Crank 1984), where the ``Stefan condition'' is provided by the pressure balance between the ambient solar wind and the MC boundary. Our approach is to keep only the main physical effects in order to derive, as well as we can, analytical results showing the effects of the main physical parameters. We assume a force-free evolution, and we will test how well the evolution could be considered as self-similar. The hypothesis and the main equations are derived in Sect. 2, with the scope to define the key variables of the evolution. The equations are solved in Sect. 3 for three very different types of flux rope. In Sect. 4, we use the known properties of the SW and our theoretical results to derive the expected and generic expansion of flux ropes in the SW. Then, we compare these results to the observed expansion of MCs and ICMEs. In Sect. 5 we analyze the evolution of flux ropes with the hypothesis of minimizing their magnetic energy with a preserved magnetic helicity; we show that under these dissipative conditions, flux ropes expand almost at the same rate as in ideal MHD. This further extends our main result, showing that the expansion of flux rope in the SW is mainly due to the decrease of the total SW pressure with solar distance. Finally, we summarize our results and conclude in Sect. 6.

 \begin{figure}
\IfFileExists{color_figures.txt}{
{\includegraphics[width=0.5\tex...
...}{
{\includegraphics[width=0.5\textwidth, clip=]{0971f1c.eps} }
}\end{figure} Figure 1:

Schematic of a flux rope ejected from the Sun (the apex of the central field line is at a distance D from the Sun). The magnetic configuration of a magnetic cloud (MC) is well represented by an elongated twisted flux tube mostly anchored to the Sun. The schematic on the left shows the locally straight approximation used in this work (A color version is available in the electronic version).

Open with DEXTER

2 Equations for a force-free expanding flux rope

In this section we present a new model to describe the evolution of an expanding flux rope in the solar wind, deriving the complete set of its equations. Since MCs typically have a low-$\beta $ plasma, we expected that their magnetic field is near to a force-free state. Below, the main hypothesis of the model is that the flux rope stays in a force-free state and conserves elementary magnetic fluxes while it evolves (i.e. we suppose an ideal-MHD evolution). The aim is to keep only the most important physical ingredients needed to understand the typical expansion rate of MCs as derived from in situ observations.

2.1 Ideal MHD and force-free field evolution

This section describes the internal evolution of the flux rope. Its interaction with the surrounding SW is described in the next section.

The magnetic configuration of a MC is typically described by a flux rope ejected from the Sun (Fig. 1). When the flux rope is moving faster than the surrounding SW, plasma and magnetic field typically accumulate in front (in the so-called sheath region). With this forcing, magnetic reconnection is expected, and indeed evidence of reconnection occurrences have been found (Gosling et al. 2005; Dasso et al. 2007). The amount of reconnected flux is highly variable in MCs, depending at least on the SW field structures overtaken, the relative orientation of the SW and MC leading fields, and the difference between the SW and MC velocities. Below we neglect this interaction that will basically ``peal'' the flux rope so that a progressively less extended flux rope is moving in the SW. The surviving flux rope is indeed expected to expand following roughly the same law as the original one since the main driver of the expansion, the decrease of SW pressure with solar distance (see below), is basically the same.

The main driver of the flux rope expansion is the quasi-balance between forces inside and outside the rope during its travel in the solar wind. Below, we retain only the Lorentz force inside the flux rope, since the plasma $\beta $ is low in MCs (typically $\beta \approx 0.1$, with values ranging from less than ${\approx} 10^{-2}$ to a few times 0.1 (e.g., Burlaga & Behannon 1982; Wu & Lepping 2007; Feng et al. 2007; Lepping et al. 2003), so the plasma pressure has a weak contribution. Other forces such as gravity are also negligible compared to the magnetic pressure gradient. In the frame of the MC, the plasma velocity is typically smaller than the Alfvén velocity (a few 100 km s-1, Burlaga & Behannon 1982), then the magnetic evolution can be described, to first approximation, by a series of force-free equilibrium ( $\vec{j} \times \vec{B} \approx \vec{0}$).

A MC typically has an elongated flux rope structure with a cross-section size much thinner than the curvature radius of its axis, so we approximate locally the magnetic field by a cylindrical field with a local invariance along its axis. We also simplify the problem by considering a circular cross section (Fig. 1). With these simplifications $\vec{B}$ is described only by its axial, Bz(r), and azimuthal, $B_{\theta}(r)$, components; they are functions of the local radius r. The internal force balance is written as

\begin{displaymath}
\frac{ {\rm d} (B_z^2+B_{\theta}^2)/2}{ {\rm d} r} + \frac{B_{\theta}^2}{r} = 0 .
\end{displaymath} (1)

The flux rope evolution should then follow approximately the above force-free equilibrium during its evolution in the SW. We define the field  $\vec{ B_{\rm o} }$ (with components $ B_{z, {\rm o}} =B_z( r_{\rm o} )$ and $ B_{\theta, {\rm o}} =B_{\theta}( r_{\rm o} )$) at a time $t= t_{\rm o} $ during the evolution. This reference field has no special requirements except to be a known field that satisfies Eq. (1), and we derive below the field at any time t as a function of the field at $ t_{\rm o} $.

In this study we suppose an ideal MHD evolution (except in Sect. 5), so at all times there is the following elementary magnetic flux conservation


 
    $\displaystyle B_{\theta} ~ {\rm d} z ~ {\rm d} r = B_{\theta, {\rm o}} ~ {\rm d} z_{\rm o} ~ {\rm d} r_{\rm o}$ (2)
    $\displaystyle B_z ~r ~ {\rm d} r = B_{z, {\rm o}} ~ r_{\rm o} ~ {\rm d} r_{\rm o}.$ (3)

These equations relate $\vec{B}$ to $\vec{ B_{\rm o} }$ through the implicit expansion function  $r( r_{\rm o} )$.

2.2 Drivers of the evolution

In this section we describe the model for the evolution of the flux rope: first, its axial expansion, then its interaction with the surrounding SW. The apex of the flux rope axis is located at the variable distance D from the Sun, and it has a length L that increases with D (Fig. 1). We suppose that the axial expansion is uniform along and across the flux tube, so the local expansion rate $ {\rm d} z / {\rm d} z_{\rm o}$ is simply provided by the relative extension of the flux rope axis

\begin{displaymath}
\tilde{L} (D) \equiv {\rm d} z / {\rm d} z_{\rm o} = L / L_{\rm o} ,
\end{displaymath} (4)

where $L_{\rm o}$ is the flux rope length at $t= t_{\rm o} $. $ \tilde{L} $ is non-dimensional, and below we add a ~ on top of all variables that are normalized at the reference state, see Table 1 for the main variables. The evolution of $ \tilde{L} $ induces an evolution of  $B_{\theta }$ (via Eq. (2)) and also of Bz (via Eq. (1)).

Table 1:   Definition of the main non-dimensional variables.

A second driver of the magnetic field evolution is the force balance with the surrounding SW. As the total SW pressure decreases with solar distance (D) the flux rope expands and its internal pressure decreases accordingly, progressively reaching new states of quasi-equilibrium. This balance of pressure, at the flux rope boundary r=R (flux rope radius), is simply written as

\begin{displaymath}
\frac{B_{\theta}^2(R) +B_z^2(R)}{2 \mu_{\rm o}} = P_{\rm t} (D) ,
\end{displaymath} (5)

where $ P_{\rm t} $ is the total SW averaged pressure around the MC.

However, when a flux rope moves faster than the SW, a dynamic pressure, typically on the order of $\rho_{\rm SW} (v-v_{\rm SW})^2/2$, is present in front of the flux rope (where $\rho_{\rm SW}$ is the SW density and $v-v_{\rm SW}$ is the relative velocity of the MC with respect to the SW). This dynamic pressure has a distribution around the flux rope that depends on the Mach number of the relative flow in pure hydrodynamic, and it is even more complex with a magnetized SW. Moreover, when the MC is moving significantly faster than the SW, a MHD shock is present in front of the flux rope. For a small fraction of MCs, typically those traveling in the fast solar wind with a comparable mean speed, both forward and reverse shocks are observed (e.g. Gosling et al. 1998). These shocks are associated with a jump in the total pressure. Each of the above complex interactions with the SW are not isotropic around the flux rope and so they contribute to its deformation as present in MHD simulations (Sect. 1).

We suppose below that the relative velocity of the flux rope with respect to the SW is sufficiently small, or equivalently that the magnetic tension is strong enough to neglect the deformation of the flux rope cross section, as well as the anisotropy of the total SW pressure around the flux rope. Then, at the boundary of the flux rope, defined by r=R, we write the pressure equilibrium as

\begin{displaymath}
\frac{B_{\theta}^2(R) +B_z^2(R)}{2 \mu_{\rm o}} = \varpi (D) ~ P_{\rm t} (D) ,
\end{displaymath} (6)

where $\varpi (D)$ is a factor that includes the internal over-pressure, the isotropic part of the dynamic pressure, and shock pressure discontinuity.

While $ P_{\rm t} (D)$ can be imposed as an external background pressure, $\varpi (D)$ can be written as an explicit function of D only a posteriori, when the flux rope evolution is solved. However, the non-dimensional factor  $\varpi (D)$ is typically in the interval between unity and several units (being larger for a stronger shock). Moreover  $\varpi (D)$ is generally not expected to change drastically with D since it would require an important modification of the flux rope velocity, or of the SW properties at some distance D. In contrast, $ P_{\rm t} (D)$ has a strong and systematic variation, typically ${\approx} D^{-3} $ (see Sect. 4.2). For example, as the flux rope moves from 0.1 to 1 AU, $ P_{\rm t} $ decreases by a factor ${\approx} 10^{3}$ well above any plausible variation of  $\varpi (D)$.

We define the non-dimensional function $ \tilde{ P_{\rm t} } $ using the reference state (field $\vec{ B_{\rm o} }$ when the flux rope is at the distance $ D_{\rm o} $)

\begin{displaymath}
\tilde{ P_{\rm t} } (D) \equiv \frac{\varpi(D)~ P_{\rm t} (...
..._{\rm o} )} = \frac{B^2(R)}{ B_{\rm o} ^2( R_{\rm o} )} \cdot
\end{displaymath} (7)

Formally, $ \tilde{ P_{\rm t} } $ incorporates the contribution of $\varpi (D)$, so $ \tilde{ P_{\rm t} } $ includes the over-pressure present in the MC with respect to its surrounding SW. Then, this over-pressure could be incorporated in the following results if $\varpi (D)$ is derived in some cases. In practice, this effect is only expected to be a correction to  $ P_{\rm t} (D)$, and the uncertainties on $ P_{\rm t} (D)$ derived presently by the observations are probably larger than the contribution of $\varpi (D)$ itself.

2.3 Main equations of the evolution

In this section we write the above equations in a minimal format, decreasing both the number of equations and variables from five (Eqs. (1), (2), (3), (4), (7)) and $B_{\theta}, B_z,B,r, {\rm d} z / {\rm d} z_{\rm o}$) to two (Eq. (10), (11)) and B,r).

Since the magnetic field strength is present both in the force-free equation and the boundary equation (Eqs. (1), (7)), we keep B as one of the main unknown function. From Eqs. ((2)-(4)), we derive $B_{\theta} = B_{\theta, {\rm o}} ( {\rm d} r_{\rm o} / {\rm d} r) / \tilde{L} $ and $B_z = B_{z, {\rm o}} ( r_{\rm o} /r)~( {\rm d} r_{\rm o} / {\rm d} r) $. Then, both $B_{\theta }$ and Bz can be written as a function of $B=\sqrt{B_{\theta}^2+B_z^2}$, r and the reference equilibrium by dividing previous equations by B. Thus, for example,

\begin{displaymath}
B_{\theta}(r) = B(r) \frac{r~ B_{\theta, {\rm o}} }{\sqrt{(...
...rm o}} )^2 + ( \tilde{L} ~ r_{\rm o} ~ B_{z, {\rm o}} )^2}} ,
\end{displaymath} (8)

where $ B_{\theta, {\rm o}} $ and $ B_{z, {\rm o}} $ are functions of $ r_{\rm o} $. Then, Eq. (1) is rewritten as

\begin{displaymath}
\frac{ {\rm d} B}{ {\rm d} r} = - r~B~ \frac{ B_{\theta, {\...
...o}} )^2 + ( \tilde{L} ~ r_{\rm o} ~ B_{z, {\rm o}} )^2} \cdot
\end{displaymath} (9)

The evolution of r is constrained by magnetic flux conservation. With Eq. (8) included in Eq. (2), we have

\begin{displaymath}
\frac{ {\rm d} r}{ {\rm d} r_{\rm o} } =
\frac{\sqrt{(r~ ...
...{L} ~ r_{\rm o} ~ B_{z, {\rm o}} )^2}}{ \tilde{L} ~r~B} \cdot
\end{displaymath} (10)

Since we suppose an ideal MHD evolution, $ r_{\rm o} $ labels a given plasma blob during the evolution (it is a Lagrangian coordinate), and  $r( r_{\rm o} )$ describes the radial evolution of this plasma blob. Then, we use $ r_{\rm o} $ as the basic variable to describe all quantities. An equation for  $B( r_{\rm o} )$ is simply obtained by multiplying Eqs. (9) by (10)

\begin{displaymath}
\frac{ {\rm d} B}{ {\rm d} r_{\rm o} } = - \frac{ B_{\theta...
...}} )^2 + ( \tilde{L} ~ r_{\rm o} ~ B_{z, {\rm o}} )^2}} \cdot
\end{displaymath} (11)

The boundary conditions of the system of Eqs. ((10), (11)) are
r = 0 $\textstyle {\rm ~~~~at}$ $\displaystyle r_{\rm o} = 0$  
$\displaystyle B( r_{\rm o} ) = \sqrt{ \tilde{ P_{\rm t} } } ~ B_{\rm o} ( R_{\rm o} )$ $\textstyle {\rm ~~~~at}$ $\displaystyle r_{\rm o} = R_{\rm o}.$              (12)

The force-free field evolution is described by a boundary-value problem coupling two first order differential Eqs. ((10), (11)); its solution provides  $r( r_{\rm o} )$ and  $B( r_{\rm o} )$, from which we can derive the magnetic profile B(r) and the field components (e.g. Eq. (8)). The magnetic field evolution is driven by the change of both the axial extension ($ \tilde{L} $) and the total pressure at the boundary ( $ \tilde{ P_{\rm t} } $).

2.4 Reformulation of the equations

Since $ B_{\theta, {\rm o}} ( r_{\rm o} =0)=0$ for any non-singular magnetic field, the system of Eqs. ((10), (11)) presents a weak singularity at $r= r_{\rm o} =0$ (in both equations the numerator and denominator vanish). We define the function  $ b_{\theta, {\rm o}} $ as

\begin{displaymath}
b_{\theta, {\rm o}} ( r_{\rm o} ) \equiv B_{\theta, {\rm o}} ( r_{\rm o} ) / r_{\rm o} .
\end{displaymath} (13)

Equations ((10), (11)) rewrite as the non-singular equations
 
$\displaystyle \frac{ {\rm d} r^2}{ {\rm d} r_{\rm o} ^2} =
\frac{\sqrt{(r~ b_{\theta, {\rm o}} )^2 + ( \tilde{L} ~ B_{z, {\rm o}} )^2}}{ \tilde{L} ~B}
,$     (14)
$\displaystyle \frac{ {\rm d} B}{ {\rm d} r_{\rm o} ^2} =
- \frac{ b_{\theta, {\...
...lde{L} \sqrt{(r~ b_{\theta, {\rm o}} )^2 + ( \tilde{L} ~ B_{z, {\rm o}} )^2}}
.$     (15)

The field strength B is a decreasing function of  $ r_{\rm o} ^2$, and since it should be finite at the boundary (Eq. (12)), B cannot vanish. The same argument applies to $ B_{\rm o} $, so $ B_{z, {\rm o}} $ and $ r_{\rm o} ~ b_{\theta, {\rm o}} $ cannot vanish together at the same value of $ r_{\rm o} $. Since r=0 only for $ r_{\rm o} =0$, $ B_{z, {\rm o}} $ and $r~ b_{\theta, {\rm o}} $ cannot also vanish together. Then, both denominators in Eqs. ((14), (15)) are well behaved for all the integration range of  $ r_{\rm o} ^2$ ( $0, R_{\rm o} ^2$).

The evolution of B is driven by the modification of both the flux rope length, $ \tilde{L} $, and the total pressure $ \tilde{ P_{\rm t} } $. $ \tilde{L} $ is present only in Eqs. ((14), (15)) and $ \tilde{ P_{\rm t} } $ is present only in the boundary conditions (Eq. (12)). In order to group them in the equations, we define a normalized field, $ \tilde{B} $, that has a fixed value at the boundary, r=R during the evolution (this fixed value is $ B_{\rm o} ( R_{\rm o} )/ B_{\rm o} (0)$). We also define a normalized radius, $ \tilde{r} $, guided by the flux conservation (in order of magnitude). Then, we normalize both quantities to the reference field,

\begin{displaymath}\fbox{$\begin{array}{rcl}
\tilde{B} & \equiv & ( \tilde{ P_{...
...\tilde{ P_{\rm t} } )^{ 1/4} ~ r/ R_{\rm o}
,
\end{array} $ } \end{displaymath} (16)

where $B_{\rm c,o} = B_{\rm o} ( r_{\rm o} =0)$ is the field magnitude at the center of the reference field.

With the normalization (16), Eqs. ((14), (15)) are rewritten as

 
$\displaystyle \frac{ {\rm d} \tilde{r} ^2}{ {\rm d} \tilde{r}_{\rm o} ^2} =
\fr...
...}} )^2 + ( \tilde{s} ~ \tilde{B}_{z, {\rm o}} )^2} }{ \tilde{s} ~ \tilde{B} }
,$     (17)
$\displaystyle \frac{ {\rm d} \tilde{B} }{ {\rm d} \tilde{r}_{\rm o} ^2} =
- \fr...
...\tilde{b}_{\theta, {\rm o}} )^2 + ( \tilde{s} ~ \tilde{B}_{z, {\rm o}} )^2} } ,$     (18)

with $ \tilde{r}_{\rm o} = r_{\rm o} / R_{\rm o} $, $ \tilde{b}_{\theta, {\rm o}} = b_{\theta, {\rm o}} ~ R_{\rm o} /B_{\rm c,o} $, $ \tilde{B}_{z, {\rm o}} = B_{z, {\rm o}} /B_{\rm c,o} $, and with the size factor, $ \tilde{s} $, defined by

\begin{displaymath}
\fbox{$\displaystyle \tilde{s} = \tilde{L} ~( \tilde{ P_{\rm t} } )^{1/4} $\space .} \end{displaymath} (19)

The boundary conditions are


$\displaystyle \tilde{r} ^2 = 0$ $\textstyle {\rm ~~at~~}$ $\displaystyle \tilde{r}_{\rm o} ^2 = 0$  
$\displaystyle \tilde{B} = B_{\rm o}( R_{\rm o} )/B_{\rm c,o}$ $\textstyle {\rm ~~at~~}$ $\displaystyle \tilde{r}_{\rm o} ^2 = 1.$            (20)

The system of Eqs. ((17)-(20)) provides the evolution of $ \tilde{B} $ and $ \tilde{r} $ (and thus the size and the internal distribution of the magnetic field in the MC) while the flux rope travels in the solar wind. This system is numerically solved with a shooting method (Press et al. 1992, chap. 17) starting the integration at $ \tilde{r}_{\rm o} ^2=0$ from $ \tilde{r} ^2=0$ and $ \tilde{B} = \tilde{B}_{\rm c} $. An iteration on $ \tilde{B}_{\rm c} $ is realized so that $ \tilde{B} ( \tilde{r}_{\rm o} ^2=1)$ converges to the boundary conditions (Eq. (20)).

The normalization defined by Eq. (16) is also applied to the field components to define $ \tilde{B} _{\theta}, \tilde{B} _z$. Then, Eq. (8) is written as

\begin{displaymath}
\tilde{B} _{\theta} = \tilde{B}
\frac{ \tilde{r} ~ \tilde...
..., {\rm o}} )^2 + ( \tilde{s} ~ \tilde{B}_{z, {\rm o}} )^2}} ,
\end{displaymath} (21)

and we have a similar equation for $ \tilde{B} _z$ $\left(=\sqrt{ \tilde{B} ^2 - \tilde{B} _{\theta}^2}\right)$.

2.5 General properties of the evolution

In the system of Eqs. ((17)-(20)), for a given reference field $( \tilde{B}_{\theta, {\rm o}} , \tilde{B}_{z, {\rm o}} )$, $ \tilde{r} $ and $ \tilde{B} $ are only functions of $ \tilde{s} $ (so only of a combination of $ \tilde{L} $ and  $ \tilde{ P_{\rm t} } $). We conclude that the normalized field $ \tilde{B} ( \tilde{r} )$, as well as its components (Eq. (21)), depends only on $ \tilde{s} $. The explicit dependance of the field B(r) on $ \tilde{ P_{\rm t} } $ is only present as a multiplicator factor: the field B(r) is deduced from $ \tilde{B} ( \tilde{r} )$ by multiplying the normalized field strength by $( \tilde{ P_{\rm t} } )^{1/2} ~ B_{\rm c,o} $ and the normalized radial distance by $( \tilde{ P_{\rm t} } )^{ -1/4} ~ R_{\rm o} $.

A typical reference magnetic field has an axial field dominating in the core of the flux rope, and an azimuthal field dominating at the periphery (recalling that a well behave field has $ b_{\theta, {\rm o}} (0)=0$). As $ \tilde{s} $ increases to large values, Eq. (18) implies a smaller gradient of $ \tilde{B} $ across the flux tube. This corresponds to a magnetic field having a smaller magnetic tension, so more dominated by its axial component (see Eq. (1)). So, as $ \tilde{s} $ increases, the core expands more than its periphery, and the azimuthal field is compressed close to the border of the flux rope (this is illustrated later for specific fields in Fig. 5). On the other hand, as $ \tilde{s} $ decreases to small values, the field strength has an increasing gradient provided by an increasing magnetic tension of the azimuthal field. So as $ \tilde{s} $ decreases the magnetic field is increasing dominated by its azimuthal component in a larger fraction of the flux rope periphery.

In the particular case where the evolution with heliodistance of the length and pressure satisfy $ \tilde{s} =1$, Eqs. ((17)-(20)) are satisfied with $ \tilde{r} = \tilde{r}_{\rm o} $ and $ \tilde{B} = \tilde{B}_{\rm o} $. In this case the only evolution is a re-scaling of the radius by L and of the field strength by L-2. So, when $ \tilde{ P_{\rm t} } =1/ \tilde{L} ^4$, the evolution is simply an isotropic self-similar evolution regardless of the reference field.

For MCs in the SW, $ \tilde{L} $ increases approximatively as the distance, $ \tilde{D} $, while $ \tilde{ P_{\rm t} } $ decreases approximatively as $ \tilde{D} ^{-n_P}$ (see Sect. 4), so $ \tilde{s} = \tilde{D} ^{1-n_P/4}$. With nP close to 4, any flux rope evolution is close to self-similar. However, in the SW, typically, $n_P \approx 3$ (see Sect. 4.2), then $ \tilde{s} \approx \tilde{D} ^{0.25}$ has a moderate but still a significant evolution, e.g. a factor ${\approx} 3$ for an evolution between 0.1 and 10 AU. Then, in the SW, is the rope evolution close to a self-similar evolution or is there an important reorganization of the field inside the flux rope? How do the flux rope radius and the field strength change with distance from the Sun? To answer these questions we need to solve the above equations for different reference fields  $\vec{ B_{\rm o} }$.

3 Evolution of the normalized magnetic field

We consider below three examples of reference fields and follow their evolution with the distance from the Sun. Generally, the spatial distribution of $ \tilde{B} _z$ and  $ \tilde{B} _{\theta}$, and their functional dependence with the radius, changes during the evolution. The first two cases considered below are exceptions to this general behavior, and they can be solved analytically. This provides a deeper understanding of the flux rope evolution than with a numerical integration of the equations. The evolution of these two extremes cases are compared to the evolution of a reference field, the Lundquist's field, as it often provides a good fit to the observed magnetic field in MCs.

 \begin{figure}
\IfFileExists{color_figures.txt}{
{\includegraphics[width=0.95\te...
...{
{\includegraphics[width=0.95\textwidth, clip=]{0971f2c.eps} }
}\end{figure} Figure 2:

Log-log plots of the normalized central magnetic field strength, $ \tilde{B}_{\rm c} $, and radius $ \tilde{R} $ defined by Eq. (16) in function of the size factor, $ \tilde{s} $ defined by Eq. (19). $ \tilde{B}_{\rm c} $ is defined as the relative evolution of the central to the border field strength normalized to the reference field ( $ \tilde{B}_{\rm c} = \frac{B(0)}{B(R)}\frac{ B_{\rm o} ( R_{\rm o} )}{ B_{\rm o} (0)}$). Three cases suppose an ideal MHD evolution with reference fields (at $ \tilde{s} =1$) being the Gold & Hoyle (GH, Eq. (22)), the ``split'' (SP, Eq. (30)), and the linear force-free field (LU, Eq. (35)). The fourth case (ME, Eq. (45)) always has the magnetic energy to keep the magnetic helicity preserved. The ratio of the azimuthal to the axial flux, $ F_{\rm ratio}$ (computed with a length $ L_{\rm o} = 2\pi R_{\rm o} $) is typically ${\approx } 2$ in MCs ( left panels), while $ F_{\rm ratio}= 5$ illustrates a very twisted flux rope ( right panels). (A color version is available in the electronic version).

Open with DEXTER

3.1 Uniformly twisted field

A uniformly twisted and non-linear force-free model (denoted GH, Gold & Hoyle 1960) is written as

                           $\displaystyle B_{z, {\rm o}}$ $\textstyle \equiv$ $\displaystyle B_{\rm c,o} ~/~ (1+Appl.~Opt.^2~ r_{\rm o} ^2) ,$  
$\displaystyle B_{\theta, {\rm o}}$ $\textstyle \equiv$ $\displaystyle B_{\rm c,o} ~Appl.~Opt.~ r_{\rm o} ~/~ (1+Appl.~Opt.^2~ r_{\rm o} ^2) .$ (22)

The parameter ApplOpt. describes both the spatial extension of the field and the amount of twist. To model a MC field, we consider this model up to the radius $ r_{\rm o} = R_{\rm o} $. In the equations below we use the normalized variables $ \tilde{a}_{\rm o} = Appl.~Opt.~ R_{\rm o} $ and $ \tilde{r}_{\rm o} = r_{\rm o} / R_{\rm o} $.

In order to compare the results of different magnetic field configurations, we compute the ratio of the azimuthal to the axial flux

\begin{displaymath}
F_{\rm ratio,GH} \equiv \frac{L \int_{0}^{ R_{\rm o} } B_{\...
... d} r_{\rm o} }
= Appl.~Opt. R_{\rm o} = \tilde{a}_{\rm o} .
\end{displaymath} (23)

In this expression, we select the length L of the flux-rope portion to be $2 \pi R_{\rm o} $ to simplify the expression (this is not an additional hypothesis since $ F_{\rm ratio}$ is just used below to compare different field distributions for a given flux rope length).

The evolution of a uniformly twisted flux rope is special since, apart from a change in magnitude, the spatial distribution of the twist cannot change due to radial expansion (nor due to axial expansion since we suppose a uniform expansion, Eq. (4)). Since the twist stays spatially uniformly distributed in the flux rope whatever is its evolution, then $B_{\theta} = a r B_z$, with ``a'' being an unknown (to be determined). Including this $B_{\theta }$ in Eq. (1) implies that the field is a GH's field, with only the freedom of modified parameters (called $B_{\rm c} ,a,R$). This result is also directly shown by writing the ratio of Eqs. (18) to (17). The terms in  $ \tilde{r}_{\rm o} ^2$ simplify when a GH reference field (Eq. (22)) is inserted, and the equation can be directly integrated in $ \tilde{r} ^2$ to write

\begin{displaymath}
\tilde{B} ( \tilde{r} ) = \tilde{B}_{\rm c} / \sqrt{1+( \tilde{r} ~ \tilde{a}_{\rm o} / \tilde{s} )^2} ,
\end{displaymath} (24)

where $ \tilde{B}_{\rm c} $ is an integration constant. This is the field strength of a GH model (field components can be easily derived using Eq. (21) and $ \tilde{B} _z=\sqrt{ \tilde{B} ^2 - \tilde{B} _{\theta}^2}$).

The solution of the evolution is completed by integration of Eq. (17) with Eqs. ((22), (24)) included. We find

\begin{displaymath}
\left(1 + \tilde{r} ^2 \tilde{a}_{\rm o} ^2 / \tilde{s} ^2\...
...lde{s} ^2} = 1+ \tilde{a}_{\rm o} ^2~ \tilde{r}_{\rm o} ^2 .
\end{displaymath} (25)

This shows that the expansion of the full flux rope is in general not self-similar as $ \tilde{r} $ is a non linear function of  $ \tilde{r}_{\rm o} $, even if the magnetic field stays a GH field (Eq. (24)). This is so because the radius $ \tilde{R} $ of the flux rope is located at a variable position within the GH profile as $ \tilde{s} $ is changed (this is illustrated later in Fig. 5 with the B components). However, a self-similar expansion is present in the region where the axial field dominated, $ \tilde{a}_{\rm o} ~ \tilde{r}_{\rm o} \ll 1$, since $ \tilde{r} \approx \tilde{r}_{\rm o} /\sqrt{ \tilde{B}_{\rm c} }$ there.

Equation (25) included in Eq. (24) provides $ \tilde{B} $ in function of $ \tilde{r}_{\rm o} $ and allows us to impose the boundary condition on $ \tilde{B} $ (Eq. (20)), defining an equation for the central field strength  $ \tilde{B}_{\rm c} $. We find

\begin{displaymath}
2 \tilde{B}_{\rm c} \ln ( \tilde{B}_{\rm c} ) + \left( \til...
...{s} ^{2}}\right) ~\ln\left(1+F_{\rm ratio,GH}^2\right) = 0 ,
\end{displaymath} (26)

where we replace $ \tilde{a}_{\rm o} $ by $F_{\rm ratio,GH}$ (Eq. (23)) for comparison with other fields. This is a key equation for the GH field evolution, since when  $ \tilde{B}_{\rm c} $ is known, the full magnetic field is known with Eq. (24). It is then worth analyzing the properties of Eq. (26).

The reference state has $ \tilde{s} =1$, then Eq. (26) is naturally satisfied with $ \tilde{B}_{\rm c} =1$ and $ \tilde{s} =1$ for all $F_{\rm ratio,GH}$ values. Since $\ln(1+F_{\rm ratio,GH}^2)>0$, Equation (26) implies

\begin{displaymath}
\tilde{B}_{\rm c} > 1/\sqrt{1 + F_{\rm ratio,GH}^2} \equiv \underline{B}_{\rm ~ GH} .
\end{displaymath} (27)

Deriving Eq. (26) with respect to $ \tilde{s} $, and with the above lower bound, $\underline{B}_{\rm ~ GH}$, one can demonstrate that $ \tilde{B}_{\rm c} $ is a decreasing function of $ \tilde{s} $. This means that the field strength of the core decreases more rapidly than at the flux rope border (where $ \tilde{B} $ is fixed due to the normalization selected in Eq. (16)). We conclude, from flux conservation, that the flux rope core expands faster than its periphery as $ \tilde{s} $ increases. As $ \tilde{s} $ becomes much larger than unity, $ \tilde{B}_{\rm c} $ tends toward $\underline{B}_{\rm ~ GH}$. This is the extreme limit where the core field has expanded so much that it fills most of the flux rope. This limit case has a nearly uniform $ \tilde{B} _z$ field component.

At the opposite limit, for decreasing $ \tilde{s} $ values $ \tilde{B}_{\rm c} $ increases slightly less rapidly than $ \tilde{s} ^{-2}$ (there is a correction depending on $- \ln \tilde{s} $). A simple approximation of $ \tilde{B}_{\rm c} $, with a relative error lower than 10%, for the ranges $0.1<F_{\rm ratio,GH}<10$, $0.1< \tilde{s} <10$ is

\begin{displaymath}
\tilde{B}_{\rm c} = \frac{(1-\underline{B}_{\rm ~ GH})}{ \tilde{s} ^{~1.65}}
+ \underline{B}_{\rm ~ GH} .
\end{displaymath} (28)

Equations ((27), (28)), as well as Fig. 2, show that  $ \tilde{B}_{\rm c} $ changes more with $ \tilde{s} $ as $F_{\rm ratio,GH}$ increases, so for flux ropes with larger twist. Indeed as $F_{\rm ratio,GH}$ increases, larger values of $ \tilde{s} $ are needed to approach the limit $\underline{B}_{\rm ~ GH}$ by the same amount. This means that starting from a more twisted flux rope a longer evolution is needed so that the core field, mainly axial, expands to fill a given spatial fraction of the flux rope.

The normalized radius, $ \tilde{R} $, of the flux rope is found by setting $ \tilde{r}_{\rm o} =1$ in Eq. (25), or equivalently by setting $ \tilde{B} $ in Eq. (24) to its boundary value (Eq. (20)) and solving for $ \tilde{r} = \tilde{R} $. The second way provides a more convenient equation for analysis,

\begin{displaymath}
\tilde{R} = \frac{ \tilde{s} }{ \tilde{a}_{\rm o} } \sqrt{(1+ \tilde{a}_{\rm o} ^2) \tilde{B}_{\rm c} ^2-1} .
\end{displaymath} (29)

A similar analysis of the behavior of $ \tilde{R} $ can be made than for $ \tilde{B}_{\rm c} $ above. However, below we limit the description of the results to the main points. $ \tilde{R} $ is also a decreasing function of $ \tilde{s} $ converging to a finite minimum value for large $ \tilde{s} $ values (Fig. 2). The amplitude of variations also increases with $F_{\rm ratio,GH}$. So $ \tilde{R} $ has a similar behavior as $ \tilde{B}_{\rm c} $, with the main difference that the range of variation is much reduced (Fig. 2). This can also be seen in Eq. (29) where the factor $ \tilde{s} $ is moderating the contribution of  $ \tilde{B}_{\rm c} $ (approximated by Eq. (28)).

3.2 ``Split'' field

We consider in this section the theoretical field where $ B_{z, {\rm o}} $ and $ B_{\theta, {\rm o}} $ are separated in two different regions


$\displaystyle %
r_{\rm o} \leq A_{\rm o} \qquad B_{z, {\rm o}}$ $\textstyle ~\equiv~$ $\displaystyle B_{\rm c,o} ~~, \quad B_{\theta, {\rm o}} \equiv 0
,$  
$\displaystyle A_{\rm o} < r_{\rm o} \leq R_{\rm o} \qquad B_{z, {\rm o}}$ $\textstyle ~\equiv~$ $\displaystyle 0~~,\quad \quad B_{\theta, {\rm o}} \equiv B_{\rm c,o} ~ A_{\rm o} / r_{\rm o}
.$ (30)

We simply refer to this magnetic field as ``SP'' below. This field is an extreme case where each component fully dominates in one region of the flux rope (in general, equilibriums have a gradual transition between a Bz dominated core to a  $B_{\theta }$ dominated periphery). The magnitude of $ B_{\theta, {\rm o}} $ for the limit $ r_{\rm o} \rightarrow A_{\rm o,+}$ is selected to have an equilibrium (magnetic pressure balance with  $ B_{z, {\rm o}} $). Below we use the normalized parameter $ \tilde{A}_{\rm o} = A_{\rm o} / R_{\rm o} $.

The ratio of the azimuthal to the axial flux, defined as in Eq. (23), for this field is

\begin{displaymath}
F_{\rm ratio,SP} = 2 R_{\rm o} / A_{\rm o} ~\ln ( R_{\rm o}...
...m o} ) = - 2 \ln ( \tilde{A}_{\rm o} )~/ \tilde{A}_{\rm o} ~.
\end{displaymath} (31)

Since the reference field (Eq. (30)) has only one non-null component in each region, the evolved field keeps this property (due to flux conservation). It implies a uniform Bz for $r\leq A$ and a 1/r dependence for $B_{\theta }$ for r>A, so that Eq. (1) is satisfied. The integration of Eqs. ((17), (18)), with Eq. (30), the continuity of $ \tilde{r} $ and $ \tilde{B} $ at $ \tilde{r} = \tilde{A} $ and the boundary condition of Eq. (20) provides the equation for $ \tilde{B}_{\rm c} $ as

\begin{displaymath}
\sqrt{ \tilde{B}_{\rm c} } \ln ( \tilde{B}_{\rm c} ) - \lef...
...frac{1}{ \tilde{s} }\right) ~\ln ( \tilde{A}_{\rm o} ) = 0 .
\end{displaymath} (32)

This equation has the same global structure than Eq. (26), then $ \tilde{B}_{\rm c} $ has a similar behavior than for the GH field: it is a decreasing function of $ \tilde{s} $, with $ \tilde{B}_{\rm c} =1$ when $ \tilde{s} =1$, and tends to a finite value ( ${=} \tilde{A}_{\rm o} $) for large $ \tilde{s} $ values. An analytical approximation, with similar precision than for the GH field, is

\begin{displaymath}
\tilde{B}_{\rm c} = \frac{1- \tilde{A}_{\rm o} }{ \tilde{s} ^{1.1}} + \tilde{A}_{\rm o} .
\end{displaymath} (33)

It implies that the dependence of $ \tilde{B}_{\rm c} $ with $ \tilde{s} $ is weaker than for the GH field (Eq. (28), see also Fig. 2). The dependence of $ \tilde{A}_{\rm o} $ with $F_{\rm ratio,SP}$ in Eq. (31) is also similar, while weaker, than for the GH field (Eq. (27)).

The radial expansion of the field is found to be 1.25pt

$\displaystyle \tilde{r}_{\rm o} \leq \tilde{A}_{\rm o} \qquad \tilde{r} = \tilde{r}_{\rm o} / \sqrt{ \tilde{B}_{\rm c} }
,$      
$\displaystyle \tilde{A}_{\rm o} < \tilde{r}_{\rm o} \leq 1 \qquad \tilde{r} = \...
...}
~ \tilde{r}_{\rm o} ^{1/\left( \tilde{s} \sqrt{ \tilde{B}_{\rm c} }\right)}
.$     (34)

implying a self-similar expansion of the core and a periphery strongly reduced in size as $ \tilde{s} $ increases. The flux rope radius is simply $ \tilde{R} = \sqrt{ \tilde{B}_{\rm c} }$. This is a much simpler link to $ \tilde{B}_{\rm c} $ than for the GH field (Eq. (29)), while it has a similar dependence (Fig. 2).

3.3 Linear force-free field

As an example of the general evolution of a flux rope field obtained by solving numerically Eqs. ((17), (18), (20)) we select the linear force-free model (Lundquist 1950) since it has shown good fits to various MC observations,

$\displaystyle B_{z, {\rm o}} \equiv B_{\rm c,o} ~J_0( \alpha_{\rm o} r_{\rm o} ) ,$      
$\displaystyle B_{\theta, {\rm o}} \equiv B_{\rm c,o} ~J_1( \alpha_{\rm o} r_{\rm o} ) ,$     (35)

where J0 and J1 are the ordinary Bessel functions of order 0 and 1, respectively. This magnetic field is simply noted LU. The ratio of the azimuthal to the axial flux, defined as in Eq. (23), is

\begin{displaymath}
F_{\rm ratio,LU} = \frac{1-J_0( \alpha_{\rm o} R_{\rm o} )}{J_1( \alpha_{\rm o} R_{\rm o} )} \cdot
\end{displaymath} (36)

Fits of Eq. (35) to observed MCs typically implies $ \alpha_{\rm o} R_{\rm o} \approx z_{J_0} \approx 2.4$, where the zJ0 is the first zero of the Bessel function J0, so that the axial field is small at the MC border. It implies that the flux ratio, as defined above, is around ${\approx } 2$ in observed MCs (more precisely $F_{\rm ratio,LU}=1.93$ for $ \alpha_{\rm o} R_{\rm o} =z_{J_0}$, and $1.26<F_{\rm ratio,LU}<3.2$ for $0.8~z_{J_0}< \alpha_{\rm o} R_{\rm o} <1.2~z_{J_0}$). The evolution of $ \tilde{B}_{\rm c} $ and $ \tilde{R} $ with $ \tilde{s} $ is shown in Fig. 2.

3.4 Comparison of the field evolutions

Below, we compare the evolution of the three fields defined in the above sections. Guided by the results on MCs, we compare the fields for $ F_{\rm ratio}=2$ (Fig. 2, left panels). We also show the case $ F_{\rm ratio}= 5$ to illustrate an extreme case of very twisted flux ropes (not observed). Such an extreme case is selected for viewing purpose since the effect of $ F_{\rm ratio}$ is weak. We do not plot cases with $ F_{\rm ratio}<2$ since they have an even lower dependence with $ \tilde{s} $, and all the fields have closer values of $ \tilde{B}_{\rm c} $ and $ \tilde{R} $ as $ F_{\rm ratio}$ decreases.

As for the GH and SP fields (Eqs. (28), (33)), $ \tilde{B}_{\rm c} $ tends towards a fixed value for large value of $ \tilde{s} $ for the LU field (Fig. 2, top panels). This corresponds to the case where the core field has expanded to fill most of the flux rope, so the central field $ \tilde{B}_{\rm c} $ is almost equal to the field strength at the boundary. For $ F_{\rm ratio}=2$, this limit is nearly the same for the three reference fields. In view of the strong differences in their twist distribution, we conclude that it is a general property for all field distributions found in MCs. Differences between the reference fields are only found for larger values of $ F_{\rm ratio}$. For $ F_{\rm ratio}= 5$, the LU and SP results are still very close together (Fig. 2, right panels).

A relatively strong variation of $ \tilde{B}_{\rm c} $ is present for small $ \tilde{s} $ values. This is because the region dominated by $B_{\theta }$ gets increasingly more extended compared to the core as $ \tilde{s} $ becomes small. This implies an increasingly stronger variation of the field strength across the flux rope, then a central field less affected by the boundary field value, and so a larger  $ \tilde{B}_{\rm c} $. For $ F_{\rm ratio}=2$, the GH and LU fields have very similar values of $ \tilde{B}_{\rm c} $, while the SP field has significant lower values (Fig. 2, top left panel).

The results for $ \tilde{R} $ are similar to the above results for $ \tilde{B}_{\rm c} $, with some minor differences that we do not describe (see Fig. 2, bottom panels). A main difference is the magnitude of the variations, with typically a factor 3 lower variation of $\log_{10} \tilde{R} $ than of $\log_{10} \tilde{B}_{\rm c} $.

While we consider three fields with very different twist profiles, the evolution of $ \tilde{B}_{\rm c} $ and $ \tilde{R} $ are remarkably similar. This is obviously the case for a low ratio of the azimuthal to axial flux, $ F_{\rm ratio}$, since all the fields are dominated by their axial component and therefore these fields are similar. More interesting, this is approximately the case with $ F_{\rm ratio}\approx 2$, a case typical of the values found in MCs. Even for the case of flux ropes much more twisted than MCs, the evolution of $ \tilde{B}_{\rm c} $ and $ \tilde{R} $ are similar despite the very important differences in the field profiles (see Sect. 4.3 and Fig. 5).

We conclude that the evolution of $ \tilde{B}_{\rm c} $ and $ \tilde{R} $ is weakly affected by the type of field profile. Moreover, the amplitude of the variation of $ \tilde{B}_{\rm c} $ and $ \tilde{R} $ versus $ \tilde{s} $ remains moderate, especially for $ F_{\rm ratio}\leq 2$. We compare this moderate evolution to the much stronger effect of $ \tilde{ P_{\rm t} } $ in the field evolution in the next section.

4 Evolution of flux ropes in the solar wind

4.1 Flux rope length

The presence of counter-streaming suprathermal electrons in the SW is one of the primary signatures used to identify ICMEs (e.g., Gosling et al. 1987). It is also generally taken as an evidence that the corresponding field line is still rooted at both ends to the Sun since the solar corona is a reservoir of fast electrons (see e.g., Crooker & Horbury 2006; Wimmer-Schweingruber et al. 2006, for a review as well as possible caveats). Shodhan et al. (2000) found that periods of counter-streaming electrons are randomly distributed in 48 MCs observed near 1 AU, with a mean frequency of 59% of the MC duration. A very similar fraction, 55%, is found for 31 MCs observed at 5 AU indicating that a large fraction of a typical MC is still linked to the Sun (Crooker et al. 2004).

When the flux rope is attached to the Sun, the flux rope length, L, is mostly proportional to the distance, D, to the Sun (with only a weak correction depending on the evolution of the shape of the flux rope axis). When the flux rope becomes detached from the Sun, by reconnecting with the SW magnetic field, the evolution of L is more involved since torsional Alfvén waves are launched in the SW field. Here we simply suppose that the axial flux rope expansion evolves proportionally to its distance to the Sun. With the normalization to the reference field (Table 1) the coefficient of proportionality simplifies and the normalized length is simply

\begin{displaymath}
\tilde{L} = \tilde{D} .
\end{displaymath} (37)

The implications of this assumption are relatively weak since $ \tilde{L} $ enters in Eqs. ((14), (15)) only with the size factor $ \tilde{s} $ (Eq. (19)), while the evolution of $ P_{\rm t} $ is dominant in the evolution of the magnetic field (see below).

4.2 Solar wind total pressure

The total pressure in the SW, $ P_{\rm t} $, has two main components: the plasma and the magnetic pressure.

Since the solar wind speed is radial and almost independent of solar distance, D, various investigations found a plasma density decreasing as D-2 (e.g., Gazis et al. 2006; Issautier et al. 1998), in agreement with mass conservation. A wider variety of results are found for the electron temperature. The data are classically fit with a power-law $T \propto D^{-n_T}$. nT is found to be typically in the range [0.3,1] (see Fig. 1 in Issautier et al. 1999). Recent estimations of nT are typically around 0.7, e.g. $n_T=0.64\pm 0.03$ (Issautier et al. 1998) and $n_T=0.78\pm 0.03$ (Gazis et al. 2006). The estimated values of nT are a function of the type of SW considered (slow/fast), of the removal or not of the ICMEs from the processed data, and of the range of distance considered. nT is typically observed to be lower in the fast SW and at larger D values, a tendency well reproduced by a kinetic collision-less SW model (Meyer-Vernet & Issautier 1998).

The mean SW magnetic field is well described by the Parker model (Parker 1958) in most of the interplanetary space between 0.3 to 5 AU (Forsyth et al. 1996, and references therein). Most of the observed deviations occur at high latitude where the Ulysses spacecraft found a slightly under wound magnetic field by ${\approx} 10{-}20^{0}$, see Zurbuchen (2007) for a review. Then the SW magnetic field strength is approximately

\begin{displaymath}
B = B_{r, {\rm o}} \sqrt{ \tilde{D} ^{-4}+( \tilde{V} ~ \tilde{D} )^{-2}} ,
\end{displaymath} (38)

where $B_{r, {\rm o}}$ is the radial field value at 1 AU, $ \tilde{D} =D/ D_{\rm o} $ is the distance normalized at $ D_{\rm o} =1$ AU, and $ \tilde{V} = V_{\rm SW} / (\Omega D_{\rm o} \cos \theta )$ with $\theta$ the latitude and $\Omega$ the angular velocity of the Sun ( $\Omega\approx 2.8\times 10^{-6}$ s-1). At 1 AU in the ecliptic plane, $ \tilde{V} \approx 1$, since typically $V_{\rm SW} \approx 400$ km s-1.

Following the above results, the total SW pressure is written in the undimensional form

\begin{displaymath}
P_{\rm t} = \frac{P_{\rm t,o}}{1+\beta_{\rm o}}
\left( \f...
...tilde{V} ^{-2}} + \beta_{\rm o} \tilde{D} ^{-2-n_T} \right) ,
\end{displaymath} (39)

where $P_{\rm t,o}$ and $\beta _{\rm o}$ are the total SW pressure and the plasma $\beta $ at $ D_{\rm o} =1$ AU.

Equation (39) provides the boundary condition, $ P_{\rm t} $, for the flux rope expansion (Eq. (6)). In the above model, $ P_{\rm t} (D)$ is a key function for the evolution of both the magnetic field strength, B, and the flux rope radius, R, (Eqs. (7), (16)). In Sect. 4.4, we compare the results of the model to studies of MCs/ICMEs that derived the evolution of both B and R with D from statistical analysis of a wide set of events observed at different D. All these studies fit the observed results with power law of D. Then, below we consider an approximation of Eq. (39) with a power law

\begin{displaymath}
P_{{\rm t},a} = P_{\rm t,o} \tilde{D} ^{-n_P} .
\end{displaymath} (40)

With this approach, we can also compare the above models to observations of expansion velocity measured at 1 AU (Sect. 4.5) if we derive the exponent nP at 1 AU.

The logarithmic derivative $ {\rm d} \ln P_{\rm t} / {\rm d} \ln D$, called n, defines the local decrease of $ P_{\rm t} $ (i.e. a local approximation of $ P_{\rm t} $ is $P_{{\rm t},a}$ with nP=n). For $ \tilde{D} \ll 1$, $n\approx 4$ since the radial component of the SW magnetic field provides the dominant pressure close to the Sun. At the opposite extreme, for $ \tilde{D} \gg 1$, $n\approx 2$. In between these two extremes, n varies between ${\approx} 3.5$ to 2.5 for D varying from 0.3 to 3 AU.

At 1 AU, the radial and azimuthal magnetic field components are comparable in the slow SW present at low latitude ( $ \tilde{V} \approx 1$). Then, the magnetic pressure alone implies $n \approx 3$ at 1 AU. The plasma pressure alone implies $n \approx 2+n_T \approx 2.7$ with the recent results summarized above. Then, around 1 AU the value of n is weakly dependent on the plasma $\beta $ since the magnetic and plasma pressure provides similar n values. With $ \tilde{V} \approx 1$, we get $n = 2.8 \pm 0.2$ for $\beta _{\rm o}$ in the range [0.3,3]. Figure 3 shows that $P_{{\rm t},a}$ (Eq. (40)), with nP = 2.8 is a good approximation of $ P_{\rm t} $ (Eq. (39)) within a relatively wide range of solar distances in the slow SW. Finally, this value of n is only increased to $3.2\pm 0.4$ for a fast SW ( ${\approx} 800$ km s-1) and at high latitude ( ${\approx} 60\hbox{$^\circ$ }$), so for $ \tilde{V} \approx 4$, since the radial component of the magnetic field is then more dominant at 1 AU.

 \begin{figure}
\IfFileExists{color_figures.txt}{
{\includegraphics[width=0.45\te...
...
{\includegraphics[width=0.45\textwidth, clip=]{0971f3bc.eps} }
}\end{figure} Figure 3:

Log-log plots showing the evolution of the SW total pressure (Eq. (39)) normalized at 1 AU, and a power-law approximation (Eq. (40) with nP=2.8) as a function of the solar distance D (in AU). $\beta _{\rm o}$ is the plasma $\beta $ in the SW at 1 AU, and nT define the SW temperature dependence with D ( $T \propto D^{-n_T}$). In both panels $ \tilde{V} =1$, a typical value for the slow SW at 1 AU and low latitude (A color version is available in the electronic version).

Open with DEXTER

4.3 Expected evolution of flux ropes

With Eqs. ((37), (40)), (19) is rewritten as

\begin{displaymath}
\tilde{s} = \tilde{D} ^{1-n_P/4} .
\end{displaymath} (41)

With nP = 2.8, $ \tilde{s} = \tilde{D} ^{0.3}$, so $ \tilde{s} $ has a moderate variation with $ \tilde{D} $. For example, for $ D_{\rm o} =1$ AU and D ranging from 0.1 to 10 AU, $\log_{10} \tilde{s} $ is in the interval [-0.3,0.3]. It corresponds to a narrow central portion of Fig. 2, so to very moderate variations of both $ \tilde{R} $ and $ \tilde{B}_{\rm c} $ ( ${\approx} \tilde{D} ^{-0.07}$ and ${\approx} \tilde{D} ^{-0.2}$, respectively, in the vicinity of $ \tilde{D} =1$). This contrasts with the major contribution of SW total pressure ( $ \tilde{ P_{\rm t} } $) in R and $B_{\rm c} $ in Eq. (16), as $ \tilde{D} ^{0.7}$ and $ \tilde{D} ^{-1.4}$, respectively.

More generally, with Eqs. (40), (16) implies


    $\displaystyle R = R_{\rm o} ~( \tilde{D} )^{ n_P/4} ~ \tilde{R} (n_P, \tilde{D} ) ,$  
    $\displaystyle B_{\rm c} = B_{\rm c,o} ~( \tilde{D} )^{-n_P/2} ~ \tilde{B}_{\rm c} (n_P, \tilde{D} ).$ (42)

Figure 4 shows typical cases of observed MCs with $ F_{\rm ratio}=2$. All the three models considered above have a central field strength decreasing only slightly faster with solar distance than the D-nP/2 scaling law. The radius of the flux rope increases with D slightly less than D nP/4 since the correction to the scaling law is even smaller than for $B_{\rm c} $. The slope difference becomes slightly more important as the flux rope becomes more twisted (larger $ F_{\rm ratio}$, for example $ F_{\rm ratio}= 5$ in the right panels of Fig. 4).

The magnetic profile evolves with solar distance according to the results of Sect. 3. The core region expands more with increasing D than the periphery (panels of Bz and $B_{\theta }$ in Fig. 5 and panels of  $r( r_{\rm o} )$ in Fig. 6). This is because the SW total pressure does not decrease fast enough (less than D-4) to allow a true self-similar expansion of the flux rope. Indeed, the peripheral region is dominated by the azimuthal component. Its magnitude is mainly fixed by the total pressure $ P_{\rm t} $ present in the surrounding SW. The extension in length, $L\propto D$, and the flux conservation restrict the radial expansion. As D increases, the field strength profile becomes flatter with the azimuthal component getting weaker than the axial component in most of the flux rope.

 \begin{figure}
\IfFileExists{color_figures.txt}{
{\includegraphics[width=0.90\te...
...{
{\includegraphics[width=0.90\textwidth, clip=]{0971f4c.eps} }
}\end{figure} Figure 4:

Log-log plots of the central magnetic field strength, $B_{\rm c} $, and radius, R, as a function of the distance from the Sun, D. We set the typical values $B_{\rm c} = 20$ nT and R = 0.1 AU at D=1 AU, and nP=2.8 (Eq. (40), Fig. 3). Three cases, the Gold and Hoyle (GH, Eq. (22)), ``split'' (SP, Eq. (30)), and linear force-free field (LU, Eq. (35)), suppose an ideal-MHD evolution with reference fields defined at $D= D_{\rm o} $. The fourth case (ME) always has the minimum magnetic energy keeping the magnetic helicity preserved. The straight lines labeled ``scaling'' have the slopes -nP/4 and nP/2 for $B_{\rm c} $ and R, respectively (Eq. (42)). The ratio of the azimuthal to the axial flux, $ F_{\rm ratio}$, is computed with a length $ L_{\rm o} = 2\pi R_{\rm o} $ (A color version is available in the electronic version).

Open with DEXTER

 \begin{figure}
\IfFileExists{color_figures.txt}{
{\includegraphics[width=\textwi...
... }
}{
{\includegraphics[width=\textwidth, clip=]{0971f5c.eps} }
}
\end{figure} Figure 5:

Evolution of the field strength, B, the axial, Bz, and azimuthal, $B_{\theta }$, field components with the distance D (D = (0.1, 1, 10) AU in the top, middle, bottom panels, respectively). The reference field is taken at $D=1~\rm AU$ with the typical value $B_{\rm c} = 20$ nT and R=0.1 AU. The same four cases as in Fig. 4 with $ F_{\rm ratio}=2$ are shown (A color version is available in the electronic version).

Open with DEXTER

4.4 Comparison to previous studies

We compare below our results with previous studies of MCs and ICMEs.

Chen (1989) developed a toroidal model of a flux rope ejected from the Sun. The governing equations for the flux rope center and radius are deduced from MHD equations by approximately averaging the magnetic and plasma force over the flux rope cross section. Then, the deduced ordinary differential equations are solved starting from an initial unstable equilibrium in the corona. Chen (1996) extended the model to study the flux rope evolution in the interplanetary space using a SW model similar to the one described in Sect. 4.2. The initially stable flux rope in the corona is forced to erupt by imposing a strong increase of the azimuthal flux. However, while important for the physics of CMEs, the precise mechanism leading to the ejection should have a weak effect on the propagation of the interplanetary flux rope (for the same flux rope having the same global velocity).

The flux rope model of Chen (1996) keeps the inertial term and the plasma pressure in the internal dynamics while it only includes an approximate description of the magnetic forces. The model finally uses a spatial average of the force balance across the flux tube. Rather, we suppose that a force-free equilibrium is achieved at all times, and we analyze the internal force balance with several spatial profiles of the magnetic field. Despite these differences, the diameter of the flux rope has a similar dependence on the solar distance D (Table 2). This is so because the inertial term, the plasma pressure and the magnetic profile have only a weak effect on the diameter evolution. Most of the evolution is imposed by the dependence of  $ P_{\rm t} (D)$ in the surrounding SW.

Kumar & Rust (1996) developed a flux rope model based on the conservation of axial magnetic flux and magnetic helicity, supposing that the flux rope field stays linear force-free (see Sect. 5). They also analyzed a set of MCs previously fit to a LU field (see Sect. 3.3) by various authors.

Bothmer & Schwenn (1998) analyzed a set of ${\approx} 50$ MCs. The size S is defined by the product of the MC duration and the mean measured velocity (so they determine the apparent size of the MC along the spacecraft trajectory). Leitner et al. (2007) extended the analysis to 130 MCs, and they fit the magnetic field observations with the LU field (see Sect. 3.3), then they found the impact parameter and the orientation of the MC axis to estimate the diameter, d, of the MCs.

In all these studies, the evolution of the MC radius (or size) is described by a power law, ${\propto} D^e$, in some range of distance, D, to the Sun. The above studies and ours agree on a larger exponent e in the internal heliosphere than outward. Indeed, the general tendency is $e\approx 1$ for D<1 AU, so that the expansion is close to self-similar, while $e \approx 0.7$ for D>1 AU (Table 2). Around 1 AU, we find a power law exponent closer to the one of Leitner et al. when the full range of D is considered ([0.3, 7] AU).

Table 2:   Summary of results on MCs and ICMEs.

We also compare our results with ICME studies (containing a subset of MCs). Liu et al. (2005) did a statistical study of about 220 ICMEs observed in situ from 0.3 to 5.4 AU. They combined two criteria to localize ICMEs: an enhanced alpha-to-proton density ratio (greater than 0.08) and a low proton temperature (lower than half the value expected in the solar wind with the same speed). Wang et al. (2005) did a similar analysis on a wider set of 600 possible ICMEs defined only by using the low temperature criteria. The size S is defined as for Bothmer & Schwenn (1998). These two similar studies have surprisingly different exponents for the power law fits (Table 2). Our result best agree with Wang et al. (2005), who found the same exponent as Leitner et al. (2007).

The same comparison is done for the magnetic field strength (with the exception of Chen (1996) who did not wrote this result). Our exponent is closer to the results of Leitner et al. (2007) within the inner heliosphere (while for the diameter, the agreement is best for the full range of D, Table 2). Our results are also compatible with both ICME studies within the uncertainties of the exponents (however, we do not compute a mean field, as for ICME studies).

We conclude that, taking into account the uncertainties, the field strength versus the distance predicted by our model is in good agreement with the observations. We also find that the relation for the flux rope diameter, and for the field strength, versus solar distance is not strictly a power-law (Fig. 4). This could be the main reason for the different results reported in Table 2 since the computed exponents dependent on the distribution of the analyzed events with solar distance in the different data sets.

 \begin{figure}
\IfFileExists{color_figures.txt}{
{\includegraphics[width=\textwi...
... }
}{
{\includegraphics[width=\textwidth, clip=]{0971f6c.eps} }
}\end{figure} Figure 6:

Evolution of a flux rope for the three ideal-MHD cases shown in Fig. 5 (D = (0.1, 1, 10) AU in the top, middle, bottom panels, respectively). r is the radial coordinate of the flux rope, with $ r_{\rm o} $ being its value in the reference field taken at $ D_{\rm o} =1$ AU (with the flux rope radius $ R_{\rm o} =0.1$ AU). v is the radial plasma velocity. v D /(r Vc) is the normalized expansion velocity (Sect. 4.5). We set $ F_{\rm ratio}=2$, nP=2.8, and $V_{\rm c} =600$ km -1 ($V_{\rm c} $ is the global velocity of the flux rope away from the Sun). The horizontal line in the right panels has an ordinate of nP/4 (A color version is available in the electronic version).

Open with DEXTER

4.5 Expansion velocity

Since the time evolution of the flux rope is ideal we compute the plasma velocity by deriving  $r( r_{\rm o} )$ with time, keeping $ r_{\rm o} $ fixed (so following the same blob of plasma). With Eqs. ((16), (40)) we find

\begin{displaymath}
v = \left. \frac{ {\rm d} r}{ {\rm d} t} \right\vert _{ r_{...
...{\rm o} }
+ \frac{n_P}{4} \right] ~\frac{V_{\rm c} ~r}{D} ,
\end{displaymath} (43)

where $V_{\rm c} = {\rm d} D/ {\rm d} t$ is the velocity of the flux rope center.

From the results of Sect. 3 and with $ \tilde{s} $ written as a function of $ \tilde{D} $ (Eq. (41)), we find that $ {\rm d} \ln \tilde{r} / {\rm d} \ln D$ is small compare to nP/4, so that the velocity inside the flux rope is approximately proportional to r (Fig. 6). A small correction to this is due to a weaker expansion of the flux-rope periphery (Sect. 4.3). This implies a smaller velocity there (than if the whole flux rope had a self-similar expansion).

The velocity profile becomes a linear function of r with increasing D since the dominated axial field region extends to a larger fraction of the flux rope (Fig. 5). Moreover, the expansion speed decreases with distance because r increases less rapidly than D. The expansion speed decreases by about a factor of two with a factor of ten increase in distance. Another effect to take into account for a wide range of D is the decrease of nP with D when the full $ P_{\rm t} (D)$ is used (Eq. (39)), but here this weaker effect is not taken into account since we use Eq. (40) with nP=2.8. This effect is at most a factor of 2 from $D\ll 1$ AU to $D\gg 1$ AU) since $n= {\rm d} \ln P_{\rm t} / {\rm d} \ln D$ evolves at most from 4 close to the Sun to 2 at very large distances (Sect. 4.2).

Next, it is remarkable that ``GH'' and ``LU'' fields have such similar velocities. Only the ``SP'' field, with its separation of field components into distinct regions, has a slightly different profile. We conclude that the velocity profiles measured in MCs are expected to be almost independent of the magnetic profile observed.

The magnitude of the expansion is directly proportional to the global velocity of the MC ($V_{\rm c} $), so that the results of Fig. 6 are scaled by $V_{\rm c} $ (Eq. (43)). In fact $V_{\rm c} $ has an even stronger effect on the observed velocity profiles since they are function of time, t, not of r, as follows. Let us just consider here the simplest case where the spacecraft crosses the center of the flux rope nose moving at a constant velocity $V_{\rm c} $, then the position observed by the spacecraft is $r=V_{\rm c} (t-t_{\rm c})$ with $t_{\rm c}$ being the time when the spacecraft crosses the flux rope center. Then, Eq. (43) rewritten in function of the observed time (t) is

\begin{displaymath}
v(t) = \left[ \left. \frac{ {\rm d} \ln \tilde{r} ( r_{\rm ...
...n_P}{4} \right] ~\frac{V_{\rm c} ^{2}~(t-t_{\rm c})}{D} \cdot
\end{displaymath} (44)

This implies that the magnitude of the observed expansion velocity scales as  $V_{\rm c} ^{2}$. This effect needs to be taken into account when comparing the difference between the front and back boundary velocities of various MCs. Fast and/or big MCs have faster apparent expansion velocity due to this effect. For example this is the case of the ``monster'' MC observed on 29 Oct. 2003 (Mandrini et al. 2007). While the observed expansion speed is as large as 200 km s-1, they found that this MC has the typical expansion rate obtained for more classical MCs when the normalized expansion rate is considered (i.e. a similar slope of $v(t)DV_{\rm c} ^{-2}$).

The generalization of the above results to all possible types of spacecraft/MC crossing can be derived following the same logic than as Démoulin et al. (2008). They consider self-similar expansion of a magnetic configurations with possible different expansion coefficients, l,m,n in three orthogonal directions (i.e., along the axial direction $z(t)=z_{\rm o} (D(t)/ D_{\rm o} )^{n}$, and similar expressions in the x,y directions). Otherwise any orientation of the flux rope, any impact parameter and any radial motion of the flux rope away from the Sun can be selected. Redoing the analysis with the above flux ropes, we get the same equations for the velocity components at the spacecraft location, except that n=1 (from Eq. (37)) and $l=m=v D /(r V_{\rm c} )= {\rm d} \ln \tilde{r} / {\rm d} \ln D \vert _{ r_{\rm o} } + n_P/4$. Then l,m are no longer constant across the flux rope. Still, nP/4 is the dominant term and only weak variations are present (Fig. 6, right panels). As expected from Figs. 2, 4, the core expands slightly more rapidly than the periphery of the flux rope.

Finally, the above theoretical results on the expansion rate are compared to the expansion rate found in 26 typical MCs observed at 1 AU (Démoulin et al. 2008). They fit the observed main velocity component (anti-solar direction) with a linear function of time. The observed slope, $ {\rm d} V_{x,GSE} / {\rm d} t$, normalized by $V_{\rm c} ^{2}/D$ provides an estimation, called $\zeta$, of the exponent l, since the MCs are observed close to their noses. They found $\zeta= 0.8 \pm 0.2$. This is indeed in close agreement with our theoretical expectations since we find $l=v D /(r V_{\rm c} )= 0.75\pm 0.1$ across all the theoretical flux ropes at 1 AU (Fig. 6, recalling that the ``SP'' case is far from observed MCs). We conclude that the expansion of MCs is dominantly controlled by the decrease of the total SW pressure with solar distance.

5 Evolution of flux ropes with energy minimization

In previous sections we have studied the evolution of flux ropes supposing an ideal MHD evolution, so preserving both the elementary axial and azimuthal magnetic fluxes. Here we investigate another extreme limit where the flux rope evolves with all the necessary reconnections to minimize its magnetic energy, keeping constant its magnetic helicity. The purpose is to check if the results of previous sections are robust enough. The true evolution of MCs is in between these two extreme limits.

5.1 Linear force-free state

In a resistive MHD evolution, for systems with low resistivity as the solar wind, the magnetic field cannot relax to its absolute minimum energy since its magnetic helicity, H, is almost conserved (on a time-scale less than the global diffusion time-scale, that is typically many orders of magnitude above the evolution time scales, Matthaeus & Goldstein 1982; Berger 1984). If there is enough dissipation (e.g., via magnetic reconnection or turbulent dissipation) to remove all the free magnetic energy, the magnetic field relaxes to a linear force-free field (Taylor 1974). We assume such minimization of the magnetic energy under the constraint of a fixed H in present section (this is just an extreme possibility, we are not claiming that this energy minimization is fully achieved in MCs). This implies that the magnetic field stays a Lundquist's field in the flux rope (r<R)

$\displaystyle B_{z} = B_{\rm c} ~J_0(\alpha r) ,$      
$\displaystyle B_{\theta } = B_{\rm c} ~J_1(\alpha r) ,$     (45)

where $B_{\rm c} $, $\alpha$ and R are evolving parameters with distance to the Sun. Such field is labeled ``ME''.

Our approach is similar to the one of Kumar & Rust (1996). A main difference is that we do not set an external curl-free magnetic field around the flux rope. Indeed, in situ observations do not show evidence of this purely azimuthal field around MCs but rather the presence of the SW (with its variable magnetic field). This difference of treatment of the region outside the flux rope is at the primary differences between Kumar & Rust (1996) and our results (the scaling with solar distance, Table 2, as well as the evolution of the magnetic field profile, Fig. 5).

5.2 Formal solution of the evolution

The magnetic helicity contained in a cylindrical flux rope of length L is (Berger 1999)

\begin{displaymath}
H = 4 \pi L~ \int_{0}^{R} A_\theta (r) ~B_\theta (r) ~r ~ {\rm d} r ,
\end{displaymath} (46)

where $A_\theta$ is the azimuthal component of the vector potential of $\vec{B}$. For the Lundquist's field, the magnetic helicity is (Dasso et al. 2003)
                                 H = $\displaystyle 2 \pi L~B_{\rm c} ^2 R^3 ~\frac{J_1^2(U) -J_0(U)~J_2(U)}{U}$  
  $\textstyle \equiv$ $\displaystyle 2 \pi L~B_{\rm c} ^2 R^3 ~f_H(U),$ (47)

with $U=\alpha R$, Jn is the ordinary Bessel function of order n, and fH(U) defines a function of U alone.

We next suppose that there is no reconnection with the solar wind magnetic field. Then, the flux rope stays attached to the Sun, and the total axial magnetic flux, Fz is conserved, where Fz can be written as

                   Fz = $\displaystyle 2 \pi ~B_{\rm c} R^2 ~J_1(U)/U$  
  $\textstyle \equiv$ $\displaystyle 2 \pi ~B_{\rm c} R^2 ~f_F(U).$ (48)

If the flux rope partly reconnects with the overtaken solar wind magnetic field, being partially peeled, as observed in some cases (Dasso et al. 2006,2007), the present derivation could be applied only to the un-reconnected central part of the flux rope.

The pressure equilibrium at the boundary (Eq. (6)) is written as

                 $\displaystyle \varpi (D) ~ P_{\rm t} (D)$ = $\displaystyle B_{\rm c} ^2 /(2 \mu_0) ~(J_0^2(U)+J_1^2(U))$  
  $\textstyle \equiv$ $\displaystyle B_{\rm c} ^2 /(2 \mu_0) ~f_P(U).$ (49)

There are three unknowns, $B_{\rm c} $, R and U (or equivalently $\alpha$), and three Eqs. ((47)-(49)) that define the evolution of the flux rope. With the normalization of B and R defined by Eq. (16), the three equations write


    $\displaystyle \tilde{s} ~ \tilde{R} ^3 ~ \tilde{B}_{\rm c} ^2 ~f_H(U) = f_H( U_{\rm o} ) ,$  
    $\displaystyle ~ \tilde{R} ^2 ~ \tilde{B}_{\rm c} ~~f_F(U) = f_F( U_{\rm o} ) ,$ (50)
    $\displaystyle \tilde{B}_{\rm c} ^2 ~f_P(U) = f_P( U_{\rm o} ) .$  

As in Sect. 2.4, the normalized field strength, $ \tilde{B}_{\rm c} $, and radius $ \tilde{R} $ only depends on $ \tilde{s} = \tilde{L} ~( \tilde{ P_{\rm t} } )^{1/4}$. Defining the function f(U) as

\begin{displaymath}
f(U) \equiv f_F^{3/2}(U)~f_P^{1/4}(U)~/f_H(U) ,
\end{displaymath} (51)

the solution of the system of Eqs. (50) is


    $\displaystyle f(U) = \tilde{s} ~f( U_{\rm o} ) ,$  
    $\displaystyle \tilde{B}_{\rm c} = \left( \frac{f_P(U)}{f_P( U_{\rm o} )} \right) ^{-1/2} ,$ (52)
    $\displaystyle \tilde{R} = \left( \frac{f_P(U)}{f_P( U_{\rm o} )} \right) ^{1/4}
~\left( \frac{f_F(U)}{f_F( U_{\rm o} )} \right) ^{-1/2} \cdot$  

$B_{\rm c} $ and R are obtained by replacing the above results for $ \tilde{B}_{\rm c} $ and $ \tilde{R} $ in Eq. (16).

 \begin{figure}
\IfFileExists{color_figures.txt} {
{\includegraphics[width=0.47\t...
...
{\includegraphics[width=0.47\textwidth, clip=]{0971f7c.eps} }
}
\end{figure} Figure 7:

The evolution of $U=\alpha R$ in function of $\log_{10}( \tilde{s} )$ (Eq. (19)) for different ratios of the azimuthal to axial flux ( $ F_{\rm ratio}$, Eq. (36)). The magnetic energy is supposed to be minimized with a preserved magnetic helicity, so the field stays a linear force-free field. The thick curves show the solutions of the first equation in Eq. (52), while the thin curves are for its approximation (Eq. (54)). (A color version is available in the electronic version).

Open with DEXTER

5.3 Evolution of flux ropes

As the flux rope gets farther from the Sun, its length L increases and the total pressure in the surrounding solar wind, $ P_{\rm t} $, decreases. Both contribute to the evolution of $ \tilde{s} $. For the case, encountered in the SW, where $ \tilde{ P_{\rm t} } $ decreases less rapidly than  $ \tilde{L} ^{-4}$, $ \tilde{s} $ increases with distance from the Sun (see Eq. (19)). This evolution of $ \tilde{s} $ determines the evolution of $U=\alpha R$ (first equation of (52)), and of the normalized profile of the field (it is function only of $\alpha r = U~r/R = U~ \tilde{r} / \tilde{R} $). Then, Eq. (45) is rewritten as
$\displaystyle B_{z} = B_{\rm c,o} ~\sqrt{ \tilde{ P_{\rm t} } } ~ \tilde{B}_{\rm c} ~J_0(U~ \tilde{r} / \tilde{R} ) ,$      
$\displaystyle B_{\theta } = B_{\rm c,o} ~\sqrt{ \tilde{ P_{\rm t} } } ~ \tilde{B}_{\rm c} ~J_1(U~ \tilde{r} / \tilde{R} ) ,$     (53)

for $0< \tilde{r} < \tilde{R} $. $ \tilde{B}_{\rm c} $, $ \tilde{R} $ are provided by the two last equations of (52).

We analyze first the evolution of U since it determines the spatial behavior of the field components. Two values of U are worth outlining. Let us call zJ0 and zJ1 the first zero of the Bessel functions J0 and J1, respectively. For U=zJ0 the axial field Bz vanishes (it frequently takes place near the observed flux rope boundary), while for U=zJ1, the axial magnetic flux Fz vanishes (Eq. (48)). From the MC data fitted with the Lundquist field, the value of U is typically found to be around zJ0, while generally U<zJ1. Then, the range $0< U_{\rm o} <z_{J_1}$ includes the observed MCs, and Fz(U) is a single-value and monotonically decreasing function of U. This implies that U stays in the same interval as $ U_{\rm o} $ to keep Fz preserved.

The evolution of U is determined by the behavior of the function f(U) defined by Eq. (51). It is a monotonically decreasing function of U for 0<U<zJ1, with the limits $f(U)\approx 2 \sqrt{2} /U$ for $U\ll 1$ and f(zJ1)=0. A simple approximation, only slightly underestimating f(U), is the function $2 \sqrt{2} (1/U-1/z_{J_1})$. Including this approximation in the first equation of (52), we get the approximate solution, plotted in Fig. 7,

\begin{displaymath}U \approx \frac{z_{J_1}}{ \tilde{s} ~(z_{J_1}/ U_{\rm o} -1)+1} \cdot
\end{displaymath} (54)

In the SW case, $ \tilde{s} $ is a growing function of the distance from the Sun, so the above equation describes also the evolution of U versus the solar distance. Starting from the reference field (with $ \tilde{s} =1$, $U= U_{\rm o} $), U decreases with an increasing $ \tilde{s} $ (Fig. 7), and the Bz profile is progressively more restricted around the origin of the function J0. This corresponds to a magnetic field more dominated by the axial field. As expected, this effect is delayed to larger $ \tilde{s} $ values for a more twisted reference field (so with a larger $ F_{\rm ratio}$, Eq. (36), then a larger $ U_{\rm o} $).

For $ \tilde{s} $ decreasing below 1, the opposite behavior is present, i.e. the flux rope is more dominated by its azimuthal field for lower $ \tilde{s} $ values. Then, we found the same behavior as with an ideal MHD evolution (Sect. 3.4). The limiting case, for $ \tilde{s} $ small enough, is $U \approx z_{J_1}$, resulting in a flux rope with a significant reverse Bz flux around the core field. Starting close to the Sun with this type of field, an evolution, that minimizes the energy, progressively removes the reverse Bz flux by reconnection to provide only an axial field, Bz, with a fixed sign within the flux rope (as often observed at and beyond 1 AU in MCs). Later, this relaxation continues and creates an increasingly flatter Bz profile, as well as a smaller $B_{\theta }$ contribution to the field strength.

The evolution of $ \tilde{B}_{\rm c} $ is very limited in this energy minimization case (Fig. 2, top panels), mostly because the field stays a Lundquist's field and so lacks any significant difference between the field strength at the center and at the border of the flux rope. Then, most of the field strength evolution at the center is included in the $( \tilde{ P_{\rm t} } )^{1/2}$ term (Eq. (16)). Indeed, from Eq. (52), we have $\vert J_0(z_{J_1}) \vert < \tilde{B}_{\rm c} < 1$, and at most only a change of a factor ${\approx} 2.5$ is possible for $ \tilde{B}_{\rm c} $ (compare to orders of magnitude for $( \tilde{ P_{\rm t} } )^{1/2}$). Then the energy minimization case follows very well the simple scaling law $B_{\rm c} \approx B_{\rm c,o} ~ \sqrt{ \tilde{ P_{\rm t} } }$ (Fig. 4, top panels).

The evolution of $ \tilde{R} $ with $ \tilde{s} $ is larger than with $ \tilde{B}_{\rm c} $ since $ \tilde{R} $ involves not only fP(U), but also fF(U) that has a broader variation range. $ \tilde{R} ( \tilde{s} )$ has a variation comparable to the ideal-MHD cases analyzed in Sect. 3 (Fig. 2, bottom panels).

As for the ideal-MHD cases, the evolutions of both $B_{\rm c} $ and R with the solar distance D are dominated by the  $ P_{\rm t} (D)$ contributions in Eq. (16). Indeed, Fig. 4 summarizes a broad range of evolutions ranging from ideal-MHD with very different twist profiles to a fully resistive evolution. The evolution of MCs is expected to be in between these cases, so their evolution of both $B_{\rm c} $ and R are expected to be mainly due to  $ P_{\rm t} (D)$. The slight differences between the above models is not expected to appear in observations since these differences can easily be covered by other effects not considered here, such as the non-circular cross section of the flux rope, the dynamic SW pressure, the MC over-pressure, the jump of pressure at shocks, the reconnection with the SW, and the interaction with SW inhomogeneities (in particular overtaking flows).

6 Summary and conclusions

MCs are expanding flux ropes in the interplanetary space with very different expansion properties than the SW, as derived from the observed profiles of plasma velocity at fixed heliodistance, but also from the plasma density and field strength evolution with solar distance. Since the magnetic field energetically dominates the plasma in MCs, the difference in expansion between MC and SW is related to the special magnetic configuration of MCs. Several questions about the main physical mechanisms involved in their expansion remain unanswered in previous studies, e.g. how does their expansion depend on their magnetic configuration? Is an internal over-pressure needed to explain the observed expansion?

We investigate these problems by focusing on the main physical properties. As a local approximation for an elongated field configuration, we consider a straight cylindrical flux rope geometry. Since the plasma $\beta $ is low and the expansion velocity frequently is below the local Alfvén speed in most of the MC volume (Burlaga et al. 1982), we analyze the flux rope evolution as a series of force-free equilibria in pressure balance with the surrounding SW. We show that the presence of shocks or/and an internal over pressure only has a weak effect compared to the variation of the SW total pressure, $ P_{\rm t} (D)$, with solar distance D. We do not include the dynamic pressure of the SW (supposing that the MC and the SW travel at sufficiently comparable velocities). These simplifications permit us to tackle the problem analytically in some cases, and more generally by solving a system of ordinary differential equations (Eqs. (14)-(15)), allowing us a broad survey of the possible magnetic configurations and their evolution with solar distance.

The main dependence of the field strength and radius is derived analytically in Sects. 3 and 5. We show that they mainly depend on the product $L~ P_{\rm t} ^{1/4}$, where L is the length of the flux rope and $ P_{\rm t} $ is the total pressure in the SW. Then, flux ropes with a length evolving proportionally to  $ P_{\rm t} ^{-1/4}$ have a self-similar evolution independently of the structure of their magnetic field (i.e. field strength and spatial scales are simply multiplied by a global factor, function of $ P_{\rm t} (D)$, as D increases).

While a flux rope still attached to the Sun departs significantly from $ P_{\rm t} \propto L^{-4} \propto D^{-4}$ with the observed  $ P_{\rm t} (D)$ dependence in the SW $ P_{\rm t} \propto D^{-2.8}$, we show that these flux ropes still have a nearly self-similar radial expansion. More precisely, their radius is mainly multiplied by $ P_{\rm t} ^{-1/4}$ and their field strength by $ P_{\rm t} ^{1/2}$ as they travel away from the Sun (Sect. 3). This result is obtained after investigating the evolution of the several spatial distributions of magnetic field inside the flux rope. First, we analyze the ideal MHD evolution of a Lundquist's field, as well as extreme twist distributions (uniform or concentrated at the periphery of the flux rope, Sect. 4). Then, we consider the minimization of the magnetic energy with constant helicity (Sect. 5). The strongest departure from these laws, that still remains moderate, is obtained for flux ropes much more twisted than actually observed in MCs (Figs. 2, 4, right panels).

For a flux rope environment given by a magnetized medium described by the Parker spiral of the SW, and values of plasma $\beta $ between 0.1 and 10, we found that the MC core expands slightly faster than its periphery with increasing D. This evolution is due to different constraints present on the axial field (dominant in the core) and on the azimuthal field (dominant at the periphery). This implies a flatter field strength profile with increasing D for both ideal and resistive evolutions. This is a plausible origin of the flattening tendency previously noticed in MCs observed at various solar distances (e.g., Burlaga et al. 1998).

The expansion velocity is found to be almost proportional to the distance from the flux rope axis. While the expansion of the flux rope core (dominated by the axial field) is almost self-similar, the periphery (dominated by the azimuthal field) is expanding at a slightly lower rate for a typical twist value found in MCs. The velocity profile is almost the same for all the cases considered (varying the spatial distribution of the magnetic field and the dissipation law permitted in the MHD system). These results explain the origin of a nearly linear velocity profile observed generically in MCs.

Moreover, the expansion velocity rate found in the models described here is in full agreement with the values observed in MCs (Démoulin et al. 2008). Simulating the crossing of a MC orthogonally to its axis by a fixed spacecraft, we find a plasma velocity equal to $V_{\rm c} - l~V_{\rm c} ^{2} t/D$, where t is the time (=0 at the MC axis), $V_{\rm c} $ is the velocity of the flux rope center, and l is an approximately constant coefficient (i.e. independent of the position within the flux rope and of the magnetic field configuration). We previously found the same functional form by analyzing a set of MCs observed at 1 AU (Démoulin et al. 2008). Furthermore, with a dependence of the SW total pressure as $ P_{\rm t} \propto D^{-n_P}$, we find $l\approx n_P/4$ in present analysis. The typical decay of the SW pressure corresponds to $n_P\approx 2.8$, providing $l\approx 0.7$, in agreement with our previous results $l\approx 0.8 \pm 0.2$ in a set of MCs.

The theoretical profile found, $V_{\rm c} - l~V_{\rm c} ^{2} t/D$, also implies that the velocity profile observed in a MC by a spacecraft, with a given position in space, is strongly dependent on the mean velocity, $V_{\rm c} $ of the observed MC. Indeed, faster MCs need to expand faster to reach an approximate pressure balance with the surrounding SW at a given distance from the Sun. Also, bigger MCs have a larger velocity difference between their front and back simply because they are larger (and the expansion is almost self-similar). Our theoretical results further justify the use of a normalized expansion rate, rather than simply the difference between the front and back velocities, in order to quantify whether a MC is over-expanding or not.

Finally, we found that the magnetic field strength evolves mostly as D-nP/2 and the MC radius as DnP/4. Both dependences predicted under general conditions by our model are in agreement with several statistical studies of MCs and ICMEs observed between 0.3 and 5 AU, taking into account the typical SW pressure dependence on solar distance.

We conclude that MCs and ICMEs are expanding magnetic structures mostly because they encounter a strongly decreasing total pressure in the surrounding SW as they travel away from the Sun. MCs are 3D expanding structures, with approximately a self-similar expansion along and across their axis, but with slightly different rate in these directions (in our model the expansion velocity across their axis is a factor ${\approx} n_P/4$ lower than along their axis). This contrasts with the almost 2D expansion of the surrounding SW. The main differences between the typical expansion profiles of MCs and the SW are due to the different force balances involved. The MC expansion is the combined result of the internal magnetic force balance, its length extension and the pressure evolution at its boundary. However, the volume exponents are not so different in MCs and the SW ( ${\approx} 1+n_P/2 \approx 2.4$ and ${\approx } 2$, respectively). This implies, as observed, only a slightly more rapid decrease of the plasma density with solar distance in MCs than in the SW.

Acknowledgements
We thank Lidia van Driel-Gesztelyi and the referee for reading carefully, and, so improving the manuscript. The authors acknowledge financial support from ECOS-Sud through their cooperative science program (No. A08U01). This work was partially supported by the Argentinean grants: UBACyT X425, PICTs 2005-03-33370 and 2007-03-00856 (ANPCyT). S.D. is member of the Carrera del Investigador Científico, CONICET.

References

All Tables

Table 1:   Definition of the main non-dimensional variables.

Table 2:   Summary of results on MCs and ICMEs.

All Figures

  \begin{figure}
\IfFileExists{color_figures.txt}{
{\includegraphics[width=0.5\tex...
...}{
{\includegraphics[width=0.5\textwidth, clip=]{0971f1c.eps} }
}\end{figure} Figure 1:

Schematic of a flux rope ejected from the Sun (the apex of the central field line is at a distance D from the Sun). The magnetic configuration of a magnetic cloud (MC) is well represented by an elongated twisted flux tube mostly anchored to the Sun. The schematic on the left shows the locally straight approximation used in this work (A color version is available in the electronic version).

Open with DEXTER
In the text

  \begin{figure}
\IfFileExists{color_figures.txt}{
{\includegraphics[width=0.95\te...
...{
{\includegraphics[width=0.95\textwidth, clip=]{0971f2c.eps} }
}\end{figure} Figure 2:

Log-log plots of the normalized central magnetic field strength, $ \tilde{B}_{\rm c} $, and radius $ \tilde{R} $ defined by Eq. (16) in function of the size factor, $ \tilde{s} $ defined by Eq. (19). $ \tilde{B}_{\rm c} $ is defined as the relative evolution of the central to the border field strength normalized to the reference field ( $ \tilde{B}_{\rm c} = \frac{B(0)}{B(R)}\frac{ B_{\rm o} ( R_{\rm o} )}{ B_{\rm o} (0)}$). Three cases suppose an ideal MHD evolution with reference fields (at $ \tilde{s} =1$) being the Gold & Hoyle (GH, Eq. (22)), the ``split'' (SP, Eq. (30)), and the linear force-free field (LU, Eq. (35)). The fourth case (ME, Eq. (45)) always has the magnetic energy to keep the magnetic helicity preserved. The ratio of the azimuthal to the axial flux, $ F_{\rm ratio}$ (computed with a length $ L_{\rm o} = 2\pi R_{\rm o} $) is typically ${\approx } 2$ in MCs ( left panels), while $ F_{\rm ratio}= 5$ illustrates a very twisted flux rope ( right panels). (A color version is available in the electronic version).

Open with DEXTER
In the text

  \begin{figure}
\IfFileExists{color_figures.txt}{
{\includegraphics[width=0.45\te...
...
{\includegraphics[width=0.45\textwidth, clip=]{0971f3bc.eps} }
}\end{figure} Figure 3:

Log-log plots showing the evolution of the SW total pressure (Eq. (39)) normalized at 1 AU, and a power-law approximation (Eq. (40) with nP=2.8) as a function of the solar distance D (in AU). $\beta _{\rm o}$ is the plasma $\beta $ in the SW at 1 AU, and nT define the SW temperature dependence with D ( $T \propto D^{-n_T}$). In both panels $ \tilde{V} =1$, a typical value for the slow SW at 1 AU and low latitude (A color version is available in the electronic version).

Open with DEXTER
In the text

  \begin{figure}
\IfFileExists{color_figures.txt}{
{\includegraphics[width=0.90\te...
...{
{\includegraphics[width=0.90\textwidth, clip=]{0971f4c.eps} }
}\end{figure} Figure 4:

Log-log plots of the central magnetic field strength, $B_{\rm c} $, and radius, R, as a function of the distance from the Sun, D. We set the typical values $B_{\rm c} = 20$ nT and R = 0.1 AU at D=1 AU, and nP=2.8 (Eq. (40), Fig. 3). Three cases, the Gold and Hoyle (GH, Eq. (22)), ``split'' (SP, Eq. (30)), and linear force-free field (LU, Eq. (35)), suppose an ideal-MHD evolution with reference fields defined at $D= D_{\rm o} $. The fourth case (ME) always has the minimum magnetic energy keeping the magnetic helicity preserved. The straight lines labeled ``scaling'' have the slopes -nP/4 and nP/2 for $B_{\rm c} $ and R, respectively (Eq. (42)). The ratio of the azimuthal to the axial flux, $ F_{\rm ratio}$, is computed with a length $ L_{\rm o} = 2\pi R_{\rm o} $ (A color version is available in the electronic version).

Open with DEXTER
In the text

  \begin{figure}
\IfFileExists{color_figures.txt}{
{\includegraphics[width=\textwi...
... }
}{
{\includegraphics[width=\textwidth, clip=]{0971f5c.eps} }
}
\end{figure} Figure 5:

Evolution of the field strength, B, the axial, Bz, and azimuthal, $B_{\theta }$, field components with the distance D (D = (0.1, 1, 10) AU in the top, middle, bottom panels, respectively). The reference field is taken at $D=1~\rm AU$ with the typical value $B_{\rm c} = 20$ nT and R=0.1 AU. The same four cases as in Fig. 4 with $ F_{\rm ratio}=2$ are shown (A color version is available in the electronic version).

Open with DEXTER
In the text

  \begin{figure}
\IfFileExists{color_figures.txt}{
{\includegraphics[width=\textwi...
... }
}{
{\includegraphics[width=\textwidth, clip=]{0971f6c.eps} }
}\end{figure} Figure 6:

Evolution of a flux rope for the three ideal-MHD cases shown in Fig. 5 (D = (0.1, 1, 10) AU in the top, middle, bottom panels, respectively). r is the radial coordinate of the flux rope, with $ r_{\rm o} $ being its value in the reference field taken at $ D_{\rm o} =1$ AU (with the flux rope radius $ R_{\rm o} =0.1$ AU). v is the radial plasma velocity. v D /(r Vc) is the normalized expansion velocity (Sect. 4.5). We set $ F_{\rm ratio}=2$, nP=2.8, and $V_{\rm c} =600$ km -1 ($V_{\rm c} $ is the global velocity of the flux rope away from the Sun). The horizontal line in the right panels has an ordinate of nP/4 (A color version is available in the electronic version).

Open with DEXTER
In the text

  \begin{figure}
\IfFileExists{color_figures.txt} {
{\includegraphics[width=0.47\t...
...
{\includegraphics[width=0.47\textwidth, clip=]{0971f7c.eps} }
}
\end{figure} Figure 7:

The evolution of $U=\alpha R$ in function of $\log_{10}( \tilde{s} )$ (Eq. (19)) for different ratios of the azimuthal to axial flux ( $ F_{\rm ratio}$, Eq. (36)). The magnetic energy is supposed to be minimized with a preserved magnetic helicity, so the field stays a linear force-free field. The thick curves show the solutions of the first equation in Eq. (52), while the thin curves are for its approximation (Eq. (54)). (A color version is available in the electronic version).

Open with DEXTER
In the text


Copyright ESO 2009

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.