ABSTRACT
The dynamics of the interstellar medium (ISM) are strongly affected by turbulence, which shows increased anisotropy in the presence of a magnetic field. We expand upon the Esquivel & Lazarian method to estimate the Alfvén Mach number using the structure function anisotropy in velocity centroid data from Position–Position–Velocity maps. We utilize three-dimensional magnetohydrodynamic simulations of fully developed turbulence, with a large range of sonic and Alfvénic Mach numbers, to produce synthetic observations of velocity centroids with observational characteristics such as thermal broadening, cloud boundaries, noise, and radiative transfer effects of carbon monoxide. In addition, we investigate how the resulting anisotropy-Alfvén Mach number dependency found in Esquivel & Lazarian might change when taking the second moment of the Position–Position–Velocity cube or when using different expressions to calculate the velocity centroids. We find that the degree of anisotropy is related primarily to the magnetic field strength (i.e., Alfvén Mach number) and the line-of-sight orientation, with a secondary effect on sonic Mach number. If the line of sight is parallel to up to ≈45 deg off of the mean field direction, the velocity centroid anisotropy is not prominent enough to distinguish different Alfvénic regimes. The observed anisotropy is not strongly affected by including radiative transfer, although future studies should include additional tests for opacity effects. These results open up the possibility of studying the magnetic nature of the ISM using statistical methods in addition to existing observational techniques.
Export citation and abstract BibTeX RIS
1. INTRODUCTION
The current picture of the interstellar medium (ISM) vitally includes magnetohydrodynamic (MHD) turbulence acting on scales ranging from kiloparsecs to sub-AU (see Armstrong et al. 1995; Elmegreen & Scalo 2004). This is in part due to the fact that MHD turbulence is of key importance for fundamental astrophysical processes, e.g., heat transport, star formation, Galactic pressure support, magnetic reconnection, and the acceleration of cosmic rays.
MHD turbulence is notoriously difficult to study both observationally and theoretically (see Elmegreen & Scalo 2004 for further discussion). In light of this, numerical simulations have tremendously influenced our understanding of the physical conditions and statistical properties of MHD turbulence (see Mac Low & Klessen 2004; Ballesteros-Paredes et al. 2007; McKee & Ostriker 2007 and references therein). Present codes can produce simulations that resemble observations in terms of structures and scaling laws, but because of their limited numerical resolution, they cannot reach the observed Reynolds4 numbers of the ISM.
Statistical studies represent the best hope to bridge the gap between simulations and observations. Thus, many techniques beyond the traditional turbulence power spectrum have been developed to study and parameterize observational magnetic turbulence. These include higher order spectra, such as the bispectrum (Burkhart et al. 2009), higher order statistical moments (Kowal et al. 2007; Burkhart et al. 2010), topological techniques (such as genus; see Chepurnov & Lazarian 2009), clump and hierarchical structure algorithms (such as dendrograms; see Goodman et al. 2009; Burkhart et al. 2013a), principle component analysis (i.e., PCA; Heyer & Schloerb 1997; Heyer et al. 2008), Tsallis function studies for ISM turbulence (Esquivel & Lazarian 2010; Tofflemire et al. 2011), Velocity Channel Analysis and Velocity Coordinate Spectrum (Lazarian & Pogosyan 2004, 2006, 2008), and structure/correlation functions as tests of intermittency and anisotropy (Cho & Lazarian 2003; Esquivel & Lazarian 2005; Kowal & Lazarian 2010).
However, the results of these statistical applications to numerics and/or observations are less insightful without being placed into the theoretical framework of turbulence. The famous Kolmogorov (1941) theory of turbulence describes the hydrodynamic counterpart of MHD turbulence. The transfer of energy from large-scale eddies to smaller scales continues without losses until the cascade reaches eddies that are small enough to dissipate energy over an eddy turnover time. In terms of the ISM, the injection scale and main energy sources are still unknown, but it is clear that turbulence in the Galaxy is driven on large scales (kiloparsec) by supernova, galactic fountain, high-velocity cloud impacts, hydrodynamical and magnetohydrodynamical instabilities, or some combination of these. At small scales, one should see the scales corresponding to sinks of energy, i.e., dissipation of energy.
The ISM is also magnetized, and therefore Alfvénic perturbations are vital to the development of an MHD cascade. Contrary to Kolmogorov turbulence, in the presence of a dynamically important magnetic field, turbulent eddies become elongated along the mean magnetic field (i.e., they become anisotropic) and Alfvénic perturbations develop an independent cascade which proceeds perpendicular to the local magnetic field and is marginally affected by the fluid compressibility (see Cho & Lazarian 2003). The dynamic influence of the magnetic field and the induced anisotropy of the eddies increases as the cascade proceeds down to smaller scales. This corresponds to the predictions of the Goldreich & Sridhar (1995, henceforth GS95) theory of Alfvénic turbulence.
It is important to stress that the above picture of magnetized turbulence is developed in the context of the local magnetic field relative to the eddies. The anisotropy of interstellar turbulence that is accessible to observations, which are averaged along the line of sight (LOS), is sampled in the global reference frame relative to the large-scale mean magnetic field. In this case, the anisotropy is determined by the largest-scale eddies. The first discussion of the possibilities of observationally measuring the large scale anisotropy induced by the magnetic field was in Lazarian et al. (2002, henceforth LPE02), who proposed to measure contours of equal correlation corresponding to data within different velocity channel thickness in H i data (i.e., in a procedure similar to the Velocity Channel Analysis). Follow-up papers by Esquivel & Lazarian (2005) and Esquivel & Lazarian (2011, henceforth EL11) showed that velocity centroids can be used for testing whether turbulence is sub- or super-Alfvénic, while Heyer et al. (2008) used PCA to recover an empirical relationship for the anisotropy found in simulations and applied this to molecular cloud observations.
EL11 developed a method to quantify the large-scale anisotropy of the turbulent cascade in velocity centroid maps of Position–Position–Velocity (PPV) data cubes. In this paper, we expand upon the EL11 method in terms of the parameter space studied and the applicability of their method to the observations. EL11 showed that this method is highly sensitive to the global Alfvénic Mach number of turbulence, defined as MA ≡ 〈VL/vA〉, as well as the direction of the magnetic field. We investigate a similar approach, but now use higher-resolution simulations and consider more realistic synthetic observations. Furthermore, we study the relation of the isotropy degree with the Alfvén Mach number (MA) and the sonic Mach number (Ms = VL/cs, where cs is the sound speed) for different LOS orientation angles, which can be useful for estimating the three-dimensional (3D) structure of the magnetic field in the ISM. The paper is organized as follows. In Section 2, we describe the GS95 model for sub- and super-Alfvénic turbulence and how this applies to the EL11 method. In Section 3, we describe our extended database of MHD simulations and describe in detail our procedure to create synthetic observations and calculate the isotropy degree. In Section 4, we describe our results followed by the discussion in Section 5 and conclusions in Section 6.
2. THE GS95 ANISOTROPY FOR SUB-ALFVÉNIC AND SUPER-ALFVÉNIC TURBULENCE
The GS95 theory assumes the injection of energy at scale L and the injection velocity equal to the Alfvén velocity in the fluid VA, i.e., the Alfvén Mach number MA ≡ (VL/VA) = 1 (i.e., trans-Alfvénic turbulence), where VL is the injection velocity. The GS95 model was later generalized for both sub-Alfvénic, i.e., MA < 1, and super-Alfvénic, i.e., MA > 1, cases (see Lazarian & Vishniac 1999 and Lazarian 2006) and thus the results of EL11 and the current work must be understood in this context.
For the eddies perpendicular to the magnetic field, the original Kolmogorov energy scaling is applicable resulting in perpendicular motions scaling as , where l⊥ denotes eddy scales measured perpendicular to the local magnetic field. Mixing motions induce Alfvénic perturbations that determine the parallel size of the magnetized eddy. This concept of critical balance, i.e., the equality of the eddy turnover time (l⊥/vl) and the period of the corresponding Alfvén wave ∼l||/VA, where l|| is the parallel eddy scale and VA is the Alfvén velocity. Making use of , one finds the scaling relation for the parallel and perpendicular eddies as . This reflects the scale-dependent anisotropy of eddies along the magnetic field lines as the energy cascades proceeds to smaller scales and has been tested using second order structure functions in the reference frame to the local magnetic field (see Cho & Lazarian 2003; Beresnyak et al. 2005; Kowal & Lazarian 2010).
The EL11 method takes advantage of the global anisotropy observed in the largest-scale eddies that can be measured via structure function analysis. The first mention of the use of the structure function anisotropy technique to study turbulence and the direction of the mean magnetic field was made by LPE02, who used synthetic spectral line emission maps obtained via MHD turbulence simulations to demonstrate the method's promise. Later studies (e.g., Esquivel et al. 2003; Vestuto et al. 2003; Heyer et al. 2008) confirmed that the anisotropy is evident from two-point statistics, i.e., the structure function of observational quantities such as velocity centroids. While EL11 studied the relation between the anisotropy and the global Alfvénic Mach number, they did so with limited resolution simulations and less attention to observational effects such as thermal broadening and unknown LOS angle relative to the mean magnetic field. We expand the parameter range of their study and include additional measures and observational considerations.
3. MHD SIMULATIONS AND STRUCTURE FUNCTIONS OF SYNTHETIC VELOCITY CENTROID MAPS
We generate 3D numerical simulations of isothermal compressible (MHD) turbulence by using the Cho & Lazarian (2003) MHD code and varying the input values for the sonic and Alfvénic Mach number. Turbulence is driven with large-scale solenoidal forcing. The magnetic field has contributions from a uniform background field and a fluctuating turbulent field: B = B0 + b. Initially b = 0. The simulations have resolutions of either 5123 or 2563 and the models are run for t ∼ 5 crossing times, to guarantee full development of the energy cascade. For more details, see Cho & Lazarian (2003), Kowal et al. (2007), Burkhart et al. (2009), and EL11.
We partition our models into three groups corresponding to their Alfvénic Mach number, which covers sub-Alfvénic (B0 = 5.0, 3.0) to trans-Alfvénic (B0 = 1.0) to super-Alfvénic (B0 = 0.1) turbulence. The initial conditions were defined with ρ = 1 and the Alfvén speed . The simulations were evolved to reach a stationary state with the rms velocity close to unity (vrms ∼ 0.7). For each group we compute several models with different values of the sonic Mach number (see Table 1, second column). The models are listed and described in Table 1, where 〈Pgas, 0〉 and B0 represent the initial gas pressure and magnetic field, respectively. The labels given by EL11 for six models (M1–M3 and M7–M9) with similar initial conditions as here are indicated in parentheses. We note that the units on these quantities are given in dimensionless code units. For a detailed discussion on how to convert code units into physical units, see Appendix A of Hill et al. (2008).
Table 1. Parameters of the MHD Simulations
Model | B0 | 〈Pgas, 0〉 | Ms | MA | Resolution |
---|---|---|---|---|---|
L1 | 0.1 | 0.0049 | ∼8.7 | ∼6.1 | 5123 |
L2(M1) | 0.1 | 0.01 | ∼5.7 | ∼5.7 | 5123 |
L3 | 0.1 | 0.025 | ∼3.5 | ∼5.5 | 5123 |
L4 | 0.1 | 0.05 | ∼2.5 | ∼5.6 | 5123 |
L5(M2) | 0.1 | 0.1 | ∼1.8 | ∼5.8 | 5123 |
L6(M3) | 0.1 | 1.0 | ∼0.6 | ∼5.8 | 5123 |
L7 | 0.1 | 2.0 | ∼0.4 | ∼5.4 | 5123 |
L8 | 1.0 | 0.0049 | ∼8.0 | ∼0.6 | 5123 |
L9 | 1.0 | 0.0077 | ∼6.3 | ∼0.6 | 5123 |
L10(M7) | 1.0 | 0.01 | ∼5.5 | ∼0.5 | 5123 |
L11 | 1.0 | 0.025 | ∼3.5 | ∼0.6 | 5123 |
L12 | 1.0 | 0.05 | ∼2.6 | ∼0.6 | 5123 |
L13(M8) | 1.0 | 0.1 | ∼1.8 | ∼0.6 | 5123 |
L14(M9) | 1.0 | 1.0 | ∼0.6 | ∼0.6 | 5123 |
L15 | 1.0 | 2.0 | ∼0.4 | ∼0.6 | 5123 |
L16 | 3.0 | 0.01 | ∼10.2 | ∼0.3 | 2563 |
L17 | 3.0 | 0.1 | ∼3.1 | ∼0.3 | 2563 |
L18 | 3.0 | 1.0 | ∼1.0 | ∼0.3 | 2563 |
L19 | 5.0 | 0.01 | ∼8.7 | ∼0.2 | 2563 |
L20 | 5.0 | 0.1 | ∼2.6 | ∼0.2 | 2563 |
L21 | 5.0 | 1.0 | ∼0.8 | ∼0.2 | 2563 |
Note. Labels in parentheses are those given by EL11 for simulations with similar initial conditions.
Download table as: ASCIITypeset image
We create synthetic PPV cubes and the corresponding velocity centroid maps assuming a fully optically thin media with the emissivity being proportional to the density. EL11 computed the velocity centroid maps directly from the simulated cubes as
This form, henceforth known as the "ideal centroid," is convenient to be applied directly on the simulation data described in Section 3, but not on observed PPV cubes. In this work, we consider the velocity centroid which is the definition applied to observed PPV cubes. We first computed a synthetic PPV cube from the intensity distribution I(X) ≡ ∫ρsdVLOS, where ρs is the density of emitters in the PPV space (i.e., the intensity values of the 3D PPV data cube), X denotes the position on the plane of the sky, and VLOS is the LOS velocity-axis. The integral is made in velocity along the entire LOS, applied at every position in the plane of the sky. For instance, in the expression in Equation (1), the LOS coincides with the z-axis, and X = (x, y).
We compute the structure function as SF(r) = 〈[f(x) − f(x + r)]2〉. We denote the structure function of velocity centroids obtained from the PPV cube with LOS along the x-axis (which, in our simulations, is parallel to the mean magnetic field) as SFC, x(R). SFC, y(R) and SFC, z(R) denote the structure functions of the velocity centroids of PPV data with LOS along y and z axes, respectively (both of which are perpendicular to the mean magnetic field). In our application, the structure function is a two-dimensional function with isocontours that are approximately circular for isotropic Kolmogorov-type turbulence and elliptical for anisotropic turbulence. We present an example case in Figure 1 comparing models L7 and L15 from Table 1. These two simulations have the same sonic Mach number but almost an order of magnitude difference in Alfvén Mach number. Model L7 (with MA = 5.8) has an isotropic (circular) structure function, while model L15 (with MA = 0.6) has and anisotropic (elliptical) structure function. As was discussed in Section 2, the anisotropy indicates the presence of a magnetic field and increases with increasing magnetic field. In order to quantify this effect observed in the two-dimensional structure function of the velocity centroid maps, we define the isotropy degree as being the ratio of the structure functions in two perpendicular directions to the LOS, intersecting at the distribution center
In what follows, we will explore the relation of the isotropy degree with MA and Ms and compare our results with EL11.
4. RESULTS
In order to illustrate the general trends, we selected 12 simulations (3 different sonic Mach number) for the four different values of magnetic field from our simulation parameter space presented in Table 1. We plot the isotropy degree versus the spatial separation (r) in Figure 2 for the different models. Vertical lines show our range of r values for obtaining the average isotropy degree (shown in Figure 3). Below a 10 grid point scale, the density from the MHD simulations is affected by numerical diffusion, and the effect of noise is more pronounced. Past 100 grid points, the simulation is dominated by the injection scale of the turbulence. We find that in nearly all cases the anisotropy is virtually scale-independent from the small scales up to separations on the order of one-fifth of the computational box (about half the size of the injection scale). We note that scale independence should not exist in the local frame of reference to the magnetic field. However, because we only sample the global frame (large-scale eddies) we observe no scale dependency in the anisotropy. We consider other lines of sight in the next section.
Download figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution imageFigure 2 shows a clear separation of the isotropy degree for simulations with different value of magnetization across a range of spatial scales. The simulations generally cluster in isotropy degree around three Alfvénic regimes: high magnetization (B = 3.0, 5.0; sub-Alfvénic turbulence), trans-Alfvénic turbulence (B = 1.0), and super-Alfvénic turbulence (B = 0.1). Sub-Alfvénic simulations with B = 3.0 and 5.0 show the lowest isotropy degrees, implying that the isocontour values of the structure function as applied to maps of the velocity centroids show a large degree of anisotropy along the mean magnetic field.
In Figure 3, we show the average degree of isotropy as a function of the sonic Mach number for all the models. The results are obtained by averaging the two cases where the LOS is perpendicular to the mean field, and over the range of separations from r = 10 cells to one-fifth of the computational box. The error bars show the maximum variation of the averaging procedure (including variation across scales). It is clear from Figure 3 that the degree of anisotropy depends primarily on the Alfvénic Mach number and is generally weakly dependent on the sonic Mach number. One can attribute such dependence to the original density field (i.e., arising from shocks in supersonic turbulence). Several additional statistical and observational methods exist to find the sonic Mach number (see Burkhart et al. 2010; Burkhart & Lazarian 2012) and these can be used to break the slight degeneracy observed in Figure 3. However, there are additional observational effects that must be considered, such as the LOS orientation of the magnetic field, the application of cloud boundaries, and the effects of noise, which we will discuss in the next subsection.
4.1. Application to Synthetic Observations
We first repeat the steps EL11 took to make our results more applicable to observations and later add to these steps. First, we (and EL11) include two different contributions to mimic observational effects: an αr−2 gradient to induce the effect of cloud boundaries, and white Gaussian noise. The white noise was produced using fractional Brownian motion (fBm) structures with a power spectrum index (Stutzki et al. 1998; Bensch et al. 2001) set to zero. As in EL11, the noise was added to the density (the mean density is 1.0) with a floor value of 0.01. Thereafter, we apply a Gaussian convolution on the Vz direction for each (x, y) cube positions to mimic the effects of thermal broadening. This is a new addition to the technique and was not performed in EL11. The FWHM of the Gaussian is estimated from the velocity dispersion, σthermal = σturb/Ms, where σturb is the turbulent velocity.
Finally, the "PPV centroid," Cz(x, y), can be computed from the PPV cube as
where Vz, ppv is the velocity axis along the PPV data cube along the z LOS in the cube. In practice, Equation (3) produces output identical to Equation (1); however, Equation (3) is an observational method for calculating the velocity centroid map while Equation (1) can only be applied to numerical simulations.
We also consider the effects of anisotropy on the maps of mean LOS velocity (which is not an observable) defined as
This is to compare the statistics of the synthetic velocity centroids with the actual average velocity of the turbulence.
In addition to creating the velocity centroid by taking the first moment with respect to V, we also investigate the second moment by instead using V2 in Equation (3), thus making the isotropy degree more sensitive to velocity and less sensitive to density. We will refer to this second moment map as "PPV2 centroid," which was not investigated in EL11.
Figure 4 shows the isotropy degree seen in the structure function versus the spatial separation (r) for a simulation with similar initial conditions as those considered by EL11, i.e., model M8 from Table 1. We plot different lines of sight relative to the mean magnetic field as columns across. The top row (panels (a), (b), and (c)) shows the structure functions applied to maps of mean velocities Vx(y, z), Vy(x, z), and Vz(x, y). Panels (d), (e), and (f) show the structure functions computed on PPV cubes which are identical to the EL11 method. Panels (g), (h), and (i) show the structure functions applied to the velocity centroids obtained by taking the first moment of a synthetic PPV cube with thermal smoothing applied (called "PPV centroids"). Finally, in the bottom row (panels (j), (k), and (l)), the structure functions were obtained from the second moment denoted PPV2.
Download figure:
Standard image High-resolution imageThe different lines in each panel denote the density field (ρ) used to obtain the centroids: the solid line corresponds to the original turbulent density field, the dashed line to the density field with an αr−2 gradient applied to it to mimic cloud boundary effects, and the dotted line to the addition of white Gaussian noise and αr−2 cloud boundaries.
The most striking result shown in Figure 4 is that the degree of isotropy is very similar in all rows. In particular, there is no noticeable difference between panels (d)–(f) and (g)–(i), respectively, showing the compatibility between the calculation of the velocity centroid map from the PPV cube and directly from the simulation velocity and density cubes. The isotropy degree of "PPV2 centroids" panels ((j)–(l)) are also comparable with panels (d)–(f) and (g)–(i), although there are notable differences, particularly with the application of a density gradient and white noise (dashed and dotted lines, respectively), due to large- or small-scale fluctuations, respectively. As expected, the X LOS (left column) shows a high degree of isotropy while the Y and Z LOS (middle and right column) show a low degree of isotropy as they are perpendicular to the mean magnetic field.
We also investigate the application of velocity centroid structure function anisotropy technique as outlined above to molecular emissions lines arising from the 13CO J2-1 transition. We apply the post-processing radiative transfer algorithm from Ossenkopf (2002) to our MHD simulations. We refer the reader to Ossenkopf (2002) and Burkhart et al. (2013b, 2013c) for a detailed description of the radiative transfer algorithm. We must scale the simulations to physical cloud parameters and choose a similar initial set up to that of Burkhart et al. (2013b, 2013c): a cloud size of 5 pc, an average density of 275 cm−3, the LOS perpendicular to the mean magnetic field, and a gas temperature of 10 K. The cube is observed at a distance of 450 pc with a beam FWHM of 18'' and a velocity resolution of 0.5 km s−1 and the CO abundance is xco = 1.5 × 10−6. The average optical depth of these simulations is slightly greater than unity.
We show the average isotropy degree versus sonic Mach number of the structure functions of the CO velocity centroids in Figure 5. From our full parameter space represented in Table 1, we choose four different sonic Mach numbers (≈ 8.5, 6.5, 3.5, 0.4) and two different initial Alfvénic number (≈0.7 and 7.0) in order to cover the bulk of the parameter space that was shown in Figure 2. We also overplot the same synthetic velocity centroid isotropy degrees without radiative transfer effects.
Download figure:
Standard image High-resolution imageFigure 5 shows very similar behavior and good agreement with the models of Figure 2 (represented as red diamond and square symbols), which do not include radiative transfer or spatial smoothing. Sub-Alfvénic CO emission creates considerable anisotropy in the isocontours of the CO velocity centroid maps, although it is slightly closer to the isotropic case then the fully optically thin emission. Super-Alfvénic CO emission remains isotropic. This effect is largely insensitive to the sonic Mach number. This gives us confidence that the method could be applied to observational CO emission cubes with success. We will test the effects of varying opacity in future works.
4.2. Anisotropy for Different LOS Orientation Angles
EL11 only considered the LOS either parallel or perpendicular to the mean magnetic field. However, the ISM has a range of LOS orientations relative to the local or global mean field. We repeat the analysis now including rotation of the LOS around an azimuthal origin (0o) and computed the parameter averages. The LOS was rotated (1) around the z-axis (with the azimuthal origin found in the y-axis) and (2) around the y-axis (with the azimuthal origin found in the z-axis). Angle zero thus represents the LOS perpendicular to the mean magnetic field and angle 90 is along the x-axis (LOS parallel to the mean magnetic field). These results are shown in Figures 6 and 7, where in the top row we considered the velocity centroids as those by EL11 ("ideal centroids"), in the middle row we obtained the centroids from PPV cubes, and in the bottom row we used the "PPV2 centroids." In Figure 6, the isotropy degree is plotted versus the LOS orientation angle for different MA values and for a fixed Ms ∼ 0.6, i.e., in the subsonic regime. In Figure 7, the plots are shown for a fixed Ms ∼ 8.0 in a supersonic regime (only considering velocity centroids as the three different centroid calculations yield similar plots).
Download figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution imageThe observed velocity centroid anisotropy is greatest for sight lines perpendicular to the mean magnetic field (i.e., at zero degrees) regardless of sonic Mach number. For super-Alfvénic turbulence, the turbulence remains isotropic regardless of observer angle relative to the mean field. For sub-Alfvénic turbulence, the observer LOS greatly alters the degree of anisotropy observed in the velocity centroid map. In general, the method provides the possibility to distinguish between sub-Alfvénic, super-Alfvénic, and trans-Alfvénic turbulence from 0 to 40 deg relative to the axis perpendicular to the mean field. Past this angle, the turbulence can begin to look isotropic despite the strength of the field. At 90 deg (i.e., parallel to the mean field) the eddies in the centroid maps look isotropic regardless of magnetic field strength. We illustrate this as an illustration in Figure 8
Download figure:
Standard image High-resolution image5. DISCUSSION
The magnetic field in the ISM can be measured from different techniques, as, for example, Zeeman Splitting (e.g., molecular OH), star light polarization, and Faraday rotation (e.g., in H i). However, these techniques often involve complex and difficult data reduction as well as large amounts of telescope time. The method presented here allows one to estimate the plane-of-sky Alfvénic Mach number in the ISM from PPV velocity centroid data, which are provided by several publicly available surveys across multiple tracers (e.g., the COMPLETE survey; see Ridge et al. 2006). In this work, we extended the original analysis of EL11 to include thermal broadening, a varying LOS relative to the mean magnetic field, and investigated the anisotropy in the second moment maps. We also tested our method on synthetic CO maps and found that the effect is preserved even when telescope beam smoothing and radiative transfer effects are included.
The statistical measurement of the Alfvénic Mach number by observing anisotropy in velocity centroid maps is not without its limitations. For example, in this work we demonstrated that we can only detect a lower limit of anisotropy due to projection effects. That is, the mean isotropy degree increases and trend to 1.0 when the LOS changes from a perpendicular direction (angle = 0o) to the mean magnetic field to a parallel direction (angle = 90o). However, dust polarization can provide complimentary plane-of-sky magnetic field directions to the method presented here. The velocity centroid method presented in this paper should be compared with estimates of the Alfvénic Mach number from the Chandrasekhar-Fermi technique when possible. The LOS field cannot be obtained using anisotropy of the velocity centroids and other methods, such as rotation measure or Zeeman splitting, should be employed. In general, we advocate that statistical techniques for studies of turbulence, including studying the Alfvénic Mach number, should be done with multiple techniques in mind. For the estimation of the Alfvénic nature of the ISM this includes techniques such as the bispectrum, PCA and phase coherence (Heyer et al. 2008; Burkhart et al. 2009).
6. CONCLUSIONS
In this work, we explored the anisotropy present in velocity centroid maps caused by an external mean magnetic field. We used a set of ideal MHD simulations to create synthetic PPV maps with a wide range of sonic and Alfvénic Mach numbers and applied the structure function to the velocity centroid maps. We tested the use of our method for observations by including several effects present in the real data such as thermal line broadening, cloud boundaries, noise, and radiative transfer effects from the 13CO J2-1 transition. We found that none of these effects greatly altered the anisotropy observed in the structure function contours. However, the LOS relative to the mean magnetic field does alter the anisotropy observed.
We conclude that anisotropy in velocity centroids is a robust means of distinguishing the Alfvénic regime of a cloud so long as the LOS relative to the global field is known. Since for a given LOS the anisotropy can only be equal to or smaller than that observed when this is perpendicular to the mean magnetic field, the Alfvén Mach number obtained with our technique for any given LOS is an upper limit. In that sense, if one finds that the anisotropy is found to be consistent with sub-Alfvénic turbulence we can be confident that it is; however, the opposite statement is not necessarily true. Our method is complimentary to optical polarization measurements of the orientation (and strength via the Chandrasekhar-Fermi technique) to obtain the plane-of-sky magnetic field and methods to obtain the LOS field, such as Faraday rotation measurements.
Research activities of the Stellar Board of the Federal University of Rio Grande do Norte are supported by continuous grants from CNPq and FAPERN brazilian agencies. J.R.M, A.L., B.B., and I.C.L. acknowledge financial support of the INCT INEspao. I.C.L. acknowledges a CNPq Post-Doctoral fellowship. A.L. is Visiting Researcher PVE/CAPES at the Federal University of Rio Grande do Norte. A.L. and B.B. acknowledge support from the Center for Magnetic Self-Organization in Astrophysical and Laboratory Plasmas and the Physics Graduate Program/UFRN, at Natal, for hospitality. A.E. acknowledges support from CONACYT grant 167611, and UNAM DGAPA grant IG100214.
Footnotes
- 4
The Reynolds number is , which is the ratio of an eddy turnover rate and the viscous dissipation rate . Therefore, large Re correspond to negligible viscous dissipation of large eddies over the cascading time τcasc which is equal to τeddy in Kolmogorov turbulence.