Ambipolar diffusion (AD) is believed to be a crucial process for redistributing magnetic flux in the dense molecular gas that occurs in regions of star formation. In this work, we carry out numerical ...simulations of this process in regions of low ionization using the heavy-ion approximation. The simulations are for regions of strong field (plasma β = 0.1) and mildly supersonic turbulence ($\mathcal{M} = 3$, corresponding to an Alfvén Mach number of 0.67). The velocity power spectrum of the neutral gas changes from an Iroshnikov–Kraichnan spectrum in the case of ideal MHD to a Burgers spectrum in the case of a shock-dominated hydrodynamic system. The magnetic power spectrum shows a similar behavior. We use a one-dimensional radiative transfer code to post-process our simulation results; the simulated emission from the CS J = 2–1 and H13CO+ J = 1–0 lines shows that the effects of AD are observable in principle. Linewidths of ions are observed to be less than those of neutrals, and we confirm previous suggestions that this is due to AD. We show that AD is unlikely to affect the Chandrasekhar–Fermi method for inferring field strengths unless the AD is stronger than generally observed. Finally, we present a study of the enhancement of AD by turbulence, finding that AD is accelerated by factor 2–4.5 for non-self-gravitating systems with the level of turbulence we consider.
Ambipolar diffusion redistributes magnetic flux in weakly ionized plasmas and plays a critical role in star formation. Simulations of ambipolar diffusion using explicit MHD codes are prohibitively ...expensive for the level of ionization observed in molecular clouds ( 10 super(-6)), since an enormous number of time steps is required to represent the dynamics of the dominant neutral component with a time step determined by the trace ion component. Here, we show that ambipolar diffusion calculations can be significantly accelerated by the "heavy-ion approximation," in which the mass density of the ions is increased and the collisional coupling constant with the neutrals decreased such that the product remains constant. In this approximation the ambipolar diffusion time and the ambipolar magnetic Reynolds number remain unchanged. We present three tests of the heavy-ion approximation: C-type shocks, the Wardle instability, and the one-dimensional collapse of a magnetized slab. We show that this approximation is quite accurate provided that (1) the square of the Alfven Mach number is small compared to the ambipolar diffusion Reynolds number for dynamical problems, and that (2) the ion mass density is negligible for quasi-static problems; a specific criterion is given for the magnetized slab problem. The first condition can be very stringent for turbulent flows with large density fluctuations.
ABSTRACT
We present the stability analysis of two regions, OMC-3 and OMC-4, in the massive and long molecular cloud complex of Orion A. We obtained 214 $\mu$m HAWC + /SOFIA polarization data, and we ...make use of archival data for the column density and C18O (1–0) emission line. We find clear depolarization in both observed regions and that the polarization fraction is anticorrelated with the column density and the polarization-angle dispersion function. We find that the filamentary cloud and dense clumps in OMC-3 are magnetically supercritical and strongly subvirial. This region should be in the gravitational collapse phase and is consistent with many young stellar objects (YSOs) forming in the region. Our histogram of relative orientation (HRO) analysis shows that the magnetic field is dynamically sub-dominant in the dense gas structures of OMC-3. We present the first polarization map of OMC-4. We find that the observed region is generally magnetically subcritical except for an elongated dense core, which could be a result of projection effect of a filamentary structure aligned close to the line of sight. The relative large velocity dispersion and the unusual positive shape parameters at high column densities in the HROs analysis suggest that our viewing angle may be close to axes of filamentary substructures in OMC-4. The dominating strong magnetic field in OMC-4 is unfavourable for star formation and is consistent with much fewer YSOs than in OMC-3.
ABSTRACT
Optical and infrared polarization mapping and recent Planck observations of the filametary cloud L1495 in Taurus show that the large-scale magnetic field is approximately perpendicular to ...the long axis of the cloud. We use the HAWC + polarimeter on SOFIA to probe the complex magnetic field in the B211 part of the cloud. Our results reveal a dispersion of polarization angles of 36°, about five times that measured on a larger scale by Planck. Applying the Davis–Chandrasekhar–Fermi (DCF) method with velocity information obtained from Institut de Radioastronomie Millimétrique 30 m C18O(1-0) observations, we find two distinct sub-regions with magnetic field strengths differing by more than a factor 3. The quieter sub-region is magnetically critical and sub-Alfv$\acute{\rm e}$nic; the field is comparable to the average field measured in molecular clumps based on Zeeman observations. The more chaotic, super-Alfv$\acute{\rm e}$nic sub-region shows at least three velocity components, indicating interaction among multiple substructures. Its field is much less than the average Zeeman field in molecular clumps, suggesting that the DCF value of the field there may be an underestimate. Numerical simulation of filamentary cloud formation shows that filamentary substructures can strongly perturb the magnetic field. DCF and true field values in the simulation are compared. Pre-stellar cores are observed in B211 and are seen in our simulation. The appendices give a derivation of the standard DCF method that allows for a dispersion in polarization angles that is not small, present an alternate derivation of the structure function version of the DCF method, and treat fragmentation of filaments.
We present calculations of the continuum ultraviolet radiation field (91.2 mn < lambda < 550 mn) penetrating both uniform and clumpy (3D turbulent supersonic magnetohydrodynamic) starless molecular ...gas layers. We find that despite the self-shielding of clumps, pristine (i.e., unreddened) radiation penetrates deeply both the cloud's volume and its mass, resulting in a brighter and bluer intracloud radiation field compared to that in an equivalent uniform cloud. Motivated by these results, we construct and test a toy model ray-tracing scheme for the radiative transfer that fits the UV-visible spectral range with a three-parameter function. We calculate the photoionization rates, Gamma , of the elements C, Na, Mg, Si, S, and Fe as functions of the visual extinction A sub(V) along lines of sight. Typically, the difference in Gamma (A sub(V)) between the clumpy and uniform clouds increases to orders of magnitude at even modest extinctions (A sub(V) similar to 2). Photoionization in the clumpy model extends 2-3 times deeper than in the uniform case, and it dominates cosmic-ray ionization throughout almost the entire volume. We encapsulate these average results in a parameterized form appropriate for when an approximate treatment of the effects of dumpiness is desired. However, the large point-to-point variance in this behavior suggests that uncertainties may arise when using mean values to model particular lines of sight in detail. Ideally, these new results would be used in conjunction with established results for homogeneous clouds in order to span a range of behavior that arises due to cloud inhomogeneities. We briefly explore the importance of the adopted dust properties, characterized by the selective extinction R sub(V) and the scattering parameter g. We find that the UV field is considerably less sensitive to these dust properties in clumpy clouds, emphasizing the preeminence of geometry.
Magnetic field plays a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained owing to the limitations in ...observations. We study the magnetic field in the massive infrared dark cloud G035.39-00.33 from dust continuum polarization observations at 850 m with SCUBA-2/POL-2 at JCMT for the first time. The magnetic field tends to be perpendicular to the densest part of the main filament (FM), whereas it has a less defined relative orientation in the rest of the structure, where it tends to be parallel to some diffuse regions. A mean plane-of-the-sky magnetic field strength of ∼50 G for FM is obtained using the Davis-Chandrasekhar-Fermi method. Based on 13CO (1-0) line observations, we suggest a formation scenario of FM due to large-scale (∼10 pc) cloud-cloud collision. Using additional NH3 line data, we estimate that FM will be gravitationally unstable if it is only supported by thermal pressure and turbulence. The northern part of FM, however, can be stabilized by a modest additional support from the local magnetic field. The middle and southern parts of FM are likely unstable even if the magnetic field support is taken into account. We claim that the clumps in FM may be supported by turbulence and magnetic fields against gravitational collapse. Finally, we identified for the first time a massive (∼200 M ), collapsing starless clump candidate, "c8," in G035.39-00.33. The magnetic field surrounding "c8" is likely pinched, hinting at an accretion flow along the filament.
ABSTRACT
Optical and infrared polarization mapping and recent Planck observations of the filametary cloud L1495 in Taurus show that the large-scale magnetic field is approximately perpendicular to ...the long axis of the cloud. We use the HAWC + polarimeter on SOFIA to probe the complex magnetic field in the B211 part of the cloud. Our results reveal a dispersion of polarization angles of 36°, about five times that measured on a larger scale by Planck. Applying the Davis–Chandrasekhar–Fermi (DCF) method with velocity information obtained from Institut de Radioastronomie Millimétrique 30 m C18O(1-0) observations, we find two distinct sub-regions with magnetic field strengths differing by more than a factor 3. The quieter sub-region is magnetically critical and sub-Alfv$\acute{\rm e}$nic; the field is comparable to the average field measured in molecular clumps based on Zeeman observations. The more chaotic, super-Alfv$\acute{\rm e}$nic sub-region shows at least three velocity components, indicating interaction among multiple substructures. Its field is much less than the average Zeeman field in molecular clumps, suggesting that the DCF value of the field there may be an underestimate. Numerical simulation of filamentary cloud formation shows that filamentary substructures can strongly perturb the magnetic field. DCF and true field values in the simulation are compared. Pre-stellar cores are observed in B211 and are seen in our simulation. The appendices give a derivation of the standard DCF method that allows for a dispersion in polarization angles that is not small, present an alternate derivation of the structure function version of the DCF method, and treat fragmentation of filaments.
ABSTRACT
We have identified 453 compact dense cores in 3 mm continuum emission maps in the ALMA Three-millimetre Observations of Massive Star-forming regions survey, and compiled three catalogues of ...high-mass star-forming cores. One catalogue, referred to as hyper/ultra compact (H/UC)-H ii catalogue, includes 89 cores that enshroud H/UC H ii regions as characterized by associated compact H40α emission. A second catalogue, referred to as pure s-cHMC, includes 32 candidate hot molecular cores (HMCs) showing rich spectra (N ≥ 20 lines) of complex organic molecules (COMs) and not associated with H/UC-H ii regions. The third catalogue, referred to as pure w-cHMC, includes 58 candidate HMCs with relatively low levels of COM richness and not associated with H/UC-H ii regions. These three catalogues of dense cores provide an important foundation for future studies of the early stages of high-mass star formation across the Milky Way. We also find that nearly half of H/UC-H ii cores are candidate HMCs. From the number counts of COM-containing and H/UC-H ii cores, we suggest that the duration of high-mass protostellar cores showing chemically rich features is at least comparable to the lifetime of H/UC-H ii regions. For cores in the H/UC-H ii catalogue, the width of the H40α line increases as the core size decreases, suggesting that the non-thermal dynamical and/or pressure line-broadening mechanisms dominate on the smaller scales of the H/UC-H ii cores.
Star formation in a filamentary infrared dark cloud (IRDC) is simulated over the dynamic range of 4.2 pc to 28 au for a period of 3.5 × 105 yr, including magnetic fields and both radiative and ...outflow feedback from the protostars. At the end of the simulation, the star formation efficiency is 4.3 per cent and the star formation rate per free-fall time is εff ≃ 0.04, within the range of observed values. The total stellar mass increases as ~t2, whereas the number of protostars increases as ~t1.5. We postulate that the density profile around most of the simulated protostars is ~ρ ∝ r-1.5. At the end of the simulation, the protostellar mass function approaches the Chabrier stellar initial mass function. We infer that the time to form a star of median mass 0.2 M⊙ is about 1.4 × 105 yr from the median mass accretion rate. We find good consensus among the protostellar luminosities observed in the large sample of Dunham et al., our simulation and a theoretical estimate, and we conclude that the classical protostellar luminosity problem is resolved. The multiplicity of the stellar systems in the simulation agrees, to within a factor of 2, with observations of Class I young stellar objects; most of the simulated multiple systems are unbound. Bipolar protostellar outflows are launched using a subgrid model, and extend up to 1 pc from their host star. The mass–velocity relation of the simulated outflows is consistent with both observation and theory.
The internal structure of the prestellar core G208.68-19.02-N2 (G208-N2) in the Orion Molecular Cloud 3 (OMC-3) region has been studied with the Atacama Large Millimeter/submillimeter Array (ALMA). ...The dust continuum emission revealed a filamentary structure with a length of $\sim$5000 au and an average H$_2$ volume density of $\sim$6 $\times$ 10$^7$ cm$^{-3}$. At the tip of this filamentary structure, there is a compact object, which we call a ``nucleus", with a radius of $\sim$150--200 au and a mass of $\sim$0.1 M$_{\odot}$. The nucleus has a central density of $\sim$2 $\times$ 10$^9$ cm$^{-3}$ with a radial density profile of $r^{-1.87{\pm}0.11}$. The density scaling of the nucleus is $\sim$3.7 times higher than that of the singular isothermal sphere. This as well as the very low virial parameter of 0.39 suggest that the gravity is dominant over the pressure everywhere in the nucleus. However, there is no sign of CO outflow localized to this nucleus. The filamentary structure is traced by the N$_2$D$^+$ 3--2 emission, but not by the C$^{18}$O 2--1 emission, implying the significant CO depletion due to high density and cold temperature. Toward the nucleus, the N$_2$D$^+$ also shows the signature of depletion. This could imply either the depletion of the parent molecule, N$_2$, or the presence of the embedded very-low luminosity central source that could sublimate the CO in the very small area. The nucleus in G208-N2 is considered to be a prestellar core on the verge of first hydrostatic core (FHSC) formation or a candidate for the FHSC.