Reflection and dissipation of Alfvén waves in interstellar clouds
^{1} LUTHObservatoire de Paris, 5 Place J. Janssen, 92190 Meudon, France
email: cecilia.pinto@obspm.fr
^{2} SolarTerrestrial Center of Excellence – SIDC, Royal Observatory of Belgium, 1180 Bruxelles, Belgium
^{3} INAF – Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 50125 Firenze, Italy
^{4} Dipartimento di Astronomia e Scienza dello Spazio, Università di Firenze, Largo E. Fermi 5, 50125 Firenze, Italy
Received: 10 February 2012
Accepted: 11 June 2012
Context. Supersonic nonthermal motions in molecular clouds are often interpreted as longlived magnetohydrodynamic (MHD) waves. The propagation and amplitude of these waves is affected by local physical characteristics, most importantly the gas density and the ionization fraction.
Aims. We study the propagation, reflection and dissipation of Alfvén waves in molecular clouds deriving the behavior of observable quantities such as the amplitudes of velocity fluctuations and the rate of energy dissipation.
Methods. We formulated the problem in terms of Elsässer variables for transverse MHD waves propagating in a onedimensional inhomogeneous medium, including the dissipation due to collisions between ions and neutrals and to a nonlinear turbulent cascade treated in a phenomenological way. We considered both steadystate and timedependent situations and solved the equations of the problem numerically with an iterative method and a LaxWendroff scheme, respectively.
Results. Alfvén waves incident on overdense regions with density profiles typical of cloud cores embedded in a diffuse gas suffer enhanced reflection in the regions of the steepest density gradient, and strong dissipation in the core’s interior. These effects are especially significant when the wavelength is intermediate between the critical wavelength for propagation and the typical scale of the density gradient. For larger wave amplitudes and/or steeper input spectra, the effects of the perpendicular turbulent cascade result in a stronger energy dissipation in the regions immediately surrounding the dense core.
Conclusions. The results may help to interpret the sharp decrease of line width observed in the environments of lowmass cloud cores in several molecular transitions.
Key words: magnetohydrodynamics (MHD) / waves / turbulence / ISM: magnetic fields / ISM: clouds / ISM: lines and bands
© ESO, 2012
1. Introduction
The supersonic nonthermal motions observed in molecular clouds have been attributed to the presence of magnetohydrodynamic (MHD) waves of subAlfvénic amplitudes (Arons & Max 1975), of which the Alfvén waves are the longestlived (Zweibel & Josafatsson 1983). Indeed, molecular clouds and the cores embedded in them are threaded by the Galactic magnetic field, as revealed by polarization maps of the thermal dust emission (see e.g. Tang et al. 2009).
Molecular cloud cores, the sites of star formation, are characterized by velocity dispersions with nonthermal motions subsonic and uniform, in contrast to the surrounding lower density gas, where nonthermal motions are supersonic and in general increasing with size (Barranco & Goodman 1998; Goodman et al. 1998; Caselli et al. 2002). The NH_{3} observations by Pineda et al. (2010) of the B5 region in Perseus revealed that the transition in velocity dispersion between the quiescent core and its more “turbulent” environment is sharp, taking place in a thin layer with a width of less than 0.04 pc.
Several groups have performed simulations and developed (semi)analytical models of the propagation of nonlinear and linear Alfvén waves in interstellar clouds, but the main focus of these studies was either on the generation of density fluctuations and velocity dispersions with properties comparable to those of molecular clouds (Heitsch & Burkert 2002; Kudoh & Basu 2003, 2006), or on the role of internally generated Alfvén waves to support clumps within molecular clouds against their own selfgravity, taking into account the ionization structure and the ionneutral drift (Martin et al. 1997; Coker et al. 2000). Our main objective is to study how density gradients, ionization fraction and nonlinear wave interaction affect the propagation and the amplitude of Alfvén waves at the small scale of the interface between a dense core and its environment.
For this goal, we considered a simple onedimensional cloud model for which the propagation of Alfvén waves is formulated in terms of Elsässer variables. This formulation allows one to treat propagation, reflection, damping, and nonlinear interactions in a simpler form compared to the original fluid equations (Elsässer 1950; Dobrowolny et al. 1980; Marsch & Mangeney 1987). Our approach is similar to that adopted for studying the propagation and reflection of Alvén waves in radially stratified stellar atmospheres and winds (Velli 1993). In this context the density gradients have the main effect of modulating the amplitude of Alfvén waves. This may affect the damping that results from the ionneutral drift and at the same time trigger the development of a (perpendicular) turbulent cascade. Here we will focus mostly on the first aspect, although we will also use a phenomenological model to describe the dissipation of turbulence (successfully applied to the solar wind, e.g. Verdini et al. 2010) to see how it affects the results.
The paper is organized as follows: in Sect. 2 we derive the wave equation and review the conditions for wave propagation; in Sect. 3 we formulate the problem in terms of Elssässer variables; in Sect. 4 we describe the method of solution; we discuss some special cases in Sect. 5 and, in Sect. 6, we consider wave propagation, damping and reflection in a onedimensional cloud model with a specific density and ionization profile; we evaluate the energy dissipation due to both ionneutral drift and nonlinear effects; finally in Sect. 7 we summarize our conclusions.
2. Wave equation and dispersion relation
We considered a onedimensional density distribution ρ(x), threaded by a magnetic field , in which the field is stationary and uniform, the latter condition resulting from the requirement of zero divergence, ∇·B = 0. The total density of the gas ρ is the sum of the densities of neutrals (ρ_{n}), ions (ρ_{i}) and electrons (ρ_{e}), and we assumed that the gas is electrically neutral. For applications to molecular clouds, we assumed that the field is frozen in the ions, which transfer momentum to the neutrals via elastic collisions. We neglected the electron’s inertia and focused on the ambipolar diffusion regime (see e.g. Pinto et al. 2008), ignoring Hall and Ohmic diffusion. In the regime of interest in this work (n(H_{2}) ~ 10^{4}–10^{5} cm^{3}, B ~ 10–100 μG), Ohmic diffusion is about 10 orders of magnitudes smaller than ambipolar diffusion and can be safely neglected; on the other hand, Hall diffusion is about one order of magnitude smaller than ambipolar diffusion, for standard grain size distribution, and about 4 orders of magnitudes smaller in the case of large (micron size) dust grains (see, e.g., Wardle & Ng 1999; Nakano et al. 2002; Pinto et al. 2008). The consequences of the Hall effect on the propagation of Alfvén waves in nonhomogeneous media have been considered by Campos & Isaeva (1993).
In our analysis, we also assumed that the magnetic force is transferred to the neutrals only by collisions with ions. Actually, in the density regime of interest, charged dust grains can contribute as much as ions in the frictional force, provided that a population of small ( ≲ μm size) grains is present (Nakano et al. 2002). If only large ( ≳ μm size) grains are present, as expected if some amount of dust coagulation has taken place (Ormel et al. 2009), the contribution of grains to the frictional force becomes much smaller. On the other hand, largesize dust grains may induce dispersive effects on the propagation of nonlinear magnetic disturbances, and even promote the formation of soliton structures, although on scales at best one order of magnitude smaller than the typical size of a molecular cloud core (Pandey et al. 2008).
We considered transverse disturbances that are characterized by magnetic field b and velocities u_{n} and u_{i} (the latter were also assumed to be incompressible). We assumed that the sources of the waves are located at sufficiently large distance from the core’s midplane to neglect their densitypressure feedback on the cloud. The equations of motion and the induction equation read \arraycolsep1.75ptwhere (4)is the collisional drag coefficient, ⟨ σw ⟩ _{in} is the momentum transfer rate coefficient, and ∇_{ ⊥ } is the gradient along the direction perpendicular to the largescale field B. On the righthand side we have grouped the nonlinear terms: note that because of the assumption of transverse fluctuations they involve only gradients in that direction.
Combining the linearized Eqs. (1)–(3), we obtain the wave equation for u_{i}, (5)where ρ = ρ_{n} + ρ_{i} and V_{Ai}, V_{An} are the Alfvén velocities in the ions and the neutrals, respectively: (6)Substituting in Eq. (5) a solution in the form u_{i}(x,t) = u_{0i}exp [i(kx − ωt)] , we obtain the wellknown dispersion relation (DR) (7)first derived by Piddington (1956) and later by Kulsrud & Pearce (1969).
The amplitudes of the fluctuations in u_{n} and u_{i} are related to the amplitude of the fluctuations of b by (8)and (9)For ρ_{i}/ρ_{n} ≪ 1, there are two separated regimes of wave propagation: long waves (for ) and short waves (for kV_{Ai} ≫ γ_{in}ρ_{n}). The physical reason is that when the timescale for wave damping becomes shorter than the wave period, wave oscillation cannot occur. Short waves are highfrequency waves propagating at the ion’s Alfvén speed, damped by collisions with neutrals. In the shortwave regime, the ratio of the amplitudes of the velocity oscillations reduces to (10)i.e. the neutrals are weakly coupled to the ions and their amplitude of oscillation is very small. Therefore, these waves cannot be responsible for the linewidths of neutral species observed in molecular clouds.
In the limit of long waves, the wave Eq. (5) becomes (11)where V_{A} is the Alfvén speed in the gas of neutrals and ions (of density ρ = ρ_{n} + ρ_{i}), and η_{AD} is the ambipolar diffusion resistivity (12)In this regime, the DR of Eq. (7) reduces to (13)For the propagation and damping of Alfvén waves in a molecular cloud a description in terms of forced waves is appropriate. Indeed, there are several possible sources of hydromagnetic waves in a molecular cloud, such as gravitational collapse of smaller parts, largescale nonradial oscillations of the cloud, radiation pressure from young stars, expansion of HII regions, stellar outflows (Arons & Max 1975). In this case, the DR should be considered as a quadratic equation for the wavenumber k (with real and imaginary part) as a function of the frequency ω (real). Defining a complex wave number k = ℜ(k) + iℑ(k), in the limit , we derive from the real and imaginary parts of the dispersion relation Eq. (13) the expressions (14)The condition for wave propagation,  ℜ(k)  ≥  ℑ(k)  , is satisfied only if the frequency of forcing is lower than a critical frequency , or, the wavelength is longer than a critical value λ_{cr} = πη_{AD}/V_{A} (Kulsrud & Pearce 1969). In the lowionization regime typical of a molecular cloud, ρ ≈ ρ_{n}, the condition for wave propagation reduces to ω ≤ 2γ_{in}ρ_{i} or, in terms of wavelength, to (15)We selected the case of depleted clouds where all metal species are frozen onto grains and the dominant ions are H^{+} and H. In this case γ_{in} = 2.38 × 10^{14} g^{1} cm^{3} s^{1} (see Pinto & Galli 2008) and the mean mass of neutrals and ions are m_{n} = 2.33m_{H} and m_{i} = 3m_{H}, respectively. Inserting numerical values, we obtain (16)where x_{i} = n_{i}/n is the ionization fraction. In the regime of long waves, the amplitudes of the velocity oscillations reduce to (17)Since  k  = ω/V_{A} to lowest order in , the two velocity amplitudes are equal in modulus  u_{0n}  =  u_{0i}  = V_{A}  b_{0}  /B. These longwavelength Alfvén waves propagate with speed V_{A} in a wellcoupled ionneutral gas, and may be responsible for the observed linewidths of ionic and neutral species observed in molecular clouds.
3. Model equations
The condition for the propagation of long waves is equivalent to assume that neutrals and ions are well coupled (u_{i} − u_{n}  /V_{An} ≪  b  /B). With this condition, and the approximation ρ_{i} ≪ ρ_{n}, the ion velocity can be obtained from Eq. (2) as (18)Upon substitution of the above equation into the induction Eq. (3) one obtains a quite complicated equation. To make it more tractable, still accounting for nonlinear interactions, we adopted the following additional assumptions: i) fluctuations have small amplitudes (b/B ≈ u_{n}/V_{a} < 1), thus, when substituting Eq. (18) into (3) we expand the nonlinear terms only to the second order (cubic terms are neglected); ii) the ambipolar diffusion coefficient varies on a scale much larger than all the other quantities do, thus we neglect terms involving its derivatives. Under these assumptions, the induction equation is of the form (19)We see that nonlinearities also appear in terms that contain the ambipolar diffusion coefficient (the last two terms), therefore we expect that the turbulent cascade is also affected by ionneutral drift. In view of the simple phenomenology that we will use (Eq. (26)) to describe the dissipation of turbulence, we additionally simplified the equation by adopting an ordering of the parallel and perpendicular length scales, namely iii) λ_{cr}/ℓ_{   } ≪ b/B ≪ ℓ_{ ⊥ }/ℓ_{   }, where ℓ_{   , ⊥ } are the typical sizes of parallel and perpendicular large scales. We thus obtain the final model equations: (20)and (21)Equations (20) and (21) can be conveniently expressed in terms of the Elsässer variables, (22)by dividing Eq. (21) by , then adding and subtracting Eqs. (20) and (21): (23)The superscripts ′ and ′′ indicate first and secondorder spatial derivatives along the x direction and P_{m} is the magnetic pressure. The first two terms on the lefthand side represent the temporal and spatial wave propagation, respectively. On the righthand side, the first two term account for nonlinear interactions; the term proportional to describes the losses related to the presence of gradients in the Alfvén velocity and wave reflection; the terms proportional to η_{AD} are associated to the dissipation by ambipolar diffusion and combined effects of density gradients and ionneutral drift. Because we assumed fluctuations to lie in the plane perpendicular to B, the gradient ∇_{ ⊥ } involves only perpendicular scales. Given the high Reynolds numbers found in molecular clouds (R_{m} ≈ 10^{8}, Myers & Khersonsky 1995), a perpendicular turbulent cascade is expected to develop. These nonlinear couplings involve only counterpropagating waves. Note however that even if only one mode is present, say z_{ + }, density gradients or ionneutral drift ensure the presence of the other mode (z_{ − }) and hence the triggering of the cascade.
The hypotheses that we have adopted in deriving Eqs. (20), (21), (23) are summarized in the following ordering: (24)Note that the first inequality is a stronger constraint than the long wavelength condition Eq. (15), but in molecular clouds one typically finds that λ_{cr}/ℓ_{} ≈ η_{AD}/V_{A}L ≲ 0.06 (see Eq. (16)) and ρ_{n}/ρ_{i} ≳ 10^{7}, therefore the above ordering is generally satisfied even for ℓ_{⊥}/ℓ_{} > 1. The condition b/B ≪ ℓ_{⊥}/ℓ_{} can be rewritten as τ_{A}/τ_{NL} < 1 where the Alfvén and nonlinear timescales are given by τ_{A} = ℓ_{}/V_{A} and τ_{NL} = ℓ_{⊥}/z_{±} (we assumed b/B ≈ u/V_{A} ≈ z_{±}/V_{a}). In other words, the equations are suited to describe weak turbulence. Our inequality is more stringent than τ_{A}/τ_{NL} ≪ B/b given in Galtier et al. (2002), which was used to build a cascade model. We will not attempt such a detailed description of the cascade for two reasons that are related to the inhomogeneity. It will be clear in the following that the density gradients couple counterpropagating waves. A first consequence is that the phase of oppositely propagating lowfrequency waves cannot be considered as random anymore, preventing us from exploiting a phase averaging method at large scales to obtain the kinetic equations. Moreover, in the inhomogeneous case typical of the cloudcore systems, the density gradients act as semireflecting boundaries, which enhance the amplitudes, and therefore the cascade rate of largescale modes. These lowfrequency waves in which energy is accumulated pose a problem with the continuity of the lowfrequency tail of the spectrum. Our nonlinear model is really only in the spirit of von Karman & Howarth (1938) (see also Mac Low 1999), where the large eddies control the dissipation rate and the specific mechanism for smallscale dissipation is not of central importance. Largescale fluctuations control the cascade, and fast, smallscale dissipation mechanisms, whether collisional or not, act as passive absorbers of cascaded energy.
Another comment concerns the possibility of nonlinear compressible interactions and processes such as parametric decay. These interactions may play a role, but in a completely turbulent situations where fluctuations have completely random phases, the decay via these resonant processes is greatly reduced (Malara & Velli 1996). Consequently, the density fluctuations that are there will play a role via processes such as phase mixing/resonant absorption, and shock heating, but will not alter, except at shocks, the dominant nature of the energy cascade given by the dimensional model described below.
4. Numerical method
We first sought solutions in form of monochromatic waves for the linearized Eq. (23). This allowed us to show the role of interference between waves as produced by reflection and ambipolar diffusion (see Sect. 5). The numerical method is presented in the next section Sect. 4.1. Then we solved the full set of Eq. (23) in their timedependent form, treating the nonlinear interactions in a phenomenological way and considering fluctuations that are characterized by a steep frequency spectrum. The details of the numerical method and the phenomenology adopted are given in Sect. 4.2.
4.1. Timeindependent scheme for linear equations
We neglected nonlinear terms and hence the perpendicular structure of Alfvén waves. We sought solutions for monochromatic waves of the form z_{ ± }(x,t) = V_{A0}y_{ ± }(x)exp( − iωt). For a cloud of typical size L and density ρ_{0}, we can define a characteristic Alfvén speed and Alfvén time t_{A} = L/V_{A0}, and the nondimensional variables ξ = x/L and . In nondimensional form, Eq. (23) can be written as (25)where W = ωL/V_{A0} and . For uniform density () and without dissipation (), Eq. (25) describes the propagation of Alfvén waves traveling to the right ( + sign) and to the left ( − sign).
Because ambipolar diffusion is a small effect in typical ISM conditions (), we adopted a perturbational approach. We initially solved Eq. (25) guessing at a functional form for the secondorder derivatives. In particular, we assumed the exponentially damped behavior predicted by analytical studies (we have verified that the numerical results are independent of the choice of the initial guess); then we used the secondorder derivatives of the solutions obtained by the guess function to integrate Eq. (25) and obtained a new set of secondorder derivatives. This procedure was iterated until convergence of the solution was attained.
Special care has to be taken when choosing boundary conditions suitable to describe waves incident on both sides of an interstellar cloud because of the interference between the incoming waves. When only ambipolar diffusion is present, the iterative procedure reduces the problem to a boundaryvalue problem for the fields z^{ ± }. For wave propagation in a density gradient without dissipation, the boundaryvalue problem is well defined only when one wave vanishes at one boundary. In the more realistic case, where both modes are injected at the boundaries, solutions depend on the relative phase between the injected z_{ ± }. The phase in fact controls the interference between the incident wave and the reflected wave inside the density profile. We will describe in detail the appropriate initial conditions in the application to interstellar clouds (Sect. 6), first dealing with density gradients, then considering both ambipolar diffusion and density gradients.
4.2. Timedependent scheme for nonlinear equations
To evaluate the effect of turbulent dissipation, which can be triggered by nonlinear wave interactions, we did not model the turbulent cascade. Instead we used a phenomenological expression (Dmitruk et al. 2001) that accounts for the damping of the waves and preserves the structure of nonlinear terms, allowing interactions only among counterpropagating waves. We therefore integrated Eq. (23) modified with the following substitution in the RHS, (26)where ℓ_{ ⊥ } represents the dimension of the largest eddies (the energy containing scale in the turbulent field).
The form of the nonlinear term may be heuristically derived in the framework of reduced MHD in which variations along the perpendicular directions are decoupled from those along the magnetic field (∇ = ∇_{ ⊥ } + ∇_{   }, with ∇_{ ⊥ } ≫ ∇_{   }). In presence of a strong magnetic field, the propagation time of the Alfvén waves τ_{A} is equal to or shorter than, the characteristic timescale for nonlinear interaction τ_{NL} over most of the Fourier space (consistent with the ordering in Eq. (24)). The nonlinear cascade becomes anisotropic, developing preferentially in the plane perpendicular to the direction of the mean field (Shebalin et al. 1983; Oughton et al. 2004; Goldreich & Sridhar 1995). One can therefore Fourierdecompose Eq. (23) only in k_{ ⊥ } (wavevector lying in plane perpendicular to the largescale magnetic field). The global effect of this perpendicular cascade can be described by means of two quantities at the large scales, namely a common similarity (correlation) scale and the average velocity difference among points belonging to the same eddy. Identifying these two quantities with the integral turbulent length (ℓ_{ ⊥ }) and the fluctuation’s amplitude of the Elsässer fields, we can construct a characteristic timescale , which accounts for nonlinear turbulent interactions in Eq. (23). This phenomenology assumes that all the interactions are resonant, in practice overestimating the heating rate for high frequency fluctuations, which instead would have a reduced nonlinear coupling because of their fast decorrelation. On the other hand, the phenomenology still accounts for the reduction of the cascade and associated damping due to the imbalance between the z_{ ± } fields.
In our approximation, parallel scales are not generated by the cascade itself, as is expected to happen for strong turbulence (Goldreich & Sridhar 1995). Therefore, here the perpendicular turbulent cascade represents a channel for energy dissipation that is alternative to ambipolar diffusion (assumed to act only on the parallel scale). We finally remark that the ordering in Eq. (24) is assumed to be valid at large scales but it breaks at smaller scales, where the cascade is expected to become strong, eventually reaching the dissipative scale. However, for the purpose of estimating the turbulent heating rate, the phenomenology is still valid and consistent with the ordering since the overall cascade rate will be controlled by the weaker cascade at the modeled large scales.
We integrated the evolution Eq. (23), with the substitution of Eq. (26), using a numerical scheme of first order in time and second order in space (a centered LaxWendroff scheme). Passing to timedependent equations allows us to better model the interference (dependence on the phase) of Alfvén waves (from incoherent sources) that are assumed to enter from both sides of the cloud with random phases. Boundary conditions are transparent, no reflection occurs there. Waves are injected into the cloud according to (27)where the forcing f(t), g(t) has a spectrum peaked at low frequencies and declining as ω^{2}. The escaping waves at the boundaries, ∂z_{ + }/∂t at x = L and ∂z_{ − }/∂t at x = −L, are not specified and obtained as part of the solution.
5. Special cases
To gain insight into the nature of the solution, we first considered the effects of ambipolar diffusion for waves propagating in a uniform medium and then the propagation of waves in a nonuniform medium without ambipolar diffusion.
5.1. Damping by ambipolar diffusion, uniform density
For a uniform density medium, , Eq. (25) reduces to (28)In the limit of small dissipation, setting , Eq. (28) has a general solution of the form (29)where x_{d} = (2ϵW)^{1} is the damping length and and a_{ ± } are constants. In dimensional units, (30)Notice that the damping length can also be expressed as ℓ_{d} = 2πλ^{2}/λ_{cr}.
In this case each wave is damped in the direction of propagation. The solution contains a primary wave traveling in one direction (first term on the RHS), whose damping produces a secondary wave traveling in the opposite direction (second term on the RHS). Notice that when two primary waves are present (a_{ ± } ≠ 0), the profile of u_{n} and b will be determined by their interference in a nontrivial way.
5.2. Wave reflection by density gradients
We then considered a nonuniform medium where density gradients imply changes in the Alfvén speed, but we neglected ambipolar diffusion. Equation (25) becomes (31)In terms of the scaled amplitude and the scaled wavevector , Eq. (31) becomes (32)The general properties of this equation have been discussed in detail by Velli (1993). Here, it is sufficient to notice that for a weak density gradient a Taylor expansion gives q(ξ) ≈ W + q_{1}ξ, which, inserted into Eq. (32), gives the approximate solution (33)where β = q_{1}/4W^{2} and a_{ ± } are constants. Thus, a density gradient produces two effects: a variation of the wave amplitude, as expressed by the ρ^{1/4} scaling (Walén 1944) valid in the WKB approximation, and the generation of secondary reflected waves (second term on RHS). Notice that this approximate solution is similar (waves traveling in opposite directions) to Eq. (29), thus, again, when two primary waves are present (a_{ ± } ≠ 0) the general solution for u_{n} and b will depend on their interference in a nontrivial way.
In the absence of damping, the wave amplitude variation is governed by the conservation of the energy flux S = S_{ + } − S_{ − } = const, where S_{ ± } =  f_{ ± }  ^{2}/8 are the fluxes carried by the individual modes. For a weak gradient, the amplitudes corresponding to the solutions (33) are (34)Thus the density gradient generates reflected waves that induce spatial oscillations in the wave amplitude, with half the wavelength of the incident wave, but of course at the same time conserves the energy flux (). For a given wavelength, the stronger the gradients, the stronger the coupling, and the reflected wave tends to have the same amplitude as the incident wave. For a weak gradient the coupling is negligible and the solution approaches the WKB solution.
Fig. 1 Large panel: normalized amplitudes of velocity fluctuations as a function of ξ = x/L for waves incident on the right side of a core with the density profile given by Eq. (35) for various values of the nondimensional wave frequency W. No dissipation is included. The dashed red line shows the WKB solution. Small panel: the red line shows the normalized density profile, and the black line the normalized Alfvén speed. 

