skip to main content

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 8:00 PM ET on Friday, March 21 until 8:00 AM ET on Saturday, March 22 due to maintenance. We apologize for the inconvenience.


Title: (Sub)millimetre dust polarization of protoplanetary discs from scattering by large millimetre-sized irregular grains
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 can 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.

 
more » « less
Award ID(s):
1910106 1910364
PAR ID:
10395179
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
520
Issue:
1
ISSN:
0035-8711
Format(s):
Medium: X Size: p. 1210-1223
Size(s):
p. 1210-1223
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT A number of young circumstellar discs show strikingly ordered (sub)millimetre polarization orientations along the minor axis, which is strong evidence for polarization due to scattering by ∼0.1 mm-sized grains. To test this mechanism further, we model the ALMA dust continuum and polarization data of HD 163296 using radmc-3d. We find that scattering by grains with a maximum size of 90  μm simultaneously reproduces the polarization observed at Band 7 and the unusually low spectral index (α ∼ 1.5) between Bands 7 and 6 in the optically thick inner disc as a result of more efficient scattering at the shorter wavelength. The low spectral index of ∼2.5 inferred for the optically thin gaps is reproduced by the same grains, as a result of telescope beam averaging of the gaps (with an intrinsic α ∼ 4) and their adjacent optically thick rings (where α ≲ 2). The tension between the grain sizes inferred from polarization and spectral index disappears because the low α values do not require large mm-sized grains. In addition, the polarization fraction has a unique azimuthal variation: higher along the major axis than the minor axis in the gaps, but vice versa in the rings. We find a rapidly declining polarization spectrum (with p ∝ λ−3 approximately) in the gaps, which becomes flattened or even inverted towards short wavelengths in the optically thick rings. These contrasting behaviours in the rings and gaps provide further tests for scattering-induced polarization via resolved multiwavelength observations. 
    more » « less
  2. ABSTRACT

    Recent (sub)millimetre polarization observations of protoplanetary discs reveal toroidally aligned, effectively prolate dust grains large enough (at least $\sim 100$$\mu$m) to efficiently scatter millimetre light. The alignment mechanism for these grains remains unclear. We explore the possibility that gas drag aligns grains through gas–dust relative motion when the grain’s centre of mass is offset from its geometric centre, analogous to a badminton birdie’s alignment in flight. A simple grain model of two non-identical spheres illustrates how a grain undergoes damped oscillations from flow-induced restoring torques which align its geometric centre in the flow direction relative to its centre of mass. Assuming specular reflection and subsonic flow, we derive an analytical equation of motion for spheroids where the centre of mass can be shifted away from the spheroid’s geometric centre. We show that a prolate or an oblate grain can be aligned with the long axis parallel to the gas flow when the centre of mass is shifted along that axis. Both scenarios can explain the required effectively prolate grains inferred from observations. Application to a simple disc model shows that the alignment time-scales are shorter than or comparable to the orbital time. The grain alignment direction in a disc depends on the disc (sub-)structure and grain Stokes number (St) with azimuthal alignment for large St grains in sub-Keplerian smooth gas discs and for small St grains near the gas pressure extrema, such as rings and gaps.

     
    more » « less
  3. null (Ed.)
    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. 
    more » « less
  4. null (Ed.)
    ABSTRACT (Sub)millimetre dust opacities are required for converting the observable dust continuum emission to the mass, but their values have long been uncertain, especially in discs around young stellar objects. We propose a method to constrain the opacity κν in edge-on discs from a characteristic optical depth τ0,ν, the density ρ0, and radius R0 at the disc outer edge through κν = τ0,ν/(ρ0R0), where τ0,ν is inferred from the shape of the observed flux along the major axis, ρ0 from gravitational stability considerations, and R0 from direct imaging. We applied the 1D semi-analytical model to the embedded, Class 0, HH 212 disc, which has high-resolution data in Atacama Large Millimetre/submillimetre Array (ALMA) bands 9, 7, 6, and 3 and Very Large Array Ka band (λ = 0.43, 0.85, 1.3, 2.9, and 9.1 mm). The modelling is extended to 2D through RADMC-3D radiative transfer calculations. We find a dust opacity of κν ≈ 1.9 × 10−2, 1.3 × 10−2, and 4.9 × 10−3 cm2 g−1 of gas and dust for ALMA bands 7, 6, and 3, respectively, with uncertainties dependent on the adopted stellar mass. The inferred opacities lend support to the widely used prescription κλ = 2.3 × 10−2(1.3mm/λ) cm2 g−1 . We inferred a temperature of ∼45 K at the disc outer edge that increases radially inwards. It is well above the sublimation temperatures of ices such as CO and N2, which supports the notion that the disc chemistry cannot be completely inherited from the protostellar envelope. 
    more » « less
  5. ABSTRACT

    Polarization is a unique tool to study the dust grains of protoplanetary discs. Polarization around HL Tau was previously imaged using the Atacama Large Millimeter/submillimeter Array (ALMA) at Bands 3 (3.1 mm), 6 (1.3 mm), and 7 (0.87 mm), showing that the polarization orientation changes across wavelength λ. Polarization at Band 7 is predominantly parallel to the disc minor axis but appears azimuthally oriented at Band 3, with the morphology at Band 6 in between the two. We present new ∼0.2 arcsec (29 au) polarization observations at Q-Band (7.0 mm) using the Karl G. Jansky Very Large Array (VLA) and at Bands 4 (2.1 mm), 5 (1.5 mm), and 7 using ALMA, consolidating HL Tau’s position as the protoplanetary disc with the most complete wavelength coverage in dust polarization. The polarization patterns at Bands 4 and 5 follow the previously identified morphological transition with wavelength. From the azimuthal variation, we decompose the polarization into contributions from scattering (s) and thermal emission (t). s decreases slowly with increasing λ, and t increases more rapidly which are expected from optical depth effects of toroidally aligned scattering prolate grains. The weak λ dependence of s is inconsistent with the simplest case of Rayleigh scattering by small grains in the optically thin limit but can be affected by factors such as optical depth, disc substructure, and dust porosity. The sparse polarization detections from the Q-band image are also consistent with toroidally aligned prolate grains.

     
    more » « less