skip to main content

Title: The transition of polarized dust thermal emission from the protostellar envelope to the disc scale
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, more » 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. « less
Authors:
; ; ; ; ; ;
Award ID(s):
1815784 1910364 1716259
Publication Date:
NSF-PAR ID:
10290454
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
507
Issue:
1
Page Range or eLocation-ID:
608 to 620
ISSN:
0035-8711
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    The size of dust grains, a, is key to the physical and chemical processes in circumstellar discs, but observational constraints of grain size remain challenging. (Sub)millimetre continuum observations often show a per cent-level polarization parallel to the disc minor axis, which is generally attributed to scattering by ${\sim}100\, \mu{\rm m}$-sized spherical grains (with a size parameter x ≡ 2$\pi$a/λ < 1, where λ is the wavelength). Larger spherical grains (with x greater than unity) would produce opposite polarization direction. However, the inferred size is in tension with the opacity index β that points to larger mm/cm-sized grains. We investigate the scattering-produced polarization by large irregular grains with a range of x greater than unity with optical properties obtained from laboratory experiments. Using the radiation transfer code, RADMC-3D, we find that large irregular grains still produce polarization parallel to the disc minor axis. If the original forsterite refractive index in the optical is adopted, then all samples can produce the typically observed level of polarization. Accounting for the more commonly adopted refractive index using the DSHARP dust model, only grains with x of several (corresponding to ∼mm-sized grains) can reach the same polarization level. Our results suggest that grains in discs canmore »have sizes in the millimetre regime, which may alleviate the tension between the grain sizes inferred from scattering and other means. Additionally, if large irregular grains are not settled to the mid-plane, their strong forward scattering can produce asymmetries between the near and far side of an inclined disc, which can be used to infer their presence.

    « less
  2. ABSTRACT We present ALMA Band 7 polarization observations of the OMC-1 region of the Orion molecular cloud. We find that the polarization pattern observed in the region is likely to have been significantly altered by the radiation field of the >104 L⊙ high-mass protostar Orion Source I. In the protostar’s optically thick disc, polarization is likely to arise from dust self-scattering. In material to the south of Source I – previously identified as a region of ‘anomalous’ polarization emission – we observe a polarization geometry concentric around Source I. We demonstrate that Source I’s extreme luminosity may be sufficient to make the radiative precession time-scale shorter than the Larmor time-scale for moderately large grains ($\gt 0.005\!-\!0.1\, \mu$m), causing them to precess around the radiation anisotropy vector (k-RATs) rather than the magnetic field direction (B-RATs). This requires relatively unobscured emission from Source I, supporting the hypothesis that emission in this region arises from the cavity wall of the Source I outflow. This is one of the first times that evidence for k-RAT alignment has been found outside of a protostellar disc or AGB star envelope. Alternatively, the grains may remain aligned by B-RATs and trace gas infall on to the Main Ridge.more »Elsewhere, we largely find the magnetic field geometry to be radial around the BN/KL explosion centre, consistent with previous observations. However, in the Main Ridge, the magnetic field geometry appears to remain consistent with the larger-scale magnetic field, perhaps indicative of the ability of the dense Ridge to resist disruption by the BN/KL explosion.« less
  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 aremore »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.

    « less
  4. ABSTRACT

    We analyse the first giant molecular cloud (GMC) simulation to follow the formation of individual stars and their feedback from jets, radiation, winds, and supernovae, using the STARFORGE framework in the GIZMO code. We evolve the GMC for $\sim 9 \rm Myr$, from initial turbulent collapse to dispersal by feedback. Protostellar jets dominate feedback momentum initially, but radiation and winds cause cloud disruption at $\sim 8{{\ \rm per\ cent}}$ star formation efficiency (SFE), and the first supernova at $8.3\, \rm Myr$ comes too late to influence star formation significantly. The per-free-fall SFE is dynamic, accelerating from 0 per cent to $\sim 18{{\ \rm per\ cent}}$ before dropping quickly to <1 per cent, but the estimate from YSO counts compresses it to a narrower range. The primary cluster forms hierarchically and condenses to a brief ($\sim 1\, \mathrm{Myr}$) compact ($\sim 1\, \rm pc$) phase, but does not virialize before the cloud disperses, and the stars end as an unbound expanding association. The initial mass function resembles the Chabrier (2005) form with a high-mass slope α = −2 and a maximum mass of 55 M⊙. Stellar accretion takes $\sim 400\, \rm kyr$ on average, but $\gtrsim 1\,\rm Myr$ for >10 M⊙ stars, so massive stars finishmore »growing latest. The fraction of stars in multiples increase as a function of primary mass, as observed. Overall, the simulation much more closely resembles reality, compared to previous versions that neglected different feedback physics entirely. But more detailed comparison with synthetic observations will be needed to constrain the theoretical uncertainties.

    « less
  5. ABSTRACT

    Satellite galaxies in the cluster environment are more likely to be quenched than galaxies in the general field. Recently, it has been reported that satellite galaxy quenching depends on the orientation relative to their central galaxies: satellites along the major axis of centrals are more likely to be quenched than those along the minor axis. In this paper, we report a detection of such anisotropic quenching up to z ∼ 1 based on a large optically selected cluster catalogue constructed from the Hyper Suprime-Cam Subaru Strategic Program. We calculate the quiescent satellite galaxy fraction as a function of orientation angle measured from the major axis of central galaxies and find that the quiescent fractions at 0.25 < z < 1 are reasonably fitted by sinusoidal functions with amplitudes of a few per cent. Anisotropy is clearer in inner regions (<r200m) of clusters and not significant in cluster outskirts (>r200m). We also confirm that the observed anisotropy cannot be explained by differences in local galaxy density or stellar mass distribution along the two axes. Quiescent fraction excesses between the two axes suggest that the quenching efficiency contributing to the anisotropy is almost independent of stellar mass, at least down to our stellarmore »mass limit of $M_{*}=1\times 10^{10}\, {\rm M}_{\odot }$. Finally, we argue that the physical origins of the observed anisotropy should have shorter quenching time-scales than $\sim 1\, \mathrm{Gyr}$, like ram-pressure stripping, because, for anisotropic quenching to be observed, satellites must be quenched before their initial orientation angles are significantly changed.

    « less