Open with DEXTER 
6. Application to molecular cloud cores
We applied our numerical schemes to the study of the propagation and damping of Alfvén waves in molecular cloud cores. We considered a onedimensional cloud consisting of a lowmass core embedded in a uniform density envelope. For the cloud’s density we assumed a modified isothermal profile, (35)with μ = 2.36, n_{env} = 10^{4} cm^{3}, n_{core} = 10^{5} cm^{3}, and ξ = x/L where L = 0.04 pc. The slab is threaded by a perpendicular, straight and uniform magnetic field of strength 70 μG and we considered the lowionization regime typical of molecular cloud, ρ ≈ ρ_{n}. We assumed a simple ionization law, (36)with κ = 1/2 and C = 4.5 × 10^{17} cm^{−3/2} g^{1/2}. With these values, and γ_{in} = 2.38 × 10^{14} cm^{3} g^{1} s^{1} (Pinto & Galli 2008), the ambipolar diffusion resistivity is (37)The effect of wave reflection is illustrated in the following example, which describes the dependence of the velocity fluctuations on the normalized wave frequency W. We neglected dissipation (assuming η_{AD} = 0 and ℓ_{ ⊥ } = ∞) and assumed that only one outward propagating wave (z_{ + }) is injected at the right boundary. Solutions are obtained integrating backward from ξ = 3 and do not depend on the phase of the injected wave. The resulting velocity amplitudes are shown in Fig. 1 (large panel) for different values of the nondimensional frequency W ∈ [0.6,6] . Clearly, for increasing frequency the solution approaches the WKB analytical solution and the highest frequency case W = 6 is indistinguishable from it. The amplitude profiles show variations on a scale 1/k = V_{A}/ω that changes according to the Alfvén speed profile (shown in the small panel of Fig. 1), being smaller in the denser region  ξ  ≲ 1. This means that small scales are naturally formed in the core region.
The inclusion of diffusive effects requires a careful treatment of the boundary conditions. The choice of a symmetric density profile, with vanishing derivative at the cloud’s center and at infinity, allows one to match the asymptotic solutions with linear combinations of the solutions obtained by retaining only the terms that involve ambipolar diffusion (see Sect. 5.1), which we require to be exponentially damped waves. In practice, we assumed that at x = ∞ (x = −∞) only leftward (rightward) propagating waves are present. The integration was carried out only in one half of the domain and we looked for solutions symmetric with respect to the core’s midplane. Since the gradients vanish at x = 0, the coupling given by reflection is absent and we can impose as boundary condition  z_{+}(0)  =  z_{−}(0)  there. We also averaged over eight phase differences (Δφ = nπ/4 with n = 0,1,...,7) between z_{+}(0) and z_{−}(0) to mimic the effect of incoherent wave sources at infinity. Results are shown in Fig. 2 for the range of normalized frequencies 0.6 < W < 6, and over a spatial interval corresponding to roughly three times the core size. The amplitudes of the waves were normalized to the values in ξ = 0. The frequency interval explored corresponds approximately to the wavelength interval L ≤ λ ≤ 10L in the core. In the highfrequency regime (W = 6), the damping length is comparable to the scale of the density gradient (x_{d} = 1.4 at ξ = 0), and the dominant effect is the quasiexponential damping of the waves on a scale ~ x_{d}. On the other hand, in the lowfrequency regime (W = 0.6) the evolution of the wave amplitude is dominated by the the density gradient on a scale ξ ~ 1, because the damping length is much longer than the core’s size (x_{d} = 140 at ξ = 0).
Fig. 2 Normalized amplitudes of z_{ − } (left panel ) and z_{ + } (right panel ) as a function of ξ = x/L in a cloud core with the density profile given by Eq. (35) for values of the nondimensional wave frequency from W = 0.6 to W = 6. 

