A&A 413, 403-414 (2004)
DOI: 10.1051/0004-6361:20034008
D. Moss^{1} - A. Shukurov^{2}
1 - Department of Mathematics, University of Manchester, Manchester M13 9PL, UK
2 - School of of Mathematics and Statistics, University of Newcastle,
Newcastle upon Tyne NE1 7RU, UK
Received 26 June 2003 / Accepted 3 October 2003
Abstract
We show that a relatively weak poloidal
external magnetic field,
imposed in the outer parts of a Keplerian accretion disc,
can open magnetic lines of a much stronger magnetic field
generated by the dynamo in the inner disc. The resulting
magnetic configuration can be favourable for launching
a centrifugally driven wind.
Even for a relatively weak external magnetic field (with energy
density in excess of about only 10% of the thermal energy density at the outer
disc radius)
the geometry of the poloidal field in the disc neigbourhood is
almost independent of the dynamo action, and is determined by
the external field.
The radial profile of the poloidal magnetic field on the disc surface
is similar to that in the self-similar solution of Blandford & Payne (1982).
We conclude that poloidal fields resulting from dynamo action can be important
for launching a magneto-centrifugally driven outflow in
accretion discs that occur in weakly magnetized environments,
i.e. where an external field is too weak to be important.
However, even a relatively weak external field
can be wound up by the
differential rotation to give a much stronger azimuthal field,
and so modify dynamo action in the disc.
Key words: accretion, accretion disks - magnetic fields - magnetohydrodynamics (MHD) - turbulence
Accretion discs and the associated winds and jets are expected to form and to evolve in a magnetic environment. As reviewed, e.g., by Königl & Ruden (1993) and Königl & Pudritz (2000), magnetic fields are essential for the accretion to occur, both as they directly brake the rotation via the Maxwell stresses and by driving an outflow (Bisnovatyi-Kogan & Ruzmaikin 1976; see, e.g., Pelletier & Pudritz 1992; Ferreira & Pelletier 1995; Casse & Ferreira 2000, and references therein), and as they generate turbulence in the disc (Balbus & Hawley 1998) to maintain viscous transport of angular momentum.
Accretion discs are believed to be a site of vigorous dynamo action (Pudritz 1981a,1981b; Stepinski & Levy 1991; Brandenburg et al. 1995), but virtually all models of magnetized accretion discs and associated outflows rely on externally supported magnetic fields, mainly because dynamos produce magnetic fields that close within the immediate vicinity of the disc, so the closed magnetic configuration cannot support an outflow extending to the observed very large distances (Pelletier & Pudritz 1992; Sect. V.B in Königl & Pudritz 2000). However, the numerical simulations of von Rekowski et al. (2003) show that a dynamo-generated magnetic field can support a vigorous outflow, provided the field can be opened up by, say, a thermally driven wind. Campbell (1999, 2000, 2001) and Campbell & Caunt (1999) discuss accretion disc models and outflows based on dynamo-generated magnetic fields using a simple, analytically solvable, one-dimensional dynamo model where the disc field matches a potential field in the disc corona; however, these authors do not discuss the global geometry of the magnetic field and how it can be opened at large distances from the disc.
It has been suggested by Reyes-Ruiz & Stepinski (1997) and Reyes-Ruiz (2000) that an external magnetic field can open magnetic lines of a dynamo-generated field to make the magnetic configuration suitable for the magneto-centrifugal wind launching (see also Moss & Shukurov 2001). We reconsider this idea here, with a fully nonlinear treatment of the dynamo action in the presence of an external magnetic field and estimate the external field required to produce the desired magnetic configuration.
A mean-field dynamo model for an accretion disc embedded into an external field is presented in Sect. 2, and we show in Sect. 3 that a rather weak external magnetic field can open up most or all the magnetic lines produced by the dynamo. We show in Sect. 3 that our results are insensitive to the details of the mean-field parameterization of the dynamo action and to the detailed distribution of the (turbulent) magnetic diffusivity in the disc and its environment. Our results are put into physical context in Sect. 4 where, in particular, we argue that external magnetic field strength required to open the dynamo field in the disc of a young stellar object is close to that of a field that can be trapped within a few tens of AU in the core of a molecular cloud.
Our interest here is with magnetic fields that can be involved in launching large-scale outflows from an accretion disc, and so their scale must be comparable to the disc size and, in particular, be much larger than the disc's turbulent scale. Therefore, a mean-field description of the accretion disc is appropriate where space-averaged quantities are involved - e.g. the standard model of Shakura & Sunyaev (1973). However, the averaged induction equation contains not only terms responsible for turbulent magnetic diffusion and induction by differential rotation, but also a term related to the vectorial nature of magnetic field, known as the -term (not to be confused with the "alpha-coefficient'' of the Shakura & Sunyaev theory which we denote as here). As similarly argued by Blackman (2001), a consistent mean-field description of thin accretion discs must also include the hydromagnetic -effect that quantifies the mean-field dynamo action.
We study a dynamo active disc embedded in a passive, conducting corona (halo),
and solve the standard
dynamo equation in the disc
(e.g., Moffatt 1978; Zeldovich et al. 1983)
The coefficient
in Eq. (1) must be an odd function
of z because it arises from the interaction of the Coriolis force
with turbulence in stratified medium,
(Parker 1979) or with
magnetic buoyancy (Brandenburg 1997).
As discussed by e.g. Ruzmaikin et al. (1988), solutions of
Eq. (1) only weakly depend on the particular choice of the
functional dependence of
on z. Solutions with independent of the magnetic field grow exponentially in time provided
the magnetic Reynolds numbers based on
and ,
as
defined in Eqs. (7) and (8) below, are large
enough. The growth will saturate as soon as the Lorentz force
becomes significant; the saturation is often modelled with an
algebraic dependence of
on magnetic field strength such that
decreases as B grows. Thus we adopt an -effect
described by
We consider only axisymmetric solutions, and so write
(3) |
We take the physical state of the disc to be that of a conventional
viscous
accretion disc ("-disc'', e.g., Frank et al. 1985), and for convenience
we reproduce the solution here, in standard notation, for the
disc scale height
We wish to impose a uniform vertical field at large distances above and below
the disc plane, and so there apply the boundary conditions
We consider solutions with fixed parity with respect to the plane z=0(in fact, we only consider dipole-like fields, see Sect. 2.5), and so calculations are performed in the region only, with suitable symmetry conditions on z=0.
The choice of the inner and outer radii of the disc is important. We need to choose the
outer radius R large enough that the important features of the
system are captured.
On the other hand, angular velocity
gradients, and so dynamo generation, are strongest in the inner regions.
But it is computationally difficult to handle a disc in which the ratio
is too small.
After some experimentation, we chose
The radial boundary conditions are
Equation (1) is nondimensionalized in terms of the length R and
time
.
Thus in dimensionless units,
,
.
Angular velocity, radial velocity and
are measured in
units of their values ,
u_{0} and
at r=R, and
is measured in units of
at r=R; the
corresponding unit of
is RB_{0}. Taking k=1 and
for definiteness, this corresponds to
in
models with
,
but
in
Sect. 3.4 where
.
We assume
to be independent
of z in
;
its value in |z|> h(r) is irrelevant, as there and
only enters via the alpha-quenching term in
Eq. (2). We take as standard values
For our initial exploratory investigations (where we took , f_{R}=1) we used the axisymmetric dynamo code written in cylindrical polar coordinates described in Moss & Shukurov (2001). For computations with , f_{R}>1 we modified this code to use as the radial variable. Having verified that odd parity solutions were strongly preferred for , we sought solutions in the half-space , and also changed to using as the independent variable rather than z, solving in , with usually. Thus the boundary conditions to select the odd parity solutions are applied at rather than at z=0; comparison with the version of the code using z as variable and showed that this inaccuracy was unimportant for our investigation.
For the revised version of the code we used a semi-implicit Dufort-Frankel scheme for the time stepping, rather than the Runge-Kutta integrator used previously, as the restrictions on the time step were then less stringent. Our standard computational grid had 101 points uniformly distributed in and 201 points in .
Table 1: of results with magnetic diffusivity independent of z(i.e., ), , , R_{u}=0.64, , , and . and are the field energies in the whole computational domain and near the disc respectively (see text), s=E_{P}/E_{T}is the ratio of energies in the poloidal and toroidal magnetic field in the whole domain, and F_{0} is the signed flux through the disc plane in . The final entry (asterisked) has , and refers to a purely advective calculation, with no dynamo action.
We describe firstly a "standard'' set of computations, with , , , giving , , R_{u}=0.64, . The magnetic diffusivity is taken to be independent of z, i.e. the form (5) holds at all z (), and increases from zero to "saturation'' values. A summary of these results is presented in Table 1. There is the total magnetic energy in the computational domain, is the magnetic energy within the disc, , is the signed flux of (vertical) magnetic field through the disc mid-plane (), and s is the ratio of poloidal to azimuthal field energies. Note that s is always small, and that the energy is predominantly in the toroidal part of the magnetic field, even though it is the poloidal field geometry that shows the most obvious changes as increases. (Of course, generation of toroidal field from poloidal by the winding action of the differential rotation is independent of the source of the poloidal field.)
Figure 1: The dependence of the near-disc energy (solid) and poloidal flux through the disc plane F_{0} (dashed) on for the dynamo solutions with independent of z(see also Table 1). The dotted lines show the same quantities for the purely advective solutions. | |
Open with DEXTER |
In Fig. 1 we plot F_{0} and as a function of the external field strength . This shows a smooth transition between pure dynamo solutions in the limit , where the disc magnetic field is independent of , and purely advective solutions for large , where the dynamo action is suppressed by the contribution of the external field to -quenching as described by Eq. (2). For example, the magnetic energy and total flux are always equal to or in excess of the greater of the values for the pure dynamo and pure advection fields. For , the solutions are essentially purely advective, and scale with the asterisked entry in Table 1 as , .
The geometry of the solutions is of interest. We first display, in Fig. 2, an intermediate case with . There is a region where field lines close within the disc, and the field generally resembles a typical disc dynamo field. Further from the axis () the imposed field is relatively unaffected by the dynamo action in the disc. Solutions for several other values of are illustrated in Fig. 3. For the pure dynamo solution, the poloidal field lines close within, or not far outside of, the disc region . Essentially no field lines get as far as the upper boundary (Fig. 3a). For , a substantial proportion of the field lines passing through z=0 reverse their curvature above the disc and continue to large values of z, see Figs. 2 and 3b. Comparison between the purely advective solution, shown in Fig. 3d, and Fig. 3c shows that the dynamo is essentially irrelevant for the poloidal field for , as the poloidal field is then dominated by the imposed (poloidal) external field. However, the azimuthal field is essentially that of the pure dynamo, unless .
We further verified that in the absence of advection (R_{u}=0), for larger values of there is no accumulation of poloidal field lines towards the inner parts of the disc, and the poloidal field is near to vertical everywhere. In other words, the bunching together of poloidal field lines when is a purely advective effect rather than a result of any dynamo action (which is suppressed for such a strong external field).
Figure 2: Poloidal field lines for the solution with , , , R_{u}=0.64, , independent of z, as represented in Table 1 Here and in other similar figures below, the solution at z<0 is obtained from that shown by the symmetry, . Labels on the curves refer to the magnitude of the poloidal flux function , and are equally spaced at intervals of 0.01 from 0.01 to 0.07. | |
Open with DEXTER |
We then sought to increase the effects of advection by taking (i.e. reducing while keeping unchanged); then R_{u}=6.4. But dynamo action is then also enhanced, with , . This was computationally more demanding, and we only calculated a few examples. For numerical reasons we here put . The essential features of the solutions are unaltered. We show in Fig. 4 the poloidal field lines with and 10. The only significant change from the models with is that the transition from dynamo dominated to external field dominated solutions appears to occur at rather larger values of (by a factor estimated as between 3 and 10), consistent with the stronger dynamo action.
We also made some exploratory calculations with a nonuniform magnetic
field imposed at the boundary z=Z,
corresponding to a concentration of the external magnetic flux
towards the axis having been produced during collapse. We took
two forms of this boundary condition:
Figure 3: As in Fig. 2, but for several values of the external magnetic field: a) - pure dynamo solution; b) ; c) . Panel d) shows the purely advective solution, . Contours are equally spaced as a function of poloidal flux , between levels 0.007 and 0.064 in a) and between 0.01 and 0.023 in b); and are logarithmically spaced (factor of 2), between levels 0.01 and 0.64 in c) and between 0.12 and 61.4 in d). | |
Open with DEXTER |
The models discussed above have a magnetic diffusivity that is independent of height. The value of the diffusivity above an accretion disc is quite uncertain. Thus we also studied solutions with the coronal magnetic diffusivity both enhanced, by a factor , and reduced, . A summary is presented in Table 2, for the same dynamo parameters as used in Table 1. Again, dynamo action can be seen to influence solutions whilst . For larger values of , the solutions approach the appropriate purely advective solutions (asterisked entries in Table 2, suitably scaled). Typical poloidal field configurations, with and 10, are shown in Fig. 5. Overall, the results are very similar to those discussed above with no z-dependence of , although the poloidal field lines are somewhat more strongly bent in a more highly conducting halo () - as naively might be expected.
As can be seen in Table 2, the results for magnetic diffusivity independent of height (, Table 1) and enhanced coronal diffusivity () are similar even quantitatively. However, a disc surrounded by a better conducting corona ( ) maintains a larger vertical magnetic flux and azimuthal field strength, presumably because of the dependence of the dynamo on boundary effects at the disc-corona interface.
Figure 4: Poloidal field lines for models with increased magnetic Prandtl number, and a range of external magnetic field strengths: a) ; b) ; c) . The other parameters are , , R_{u}=6.4, and . Contours are equally spaced as a function of poloidal flux . | |
Open with DEXTER |
We also found for standard choices of parameters, solutions for which the external field is excluded from the inner part of the disc, that is the external field is in the opposite sense to the dynamo dominated field close to the axis; of course, such configurations occur only when the external field and the seed field for the dynamo are independent. In this case magnetic fields generated by the dynamo cannot merge with the external field and the vicinity of the disc is magnetically isolated from the surroundings. The external field bends to "flow around'' the disc region so that the disc is surrounded by a closed current sheet. The net signed magnetic flux of B_{z} within radius r in the disc plane changes sign at some radius , i.e. is the smallest radius at which the external field passes through the plane z=0. In other words, is the radius in the plane z=0 of the magnetic cocoon surrounding the disc. Figure 6 shows the poloidal field geometry for solutions of this type for several values of . As increases, decreases (e.g. for respectively), and for the solutions are almost indistinguishable (apart from sign of the field) from the corresponding solutions of Table 1. (There remains a small reduction in the value of .)
Clearly, neither in configurations with a closed magnetosphere (as Fig. 6a), nor in those with an externally dominated field, that bends away from the z-axis in the vicinity of the disc plane (Fig. 6b), will wind launching occur.
In the limit of large external field, the field configuration is identical with that in the "standard'' case, being now independent of the relative directions of the dynamo and external fields, cf. Figs. 3d and 6c. We note that the influence of the external field dominates over that of the dynamo generated field at approximately the same strength as in the previous set of simulations where it is directed in the same sense as the dynamo field.
Figure 5: Poloidal field lines for models with magnetic diffusivity either increased () or reduced () in the corona, for various values of : a) , ; b) , ; c) , ; and d) , . The other parameters are , , R_{u}=0.64 and . Contours are equally spaced as a function of poloidal flux . | |
Open with DEXTER |
Table 2: Table 1, but with enhanced and reduced coronal diffusivities ( and 0.3). The asterisked entries have , and refer to purely advective calculations, with no dynamo action.
We recognize the possibility that the disc at large radii may be sufficiently
dense and cool that ionization even by cosmic rays would be insufficient
to justify treating it as a conducting medium.
We can expect that this would affect the nature of
the dynamo solutions, since the local dynamo number
When presented in dimensionless form, quantitative results obtained for the smaller disc, as shown in Table 3, are not dissimilar to those found with the standard larger disc, as given in Table 1. More qualitatively, the shapes of the magnetic field lines are very similar for the same dimensionless , so we do not show them for the smaller disc. In terms of dimensional quantities, this means that a stronger external field is required to open up the dynamo-generated field in a smaller disc. The difference is a factor of 420 if the disc radii differ by a factor of 100. As we argue in Sect. 4, this feature is related to the global nature of the dynamo solutions.
A peculiar feature of the standard accretion disc model is that the dynamo number does not vary much with radius because a decrease in the rotational shear is compensated by the disc flaring (e.g. Pudritz 1981b). As discussed in Sects. 2.2 and 2.4, the dynamo number is estimated as . The dynamo action can maintain a non-decaying dipolar magnetic field in a disc surrounded by vacuum if -70 (Parker 1979; Soward 2003), i.e., for for . Therefore, the growth rate of the magnetic field remains positive at all radii if is small enough. In this sense, the dynamo region in an accretion disc has no natural outer radius. An effective outer radius can be identified with that at which the steady-state magnetic field becomes negligible or comparable with an external magnetic field. Thus, the effective radius of the dynamo region may depend on the external magnetic field - a rather unusual situation.
Figure 6: Poloidal field lines for models with excluded external field, with various values of the latter: a) ; b) ; c) . Solid contours indicate field lines which originate in the dynamo-dominated part of the disc (), and do not connect with the external field, broken contours show field lines that connect with the external field (). The other parameters are , , R_{u}=0.64, and . Contours are equally spaced as a function of poloidal flux . | |
Open with DEXTER |
Table 3: of results with a smaller disc ( ), else as for Table 1. F_{0} is the signed flux through the disc plane in (corresponding to 10^{14} cm) and is the dimensionless strength of the external field, but now measured in the units of at that radius.
As discussed in Sect. 3.4, a stronger external magnetic field is required to open up poloidal magnetic lines produced by the dynamo in a disc of smaller radius. This can be understood as follows. Large-scale magnetic fields generated by the dynamo have a global nature: since most of magnetic lines are closed within the disc, a change in the size of the dynamo system leads to a change in the magnetic field throughout the whole region. The dynamo action is most efficient near the axis, but magnetic lines that pass through that region cross the disc midplane again near the outer boundary of the disc. This explains why a relatively weak external magnetic field can open up the strong field in the inner parts of the disc: magnetic lines that reconnect with the external field at large radii are anchored near the disc axis; automatically the lever arm of such open field lines will be large. This feature makes dynamo solutions slightly modified by an external field especially favourable for wind launching. We show in Fig. 7 the opening angle of magnetic lines at the disc surface, for the model presented in Figs. 2 and 3. For a cold flow, the magneto-centrifugal wind can be launched if (Blandford & Payne 1982); pressure effects (Pelletier & Pudritz 1992) and finite disc thickness (Campbell 1999) result in slightly larger critical values of the launching angle.
Figure 7: The opening angle of magnetic lines, at the disc surface z=h(r)(the "launch angle'') for the model presented in Figs. 2 and 3, with a range of external field strengths, (solid), 0.1 (dashed) and 1.0 (dotted). The values favourable for launching a magneto-centrifugal wind are in the range . | |
Open with DEXTER |
Figure 8: The variation of the function F(r) (the flux passing through the disc plane within radius r) with radius for the model presented in Figs. 2 and 3, with external field strengths (solid), 0.1 (long-dashed), 1.0 (short-dashed) and 10.0 (no dynamo action, dotted). The approximately linear parts of these plots through the bulk of the disc region correspond approximately to power laws r^{0.66}, r^{0.73}, r^{0.92}, r^{1.0}, respectively. | |
Open with DEXTER |
Apart from having an appropriate opening angle, the magnetic field of the disc must satisfy several other requirements in order to support a magnetically driven outflow. For example, the vertical current must not diverge at in order to avoid instabilities that would disrupt the outflow (Pelletier & Pudritz 1992). The radial variation of the azimuthal magnetic field far from the axis can be expected to (and, in our model, does) follow closely that of the equipartition field - this expectation arises from the form of the nonlinearity adopted, Eq. (2). This implies that away from the axis, and so there is a physically reasonable reduction of the vertical current at large r. The axisymmetric dynamo solutions discussed here of course have for , keeping the current finite at the axis.
A constraint on the radial distribution of poloidal magnetic flux in a steady-state accretion disc, with angular momentum transport by the outflow, is discussed in detail by Pelletier & Pudritz (1992) whose solutions have (with ) where is a constant and all quantities are taken at z=0. Depending on the value of , the flow at large z collimates to a cylindrical shell if or diverges if ; the case corresponds to the self-similar solution of Blandford & Payne (1982). From the requirement that the terminal outflow speed decreases with radius, Pelletier & Pudritz (1992) obtain ; these authors favour the value because it produces and a vertical current that remains independent of radius.
Figure 9: The variation of the toroidal magnetic field along open poloidal magnetic field lines as a function of height z. a) The model with along field lines with values of the flux function (solid), 0.02 (broken) and 0.04 (dotted) - see Fig. 2 where magnetic field lines are labelled with the values of ; these field lines intercept the disc plane z=0 at radii r = 0.0059, 0.13, 0.25 respectively. b) The model with , see Fig. 3c, along the field lines with (solid), 0.04 (broken) and 0.10 (dotted), which meet the plane z = 0 at r = 0.049, 0.13, 0.43respectively. In each case the position at which the field line meets the disc surface is indicated by a square. | |
Open with DEXTER |
Kinematic (linear) dynamo solutions of dipolar parity in a thin disc with aspect ratio have the asymptotic scaling for (Soward 1978, 1992, 2003), where a(r) is a certain slowly varying function of radius. Assuming that this scaling applies to the (nonlinear) steady-state solutions as well, and that all the variables can be replaced by their values at the local radius r(in particular, with replaced by ), we obtain for the scalings adopted in Sect. 2.2, which gives . We show in Fig. 8 the radial distribution of the the flux function F(r) for models illustrated in Figs. 2 and 3. The solutions where the external field is not dominant, i.e. , are in excellent agreement with the asymptotic results, given that a(r)varies insignificantly with r away from the boundaries. The corresponding value of , in the radial range where the power-law approximation is meaningful, is , a value remarkably similar to that in the model of Blandford & Payne (1982). The agreement is especially close for external fields of moderate strength, -0.1. Figure 10 displays the dependence on radius of the field components B_{r}, B_{z} and at the disc surface z=h(r).
A notable difference of the dynamo solutions with a conducting corona around the disc from the magnetic configurations employed in modelling disc outflows is that not only within the dynamo active disc but also in its immediate vicinity (see Fig. 10, and below in this section). We show in Fig. 9 the variation of the toroidal magnetic field along selected open poloidal field lines. A strong gradient of magnetic pressure that can contribute to driving an outflow is evident in these solutions. If efficient, such magnetic pressure driving can relax the limitation on the launch angle discussed above. We note, however, that pressure-driven winds require that the total pressure gradient is maintained by effects whose time scale is shorter than the sound crossing time across the disc, . But the dynamo regeneration time can hardly be shorter than . Indeed, must be longer than both the rotation period and the turbulent eddy turnover time, , where is the turbulent scale and . Therefore, the total pressure will remain equilibrated in spite of the growth of the magnetic field due to dynamo action.
According to our numerical solutions,
the external field affects the disc dynamo significantly if
The strength of the external magnetic field around an accretion disc is not well known from either theory or observations. An extrapolation of magnetic field strengths observed in molecular clouds (Crutcher 1999) to a number density of yields a field of about 10^{-2} G which would be sufficient to affect the dynamo in a disc of radius about cm or larger. However, only an uncertain fraction of this field will be retained in the collapsing gas that eventually forms an accretion disc (e.g. Desch & Mouschouvias 2001; Nakano et al. 2002). Most discussions of wind launching from accretion discs use external magnetic fields of strength equal to without any detailed justification.
A useful estimate of the strength of the dipolar poloidal magnetic field
within a
dynamo active disc and in its vicinity can be obtained from asymptotic dynamo
solutions of Soward (1978, 1992, 2003) for a disc surrounded
by vacuum. From the requirement that the production rate of the poloidal
magnetic field by the -effect balances turbulent diffusion, we obtain
,
which yields an
order of magnitude estimate
,
where the factor
arises because the lowest-order term in the
expansion of
in
is independent of z for solutions of dipolar symmetry (see Soward 2003
for a review). With
,
where Lis the radial scale of the eigenmode, we obtain (Soward, private
communication)
To obtain a
tentative estimate for a steady-state, nonlinear solution, it is often assumed
that nonlinear magnetic effects reduce the dynamo number to its critical value
to saturate the dynamo action, but the structure of the solution remains
similar to that of the marginally stable eigenmode. Then one can use
in the steady state, assuming that
magnetic field is
still too weak to affect the differential rotation, and so
remains
unchanged. With
,
we obtain
Figure 10: The variation as a function of radius of the unsigned field components B_{r} (dotted), B_{z} (solid) and , dashed) at the disc boundary, for the model with , compare Fig. 2. In this solution, B_{r}, B_{z}>0, . The broken straight line illustrates the dependence r^{-5/4}, corresponding to that of the fields in the Blandford & Payne (1982) model - see discussion in the text. | |
Open with DEXTER |
The estimates Eqs (12) and (13) have certain
limitations, and
Fig. 10 shows that
at the disc surface in our
solutions. Firstly, Eqs. (12) and (13) are valid for a
disc surrounded by vacuum. Our numerical solutions
with a conducting corona indicate that the critical value of the dynamo number
is significantly (by a factor of about 3) smaller than for a disc surrounded by vacuum.
Correspondingly, the ratio
is smaller by the same factor.
Secondly, and perhaps more importantly, nonlinear effects can affect the
solution more strongly than is assumed in Eq. (13).
It can reasonably be expected
that the radial scale of the nonlinear solution is comparable to the disc radius,
i.e.
,
which would yield, instead of Eq. (13),
Altogether, it appears that magnetic field generated by the mean-field dynamo in accretions discs have properties favourable for launching MHD outflows. However, our dynamo solutions are based on an accretion disc structure in which angular momentum is carried away by viscous stress. Magnetic braking is plausibly a stronger mechanism of angular momentum transfer (e.g. Pelletier & Pudritz 1992). Therefore, the radial velocity can be larger than that used here; moreover, it can be a function of magnetic field. As discussed by Moss et al. (2000), stronger inflow inhibits the dynamo action and enhances the dragging of the external magnetic field. However, this effect should be rather weak under realistic conditions where the turbulent magnetic Reynolds number based on the inflow speed is less than or of order unity. Furthermore, we hope that the qualitative features of the magnetic field are sufficiently robust - mean-field dynamos are in general remarkably insensitive to the details of the model.
If the magnetic field in the disc is produced by the dynamo without being significantly affected by the external field, wind launching becomes problematic. However, pressure-driven winds can open the magnetic lines produced by the dynamo and thus help to launch a magneto-centrifugal wind (cf. von Rekowski et al. 2003). On the other hand, a weak external magnetic field opens up the magnetic configuration produced by the dynamo to create conditions favourable for magneto-centrifugal wind launching. An even stronger external magnetic field, capable of supporting strong enough accretion by contributing to angular momentum transfer, may not need any help from the dynamo action if it is bent strongly enough by the inflow.
We have shown that a relatively weak external poloidal magnetic field can open up magnetic lines generated by the mean-field dynamo in a standard accretion disc. The resulting geometry of poloidal magnetic lines is favourable for magneto-centrifugally driven outflows, with magnetic lines starting near the axis opened first (i.e., for smaller ). The structure of the poloidal field lines is very sensitive to the external magnetic field strength , and becomes practically independent of the dynamo if exceeds about 10% of the value corresponding to equipartition with thermal energy at the outer disc radius, denoted B_{0}. Even a relatively weak external poloidal magnetic field can be wound up by the differential rotation in the disc to give a strong azimuthal field, which in turn can dominate the dynamo generated field. There is a small range of , (somewhat depending on the parameters chosen - see Sect. 3) where the field is a true hybrid of the dynamo and advected fields; these fields have both suitable launching angles and field lines extending to large z. The radial scaling of the magnetic field for is very close to that of Blandford & Payne's (1982) self-similar solution. However, the range of values of in physical units where these hybrid solutions can be found depends crucially on the effective radius of the accretion disc, which is not a priori a well known quantity. Nevertheless, for plausible values of the disc radius, these fields are not exceptionally strong in the context of current models of the star formation process.
Acknowledgements
We are grateful to M. Camenzind, R. Pudritz and A. M. Soward for useful discussions and to an anonymous referee for careful reading of the manuscript and very helpful suggestions. We acknowledge support from PPARC (Grant PPA/G/S/2000/00528) and NATO (Grant PST.CLG 974737).