Search for: All records

We use the polaris radiative transfer code to produce simulated circular polarization Zeeman emission maps of the cyanide (CN) J = 1–0 molecular line transition for two types of protostellar envelope magnetohydrodynamic simulations. Our first model is a low-mass disc envelope system (box length L = 200 au), and our second model is the envelope of a massive protostar (L = 104 au) with a protostellar wind and a CN-enhanced outflow shell. We compute the velocity-integrated Stokes I and V, as well as the implied V/I polarization percentage, for each detector pixel location in our simulated emission maps. Our results show that both types of protostellar environments are in principle accessible with current circular polarization instruments, with each containing swaths of envelope area that yield percentage polarizations that exceed the 1.8 per cent nominal sensitivity limit for circular polarization experiments with the Atacama Large Millimeter/submillimeter Array. In both systems, high-polarization (≳1.8 per cent) pixels tend to lie at an intermediate distance away from the central star and where the line-centre opacity of the CN emission is moderately optically thin (τLC ∼ 0.1–1). Furthermore, our computed V/I values scale roughly with the density-weighted mean line-of-sight magnetic field strength, indicating that Zeeman observations can effectively diagnose the strength of envelope-scale magnetic fields. We also find that pixels with large V/I are preferentially co-located where the absolute value of the velocity-integrated V is also large, suggesting that locations with favourable percentage polarization are also favourable in terms of raw signal.

 

The internal velocity structure within dense gaseous cores plays a crucial role in providing the initial conditions for star formation in molecular clouds. However, the kinematic properties of dense gas at core scales (∼0.01−0.1 pc) has not been extensively characterized because of instrument limitations until the unique capabilities of GBT-Argus became available. The ongoing GBT-Argus Large Program, Dynamics in Star-forming Cores (DiSCo) thus aims to investigate the origin and distribution of angular momentum of star-forming cores. DiSCo will survey all starless cores and Class 0 protostellar cores in the Perseus molecular complex down to ∼0.01 pc scales with <0.05 km s−1 velocity resolution using the dense gas tracer N2H+. Here, we present the first data sets from DiSCo towards the B1 and NGC 1333 regions in Perseus. Our results suggest that a dense core’s internal velocity structure has little correlation with other core-scale properties, indicating these gas motions may be originated externally from cloud-scale turbulence. These first data sets also reaffirm the ability of GBT-Argus for studying dense core velocity structure and provided an empirical basis for future studies that address the angular momentum problem with a statistically broad sample.

 

We compare the structure of synthetic dust polarization with synthetic molecular line emission from radiative transfer calculations using a three-dimensional, turbulent collapsing-cloud magnetohydrodynamics simulation. The histogram of relative orientation (HRO) technique and the projected Rayleigh statistic (PRS) are considered. In our trans-Alfvénic (more strongly magnetized) simulation, there is a transition to perpendicular alignment at densities above ∼4 × 103 cm−3. This transition is recovered in most of our synthetic observations of optically thin molecular tracers; however, for 12CO it does not occur and the PRS remains in parallel alignment across the whole observer space. We calculate the physical depth of the optical depth τ = 1 surface and find that for 12CO it is largely located in front of the cloud mid-plane, suggesting that 12CO is too optically thick and instead mainly probes low-volume density gas. In our super-Alfvénic simulation, the magnetic field becomes significantly more tangled, and all observed tracers tend towards no preference for perpendicular or parallel alignment. An observable difference in alignment between optically thin and optically thick tracers may indicate the presence of a dynamically important magnetic field, though there is some degeneracy with viewing angle. We convolve our data with a Gaussian beam and compare it with HRO results of the Vela C molecular cloud. We find good agreement between these results and our sub-Alfvénic simulations when viewed with the magnetic field in the plane of the sky (especially when sensitivity limitations are considered), though the observations are also consistent with an intermediately inclined magnetic field.

 

ABSTRACT Despite the rich observational results on interstellar magnetic fields in star-forming regions, it is still unclear how dynamically significant the magnetic fields are at varying physical scales, because direct measurement of the field strength is observationally difficult. The Davis–Chandrasekhar–Fermi (DCF) method has been the most commonly used method to estimate the magnetic field strength from polarization data. It is based on the assumption that gas turbulent motion is the driving source of field distortion via linear Alfvén waves. In this work, using MHD simulations of star-forming clouds, we test the validity of the assumption underlying the DCF method by examining its accuracy in the real 3D space. Our results suggest that the DCF relation between turbulent kinetic energy and magnetic energy fluctuation should be treated as a statistical result instead of a local property. We then develop and investigate several modifications to the original DCF method using synthetic observations, and propose new recipes to improve the accuracy of DCF-derived magnetic field strength. We further note that the biggest uncertainty in the DCF analysis may come from the linewidth measurement instead of the polarization observation, especially since the line-of-sight gas velocity can be used to estimate the gas volume density, another critical parameter in the DCF method. 