Open with DEXTER 
Figure 3 shows the corresponding profiles of velocity fluctuations for the same range of frequencies as in Fig. 2. Again an average procedure was performed to take into account relative phases of uncorrelated waves moving in opposite directions from distant sources. As in the case of the amplitude of the Elsässer variables, the velocity amplitudes are mostly sensitive to the density gradient at lower frequencies W ~ 0.6 and mostly sensitive to the ambipolar damping at the highest frequencies (W ~ 6). This is easy to understand given the dependence of the ambipolar damping length on the inverse square of the frequency (Eq. (30)). Despite the simplicity and idealization of the model, it is clear that the sharp gradients of velocity dispersion observed in the environments of dense cores (Barranco & Goodman 1998; Caselli et al. 2002; Pineda et al. 2010) cannot be produced by a density gradient alone, but require the combined effects of a density gradient and ambipolar diffusion. In particular, sharp gradients in velocity amplitude are only possible when x_{d} ≈ 1, which implies waves with a frequency such as to give ℓ_{d} ≈ L. These waves have wavelength λ ~ (λ_{cr}L)^{1/2}, intermediate between the critical wavelength λ_{cr} and the typical scale of the density gradient L.
Fig. 3 Normalized amplitude of velocity fluctuations as a function of ξ = x/L in a cloud core with the density profile given by Eq. (35) for values of the nondimensional wave frequency from W = 0.6 to W = 6. 

