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
; ; ; ; ; ;
Award ID(s):
1815784 1910364 1716259
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
608 to 620
Sponsoring Org:
National Science Foundation
More Like this
  1. 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

    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

    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

    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
  5. ABSTRACT FU Ori is the prototype of FU Orionis systems that are outbursting protoplanetary discs. Magnetic fields in FU Ori’s accretion discs have previously been detected using spectropolarimetry observations for Zeeman effects. We carry out global radiation ideal MHD simulations to study FU Ori’s inner accretion disc. We find that (1) when the disc is threaded by vertical magnetic fields, most accretion occurs in the magnetically dominated atmosphere at z ∼ R, similar to the ‘surface accretion’ mechanism in previous locally isothermal MHD simulations. (2) A moderate disc wind is launched in the vertical field simulations with a terminal speed of ∼300–500 km s−1 and a mass-loss rate of 1–10 per cent the disc accretion rate, which is consistent with observations. Disc wind fails to be launched in simulations with net toroidal magnetic fields. (3) The disc photosphere at the unit optical depth can be either in the wind launching region or the accreting surface region. Magnetic fields have drastically different directions and magnitudes between these two regions. Our fiducial model agrees with previous optical Zeeman observations regarding both the field directions and magnitudes. On the other hand, simulations indicate that future Zeeman observations at near-IR wavelengths or towards other FU Orionis systems maymore »reveal very different magnetic field structures. (4) Due to energy loss by the disc wind, the disc photosphere temperature is lower than that predicted by the thin disc theory, and the previously inferred disc accretion rate may be lower than the real accretion rate by a factor of ∼2–3.« less