ABSTRACT Polarized dust continuum emission has been observed with Atacama Large Millimeter/submillimeter Array in an increasing number of deeply embedded protostellar systems. It generally shows a sharp transition going from the protostellar envelope to the disc scale, with the polarization fraction typically dropping from ${\sim } 5{{\ \rm per\ cent}}$ to ${\sim } 1{{\ \rm per\ cent}}$ and the inferred magnetic field orientations becoming more aligned with the major axis of the system. We quantitatively investigate these observational trends using a sample of protostars in the Perseus molecular cloud and compare these features with a non-ideal magnetohydrodynamic disc formation simulation. We find that the gas density increases faster than the magnetic field strength in the transition from the envelope to the disc scale, which makes it more difficult to magnetically align the grains on the disc scale. Specifically, to produce the observed ${\sim } 1{{\ \rm per\ cent}}$ polarization at ${\sim } 100\, \mathrm{au}$ scale via grains aligned with the B-field, even relatively small grains of $1\, \mathrm{\mu m}$ in size need to have their magnetic susceptibilities significantly enhanced (by a factor of ∼20) over the standard value, potentially through superparamagnetic inclusions. This requirement is more stringent for larger grains, with the enhancement factor increasing linearly with the grain size, reaching ∼2 × 104 for millimetre-sized grains. Even if the required enhancement can be achieved, the resulting inferred magnetic field orientation in the simulation does not show a preference for the major axis, which is inconsistent with the observed pattern. We thus conclude that the observed trends are best described by the model where the polarization on the envelope scale is dominated by magnetically aligned grains and that on the disc scale by scattering. 

ABSTRACT In typical environments of star-forming clouds, converging supersonic turbulence generates shock-compressed regions, and can create strongly magnetized sheet-like layers. Numerical magnetohydrodynamic simulations show that within these post-shock layers, dense filaments and embedded self-gravitating cores form via gathering material along the magnetic field lines. As a result of the preferred-direction mass collection, a velocity gradient perpendicular to the filament major axis is a common feature seen in simulations. We show that this prediction is in good agreement with recent observations from the CARMA Large Area Star Formation Survey (CLASSy), from which we identified several filaments with prominent velocity gradients perpendicular to their major axes. Highlighting a filament from the north-west part of Serpens South, we provide both qualitative and quantitative comparisons between simulation results and observational data. In particular, we show that the dimensionless ratio Cv ≡ Δvh2/(GM/L), where Δvh is half of the observed perpendicular velocity difference across a filament, and M/L is the filament’s mass per unit length, can distinguish between filaments formed purely due to turbulent compression and those formed due to gravity-induced accretion. We conclude that the perpendicular velocity gradient observed in the Serpens South north-west filament can be caused by gravity-induced anisotropic accretion of material from a flattened layer. Using synthetic observations of our simulated filaments, we also propose that a density-selection effect may explain observed subfilaments (one filament breaking into two components in velocity space) as reported in recent observations. 

ABSTRACT It is well known that the polarized continuum emission from magnetically aligned dust grains is determined to a large extent by local magnetic field structure. However, the observed significant anticorrelation between polarization fraction and column density may be strongly affected, perhaps even dominated by variations in grain alignment efficiency with local conditions, in contrast to standard assumptions of a spatially homogeneous grain alignment efficiency. Here we introduce a generic way to incorporate heterogeneous grain alignment into synthetic polarization observations of molecular clouds (MCs), through a simple model where the grain alignment efficiency depends on the local gas density as a power law. We justify the model using results derived from radiative torque alignment theory. The effects of power-law heterogeneous alignment models on synthetic observations of simulated MCs are presented. We find that the polarization fraction-column density correlation can be brought into agreement with observationally determined values through heterogeneous alignment, though there remains degeneracy with the relative strength of cloud-scale magnetized turbulence and the mean magnetic field orientation relative to the observer. We also find that the dispersion in polarization angles-polarization fraction correlation remains robustly correlated despite the simultaneous changes to both observables in the presence of heterogeneous alignment. 

ABSTRACT Discs are essential to the formation of both stars and planets, but how they form in magnetized molecular cloud cores remains debated. This work focuses on how the disc formation is affected by turbulence and ambipolar diffusion (AD), both separately and in combination, with an emphasis on the protostellar mass accretion phase of star formation. We find that a relatively strong, sonic turbulence on the core scale strongly warps but does not completely disrupt the well-known magnetically induced flattened pseudo-disc that dominates the inner protostellar accretion flow in the laminar case, in agreement with previous work. The turbulence enables the formation of a relatively large disc at early times with or without AD, but such a disc remains strongly magnetized and does not persist to the end of our simulation unless a relatively strong AD is also present. The AD-enabled discs in laminar simulations tend to fragment gravitationally. The disc fragmentation is suppressed by initial turbulence. The AD facilitates the disc formation and survival by reducing the field strength in the circumstellar region through magnetic flux redistribution and by making the field lines there less pinched azimuthally, especially at late times. We conclude that turbulence and AD complement each other in promoting disc formation. The discs formed in our simulations inherit a rather strong magnetic field from its parental core, with a typical plasma-β of order a few tens or smaller, which is 2–3 orders of magnitude lower than the values commonly adopted in magnetohydrodynamic simulations of protoplanetary discs. To resolve this potential tension, longer term simulations of disc formation and evolution with increasingly more realistic physics are needed.