Free Access
Issue
A&A
Volume 560, December 2013
Article Number A68
Number of page(s) 6
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201321761
Published online 06 December 2013

© ESO, 2013

1. Introduction

With the recent observations performed with Herschel (e.g. André et al. 2010; Molinari et al. 2010), it has become clearer that filaments are ubiquitous in molecular clouds and that they likely play a central role in star formation. While the exact influence filaments may have on the star formation process remains to be clarified, it is important to understand the properties of these filaments since they are a direct consequence of the physics at play within molecular clouds. In this respect, a particularly intriguing observational result was found by Arzoumanian et al. (2011) who showed that the central widths of the interstellar filaments have a narrow distribution that is peaked around a value of 0.1 pc. Moreover, this characteristic width does not depend on the column density within the filament. Since it is believed that turbulence is important in molecular clouds and largely triggers their evolution, this result is counterintuitive at first sight since turbulence is responsible in a great variety of contexts for producing scale-free power-law distributions. This suggests that there is probably a physical process involved in setting this distribution, which unlike turbulence, presents a characteristic scale.

A recent proposal made by Fischera & Martins (2012, see also Heitsch 2013; Gomez & Vázquez-Semadeni 2013) is that it may result from self-gravitating equilibrium. By solving the hydrostatic equilibrium in an isothermal filament, Fischera & Martins (2012) show that the filament width does not vary significantly and remains on a scale close to the observed 0.1 pc. While this explanation is appealing, a few questions arise. First, it assumes the existence of some confining pressure outside the filament whose nature remains to be specified. Second, it fails to explain why the gravitationally unstable filaments that are collapsing also present this typical width.

In this paper, we explore the idea that the typical width of a self-gravitating filament is due to the combination of accretion-driven turbulence onto the filament, as suggested by the recent velocity dispersion measurements of Arzoumanian et al. (2013), and to the dissipation of this turbulence by ion-neutral friction which does have a characteristic scale. We stress that although ambipolar diffusion is considered here, the underlying idea is totally different from the classical magnetically regulated star formation (e.g. Shu et al. 1987) in which the magnetic field is envisaged as the dominant support and ambipolar diffusion as the process through which the support can be circumvented. In the present picture clouds are typically supercritical. We also note that we do not address the reason why filaments form. As discussed in Hennebelle (2013) this may be due to the shear of the turbulence with a possible further amplification of the anisotropy by gravity for the most massive of them.

The plan of the paper is as follows. In Sect. 2, we present the various assumptions and physical processes used in our simple model. Section 3 describes the results and Sect. 4 concludes the paper.

2. Model and assumptions

2.1. Characteristics of the filament

Herschel observations (e.g. Palmeirim et al. 2013) suggest that the typical average structure of self-gravitating filaments is constituted by i) a central cylinder of nearly uniform density ρf and radius rf; and ii) an envelope whose density profile is ∝r-2. This radial structure is reminiscent of many self-gravitating objects such as Bonnor-Ebert spheres. More precisely, filaments that are collapsing in a self-similar manner are expected to present an envelope with a profile ∝r− 2/(2 − γ) where γ is the adiabatic index of the gas (Kawachi & Hanawa 1998). As it is likely that self-gravitating filaments are collapsing in a way not too different from, although not identical to, a self-similar collapse, assuming an r-2 profile is thus a well-motivated assumption, both from observations and theory. One important difference with such self-similar solutions, however, is that the central density plateau does not seem to be shrinking with time.

If L is the length of the filament, the mass of the central part is obviously Mf=πρfrf2L=mfL\hbox{$M_{\rm f} = \pi \rho_{\rm f} r_{\rm f}^2 L = m_{\rm f} L$}, where Mf is the mass and mf the mass per unit length. The total mass of the central part plus the surrounding envelope is Mtot=Mf[1+2ln(rext/rf)]Mfμ(rext/rf),\begin{eqnarray} M_{\rm tot} = M _{\rm f} [ 1 + 2 \ln ({r_{\rm ext}/r_{\rm f}}) ] \equiv M _{\rm f} \mu ({r_{\rm ext}/r_{\rm f}}), \label{mass_fil} \end{eqnarray}(1)where μ(x) = 1 + 2ln(x). It has been assumed that the density outside the filament is ∝1/r2 and that the filament stops at some radius rext. Below we assume rext ≃ L/2.

2.2. Gravitational potential within the filament