Open with DEXTER 
6.1. Heating rate
The energy generated by the friction between streaming particles is a significant source of heating for molecular clouds (Scalo 1977; Lizano & Shu 1987). The energy produced heats the bulk of the gas, increasing its pressure and allowing chemical reactions to proceed at a faster rate (Flower et al. 1985; Pineau des Forêts et al. 1992). Therefore, it is important to evaluate the heating rate associated to the dissipation of the waves considered in the previous sections.
In a weakly ionized medium, the rate of energy generation by ambipolar diffusion (per unit time and unit volume) is (38)(Braginskii 1965; see also Pinto et al. 2008). In the special case of lowamplitude Alfvén waves propagating in a onedimensional medium, Eq. (38) reduces to (39)where the brackets indicate the average over time and the prime a space derivative.
The heating rate as a function of position for waves propagating in the density profile given by Eq. (35) is shown in Fig. 4, as before for waves in the frequency range 0.6 ≤ W ≤ 6. To show the results in physical units, we have assumed b = B. Therefore, the values shown in the figure should be considered as upper limits to the actual value of the heating rate. For comparison, the red curve shows the cosmicray heating rate computed assuming a cosmicray ionization rate 3 × 10^{17} s^{1} and a mean energy deposited per ionization ΔQ = 20 eV (Goldsmith 2001). The ambipolar diffusion heating is clearly important both in the core and in the envelope, and competes with cosmic rays as a heating source for the cloud, especially for the highest values of the frequency. The behavior of the heating rate deviates significantly from the profile obtained assuming b′ ~ kb (see, e.g., Hennebelle & Inutsuka 2006), because of the modulation of the wave amplitude by the density gradient, and the space dependence of the damping length.
Fig. 4 Ambipolar diffusion heating as a function of the normalized distance from the core center. For comparison, the red line shows the cosmicray heating rate for a cosmicray ionization rate 3 × 10^{17} s^{1} and energy deposition ΔQ = 20 eV. 

