A&A 477, 25-34 (2008)
P. Hennebelle1 - R. Teyssier2
1 - Laboratoire de radioastronomie millimétrique, UMR 8112 du CNRS, École normale supérieure et Observatoire de Paris, 24 rue Lhomond, 75231 Paris Cedex 05, France
2 - Service d'Astrophysique, CEA/DSM/DAPNIA/SAp, Centres d'Études de Saclay, 91191 Gif-sur-Yvette Cedex, France
Received 18 July 2007 / Accepted 17 September 2007
Context. A large fraction of stars are found in binary systems. It is therefore important for our understanding of the star formation process, to investigate the fragmentation of dense molecular cores.
Aims. We study the influence of the magnetic field, ideally coupled to the gas, on the fragmentation in multiple systems of collapsing cores.
Methods. We present high resolution numerical simulations performed with the RAMSES MHD code starting with a uniform sphere in solid body rotation and a uniform magnetic field parallel to the rotation axis. We pay particular attention to the strength of the magnetic field and interpret the results using the analysis presented in a companion paper.
Results. The results depend much on the amplitude, A, of the perturbations seeded initially. For a low amplitude, A=0.1, we find that for values of the mass-to-flux over critical mass-to-flux ratio, , as high as , the centrifugally supported disk which fragments in the hydrodynamical case is stabilized and remains axisymmetric. Detailed investigations reveal that this is due to the rapid growth of the toroidal magnetic field induced by the differential motions within the disk. For values of smaller than 5, corresponding to higher magnetic intensities, there is no centrifugally supported disk because of magnetic braking. When the amplitude of the perturbation is equal to A=0.5, each initial peak develops independently and the core fragments for a large range of . Only for values of close to 1 is the magnetic field able to prevent the fragmentation.
Conclusions. Since a large fraction of stars are binaries, the results of low magnetic intensities preventing the fragmentation in the case of weak perturbations is problematic. We discuss three possible mechanisms which could lead to the formation of binary systems, namely the presence of high amplitude fluctuations in the core initially, ambipolar diffusion and fragmentation during the second collapse.
Key words: magnetohydrodynamics (MHD) - instabilities - ISM: kinematics and dynamics - ISM: structure - ISM: clouds
Understanding the fragmentation of collapsing prestellar dense core, is of great importance to our understanding of the star formation process. In particular, determining the number of fragments, their masses and their orbital characteristics are fundamental and rather challenging problems. A significant fraction of stars are found in binary or multiple systems (Duquenoy & Mayor 1991; although see Lada 2006). Many studies have investigated this issue in the context of hydrodynamical calculations (e.g. Miyama 1992; Boss 1993; Bonnell 1994; Truelove et al. 1998; Bate & Burkert 1997; Bodenheimer et al. 2000; Matsumoto & Hanawa 2003; Hennebelle et al. 2004; Goodwin et al. 2004; Banerjee et al. 2004). Although detailed conclusions appear to be very sensitive to the initial conditions (i.e. core shape, thermal, rotational and turbulent energy as well as to the equations of state), some trends can nevertheless be inferred. It is widely accepted that, under realistic initial conditions, a collapsing dense core fragments into a few objects, the exact number depending on the specific conditions. Therefore, fragmentation of a rotating collapsing dense core appears to be the most widely accepted mechanism to explain the formation of multiple systems. Indeed, it is today the only viable mechanism since other possibilities such as fission of a rotating protostar or capture of a companion fail in realistic conditions to produce a large fraction of binaries (Bodenheimer et al. 2000).
Due to the difficulty of the problem and despite its importance, the question of the role of the magnetic field in this process has remained little addressed. With the recent progress achieved in numerical techniques as well as increasing computing power, several studies have investigated this issue. Hosking & Whitworth (2003) using an SPH two fluid code conclude that magnetically subcritical cores do not fragment. Machida et al. (2005) using a nested grid code performed an extensive number of calculations varying the initial core rotation and magnetic field. They find that fragmentation is possible if the rotation is sufficiently large and the magnetic field strength sufficiently small. Fragmentation of a collapsing magnetized cloud using adaptive mesh refinement techniques have also been studied by Ziegler (2005), Banerjee & Pudritz (2006) and Fromang et al. (2006) which all find that the magnetic field has a strong influence. Price & Bate (2007), using magnetized SPH techniques which insure the nullity of the divergence of the magnetic field, find that with large initial perturbations, even for large values of the magnetic strength, fragmentation is possible. In all these studies, it has been found that the magnetic field has a strong impact on the fragmentation of the collapsing dense core.
In this paper, we further study the fragmentation of collapsing cores, focusing mainly on the influence of the magnetic field strength and the amplitude of the initial perturbation. Using the analysis developed in the companion paper by Hennebelle & Fromang (2008, hereafter Paper I), we try to understand, when possible quantitatively, the physical reasons for the numerical results. As in Paper I, we consider initially a uniform density cloud in solid body rotation, threaded by a uniform magnetic field. At this stage, we restrict the problem to the case where the magnetic field and rotation axis are aligned. We note that Price & Bate (2007) have investigated the influence of initially perpendicular rotation axis and magnetic field. More realistic initial conditions, including a non uniform density field, turbulent velocity field and rotation axis non aligned with the magnetic field, will be considered in future studies. In particular, we determine the lowest value of , the mass-to-flux over critical mass-to-flux ratio for which fragmentation is suppressed. It turns out that the results depend strongly on the amplitude of the initial perturbation.
The paper is organized as follow. In the second section, we consider the same initial conditions as used in Paper I, namely a spherical uniform one solar mass dense core. The thermal over gravitational energy of the core is equal to 0.37 and the rotation over gravitational energy is 0.045. Note that such rotations are typical of (perhaps slightly larger than) values observed in dense cores (Goodman et al. 1993). The initial density is about g cm-3and the cloud radius pc. The freefall time is thus of the order of years. We use a barotropic equation of state: , where km s-1 is the sound speed and g cm-3. This set of cloud parameters is known to give rise to a disk that unambiguously fragments in the pure hydrodynamics case. The only difference to Paper I is that an m=2 perturbation of amplitude A=0.1 in the density field, , as well as in the Bz field, is added. In the third section, we further discuss the physical interpretation of the results obtained in our simulations and we estimate analytically the critical value of for which it is expected that the magnetic field stabilizes the disk. In the fourth section, we consider initial density perturbations of amplitude 0.5 and show that in this case, fragmentation can be obtained for a much wider range of . The fifth section provides a discussion of possible mechanisms leading to disk fragmentation and the formation of binaries, even in the presence of magnetic fields. The sixth section concludes the paper.
|Figure 1: Equatorial density and velocity field for mass-to-flux over critical mass-to-flux ratio and perturbation amplitude A=0.1. The size of the figure corresponds to about one ten of the cloud initial radius. In physical units, this is about 300 AU. The freefall time is about years.|
|Open with DEXTER|
We investigate fragmentation in the case where the amplitude of the perturbations is initially A=0.1. In this case it is found that without rotation no fragmentation occurs. Therefore, we call this type of fragmentation, rotationally driven fragmentation. We next consider various values of the magnetic mass-to-flux over critical mass-to-flux ratio, namely (quasi hydrodynamical case), , 20, 5, 2 and 1.25. Recalling that corresponds to the case of a cloud supported by the magnetic field, the last value corresponds to a strongly magnetized supercritical cloud.
In the following, we display the inner part of the collapsing cloud. The size of each plot is about R0/10. This corresponds to a physical size of pc or about 300 AU.
We present results for the weak field cases, and . Figure 2 shows results for . Due to the difference in the initial magnetic strengths, it is not possible to compare the simulations at exactly the same times. The first snapshot shows a spiral pattern which is similar to the one seen in the previous case. A difference however, is that the central bar-like structure appears to be much shorter. This is likely due to the central magnetic field which is strongly amplified by the rapid twisting of field lines. As a result, the cloud first fragments into three fragments instead of two. This is due to the fact that since the central bar-like structure does not fragment in two objects as in the previous case, more material and angular momentum is available in its vicinity to produce one object on each side. The next snapshot shows that a symmetry breaking occurs due to one of the two satellites having merged with the central more massive fragment. The structure in the outer part ( x>0.03 R0) is still quite symmetric. Further fragments are forming in the outer part of the spiral pattern at .
From these results, one concludes that even for values of as large as 50, the magnetic field has a significant impact on the evolution of the centrifugally supported inner structure, particularly on its fragmentation. This is due to the strong amplification of the toroidal and radial magnetic field generated by the differential motions in the collapsing core and in the centrifugally supported structure (see Paper I and the following sections).
Figure 3 shows three snapshots for . The first snapshot shows a weak spiral pattern in the inner part. The second and third panels shows that the centrifugally supported structure grows as in the previous cases but it remains much more uniform and the spiral pattern is much less pronounced. As shown in Paper I, the angular momentum profile is similar to the hydrodynamical case since the magnetic braking is very weak during the collapse phase. Therefore, the difference to the case , is not less angular momentum but rather a stronger magnetic field. Indeed, with a magnetic toroidal component, the velocity of the fast MHD waves, which can be loosely seen as an effective sound speed, is: where (see Sect. 3). Therefore, the disk is stabilized against gravitational fragmentation. As a result, no fragmentation is obtained, only one central star forms and grows by accretion. We stress that the value of the magnetic field corresponding to is very modest, far below the values inferred from observations which indicates that (Crutcher 1999).
Finally, to explore whether the suppression of the fragmentation depends on initial conditions, such as the values of and (the ratios of thermal and rotational to gravitational energies), we have performed runs with different and values. First, we explore the influence of stronger and weaker rotation choosing and respectively, keeping the value of constant. The behaviour of these runs is very similar to the one with presented here. A large disk forms but it remains well axisymmetric and does not fragment. We have also explored the effect of a lower thermal energy, taking and keeping . The disk is slightly more axisymmetric than in the case but does not fragment either.
We conclude that, for small initial density perturbations, a magnetized dense core ideally coupled to the magnetic field and having does not fragment for a large range of initial conditions.
We compare our results with the study of Machida et al. (2005) although it is not straightforward because their initial conditions consist of a filament and therefore are different from ours. Also they used different definitions to quantify the amount of rotation and magnetic energies in their simulations. Comparison with their Fig. 10 is worthwhile nevertheless. Their parameter turns out to be equal to . However, their initial state has a peak density of about 104 cm-3 which is less dense than ours by a factor of roughly 100. Since the ratio of rotational over gravitational energies increases during the collapse, we can say that our correspond to in Machida et al. study (for a homologous contraction , this would indicate that should be divided by roughly 4-5 to be compared with the Machida et al.'s results). For values of , they found that for , they have no fragmentation. Comparing this value with our parameters is again not straightforward but is about . For and , . Therefore, our results broadly agree with the results of Machida et al. (2005).
Comparison with the study of Price & Bate (2007) is not possible at this stage, since they use a barotropic equation of state which becomes adiabatic at 10-14 g cm-3. Thus, even the hydrodynamical calculations they present do not fragment when the perturbation is weak (their Fig. 3).
We present results for the intermediate magnetized cases, namely and . Figure 4 shows two snapshots for . As shown in Paper I, no centrifugally supported structure forms, instead a magnetized pseudo-disk develops. Pseudo-disks arise when disk-like structures form, say oblate ellipsoids, which are not supported by rotation and are produced by the magnetic field. According to the analysis presented in Paper I, the angular momentum is lower for these values of primarily because the collapse occurs first along the field lines. Thus, the material within the pseudo-disk and central object was initially located along the pole and has less angular momentum. Some magnetic braking also occurs, reducing the angular momentum further. As a consequence, no centrifugally supported disk is observed and no fragmentation occurs.
Figure 5 shows two snapshots for . The first panel shows that a filamentary structure has developed after t=1.51 freefall times. We believe that this filament is due to the non-linear evolution of the initial m=2 perturbation. The reason why it has a shape different to the previous case is that since the cloud is more supported by the magnetic field, the collapse lasts 30 longer. Therefore the perturbation has more time to develop and to become non-linear. Indeed, if no perturbation is included initially, an axisymmetric pseudo-disk develops.
One of the important conclusions reached in the previous section is that fragmentation can be suppressed even for large values of , i.e. weak magnetic fields. Here we discuss the physical reason for the disk stabilization arising at low magnetic strength. We stress that since the angular momentum is nearly identical for to the hydrodynamical case (see Paper I), the reason for this stabilization is not magnetic braking.
Elmegreen (1987) and Gammie (1996) both consider the influence of the magnetic field on the stability of a differentially rotating system. In that case, the shear drives the growth of a toroidal magnetic component which stabilizes the disk. They both conclude that whereas in the absence of shear the magnetic field is strongly destabilizing, it has a stabilizing influence when significant shear is present. More precisely, Gammie (1996) computed the response of the disk to nonaxisymmetric perturbations for various value of Q, the Toomre parameter. He shows that the response is much weaker when a substantial magnetic field is present and concludes that the stabilizing effect of the field becomes significant once the magnetic pressure is comparable to the gas pressure. This is qualitatively in good agreement with our results. In particular, we see that the disk in the case remains much more axisymmetric than in the case and .
|Figure 7: within the disk.|
|Open with DEXTER|
For simplicity, the works mentioned above have restricted the analysis to a thin disk. However, the expansion of the magnetic tower triggered by the growth of the toroidal magnetic field removes some material from the disk and therefore reduces the disk surface density. This makes the disk even less prone to fragmentation. This effect, which is not taken into account in the thin disk analysis, certainly contributes to further stabilize the centrifugally supported structure against fragmentation. Indeed, the mass of the disk and the mass within the tower, are roughly comparable.
We estimate analytically the value of at which one expects to find a significant influence of the magnetic field on the fragmentation. Following Gammie (1996), we will consider that a strong stabilizing influence is achieved when the Alfvén speed associated with the toroidal component is comparable to the sound speed.
The principle of the analysis is as follow: we compute the growth rate of the toroidal component inside the centrifugally supported structure, so that we can estimate the time, , needed for the Alfvén speed to become comparable to the sound speed. We then compute the dynamical time of the disk, , over which fragmentation will occur. The critical value of , below which fragmentation is quenched, is obtained when is equal to .
The growth time of the toroidal magnetic component within the disk has two contributions. First the twisting of the radial component due to differential rotation proportional to and second the wrapping of the vertical magnetic component due to the vertical gradient of Bz proportional to , where z is the disk height. The divergence constraint shows that , indicating that the 2 contributions are comparable.
Thus, the growth of can be estimated by writing where h is the disk height. Since we are interested in computing the growth of the magnetic field from initially small values, we neglect the influence of the magnetic tower on the disk vertical expansion and write: . With Eq. (1) of Paper I, we obtain , where d is the ratio of the density over singular isothermal density.
If we assume that the disk is roughly Keplerian, we get
being the mass of
the star and the disk.
Taking the expression of Bz stated by Eqs. (1) of Paper I, we obtain:
For , we obtain . This analytical estimate is therefore in reasonable agreement with our numerical simulation, since for , fragmentation is suppressed. In this case, the toroidal magnetic pressure becomes comparable to the sound speed in about half a rotation period. Since typically fragmentation occurs over a few rotation period, such a fast growth time appears to be sufficient to prevent the disk from fragmenting. On the other hand for the disk fragments, confirming that this case is slightly above the critical value.
For , one has . Therefore in such a case, the growth of the toroidal component is too slow to significantly influence the disk evolution.
|Figure 8: Same as Fig. 1 for and A=0.5. The figure corresponds to a length of 1200 AU, four times bigger than the length displayed in Fig. 1.|
|Open with DEXTER|
We present results for weak and intermediate values of magnetic strengths, namely and . Figure 8 shows results for . The first panel shows that a filamentary structure develops in which two protostars have formed. The two protostars are well separated by roughly as opposed to in the case of weaker perturbations. As explained previously, the fragmentation in this case has a purely thermal origin and is simply due to the development of the initial perturbations (in the case of weak perturbations, rotationally driven fragmentation does not strongly depend on the perturbation amplitude). The second panel shows that the two fragments are approaching each other falling along the filament. The final panel shows that the two fragments have merged and that a protostellar disk forms. However, the second collapse is not treated in the simulation, thus only the first Larson core is considered. Therefore the merging should be considered with great care. Careful thermodynamic treatment of the first Larson core is required. Indeed the first core is known to have a short life time of roughly 103 years. Therefore since the timeshift between the first and second panel is larger than 104 years, it is expected that second collapse would have already occurred by the time the two fragments approach each other. Since stars have very small radii, it is unlikely that the binary will merge. Therefore the final product could be a binary rather than a single protostar.
Figure 9 shows results for . Both panels show that two fragments develop in a similar way as for . Unfortunately, the relative strong value of the magnetic intensity makes the time steps much shorter than in the previous case and it becomes computationally expensive to follow the long range evolution. However, it seems likely that a similar evolution is expected leading to the merging and the two first cores.
Figure 10 shows results for and A=0.5. As revealed by the first panel, the m=2 density perturbation has only weakly developed and only two shallow density maxima are visible. The second panel shows a single object in the centre surrounded by a pseudo-disk. Since the timeshift between the two panels is about 0.015 freefall time ( years) there is not sufficient time for the second collapse to arise. It seems plausible that in this particular case, a single star will form.
Therefore, it appears that a nearly critically magnetized dense core, ideally coupled to the magnetic field, is difficult to fragment into a few objects even if relatively large fluctuations are initially present. This conclusion is in good agreement with the conclusion reached by Hosking & Whitworth (2004), that a magnetically subcritical cloud, which collapses after sufficient magnetic flux has been lost because of ambipolar diffusion, does not fragment.
We conclude that, with large amplitude density fluctuations, a collapsing dense core fragments over a large range of initial conditions as long as it is not too strongly magnetized. The fragments are due to each initial seed collapsing individually. The question of the survival of the fragments should however be carefully investigated. We note that Price & Bate (2007) reach very similar conclusions. They also propose that in some cases, such as when the magnetic field is perpendicular to the initial rotation axis, the magnetic field may help the fragments to survive.
We start with the possibility of having initial sufficiently large fluctuations. As demonstrated in this paper and in Price & Bate (2007), fragmentation is possible in such a situation. The question is then how likely are large initial perturbations? Observationally the question is not easy to address. It is nevertheless known that the dense cores usually are not very uniform and that typically sonic velocity dispersion is observed. This is broadly compatible with the presence of density perturbations although perhaps not as large as 50. On the other hand, cores with a high aspect ratio are observed. If some of them are elongated, i.e. prolate objects, this may be equivalent to a substantial m=2perturbation. Indeed, simulations of quiescent cores having transonic turbulence present significant fluctuations (Goodwin et al. 2004). It seems however difficult at this stage to infer quantitatively whether the perturbations generated by this weak turbulence are sufficient to produce multiple systems.
Another possibility is that the collapse may be induced by an external agent such as a supernova remnant or protostellar jet. These are likely to generate large perturbations. Indeed, all simulations considering clouds evolving far from equilibrium (Bate et al. 2003; Ballesteros-Paredes et al. 2003; Hennebelle et al. 2006; Peretto et al. 2007) do find dense cores with initially strong perturbations.
Future works should specify under what conditions such initial fluctuations are likely to be produced. This may depend on the physical properties of the molecular clouds in which the dense core is embedded. For example, the answer could be different in the Taurus molecular cloud in which star formation is relatively quiescent and in the Orion molecular cloud where star formation is more active.
We stress that this mechanism constitutes a change of paradigm with respect to the standard hydrodynamical scenario (corresponding to Fig. 1). Indeed, with this mechanism the fragmentation appears to be the result of perturbations seeded at large scales rather than an intrinsic properties of the collapsing dense core.
An important prediction of this model is that the binary should not necessarily be located in the equatorial plane of the core since the initial perturbations likely should be randomly distributed.
However, the possibility remains that ambipolar diffusion could occur in the disk allowing
the toroidal magnetic field to be diffused out. To investigate this issue, we estimate
the ratio of ambipolar diffusion time in the disk along the vertical direction and the
growth time of the toroidal magnetic field given by Eq. (2).
Here we write it as
On the other hand, Machida et al. (2007) recently explored magnetized second collapse taking into account the large Ohmic dissipation (we note that this is not rigorously Ohmic dissipation) that is predicted to occur in the first Larson core by various models (e.g. Nakano et al. 2002). As a result, most of the magnetic flux is lost making fragmentation easier. It therefore appears possible that fragmentation could occur during this phase, possibly driven by rotation. In this case, the binary should gain sufficient angular momentum from the accretion to increase the separation between the two stars (see e.g. Goodwin et al. 2004). Interestingly, Banerjee & Pudritz (2006) report the formation of a very close binary in their second collapse calculations despite the ideal MHD assumption.
The prediction of this model is that the binaries should be in the equatorial plane since the orbital angular momentum of the binary is due to the angular momentum of the accreting gas. Also the separation between the 2 stars should increase with time implying that closer binaries should be observed on average in younger cores. Note that this mechanism could possibly work even if the magnetic field is initially very strong in the core.
In order to test each of these scenarios, high resolution observations would be necessary. In this respect, ALMA will certainly be a very powerful tool, as demonstrated by the synthetic observations performed by André et al. (2007) using the present calculations.
With our choice of initial conditions, and for perturbations of low amplitude, unmagnetized cores do fragment, while magnetized cores having values of as large as 20 do not fragment. Based on an analytic estimate, we suggest that for cores having larger than , the growth rate of the toroidal component is too slow to stabilize the disk. We stress that the suppression of fragmentation in this range of parameter, is not due to magnetic braking but to the rapid growth of the toroidal component of the magnetic field induced by the differential rotation within the disk. The Alfvén speed associated with this toroidal component adds up to the sound speed of the disk and stabilize it.
For values of smaller than 5, no big centrifugal disk forms because first the collapse occurs mainly along the field lines bringing less angular momentum, and second, magnetic braking removes angular momentum. This makes these cores even less prone to fragment.
The situation is different if large amplitude perturbations are initially seeded. In this case, each perturbation develops independently even without rotation. Since the strong field amplification is primarily due to differential motions in the disk, the magnetic field is unable to suppress fragmentation except if the core is almost critical, i.e. the field is initially very strong. The following evolution of these fragments requires a careful treatment of the thermodynamics of the first Larson core. It is indeed likely that if the second collapse has not occurred by the time these fragments approach each other, they will merge. On the other hand, if the protostars have already formed, the binary system is likely to survive.
For dense cores having rotation and magnetic strength typical of values inferred from observations, we find that fragmentation is suppressed by the magnetic field if the initial density perturbations are too small. This constitutes a severe problem, since there is no question that a significant fraction of stars are binaries. In view of this, we discuss the likelihood of having sufficient perturbations within the cores initially as well as the impact of ambipolar diffusion and the possibility of fragmentation during the second collapse phase. We speculate that the first may depend on the physical characteristic of the molecular cloud in which the dense core is embedded whereas the latter should berelatively independent of the large scales and could work even if the magnetic field is initially very strong.
Some of the simulations presented in this paper were performed at the IDRIS supercomputing center and on the CEMAG computing facility supported by the French Ministry of Research and Education through a Chaire d'excellence awarded to Steven Balbus. We thank Sébastien Fromang, Doug Johnstone and Philippe André for a critical reading of the manuscript. We thank Frank Shu, the referee, for helpful comments. P.H. thanks Masahiro Machida for related discussions.