The gravitational potential, in the radial direction is obtained by the Gauss’s theorem r<rf(r)=πGρfr2=Gmf.r>rf(r)=Gmf(1+2log(r/rf)+2(log(r/rf))2),Gmf(1+𝒢(r/rf)).\begin{eqnarray} \label{grav_pot} r<r_{\rm f}, \; \phi(r) &=& \pi G \rho_{\rm f} r^2 = G m_{\rm f} . \\ \nonumber r>r_{\rm f}, \; \phi (r) &=& G m_{\rm f} \left( 1 + 2 \log(r/r_{\rm f}) + 2 (\log(r/r_{\rm f}))^2 \right), \\ & \equiv & G m_{\rm f} ( 1+ {\cal G} (r/r_{\rm f}) ). \nonumber \end{eqnarray}(2)where \hbox{${\cal G}(x)= 2 \ln(x) + 2 (\ln(x))^2$}.

2.3. Magnetic field

As the ion-neutral friction dissipation depends on the magnetic field, it is necessary to know its dependence.

We proceed in two steps. First, we discuss the expected value of the magnetic field in the parent clump Bc. Second, we infer the value of the magnetic field in the filament Bf from the value of Bc.

2.3.1. Magnetic field in the parent clump

The magnetic field in the clump is assumed to be proportional to the square root of the density Bc=B0ρc/ρ0\hbox{$B_{\rm c} = B_0 \sqrt{\rho_{\rm c} / \rho_0}$} as has been observed (e.g. Crutcher 1999). Typical values are B0 ≃ 25  μG and n0 ≃ 103 cm-3, where n0 = ρ0/mp. In the following we will use the vA,0=B0/4πρ01\hbox{$v_{\rm A,0} = B_0 / \sqrt{4 \pi \rho_0} \simeq 1$} km s-1 as a fiducial value.

We note that the magnetic field dependence is still under debate and that there are alternative choices. First, as suggested by Basu (2000), the magnetic field could scale as σρ\hbox{$\sigma \sqrt{\rho}$}, where σ is the velocity dispersion, rather than just as ρ\hbox{$\sqrt{\rho}$}. Second, Crutcher et al. (2010) now favor B ∝ ρ2/3. These relations do not represent large variations and would therefore not affect our results very significantly.

2.3.2. Magnetic field in the filament

To link the magnetic field in the filament to the magnetic field in the parent clump, we proceed as follows. First we assume that the magnetic field is perpendicular to the filament. This configuration is well supported by observations in massive filaments such as Taurus (Heyer et al. 2008; Palmeirin et al. 2013) or DR21 (e.g. Kirby 2009), and is also natural on physical grounds as the gas is expected to accumulate preferentially along the field lines. This implies that at least a fraction of the gas accreted by the filament, is not impeded by the magnetic Lorentz force. There may also be gas accreted perpendicularly to the field lines, which is therefore probably slowed down by magnetic pressure. However, if the field is strong, then most of the material is presumably channeled along the field lines, while if the field is weak it also has a weak influence.

Second, we assume flux freezing which at these scales is a reasonable assumption. This implies that the magnetic field in the filament is simply the magnetic field in the clump compressed along the direction perpendicular to the field and the filament axis. Thus Bf ≃ Bc × ηL/rf since the matter that is inside the radius rf comes from a distance comparable to the clump’s size, ηL, where η typically varies with time between 0 and 1/2. A similar reasoning can be applied to get a relation between ρc and ρf since Mfρfrf2Lρcη2L3\hbox{$M_{\rm f} \propto \rho_{\rm f} r_{\rm f} ^2 L \simeq \rho_{\rm c} \eta^2 L^3$}. We thus obtain ρf = ρc(ηL/rf)2. Combining the expression for Bc obtained above with the last expression, we get Bf=B0ρfρ0\hbox{$B_{\rm f} = B_0 \sqrt{ \rho_{\rm f} \over \rho_0}$} that is to say, the magnetic field in the filament is also expected to be nearly proportional to ρf\hbox{$\sqrt{\rho_{\rm f}}$} implying that the Alfvén velocity should remain nearly constant.

This relation is valid as long as flux freezing can be assumed. While this is a reasonable assumption in the collapsing envelope, which is not magnetically supported, it is not the case in the central part, which is presumably close to equilibrium. The typical ambipolar diffusion timescale is given by Eq. (13) below. For densities on the order of 104 − 7 cm-3, the field is diffused in about 0.1−1 Myr, which is shorter than or comparable to, the accretion time of the filaments. Moreover, as emphasized by Santos-Lima et al. (2012) and Joos et al. (2013), turbulence also tends to diffuse out the field. However since self-gravitating filaments are probably accreting, the magnetic flux cannot leak out far away since it is confined by the infalling gas. Therefore, while the mean magnetic intensity within the central part of these filaments is probably on the order of Bf=B0ρfρ0\hbox{$B_{\rm f} = B_0 \sqrt{ \rho_{\rm f} \over \rho_0}$}, it is likely that the magnetic field gradient is greatly reduced with respect to the ideal MHD case.