Open with DEXTER 
Nondimensional parameters for the runs: input amplitude (b/B) and ratio of the perpendicular to parallel scales of fluctuations (ℓ_{⊥}/ℓ_{}).
6.2. Energy dissipation by turbulent cascade
In the previous section we have considered the dissipation of linear Alfvén waves caused by ionneutral collisions. Here we estimate the effect of turbulent dissipation as given by the phenomenological expression (26). We decided to switch off ambipolar diffusion to better explore the effect of turbulence, varying its strength through the control parameter τ_{A}/τ_{NL}. Setting η_{AD} = 0 we set λ_{cr} = 0 and we did not have to worry about the first inequality in Eq. (24). We imposed at both sides of the cloud incoming fluctuations characterized by a frequency spectrum of width 1/t_{A} and then falling as ω^{2}: energy is at long wavelengths and ℓ_{} ≈ V_{A0}t_{A} = L. We here discuss four runs in which we varied the turbulent strength τ_{A}/τ_{NL} by changing the input amplitude and perpendicular scale as listed in Table 1. The strong turbulence cases are runs A and D, the last even though it has a smaller input amplitude. The timeaveraged, turbulent heating rate per unit volume follows directly from Eq. (26) and is given by (40)where (41)is the (normalized) cross helicity, E =  z_{ + }  ^{2} +  z_{ − }  ^{2} is the total turbulent energy, and is a symmetric function that vanishes at  σ_{c}  = 1 and increases monotonically, reaching a maximum f(0) = 1.
The small amount of energy at high frequency justifies neglecting ambipolar diffusion, we verified that ⟨ H ⟩ ≪ H_{T} for .
Fig. 5 Timeaveraged normalized energy fluxes entering the cloud from both sides, S_{+}/S_{+}(ξ = −5) and S_{−}/S_{−}(ξ = 5), for run B. The normalized density profile is shown for comparison by the dotted line. 