2.4. Virial equilibrium

The virial theorem is applied to the filament inner part of radius rf. The expression for a filament is (e.g. Fiege & Pudritz 2000) 2PdV+2Ekincyl=Wgrav+2PextV,\begin{eqnarray} 2 \int P {\rm d}V + 2 E_{\rm kin}^{\rm cyl} = W_{\rm grav} + 2 P_{\rm ext} V, \label{virial_pasgen} \end{eqnarray}(3)where Wgrav is the gravitational term, P and Pext are the internal and external pressure, and V the volume; Ekincyl=0.5Mf(2σ1D2)\hbox{$E_{\rm kin}^{\rm cyl} = 0.5 M_{\rm f} (2 \sigma_{\rm 1D} ^2)$} is the kinetic energy in the direction perpendicular to the filament axis and σ1D is the non-thermal one-dimensional velocity dispersion. This expression is, however, not strictly valid since our model filament is accreting. Additional terms should be taken into account (see Goldbaum et al. 2011 and Hennebelle 2012) which corresponds to the surface terms that do not cancel out as it is usually the case. When the surface terms are taken into account, the virial expression becomes 122drf2dM+2PdV+2Ekincyl=Wgrav+2(Pext+Pram)V.\begin{eqnarray} {1 \over 2 } {\dot{M}^2} {{\rm d}r_{\rm f}^2 \over {\rm d}M} + 2 \int P {\rm d}V + 2 E_{\rm kin}^{\rm cyl} = W_{\rm grav} + 2 (P_{\rm ext} + P_{\rm ram}) V .\label{virial_gen} \end{eqnarray}(4)At this stage we do not consider the influence that the magnetic field may have on the equilibrium because it is not expected to change our conclusions qualitatively. In particular, as discussed in the previous section, the magnetic gradient within the central part of the filament is probably smoothed because of ambipolar diffusion. Another complication arises because the anisotropy introduced by the magnetic field is perpendicular to the filament axis, which would require a bi-dimensional analysis.

Using the different expressions obtained above we get 122fmfdrf2dmf+2cs2+2σ1D2=Gmf+2Pramρf,\begin{eqnarray} {1 \over 2 } {\dot{m}_{\rm f}^2 \over m_{\rm f}} { {\rm d}r_{\rm f}^2 \over {\rm d}m_{\rm f}} + 2 c_{\rm s}^2 + 2 \sigma_{\rm 1D}^2 = G m _{\rm f} + 2 {P_{\rm ram} \over \rho_{\rm f}},\label{virial_comp} \end{eqnarray}(5)where Pram is the ram pressure exerted by the incoming flow and where the pressure of the external medium has not been taken into account. The ram pressure will be estimated below. For the sake of simplicity, we will also use the simplified form of the virial equilibrium 2σ1D2Gmf.\begin{eqnarray} 2 \sigma_{\rm 1D} ^2 \simeq G m_{\rm f}.\label{virial} \end{eqnarray}(6)

2.5. Mechanical energy balance

The mechanical energy balance within the cylinder of radius rf leads to Mfσ3D2(rf)2τdissϵeff(φ(rext)φ(rf)).\begin{eqnarray} M_{\rm f} { \sigma^2_{\rm 3D} (r_{\rm f}) \over 2 \tau_{\rm diss} } \simeq \epsilon _{\rm eff} \dot{M} \left( \phi(r_{\rm ext}) - \phi(r _{\rm f}) \right)\!.\label{ener_balance} \end{eqnarray}(7)Obviously, the left hand-side is the dissipation which can be due either to the turbulence cascade or to the friction between ions and neutrals as described below. We note that it is assumed that the turbulence is isotropic which is why σ3D is used. The right-hand side describes the source of turbulence which is due to the accretion onto the central part of the filament (Klessen & Hennebelle 2010). The efficiency ϵeff is not well known. Klessen & Hennebelle (2010) proposed that it can be related to the density contrast between the density of the incoming flow and the density of the actual gas in which energy is injected. In the present case, the accretion shock may not be clearly defined because of the turbulent nature of the flow. In any case, the present calculation remains largely indicative at this stage. Below, the value ϵeff = 0.5 is used because it leads to good agreement with the data. We stress that since our model remains indicative, ϵeff could also take into account various other uncertainties. For the sake of simplicity, we have also ignored terms associated to the volume variation and the external pressure (e.g. Goldbaum et al. 2011) as they do not modify the results substantially, but make the mechanical energy balance much more complex.