Open with DEXTER 
Let us examine the case of intermediate turbulence strength, (run B). In the linear case with η_{AD} = 0, the coupling between z_{±} in proximity of the density gradients causes an enhancement of the respective energy fluxes S_{ ± }, which have a bump toward the cloud’s center. This ensures the conservation of the net energy flux (S_{+} − S_{−}), which is uniform. When nonlinear dissipation is taken into account, the profiles of S_{+} and S_{−} are remarkably different, as shown in Fig. 5 where we also plot the density profile for comparison (dotted line). Both energy fluxes decrease toward the cloud’s center, indicating the damping of the wave (heating). In contrast to the linear case, there is a sharp decrease on the core size due to the turbulent heating H_{T}, which is stronger in the whole core (being approximately proportional to ρ_{n}). From Fig. 5 one can also deduce the behavior of σ_{c}, which is also equal to (S_{ + } − S_{ − })/(S_{ + } + S_{ − }) (the fluxes are normalized to their respective incident value, since S_{ + }(ξ = −5) ≈ S_{ − }(ξ = 5) the normalization factors can be considered equal). The normalized cross helicity vanishes at the core’s center and increases rapidly in the core’s envelope, maintaing a value  σ_{c}  ≳ 0.8 at the edges of the core. The cross helicity and the total turbulent energy affect the average heating rate per unit mass (H_{T}/ρ_{n}), which is plotted in Fig. 6 for the four runs (run B corresponding to the dashed line). For weak turbulence (runs B and C) the profile is basically flat, varying by 20% at maximum, and thus reflecting H_{T} ∝ ρ_{n}. Interestingly, it is flat on a region entirely surrounding the core and peaks outside it, showing that the wave damping is maximal when passing from the cloud to the core. The relatively flat profile can be explained in terms of a compensation of the terms related to the crosshelicity and total turbulent energy in Eq. (40). Indeed, the cascade is more balanced (σ_{c} ≈ 0) at the center, but less energetic, while the turbulent energy grows outside the core, but the cascade is less balanced. For stronger turbulence (e.g. τ_{A}/τ_{NL} = 1) these two effects do not compensate each other anymore and the minima and maxima in the heating rate are more pronounced (H_{T} follows the scaling with ρ_{n} less closely). Note that keeping τ_{A}/τ_{NL} = 1 but varying the input amplitudes b/B = ℓ_{ ⊥ }/ℓ_{   } = 1/10 the profile does not vary, but the overall dissipation decreases (compare runs A and D). One can verify that the average dissipation per unit mass scales as , while the relative excursion at the core edges ΔH_{T}/H_{T} grows with τ_{NL}/τ_{A}.
Fig. 6 Profile of the timeaveraged heating rate per unit mass (in nondimensional units) for different turbulent strength τ_{A}/τ_{NL} = 1, 1/, corresponding to runs A, B, C, D as labeled on each curve (solid, dashed, dotdashed and longdashed lines, respectively). The normalized density profile is shown for comparison as a dotted line. 

Open with DEXTER 
7. Conclusions
We have considered the propagation of Alfvén waves in a density profile typical of a molecular cloud core embedded in a lowdensity medium. The formulation of the problem in terms of Elsässer variables was found particularly appropriate to treat the effects of wave reflection and dissipation. We have studied with timedependent as well as timeindependent numerical codes the interaction of incoherent wave trains incident on the opposite edges of a onedimensional model cloud, determining the amplitude variations associated to the reflection and damping of the waves.
We found that for waves with a wavelength longer than the critical wavelength, of the order of a few 10^{3} pc in a molecular cloud core, but shorter or of the order of the size of the core, a significant decrease in the velocity amplitude is produced by the combined effects of reflection at the core’s boundaries and dissipation due to ambipolar diffusion in the core’s interior. For these waves, the damping length associated to ionneutral collisions is of the same order as the size of the core. The resulting behavior of the amplitude of velocity fluctuations may explain, at least in part, the sharp decrease of line widths observed in the environments of lowmass cores (see Sect. 1). Moreover, for these waves the energy generated by ionneutral drift is a significant source of heating both in the core and its envelope (see Sect. 6.1).
We also considered the turbulent dissipation and found that the average heating per unit mass scales as , thus becoming significant for high enough wave amplitudes (0.1 < b/B < 1). This regime is particularly relevant for molecular clouds and starforming regions where a statistical analysis of dust polarization maps suggests that the ratio of the strength of turbulent magnetic field to the largescale mean field is of the order of 0.1–0.5 (Hildebrand et al. 2009; Houde et al. 2009). We used a phenomenological term to mimic the effect of a perpendicular turbulent cascade. The results of our calculations show that the turbulent cascade affects the amplitude of fluctuations both at the cloud’s edge and in the cloud’s interior. In particular, the energy dissipation (per unit mass) associated to this process is peaked in the regions immediately surrounding the dense central core, indicating a stronger damping of the wave amplitude. The relative excursion between maxima and minima ΔH_{T}/H_{t} depends on the turbulence strength τ_{NL}/τ_{A}, so it can be significant even for small amplitudes, provided the turbulent correlation length is of the order of a few times the size of the core. In a more realistic situation, the turbulent cascade, although it is stronger in the perpendicular direction, also generates small parallel scales. Therefore we expect that turbulence and ambipolar diffusion will be at work at the same time, enhancing the sharp variation of velocity amplitude at the interface between the cloud and the core. In the present work the investigation of the parameter space was limited by the several assumptions (Eq. (24)) we made to obtain the model equations. Interesting features appear when we are close to the limit of the allowed parameter space. Additional work is needed to address these issues with a more realistic and detailed modeling of turbulence that includes the effects of ambipolar diffusion.
Acknowledgments
We would like to thank Stella Galligani for carrying out some of the nonlinear simulations with the phenomenological model of turbulence.
References
 Arons, J., & Max, C. E. 1975, ApJ, 196, L 77 (In the text)
 Barranco, J. A., & Goodman, A. A. 1998, ApJ, 504, 207 [NASA ADS] [CrossRef] (In the text)
 Braginskii, S. 1965, in Rev. Plasma Physics, ed. M. Leontovich (New York: Consultants Bureau), 1, 205 (In the text)
 Campos, L. M. B. C., & Isaeva, N. L. 1993, Adv. Space Res., 13, 123 [NASA ADS] [CrossRef] (In the text)
 Caselli, P., Benson, P. J., Myers, P. C., & Tafalla, M. 2002, ApJ, 572, 238 [NASA ADS] [CrossRef] (In the text)
 Coker, R. F., Rae, J. G. L., & Hartquist, T. W. 2000, A&A, 360, 290 [NASA ADS] (In the text)
 Dobrowolny, M., Mangeney, A., & Veltri, P. 1980, Phys. Rev. Lett., 45, 144 [NASA ADS] [CrossRef] (In the text)
 Dmitruk, P., Milano, L. J., & Matthaeus, W. H. 2001, ApJ, 548, 482 [NASA ADS] [CrossRef] (In the text)
 Elsässer, W. M. 1950, Phys. Rev., 79, 153 [NASA ADS] [CrossRef] (In the text)
 Flower, D. R., Pineau des Forêts, G., & Hartquist, T. W. 1985, MNRAS, 216, 775 [NASA ADS] [CrossRef] (In the text)
 Galtier, S., Nazarenko, S. V., & Newell, A. C. 2002, ApJ, 564, 49 [NASA ADS] [CrossRef] (In the text)
 Goldreich, P., & Sridhar, S. 1995, ApJ, 438, 763 [NASA ADS] [CrossRef] (In the text)
 Goldsmith, P. F. 2001, ApJ, 557, 763 [NASA ADS] [CrossRef] [MathSciNet] (In the text)
 Goodman, A. A., Barranco, J. A., Wilner, D. J., & Heyer, M. H. 1998, ApJ, 504, 223 [NASA ADS] [CrossRef] (In the text)
 Hennebelle, P., & Inutsuka, S.I. 2006, ApJ, 647, 404 [NASA ADS] [CrossRef] (In the text)
 Heitsch, F., & Burkert, A. 2002, in Modes of Star Formation and the Origin of Field Populations, eds. E. K. Grebel, & W. Brandner, ASP Conf. Ser., 285, 13 (In the text)
 Hildebrand, R. H., Kirby, L., Dotson, J. L., Houde, M., & Vaillancourt, J. E. 2009, ApJ, 696, 567 [NASA ADS] [CrossRef] (In the text)
 Hollweg, J. V. 1981, Sol. Phys., 70, 25 [NASA ADS] [CrossRef] (In the text)
 Houde, M., Vaillancourt, J. E., Hildebrand, R. H., Chitsazzadeh, S., & Kirby, L. 2009, ApJ, 706, 1504 [NASA ADS] [CrossRef] (In the text)
 Kudoh, T., & Basu, S. 2003, ApJ, 595, 842 [NASA ADS] [CrossRef] (In the text)
 Kudoh, T., & Basu, S. 2006, ApJ, 642, 270 [NASA ADS] [CrossRef] (In the text)
 Kulsrud, R., & Pearce, W. P. 1969, ApJ, 156, 445 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] (In the text)
 Lizano, S., & Shu, F. H. 1987, in Physical Processes in Interstellar Clouds, NATO ASIC Proc., 210, 173 (In the text)
 Mac Low, M.M. 1999, ApJ, 169, 178 (In the text)
 Malara, F., & Velli, M. 1996, J. Geophys. Res., 92, 7363 (In the text)
 Marsch, E., & Mangeney, A. 1987, Phys. Plasmas, 3, 4427 (In the text)
 Martin, C. E., Heyvaerts, J., & Priest, E. R. 1997, A&A, 326, 1176 [NASA ADS] (In the text)
 Myers, P. C., & Khersonsky, V. K. 1995, ApJ, 442, 186 [NASA ADS] [CrossRef] (In the text)
 Nakano, T., Nishi, R., & Umebayashi, T. 2002, ApJ, 573, 199 [NASA ADS] [CrossRef] (In the text)
 Ormel, C. W., Paszun, D., Dominik, C., & Tielens, A. G. G. M. 2009, A&A, 502, 845 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Oughton, S., Dmitruk, P., & Matthaeus, W. H. 2004, Phys. Plasmas, 11, 2214 [NASA ADS] [CrossRef] (In the text)
 Pandey, B. P., Vladimirov, S. V., & Samarian, A. 2008, Phys. Plasmas, 15, 053705 [NASA ADS] [CrossRef] (In the text)
 Piddington, J. H. 1956, MNRAS, 116, 314. [NASA ADS] [CrossRef] (In the text)
 Pineda, J. E., Goodman, A. A., Arce, H. G., et al. 2010, ApJ, 712, L116 [NASA ADS] [CrossRef] (In the text)
 Pineau des Forêts, G., Roueff, E., & Flower, D. R. 1992, MNRAS, 258, 45 [NASA ADS] [CrossRef] (In the text)
 Pinto, C., & Galli, D. 2008, A&A, 484, 17 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Pinto, C., Galli, D., & Bacciotti, F. 2008, A&A, 484, 1 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Scalo, J. M. 1977, ApJ, 213, 705 [NASA ADS] [CrossRef] (In the text)
 Shebalin, J. V., Matthaeus, W. H., & Montgomery, D. 1983, J. Plasma Phys., 29, 525 [NASA ADS] [CrossRef] (In the text)
 Tang, Y.W., Ho, P. T. P., Koch, P. M., et al. 2009, ApJ, 700, 251 [NASA ADS] [CrossRef] (In the text)
 Velli, M. 1993, A&A, 270, 304 [NASA ADS] (In the text)
 Verdini, A., Velli, M., Matthaeus, W.H., et al. 2010, ApJ, 708, 116 [NASA ADS] [CrossRef] (In the text)
 von Kármán, T., & Howarth, L. 1938, Proc. Roy. Soc. London Ser. A, 164, 192 [NASA ADS] [CrossRef] (In the text)
 Walén, C. 1944, Arkiv. Mat. Astr. Fys., 30A, 15 (In the text)
 Wardle, M., & Ng, C. 1999, MNRAS, 303, 239 [NASA ADS] [CrossRef] (In the text)
 Zweibel, E. G., & Josafatsson, K. 1983, ApJ, 270, 511 [NASA ADS] [CrossRef] (In the text)