2.6. Accretion rate

The accretion rate remains uncertain since it is difficult to infer observationally (see however Palmeirim et al. 2013 for an estimate in the case of the Taurus B211/3 filament). Here, we consider two different possibilities. This will allow us to test the robustness of our conclusion.

2.6.1. Gravitational accretion rate

To estimate the accretion rate, we assume that it can be computed from the density within the parent clump and the infall velocity. Assuming that the parent clump has a cylindrical radius of about L/2 and a length equal to L, we get  = πL2ρcVinf.

The infall velocity is due to the gravitational field of the filament which is given by Eq. (2). We assume that the material which enters the clump at radius r ≃ L/2 has no initial velocity and we estimate the infall velocity at r = L/4 leading for Vinf to Vinf2Gmf(𝒢(L2rf)𝒢(L4rf))1/2·\begin{eqnarray} V_{\rm inf} \simeq \sqrt{2 G m_{\rm f}} \left( {\cal G} \left({L \over 2 r_{\rm f}}\right) - {\cal G} \left( {L \over 4 r_{\rm f}}\right) \right)^{1/2}\cdot\label{vinf} \end{eqnarray}(8)We note that this is again a rough estimate, but since the gravitational potential varies logarithmically with r, this estimate does not depend severely on these assumptions.

The density within the clump is also needed to get the accretion rate and we assume that ρc = mpnc follows the Larson relations (Larson 1981; Falgarone et al. 2009; Hennebelle & Falgarone 2012), mp being the mass per particle nc=n0(Rc1pc)ηd,σ3D=σ0(Rc1pc)η,\begin{eqnarray} n_{\rm c} = n_0 \left( { R_{\rm c} \over 1~{\rm pc}} \right)^{-\eta_{\rm d}},\,\,\sigma_{\rm 3D} = \sigma_0 \left( { R_{\rm c} \over 1~{\rm pc}} \right)^{\eta},\label{larson} \end{eqnarray}(9)where nc is the clump gas density and σ3D the internal rms velocity. The exact values of the various coefficients remain somewhat uncertain. Originally, Larson (1981) estimated ηd ≃ 1.1 and η ≃ 0.38, but more recent estimates (Falgarone et al. 2009) using larger sets of data suggest that ηd ≃ 0.7 and η ≃ 0.45 − 0.5. For the sake of simplicity, we use ηd = 1 and take n0 = 1000 cm-3. πL2ηdρ0(1pc)ηd2Gmf(𝒢(L2rf)𝒢(L4rf))1/2·\begin{eqnarray} \dot{M} \simeq \pi L^{2-\eta_{\rm d}} \rho_0 (1 {\rm pc})^{\eta_{\rm d}} \sqrt{2 G m_{\rm f}} \left( {\cal G} \left( {L \over 2 r_{\rm f}} \right) - {\cal G} \left( {L \over 4 r_{\rm f}} \right) \right)^{1/2}\cdot\label{accret_rate_grav0} \end{eqnarray}(10)The typical accretion timescale is simply given by τaccret = M/. With Eq. (10), it is easy to show that τaccretρf\hbox{$\tau_{\rm accret} \propto \sqrt{\rho_{\rm f}}$}.

2.6.2. Turbulent accretion rate

As it is not clear what controls the accretion rates onto interstellar filaments, we also consider the turbulent accretion rate constructed from the Larson relations described in Hennebelle (2012), MtotτcMtot2Rc/(σ3D/3)0(MtotM0)ηacc,\begin{eqnarray} \label{dM_dt} \dot{M} &\simeq& {M_{\rm tot} \over \tau_{\rm c}} \simeq {M_{\rm tot} \over 2 R_{\rm c}/(\sigma_{\rm 3D}/\!\!\sqrt{3})} \simeq \dot{M}_0 \left( { M_{\rm tot} \over M_0} \right)^{\eta_{\rm acc}} , \end{eqnarray}(11)where ηacc ≃ 0.7 − 0.8 depending on the exact choice of the parameter η and ηd that is retained. We will adopt ηacc ≃ 0.75 as a fiducial value. We typically have 0 = 10-3  M yr-1 for M0 = 104  M.

2.6.3. Ram pressure

The ram pressure which appears in Eq. (5) can be estimated as follows. It is equal to the product of Vinf(rf)2 and ρin = /(2πrfLVinf(rf)), where ρin is the density that is obtained assuming a constant accretion rate. We note that ρin < ρf, which implies that an accretion shock is connecting the infalling envelope and the central part of the filament. Since the above expression assumes that the flow is isotropic, and since we are assuming the structure of the flow rather than inferring an exact solution, this value remains hampered by large uncertainties and is certainly valid within a factor of a few. We have tested the influence of vaying the ram pressure by a factor of a few and found that it has a limited influence on the solution at low density, while it has no influence at high density.

2.7. Dissipation timescales

The dissipation timescale to be used in Eq. (7) is a crucial issue. Here we emphasize two dissipation mechanisms, the turbulent cascade time and the ambipolar diffusion time. Quantitative estimates of these two timescales have been recently estimated by Li et al. (2012) in turbulent two fluid MHD simulations. They found that under typical molecular cloud conditions both contribute, but the second dominates over the first.

2.7.1. Dissipation by turbulent cascade

First, we consider the standard turbulent cascade timescale, which is the crossing time of the system τdiss,c2rfσ1D·\begin{eqnarray} \tau_{\rm diss,c} \simeq {2 r_{\rm f} \over \sigma _{\rm 1D}}\label{tau_diss_f}\cdot \end{eqnarray}(12)The energy cascades to smaller and smaller scales until the size of the eddies reaches the viscous scale.

2.7.2. Dissipation by ion-neutral friction

Second, we investigate the dissipation induced by the ion-neutral friction. Its expression was first inferred by Kulsrud & Pearce (1969, see also Lequeux 2005) and is given by τdiss,amb=2γdampρivA2(2π/λ)2=2νnivA2(2π/λ)2,\begin{eqnarray} \tau_{\rm diss, amb} = { 2 \gamma_{\rm damp} \rho_i \over v_{\rm A}^2 (2 \pi / \lambda)^2} = { 2 \nu_ {\rm ni} \over v_{\rm A}^2 (2 \pi / \lambda)^2}, \label{tau_diss_amb} \end{eqnarray}(13)where vA is the Alfvén speed, λ is the wavelength assumed to be equal to rf and νni is the ion-neutral coupling coefficient. The reason we choose λ ≃ rf is that if most of the energy is dissipated at a scale much smaller than rf, then the relevant time would be the crossing or cascading time that would be necessary for the energy to cascade from rf. In this case, the timescale would thus be similar to the turbulent cascade timescale discussed above. The coefficient γdamp = 3.5 × 1013 cm3 g-1 s-1 is the damping rate. The ion density, ρi is assumed to be ρi=Cρn\hbox{$\rho_i = C \sqrt{\rho_n}$} where C = 3 × 10-16 cm−3/2 g1/2. For wavelengths λ < λc = πvA/(γdampρi), the critical wavelength, the Alfvén waves do not even propagate except at very small wavelengths when the ions and the neutrals are entirely decoupled. With n ≃ 104 cm-3 and vA ≃ 1 km s-1, λc ≃ 5 × 10-2 pc.

We note that the expression stated in Eq. (13) is strictly valid only for Alfvén waves. The corresponding expression for the compressible modes has been inferred by Ferrière et al. (1988). They are slightly more complex as it entails the angle between the field and the direction of propagation, but the order of magnitude is not different.

3. Results

To infer the radius of the filament as a function of the central density ρf we have to combine Eqs. (2), (6), and (7) together with the accretion rate given by Eq. (10) or Eq. (11). With the gravitational accretion rate (Eq. (10)), we obtain ρf1/2rfτdiss=432πϵeffρ0(1pc)G𝒢1(L/rf),\begin{eqnarray} { \rho _{\rm f} ^{1/2} r_{\rm f} \over \tau_{\rm diss} } = {4 \over 3} \sqrt{2 \pi} \epsilon_{\rm eff} \rho_0 (1 {\rm pc}) \sqrt{G} {\cal G}_1 (L/r_{\rm f}),\label{combi_grav} \end{eqnarray}(14)where \hbox{${\cal G}_1(u) = \left( {\cal G} (u/2) - {\cal G}(u/4) \right)^{1/2} {\cal G}(u/2)$}. With the turbulent accretion rate (Eq. (11)), we get the relation (πρfrf3)1ηaccτdiss=4ϵeff0𝒢2(L/rf)3M0ηacc(L/rf)1ηacc\begin{eqnarray} { (\pi \rho_{\rm f} r_{\rm f}^3) ^{1-\eta_{\rm acc}} \over \tau_{\rm diss}} = {4 \epsilon_{\rm eff} \dot{M}_0 {\cal G}_2(L / r_{\rm f}) \over 3 M_0^{\eta_{\rm acc}} (L/r_{\rm f}) ^{1-\eta_{\rm acc}}}\label{combi} \end{eqnarray}(15)where \hbox{${\cal G}_2 (u) = \mu (u/2)^{\eta_{\rm acc}} {\cal G} (u/2)$}.