All Tables
Nondimensional parameters for the runs: input amplitude (b/B) and ratio of the perpendicular to parallel scales of fluctuations (ℓ_{⊥}/ℓ_{}).
All Figures
Fig. 1 Large panel: normalized amplitudes of velocity fluctuations as a function of ξ = x/L for waves incident on the right side of a core with the density profile given by Eq. (35) for various values of the nondimensional wave frequency W. No dissipation is included. The dashed red line shows the WKB solution. Small panel: the red line shows the normalized density profile, and the black line the normalized Alfvén speed. 

Open with DEXTER  
In the text 
Fig. 2 Normalized amplitudes of z_{ − } (left panel ) and z_{ + } (right panel ) as a function of ξ = x/L in a cloud core with the density profile given by Eq. (35) for values of the nondimensional wave frequency from W = 0.6 to W = 6. 

Open with DEXTER  
In the text 
Fig. 3 Normalized amplitude of velocity fluctuations as a function of ξ = x/L in a cloud core with the density profile given by Eq. (35) for values of the nondimensional wave frequency from W = 0.6 to W = 6. 

Open with DEXTER  
In the text 
Fig. 4 Ambipolar diffusion heating as a function of the normalized distance from the core center. For comparison, the red line shows the cosmicray heating rate for a cosmicray ionization rate 3 × 10^{17} s^{1} and energy deposition ΔQ = 20 eV. 

Open with DEXTER  
In the text 
Fig. 5 Timeaveraged normalized energy fluxes entering the cloud from both sides, S_{+}/S_{+}(ξ = −5) and S_{−}/S_{−}(ξ = 5), for run B. The normalized density profile is shown for comparison by the dotted line. 

Open with DEXTER  
In the text 
Fig. 6 Profile of the timeaveraged heating rate per unit mass (in nondimensional units) for different turbulent strength τ_{A}/τ_{NL} = 1, 1/, corresponding to runs A, B, C, D as labeled on each curve (solid, dashed, dotdashed and longdashed lines, respectively). The normalized density profile is shown for comparison as a dotted line. 

Open with DEXTER  
In the text 