3.1. Dynamical equilibrium with turbulent dissipation

Combining Eq. (12) with Eq. (6), we get τdiss=22πGρf\hbox{$\tau_{\rm diss} = {2 \sqrt{2} \over \sqrt{\pi G \rho_{\rm f}}}$} which in turn together with Eq. (15) leads to the expression rf=ρfηacc3/23(1ηacc)(823π3/2ηaccGϵeff0𝒢2(L/rf)M0ηacc(L/rf)1ηacc)13(1ηacc)·\begin{eqnarray} r_{\rm f} = \rho_{\rm f} ^{ \eta_{\rm acc} -3/2 \over 3 (1-\eta_{\rm acc}) } \left( { 8 \sqrt{2} \over 3 \pi ^{3/2-\eta_{\rm acc}} \sqrt{G}} {\epsilon_{\rm eff} \dot{M}_0 {\cal G}_2 (L / r_{\rm f}) \over M_0^{\eta_{\rm acc}} (L/r_{\rm f}) ^{1-\eta_{\rm acc}} } \right)^{1 \over 3(1-\eta_{\rm acc})}\cdot\label{radius_turb} \end{eqnarray}(16)For the canonical value ηacc = 0.75, we find that rf ∝ 1/ρf. That is to say, the typical filament radius decreases with the central density as displayed in Fig. 1 (see the line labeled turbulence). The value ϵeff = 0.5 has been used in this calculation. For larger values of ηacc we still get significant variations of rf with ρf. For ηacc = 1 we would even predict that the central density is independent of rf. A filament radius independent of the central density is obtained only for ηacc ≃ 3/2. The gravitational accretion rate expression leads to a very similar expression with rf ∝ 1/ρf and the corresponding expression is not given here for conciseness as the obvious conclusion is that this behaviour is incompatible with a nearly constant filament width.

thumbnail Fig. 1

Filament radius as a function of filament density for four models. The curve labeled turbulence shows the result for a turbulent crossing time (Eq. (16)); the three curves labeled friction show results when ion-neutral friction time is assumed to be the energy dissipation time (Eqs. (17)−(18)). Frictions A and C used a gravitational accretion rate and Friction B a turbulent rate.

3.2. Dynamical equilibirum with ion-neutral friction

Combining Eqs. (13) and (6), with (14), we infer rf=3π3/2va,0222ϵeffγdampCρ0(1pc)G𝒢1(L/rf),\begin{eqnarray} r_{\rm f} = { 3 \pi ^{3/2} v_{\rm a,0}^2 \over 2 \sqrt{2} \epsilon_{\rm eff} \gamma_{\rm damp} C \rho_0 (1 {\rm pc}) \sqrt{G} {\cal G}_1 (L/r_{\rm f})},\label{ambi_grav} \end{eqnarray}(17)where we see that rf does not depend on ρf. That is to say, the width of the filament does not change with its column density as suggested from the results of Arzoumanian et al. (2011). To test the robustness of this result, it is worth investigating what the turbulent accretion rate stated by Eq. (11) is predicting. The corresponding expression is rf=ρfηacc1/213ηacc(2γdampC3π3ηaccva,02ϵeff0𝒢2(L/rf)M0ηacc(L/rf)1ηacc)113ηacc·\begin{eqnarray} r_{\rm f} = \rho_{\rm f} ^{ \eta_{\rm acc} -1/2 \over 1 - 3 \eta_{\rm acc}} \left( { 2 \gamma_{\rm damp} C \over 3 \pi^{3-\eta_{\rm acc}} v_{\rm a,0} ^2 } {\epsilon_{\rm eff} \dot{M}_0 {\cal G}_2 (L / r_{\rm f}) \over M_0^{\eta_{\rm acc}} (L/r_{\rm f}) ^{1-\eta_{\rm acc}} } \right)^{1 \over 1-3 \eta_{\rm acc} }\cdot \label{ambi_turb} \end{eqnarray}(18)As can be seen for an accretion rate exponent ηacc of the order of 0.75, we find that rfρf-0.2\hbox{$r_{\rm f} \propto \rho_{\rm f} ^{-0.2}$}, which implies a very weak dependence on the filament radius rf. For a value of ηacc = 1, we have rfρf1/4\hbox{$r_{\rm f} \propto \rho_{\rm f} ^{-1/4}$} which is still a shallow dependence as shown in Fig. 1 where the two expressions stated by Eqs. (17) and (18) are displayed (labeled as friction A and B, respectively).

thumbnail Fig. 2

Amplitude of the various terms which appear in the virial equilibrium (Eq. (5)). While at high density the filament equilibrium is due to the balance between gravity and velocity dispersion, it is due to the balance between ram and thermal pressure at low density.

3.3. A more accurate model

Finally, we investigate the solutions where the thermal support and ram pressure are considered as stated in Eq. (5). The corresponding curve is labeled as friction C. Equation (5) is an ordinary differential equation in rf, which can be solved using a standard Runge-Kutta method. Since it is necessary to specify a radius and a density to start the integration, we have explored various cases. We found that for a large range of radii (rstart ≃ 0.01 − 0.1 pc) at low density, the solutions quickly converge towards the one that is presented here and for which the radius at n = 103 cm-3 is equal to about 0.05 pc.

As can be seen, more variability is introduced, particularly at low density, where we see that the filament radius decreases at low density. In order to better understand the physical meaning of this solution we plot in Fig. 2, the values of the various virial terms as a function of density.

While the equilibrium between gravity and turbulent support at high density is robust and independent of the choice of the boundary condition, rstart, the behaviour at low density is less robust and varies with it.

It is important to stress a few points. First, the ram pressure term which causes most of the variation remains uncertain since our model is not fully self-consistent in the sense that the density and velocity fields, although reasonable, are not proper solutions of the problem. Second, in the low density regime, the filament is not gravitationally accreting and it is likely that the validity of the model is questionable.

3.4. Comparison of the two dissipation timescales

It is enlighting to compare the values of the dissipation timescales as a function of density for the filament radius corresponding to model B. As expected, the turbulent dissipation timescale is much longer than the ion-neutral friction timescale for densities lower than 106 cm-3. It also increases with density while the turbulent timescale decreases with density. This behaviour is the very reason which explains the nearly constant width of the filaments in our model because as can be seen in Fig. 3, the accretion time presents the same dependence.

thumbnail Fig. 3

Comparison between the turbulent dissipation and ion-neutral friction times as a function of density. Also shown is the accretion timescale. It scales exactly as the ion-neutral friction time does.

To summarize, assuming that the relevant timescale for energy dissipation within the central part of the filament is the turbulent crossing time, we find under reasonable assumptions for the accretion rate that the width changes significantly with density. This is because τdiss,c ∝ 1/ρf\hbox{$\!\sqrt{\rho_{\rm f}}$}. On the other hand, when we assume that the relevant timescale for energy dissipation is the ion-neutral friction time, we find that the width varies much less with the filament density. This is because, τdiss,ambρiρfτaccret\hbox{$\tau_{\rm diss,amb} \propto \rho_i \propto \sqrt{\rho_{\rm f}} \propto \tau _{\rm accret}$}. Since the relevant timescale is the shortest one, which corresponds to the smallest value of rf, one expects ion-neutral friction to be the dominant mechanism for energy dissipation up to densities of a few 105 − 106 cm-3 (see Fig. 1).

thumbnail Fig. 4

Comparison between the models and the filaments width distribution (adapted from Arzoumanian et al. 2011).

3.5. Comparison with observations

Finally, we confront the present models with the Herschel observational result of Arzoumanian et al. (2011). Figure 4 shows filament width as a function of filament column density. As can be seen, the model based on ion-neutral friction and gravitational accretion (friction B and C) work very well. We note that the model based on turbulent dissipation predicts a constant column density and a variable radius.

4. Conclusion

We have presented a simple model to describe the evolution of accreting self-gravitating filaments within molecular clouds. It assumes virial equilibrium between gravity and turbulence and mechanical energy balance between accretion which drives turbulence in the filament and the dissipation of this energy. We show that while dissipation based on turbulent cascade fails to reproduce the narrow range of radius inferred from Herschel observations, dissipation based on ion-neutral friction leads to a filament width that depends only weakly on the filament density and is very close to the ≃0.1 pc value, although our analytical approach is hampered by significant uncertainties. We conclude that the combination of accretion-driven turbulence and ion-neutral friction is a promising mechanism to explain the structure of self-gravitating filaments and deserves further investigation.

Acknowledgments

We thank the anonymous referee for a constructive and helpful report. We thank Doris Arzoumanian and Evangelia Ntormousi for discussions on this topic. P.H. acknowledge the financial support of the Agence National pour la Recherche through the COSMIS project. This research has received funding from the European Research Council under the European Community’s Seventh Framework Programme (FP7/2007−2013 Grant Agreements No. 306483 and No. 291294).

References

  1. André, P., Men’shchikov, A., Bontemps, S., et al. 2010, A&A, 518, L102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Arzoumanian, D., André, Ph., Didelon, P., et al. 2011, A&A, 529, L6 [Google Scholar]
  3. Arzoumanian, D., André, Ph., Peretto, N., & Konyves, V. 2013, A&A, 553, A119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Basu, S. 2000, ApJ, 540, L103 [NASA ADS] [CrossRef] [Google Scholar]
  5. Crutcher, R. 1999, ApJ, 520, 706 [NASA ADS] [CrossRef] [Google Scholar]
  6. Crutcher, R., Wandelt, B., Heiles, C., Falgarone, E., & Troland, T. 2010, ApJ, 725, 466 [NASA ADS] [CrossRef] [Google Scholar]
  7. Falgarone, E., Pety, J., & Hily-Blant, P. 2009, A&A, 507, 355 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Ferrière, K., Zweibel, E., & Shull, J. 1988, ApJ, 332, 984 [NASA ADS] [CrossRef] [Google Scholar]
  9. Fiege, J., & Pudritz, R. 2000, MNRAS, 311, 85 [NASA ADS] [CrossRef] [Google Scholar]
  10. Fischera, J., & Martin, P. 2012, A&A, 542, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Goldbaum, N., Krumholz, M., Matzner, C., & McKee, C. 2011, ApJ, 738, 101 [NASA ADS] [CrossRef] [Google Scholar]
  12. Gomez, G., & Vàquez-Semadeni, E. 2013, ApJ, submitted [arXiv:1308.6298] [Google Scholar]
  13. Heitsch, F. 2013, ApJ, 769, 115 [NASA ADS] [CrossRef] [Google Scholar]
  14. Hennebelle, P. 2012, A&A, 545, A147 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Hennebelle, P. 2013, A&A, 556, A153 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Hennebelle, P., & Falgarone, E. 2012, A&ARv, 20, 55 [NASA ADS] [CrossRef] [Google Scholar]
  17. Heyer, M., Gong, H., Ostriker, E., & Brunt, C. 2008, ApJ, 680, 420 [NASA ADS] [CrossRef] [Google Scholar]
  18. Joos, M., Hennebelle, P., Ciardi, A., & Fromang, S. 2013, A&A, 554, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Kawachi, T., & Hanawa, T. 1998, PASP, 50, 577 [Google Scholar]
  20. Kirby, L. 2009, ApJ, 694, 1056 [NASA ADS] [CrossRef] [Google Scholar]
  21. Klessen, R., & Hennebelle, P. 2010, A&A, 520, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Kulsrud, R., & Pearce, W. 1969, ApJ, 156, 445 [Google Scholar]
  23. Larson, R. 1981, MNRAS, 194, 809 [CrossRef] [Google Scholar]
  24. Lequeux, J. 2005, The interstellar medium (Springer) [Google Scholar]
  25. Li, P.-S., Myers, A., & McKee, C. 2012, ApJ, 760, 33 [NASA ADS] [CrossRef] [Google Scholar]
  26. Molinari, S., Swinyard, B., & Bally, J. 2010, A&A, 518, L100 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Palmeirim, P., André, Ph., Kirk, J., et al. 2013, A&A, 550, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Santos-Lima, R., de Gouveia Dal Pino, E., & Lazarian, A. 2012, ApJ, 747, 21 [NASA ADS] [CrossRef] [Google Scholar]
  29. Shu, F. H., Adams, F. C., & Lizano, S. 1987, ARA&A, 25, 23 [Google Scholar]

All Figures

thumbnail Fig. 1

Filament radius as a function of filament density for four models. The curve labeled turbulence shows the result for a turbulent crossing time (Eq. (16)); the three curves labeled friction show results when ion-neutral friction time is assumed to be the energy dissipation time (Eqs. (17)−(18)). Frictions A and C used a gravitational accretion rate and Friction B a turbulent rate.

In the text
thumbnail Fig. 2

Amplitude of the various terms which appear in the virial equilibrium (Eq. (5)). While at high density the filament equilibrium is due to the balance between gravity and velocity dispersion, it is due to the balance between ram and thermal pressure at low density.

In the text
thumbnail Fig. 3

Comparison between the turbulent dissipation and ion-neutral friction times as a function of density. Also shown is the accretion timescale. It scales exactly as the ion-neutral friction time does.

In the text
thumbnail Fig. 4

Comparison between the models and the filaments width distribution (adapted from Arzoumanian et al. 2011).

In the text

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.