Showing 1–50 of 69 results for author: Kataoka, A

Forming localized dust concentrations in a dust ring: DM Tau case study

Authors: Hauyu Baobab Liu, Takayuki Muto, Mihoko Konishi, Chia-Ying Chung, Jun Hashimoto, Kiyoaki Doi, Ruobing Dong, Tomoyuki Kudo, Yasuhiro Hasegawa, Yuka Terada, Akimasa Kataoka

Abstract: The previous, high angular resolution 225 GHz ($\sim$1.3 mm) continuum observations on the transitional disk DM Tau have resolved an outer ring at 20-120 au radii that is weakly azimuthally asymmetric. We aimed to examine dust growth and filtration in the outer ring. We performed the $\sim$0$''$.06 ($\sim$8.7 au) resolution Karl G. Jansky Very Large Array (JVLA) 40-48 GHz ($\sim$7 mm; Q band) cont… ▽ More The previous, high angular resolution 225 GHz ($\sim$1.3 mm) continuum observations on the transitional disk DM Tau have resolved an outer ring at 20-120 au radii that is weakly azimuthally asymmetric. We aimed to examine dust growth and filtration in the outer ring. We performed the $\sim$0$''$.06 ($\sim$8.7 au) resolution Karl G. Jansky Very Large Array (JVLA) 40-48 GHz ($\sim$7 mm; Q band) continuum observations and the complementary observations at lower frequencies. In addition, we analyzed the archival JVLA observations that were taken since 2010. Intriguingly, the Q band image resolved the azimuthally highly asymmetric, knotty dust emission sources close to the inner edge of the outer ring. Fitting the 8-700 GHz spectral energy distribution (SED) with two dust components indicates that the maximum grain size in these knotty dust emission sources is likely $\gtrsim$300 $μ$m while it is $\lesssim$50 $μ$m in the rest of the ring. These results may be explained by trapping of inward migrating grown dust close to the ring inner edge. The exact mechanism for developing the azimuthal asymmetry has not yet been identified, which may be due to planet-disk interaction that might also be responsible for the creation of the dust cavity and pressure bump, or the fluid instabilities and vortex formation due to shear motions. Finally, we remark that the asymmetries in DM Tau are hard to diagnose from the $\gtrsim$225 GHz observations owing to a high optical depth at the ring. In other words, the apparent symmetric or asymmetric morphology of the transitional disks may be related to the optical depths of those disks at the observing frequency. △ Less

Submitted 13 February, 2024; v1 submitted 5 February, 2024; originally announced February 2024.

Comments: 16 pages, 7 figures; accepted to A&A

Filamentary Network and Magnetic Field Structures Revealed with BISTRO in the High-Mass Star-Forming Region NGC2264 : Global Properties and Local Magnetogravitational Configurations

Authors: Jia-Wei Wang, Patrick M. Koch, Seamus D. Clarke, Gary Fuller, Nicolas Peretto, Ya-Wen Tang, Hsi-Wei Yen, Shih-Ping Lai, Nagayoshi Ohashi, Doris Arzoumanian, Doug Johnstone, Ray Furuya, Shu-ichiro Inutsuka, Chang Won Lee, Derek Ward-Thompson, Valentin J. M. Le Gouellec, Hong-Li Liu, Lapo Fanciullo, Jihye Hwang, Kate Pattle, Frédérick Poidevin, Mehrnoosh Tahani, Takashi Onaka, Mark G. Rawlings, Eun Jung Chung , et al. (132 additional authors not shown)

Abstract: We report 850 $μ$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from… ▽ More We report 850 $μ$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from north to east. Field strengths estimates and a virial analysis for the major clumps indicate that NGC 2264C is globally dominated by gravity while in 2264D magnetic, gravitational, and kinetic energies are roughly balanced. We present an analysis scheme that utilizes the locally resolved magnetic field structures, together with the locally measured gravitational vector field and the extracted filamentary network. From this, we infer statistical trends showing that this network consists of two main groups of filaments oriented approximately perpendicular to one another. Additionally, gravity shows one dominating converging direction that is roughly perpendicular to one of the filament orientations, which is suggestive of mass accretion along this direction. Beyond these statistical trends, we identify two types of filaments. The type-I filament is perpendicular to the magnetic field with local gravity transitioning from parallel to perpendicular to the magnetic field from the outside to the filament ridge. The type-II filament is parallel to the magnetic field and local gravity. We interpret these two types of filaments as originating from the competition between radial collapsing, driven by filament self-gravity, and the longitudinal collapsing, driven by the region's global gravity. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: Accepted for publication in the Astrophysical Journal. 43 pages, 32 figures, and 4 tables (including Appendix)

Aligned Grains and Scattered Light Found in Gaps of Planet-Forming Disk

Authors: Ian W. Stephens, Zhe-Yu Daniel Lin, Manuel Fernandez-Lopez, Zhi-Yun Li, Leslie W. Looney, Haifeng Yang, Rachel Harrison, Akimasa Kataoka, Carlos Carrasco-Gonzalez, Satoshi Okuzumi, Ryo Tazaki

Abstract: Polarized (sub)millimeter emission from dust grains in circumstellar disks was initially thought to be due to grains aligned with the magnetic field. However, higher resolution multi-wavelength observations along with improved models found that this polarization is dominated by self-scattering at shorter wavelengths (e.g., 870 $μ$m) and by grains aligned with something other than magnetic fields a… ▽ More Polarized (sub)millimeter emission from dust grains in circumstellar disks was initially thought to be due to grains aligned with the magnetic field. However, higher resolution multi-wavelength observations along with improved models found that this polarization is dominated by self-scattering at shorter wavelengths (e.g., 870 $μ$m) and by grains aligned with something other than magnetic fields at longer wavelengths (e.g., 3 mm). Nevertheless, the polarization signal is expected to depend on the underlying substructure, and observations hitherto have been unable to resolve polarization in multiple rings and gaps. HL Tau, a protoplanetary disk located 147.3 $\pm$ 0.5 pc away, is the brightest Class I or Class II disk at millimeter/submillimeter wavelengths. Here we show deep, high-resolution 870 $μ$m polarization observations of HL Tau, resolving polarization in both the rings and gaps. We find that the gaps have polarization angles with a significant azimuthal component and a higher polarization fraction than the rings. Our models show that the disk polarization is due to both scattering and emission from aligned effectively prolate grains. The intrinsic polarization of aligned dust grains is likely over 10%, which is much higher than what was expected in low resolution observations (~1%). Asymmetries and dust features are revealed in the polarization observations that are not seen in non-polarimetric observations. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: Published in Nature

Panchromatic (Sub)millimeter Polarization Observations of HL Tau Unveil Aligned Scattering Grains

Authors: Zhe-Yu Daniel Lin, Zhi-Yun Li, Ian W. Stephens, Manuel Fernández-López, Carlos Carrasco-González, Claire J. Chandler, Alice Pasetto, Leslie W. Looney, Haifeng Yang, Rachel E. Harrison, Sarah I. Sadavoy, Thomas Henning, A. Meredith Hughes, Akimasa Kataoka, Woojin Kwon, Takayuki Muto, Dominique Segura-Cox

Abstract: Polarization is a unique tool to study the properties of dust grains of protoplanetary disks and detail the initial conditions of planet formation. 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 $λ$. The polarization m… ▽ More Polarization is a unique tool to study the properties of dust grains of protoplanetary disks and detail the initial conditions of planet formation. 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 $λ$. The polarization morphology at Band 7 is predominantly parallel to the disk minor axis but appears azimuthally oriented at Band 3, with the morphology at Band 6 in between the two. We present new ~0.2" (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 disk with the most complete wavelength coverage in dust polarization. The polarization patterns at Bands 4 and 5 continue to follow the morphological transition with wavelength previously identified in Bands 3, 6, and 7. Based on the azimuthal variation, we decompose the polarization into contributions from scattering ($s$) and thermal emission ($t$). We find that $s$ decreases slowly with increasing $λ$, and $t$ increases more rapidly with $λ$ which are expected from optical depth effects of toroidally aligned, scattering prolate grains. The relatively weak $λ$ dependence of $s$ is consistent with large, porous grains. The sparse polarization detections from the Q-band image are also consistent with toroidally aligned prolate grains. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: 22 pages, 13 figures, submitted to MNRAS

Constraints on the dust size distributions in the HD 163296 disk from the difference of the apparent dust ring widths between two ALMA Bands

Authors: Kiyoaki Doi, Akimasa Kataoka

Abstract: The dust size in protoplanetary disks is a crucial parameter for understanding planet formation, while the observational constraints on dust size distribution have large uncertainties. In this study, we present a new method to constrain the dust size distribution from the dust spatial distribution, utilizing the fact that larger dust grains are more spatially localized. We analyze the ALMA Band 6… ▽ More The dust size in protoplanetary disks is a crucial parameter for understanding planet formation, while the observational constraints on dust size distribution have large uncertainties. In this study, we present a new method to constrain the dust size distribution from the dust spatial distribution, utilizing the fact that larger dust grains are more spatially localized. We analyze the ALMA Band 6 (1.25 mm) and Band 4 (2.14 mm) high-resolution images and constrain the dust size distribution in the two rings of the HD 163296 disk. We find that the outer ring at 100 au appears narrower at the longer wavelengths, while the inner ring at 67 au appears to have similar widths across the two wavelengths. We model dust rings trapped at gas pressure maxima, where the dust grains follow a power-law size distribution, and the dust grains of a specific size follow a Gaussian spatial distribution with the width depending on the grain size. By comparing the observations with the models, we constrain the maximum dust size $a_{\mathrm{max}}$ and the exponent of the dust size distribution $p$. We constrain that $0.9 \ \mathrm{mm} < a_{\mathrm{max}} < 5 \ \mathrm{mm}$ and $p < 3.3$ in the inner ring, and $a_{\mathrm{max}} > 3 \times 10^1 \ \mathrm{mm}$ and $3.4 < p < 3.7$ in the outer ring. The larger maximum dust size in the outer ring implies a spatial dependency in dust growth, potentially influencing the formation location of the planetesimals. We further discuss the turbulence strength $α$ derived from the constrained dust spatial distribution, assuming equilibrium between turbulent diffusion and accumulation of dust grains. △ Less

Submitted 31 August, 2023; originally announced August 2023.

Comments: 29 pages, 16 figures. Accepted for publication in ApJ

Dust enrichment and grain growth in a smooth disk around the DG Tau protostar revealed by ALMA triple bands frequency observations

Authors: Satoshi Ohashi, Munetake Momose, Akimasa Kataoka, Aya E Higuchi, Takashi Tsukagoshi, Takahiro Ueda, Claudio Codella, Linda Podio, Tomoyuki Hanawa, Nami Sakai, Hiroshi Kobayashi, Satoshi Okuzumi, Hidekazu Tanaka

Abstract: Characterizing the physical properties of dust grains in a protoplanetary disk is critical to comprehending the planet formation process. Our study presents ALMA high-resolution observations of the young protoplanetary disk around DG Tau at a 1.3 mm dust continuum. The observations, with a spatial resolution of $\approx 0.04''$, or $\approx5$ au, revealed a geometrically thin and smooth disk witho… ▽ More Characterizing the physical properties of dust grains in a protoplanetary disk is critical to comprehending the planet formation process. Our study presents ALMA high-resolution observations of the young protoplanetary disk around DG Tau at a 1.3 mm dust continuum. The observations, with a spatial resolution of $\approx 0.04''$, or $\approx5$ au, revealed a geometrically thin and smooth disk without substantial substructures, suggesting that the disk retains the initial conditions of the planet formation. To further analyze the distributions of dust surface density, temperature, and grain size, we conducted a multi-band analysis with several dust models, incorporating ALMA archival data of the 0.87 mm and 3.1 mm dust polarization. The results showed that the Toomre $Q$ parameter is $\lesssim2$ at a 20 au radius, assuming a dust-to-gas mass ratio of 0.01. This implies that a higher dust-to-gas mass ratio is necessary to stabilize the disk. The grain sizes depend on the dust models, and for the DSHARP compact dust, they were found to be smaller than $\sim400$ $μ$m in the inner region ($r\lesssim20$ au), while exceeding larger than 3 mm in the outer part. Radiative transfer calculations show that the dust scale height is lower than at least one-third of the gas scale height. These distributions of dust enrichment, grain sizes, and weak turbulence strength may have significant implications for the formation of planetesimals through mechanisms such as streaming instability. We also discuss the CO snowline effect and collisional fragmentation in dust coagulation for the origin of the dust size distribution. △ Less

Submitted 26 July, 2023; originally announced July 2023.

Comments: 29 pages, 21 figures, 5 tables. Accepted for publication in ApJ

Formulating Compressive Strength of Dust Aggregates from Low to High Volume Filling Factors with Numerical Simulations

Authors: Misako Tatsuuma, Akimasa Kataoka, Satoshi Okuzumi, Hidekazu Tanaka

Abstract: Compressive strength is a key to understanding the internal structure of dust aggregates in protoplanetary disks and their resultant bodies, such as comets and asteroids in the Solar System. Previous work has modeled the compressive strength of highly-porous dust aggregates with volume filling factors lower than 0.1. However, a comprehensive understanding of the compressive strength from low (… ▽ More Compressive strength is a key to understanding the internal structure of dust aggregates in protoplanetary disks and their resultant bodies, such as comets and asteroids in the Solar System. Previous work has modeled the compressive strength of highly-porous dust aggregates with volume filling factors lower than 0.1. However, a comprehensive understanding of the compressive strength from low ($<0.1$) to high ($>0.1$) volume filling factors is lacking. In this paper, we investigate the compressive strength of dust aggregates by using aggregate compression simulations resolving constituent grains based on JKR theory to formulate the compressive strength comprehensively. We perform a series of numerical simulations with moving periodic boundaries mimicking the compression behavior. As a result, we find that the compressive strength becomes sharply harder when the volume filling factor exceeds 0.1. We succeed in formulating the compressive strength comprehensively by taking into account the rolling motion of aggregates for low volume filling factors and the closest packing of aggregates for high volume filling factors. We also find that the dominant compression mechanisms for high volume filling factors are sliding and twisting motions, while rolling motion dominates for low volume filling factors. We confirm that our results are in good agreement with previous numerical studies. We suggest that our analytical formula is consistent with the previous experimental results if we assume the surface energy of silicate is $\simeq210\pm90\mathrm{\ mJ\ m^{-2}}$. Now, we can apply our results to properties of small compact bodies, such as comets, asteroids, and pebbles. △ Less

Submitted 15 June, 2023; originally announced June 2023.

Comments: 15 pages, 12 figures, accepted for publication in ApJ

Porous Dust Particles in Protoplanetary Disks: Application to the HL Tau Disk

Authors: Shangjia Zhang, Zhaohuan Zhu, Takahiro Ueda, Akimasa Kataoka, Anibal Sierra, Carlos Carrasco-González, Enrique Macías

Abstract: Dust particle sizes constrained from dust continuum and polarization observations by radio interferometry are inconsistent by at least an order of magnitude. Motivated by porous dust observed in small Solar System bodies (e.g., from the Rosetta mission), we explore how the dust particle's porosity affects the estimated particle sizes from these two methods. Porous particles have lower refractive i… ▽ More Dust particle sizes constrained from dust continuum and polarization observations by radio interferometry are inconsistent by at least an order of magnitude. Motivated by porous dust observed in small Solar System bodies (e.g., from the Rosetta mission), we explore how the dust particle's porosity affects the estimated particle sizes from these two methods. Porous particles have lower refractive indices, which affect both opacity and polarization fraction. With weaker Mie interference patterns, the porous particles have lower opacity at mm wavelengths than the compact particles if the particle size exceeds several hundred microns. Consequently, the inferred dust mass using porous particles can be up to a factor of six higher. The most significant difference between compact and porous particles is their scattering properties. The porous particles have a wider range of particle sizes with high linear polarization from dust self-scattering, allowing mm-cm-sized particles to explain polarization observations. With a Bayesian approach, we use porous particles to fit HL Tau disk's multi-wavelength continuum and mm-polarization observations from ALMA and VLA. The moderately porous particles with sizes from 1 mm-1 m can explain both continuum and polarization observations, especially in the region between 20-60 au. If the particles in HL Tau are porous, the porosity should be from 70% to 97% from current polarization observations. We also predict that future observations of the self-scattering linear polarization at longer wavelengths (e.g., ALMA B1 and ngVLA) have the potential to further constrain the particle's porosity and size. △ Less

Submitted 31 May, 2023; originally announced June 2023.

Comments: 39 pages, 18 figures, 1 table, accepted to be published in The Astrophysical Journal

Probing the Temperature Structure of the Inner Region of a Protoplanetary Disk

Authors: Takahiro Ueda, Satoshi Okuzumi, Akimasa Kataoka, Mario Flock

Abstract: Midplane heating induced by disk accretion plays a key role in determining the disk temperature particularly at the inner disk midplane where planets form. However, the efficiency of accretion heating has been not well constrained by observations. We construct two-dimensional models of the Class II disk around CW Tau, taking into account the midplane heating. The models are compared with the ALMA… ▽ More Midplane heating induced by disk accretion plays a key role in determining the disk temperature particularly at the inner disk midplane where planets form. However, the efficiency of accretion heating has been not well constrained by observations. We construct two-dimensional models of the Class II disk around CW Tau, taking into account the midplane heating. The models are compared with the ALMA dust continuum observations at Bands 4, 6, 7 and 8, with an angular resolution of 0.1 arcsec. The observed brightness temperatures are almost wavelength-indenpendent at $\lesssim$10 au. We find that if the maximum dust size $a_{\rm max}$ is $\lesssim100~{\rm μm}$, the brightness temperatures predicted by the model exceed the observed values, regardless of the efficiency of accretion heating. The low observed brightness temperatures can be explained if millimeter scattering reduces the intensity. If the disk is passive, $a_{\rm max}$ needs to be either $\sim150~{\rm μm}$ or $\gtrsim$ few ${\rm cm}$. The accretion heating significantly increases the brightness temperature particularly when $a_{\rm max}\lesssim300~{\rm μm}$, and hence $a_{\rm max}$ needs to be either $\sim300~{\rm μm}$ or $\gtrsim$ few ${\rm cm}$. The midplane temperature is expected to be $\sim$1.5-3 times higher than the observed brightness temperatures, depending on the models. The dust settling effectively increases the temperature of the dust responsible for the millimeter emission in the active disk, which makes the model with $300~{\rm μm}$-sized dust overpredicts the brightness temperatures when strong turbulence is absent. Porous dust (porosity of 0.9) makes the accretion heating more efficient so that some sort of reduction in accretion heating is required. Future longer wavelength and higher angular resolution observations will help us constrain the heating mechanisms of the inner protoplanetary disks. △ Less

Submitted 21 May, 2023; originally announced May 2023.

Comments: 19 pages, 25 figures, accepted for publication in A&A

Journal ref: A&A 675, A176 (2023)

The JCMT BISTRO Survey: Studying the Complex Magnetic Field of L43

Authors: Janik Karoly, Derek Ward-Thompson, Kate Pattle, David Berry, Anthony Whitworth, Jason Kirk, Pierre Bastien, Tao-Chung Ching, Simon Coude, Jihye Hwang, Woojin Kwon, Archana Soam, Jia-Wei Wang, Tetsuo Hasegawa, Shih-Ping Lai, Keping Qiu, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun, Huei-Ru Vivien Chen, Wen Ping Chen, Mike Chen, Zhiwei Chen, Jungyeon Cho, Minho Choi , et al. (133 additional authors not shown)

Abstract: We present observations of polarized dust emission at 850 $μ$m from the L43 molecular cloud which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ($N_{\rm H_2}\sim 10^{22}$-10$^{23}$ cm$^{-2}$) complex molecular cloud with a submillimetre-bright starless core and two protostellar… ▽ More We present observations of polarized dust emission at 850 $μ$m from the L43 molecular cloud which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ($N_{\rm H_2}\sim 10^{22}$-10$^{23}$ cm$^{-2}$) complex molecular cloud with a submillimetre-bright starless core and two protostellar sources. There appears to be an evolutionary gradient along the isolated filament that L43 is embedded within, with the most evolved source closest to the Sco OB2 association. One of the protostars drives a CO outflow that has created a cavity to the southeast. We see a magnetic field that appears to be aligned with the cavity walls of the outflow, suggesting interaction with the outflow. We also find a magnetic field strength of up to $\sim$160$\pm$30 $μ$G in the main starless core and up to $\sim$90$\pm$40 $μ$G in the more diffuse, extended region. These field strengths give magnetically super- and sub-critical values respectively and both are found to be roughly trans-Alfvénic. We also present a new method of data reduction for these denser but fainter objects like starless cores. △ Less

Submitted 22 May, 2023; v1 submitted 18 May, 2023; originally announced May 2023.

Comments: Accepted for publication in ApJ. 23 pages, 9 figures (7 main text, 2 appendix)

Thermal Tomography of the Inner Regions of Protoplanetary Disks with the ngVLA and ALMA

Authors: Satoshi Okuzumi, Munetake Momose, Akimasa Kataoka

Abstract: Understanding the temperature structure of protoplanetary disks is crucial for answering the fundamental question of when and where in the disks rocky planets like our own form. However, the thermal structure of the inner few au of the disks is poorly understood not only because of lack of observational constraints but also because of the uncertainty of accretion heating processes. Here, we propos… ▽ More Understanding the temperature structure of protoplanetary disks is crucial for answering the fundamental question of when and where in the disks rocky planets like our own form. However, the thermal structure of the inner few au of the disks is poorly understood not only because of lack of observational constraints but also because of the uncertainty of accretion heating processes. Here, we propose thermal tomography of the inner regions of protoplanetary disks with the ngVLA and ALMA. The proposed approach is based on the assumption that the inner disk regions are optically thick at submillimeter wavelengths but are marginally optically thin at longer millimeter wavelengths. By combining high-resolution millimeter continuum images from the ngVLA with submillimeter images at comparable resolutions from ALMA, we will be able to reconstruct the radial and vertical structure of the inner few au disk regions. We demonstrate that the thermal tomography we propose can be used to constrain the efficiency of midplane accretion heating, a process that controls the timing of snow-line migration to the rocky planet-forming region, in the few au regions of protoplanetary disks at a distance of 140 pc. △ Less

Submitted 27 April, 2023; originally announced April 2023.

Comments: 4 pages, 3 figures, ngVLA-Japan Memo Series P006 (2021), https://ngvla.nao.ac.jp/researcher/memo/

First BISTRO observations of the dark cloud Taurus L1495A-B10: the role of the magnetic field in the earliest stages of low-mass star formation

Authors: Derek Ward-Thompson, Janik Karoly, Kate Pattle, Anthony Whitworth, Jason Kirk, David Berry, Pierre Bastien, Tao-Chung Ching, Simon Coude, Jihye Hwang, Woojin Kwon, Archana Soam, Jia-Wei Wang, Tetsuo Hasegawa, Shih-Ping Lai, Keping Qiu, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun, Huei-Ru Vivien Chen, Wen Ping Chen, Mike Chen, Zhiwei Chen, Jungyeon Cho, Minho Choi , et al. (133 additional authors not shown)

Abstract: We present BISTRO Survey 850 μm dust emission polarisation observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the JCMT. We observe a roughly triangular network of dense filaments. We detect 9 of the dense starless cores embedded within these filaments in polarisation, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to… ▽ More We present BISTRO Survey 850 μm dust emission polarisation observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the JCMT. We observe a roughly triangular network of dense filaments. We detect 9 of the dense starless cores embedded within these filaments in polarisation, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to the filaments in almost all cases. We also find that the large-scale magnetic field orientation measured by Planck is not correlated with any of the core or filament structures, except in the case of the lowest-density core. We propose a scenario for early prestellar evolution that is both an extension to, and consistent with, previous models, introducing an additional evolutionary transitional stage between field-dominated and matter-dominated evolution, observed here for the first time. In this scenario, the cloud collapses first to a sheet-like structure. Uniquely, we appear to be seeing this sheet almost face-on. The sheet fragments into filaments, which in turn form cores. However, the material must reach a certain critical density before the evolution changes from being field-dominated to being matter-dominated. We measure the sheet surface density and the magnetic field strength at that transition for the first time and show consistency with an analytical prediction that had previously gone untested for over 50 years (Mestel 1965). △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 14 pages, 5 figures. ApJ accepted

JCMT BISTRO Observations: Magnetic Field Morphology of Bubbles Associated with NGC 6334

Authors: Mehrnoosh Tahani, Pierre Bastien, Ray S. Furuya, Kate Pattle, Doug Johnstone, Doris Arzoumanian, Yasuo Doi, Tetsuo Hasegawa, Shu-ichiro Inutsuka, Simon Coudé, Laura Fissel, Michael Chun-Yuan Chen, Frédérick Poidevin, Sarah Sadavoy, Rachel Friesen, Patrick M. Koch, James Di Francesco, Gerald H. Moriarty-Schieven, Zhiwei Chen, Eun Jung Chung, Chakali Eswaraiah, Lapo Fanciullo, Tim Gledhill, Valentin J. M. Le Gouellec, Thiem Hoang , et al. (120 additional authors not shown)

Abstract: We study the HII regions associated with the NGC 6334 molecular cloud observed in the sub-millimeter and taken as part of the B-fields In STar-forming Region Observations (BISTRO) Survey. In particular, we investigate the polarization patterns and magnetic field morphologies associated with these HII regions. Through polarization pattern and pressure calculation analyses, several of these bubbles… ▽ More We study the HII regions associated with the NGC 6334 molecular cloud observed in the sub-millimeter and taken as part of the B-fields In STar-forming Region Observations (BISTRO) Survey. In particular, we investigate the polarization patterns and magnetic field morphologies associated with these HII regions. Through polarization pattern and pressure calculation analyses, several of these bubbles indicate that the gas and magnetic field lines have been pushed away from the bubble, toward an almost tangential (to the bubble) magnetic field morphology. In the densest part of NGC 6334, where the magnetic field morphology is similar to an hourglass, the polarization observations do not exhibit observable impact from HII regions. We detect two nested radial polarization patterns in a bubble to the south of NGC 6334 that correspond to the previously observed bipolar structure in this bubble. Finally, using the results of this study, we present steps (incorporating computer vision; circular Hough Transform) that can be used in future studies to identify bubbles that have physically impacted magnetic field lines. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: Accepted for publication in Astrophysical Journal (ApJ)

The JCMT BISTRO-2 Survey: Magnetic Fields of the Massive DR21 Filament

Authors: Tao-Chung Ching, Keping Qiu, Di Li, Zhiyuan Ren, Shih-Ping Lai, David Berry, Kate Pattle, Ray Furuya, Derek Ward-Thompson, Doug Johnstone, Patrick M. Koch, Chang Won Lee, Thiem Hoang, Tetsuo Hasegawa, Woojin Kwon, Pierre Bastien, Chakali Eswaraiah, Jia-Wei Wang, Kyoung Hee Kim, Jihye Hwang, Archana Soam, A-Ran Lyo, Junhao Liu, Valentin J. M. Le Gouellec, Doris Arzoumanian , et al. (132 additional authors not shown)

Abstract: We present 850 $μ$m dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the sub-filaments, the magnetic fields are mainly parall… ▽ More We present 850 $μ$m dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the sub-filaments, the magnetic fields are mainly parallel to the filamentary structures and smoothly connect to the magnetic fields of the main filament. We compare the POL-2 and Planck dust polarization observations to study the magnetic field structures of the DR21 filament on 0.1--10 pc scales. The magnetic fields revealed in the Planck data are well aligned with those of the POL-2 data, indicating a smooth variation of magnetic fields from large to small scales. The plane-of-sky magnetic field strengths derived from angular dispersion functions of dust polarization are 0.6--1.0 mG in the DR21 filament and $\sim$ 0.1 mG in the surrounding ambient gas. The mass-to-flux ratios are found to be magnetically supercritical in the filament and slightly subcritical to nearly critical in the ambient gas. The alignment between column density structures and magnetic fields changes from random alignment in the low-density ambient gas probed by Planck to mostly perpendicular in the high-density main filament probed by JCMT. The magnetic field structures of the DR21 filament are in agreement with MHD simulations of a strongly magnetized medium, suggesting that magnetic fields play an important role in shaping the DR21 main filament and sub-filaments. △ Less

Submitted 4 December, 2022; originally announced December 2022.

Comments: 26 pages, 13 figures, ApJ accepted

The JCMT BISTRO Survey: A Spiral Magnetic Field in a Hub-filament Structure, Monoceros R2

Authors: Jihye Hwang, Jongsoo Kim, Kate Pattle, Chang Won Lee, Patrick M. Koch, Doug Johnstone, Kohji Tomisaka, Anthony Whitworth, Ray S. Furuya, Ji-hyun Kang, A-Ran Lyo, Eun Jung Chung, Doris Arzoumanian, Geumsook Park, Woojin Kwon, Shinyoung Kim, Motohide Tamura, Jungmi Kwon, Archana Soam, Ilseung Han, Thiem Hoang, Kyoung Hee Kim, Takashi Onaka, Eswaraiah Chakali, Derek Ward-Thompson , et al. (135 additional authors not shown)

Abstract: We present and analyze observations of polarized dust emission at 850 $μ$m towards the central 1 pc $\times$ 1 pc hub-filament structure of Monoceros R2 (Mon R2). The data are obtained with SCUBA-2/POL-2 on the James Clerk Maxwell Telescope (JCMT) as part of the BISTRO (B-fields in Star-forming Region Observations) survey. The orientations of the magnetic field follow the spiral structure of Mon R… ▽ More We present and analyze observations of polarized dust emission at 850 $μ$m towards the central 1 pc $\times$ 1 pc hub-filament structure of Monoceros R2 (Mon R2). The data are obtained with SCUBA-2/POL-2 on the James Clerk Maxwell Telescope (JCMT) as part of the BISTRO (B-fields in Star-forming Region Observations) survey. The orientations of the magnetic field follow the spiral structure of Mon R2, which are well-described by an axisymmetric magnetic field model. We estimate the turbulent component of the magnetic field using the angle difference between our observations and the best-fit model of the underlying large-scale mean magnetic field. This estimate is used to calculate the magnetic field strength using the Davis-Chandrasekhar-Fermi method, for which we also obtain the distribution of volume density and velocity dispersion using a column density map derived from $Herschel$ data and the C$^{18}$O ($J$ = 3-2) data taken with HARP on the JCMT, respectively. We make maps of magnetic field strengths and mass-to-flux ratios, finding that magnetic field strengths vary from 0.02 to 3.64 mG with a mean value of 1.0 $\pm$ 0.06 mG, and the mean critical mass-to-flux ratio is 0.47 $\pm$ 0.02. Additionally, the mean Alfvén Mach number is 0.35 $\pm$ 0.01. This suggests that in Mon R2, magnetic fields provide resistance against large-scale gravitational collapse, and magnetic pressure exceeds turbulent pressure. We also investigate the properties of each filament in Mon R2. Most of the filaments are aligned along the magnetic field direction and are magnetically sub-critical. △ Less

Submitted 13 December, 2022; v1 submitted 12 October, 2022; originally announced October 2022.

Comments: This paper is accepted to the ApJ

Massive compact dust disk with a gap around CW Tau revealed by ALMA multi-band observations

Authors: Takahiro Ueda, Akimasa Kataoka, Takashi Tsukagoshi

Abstract: Compact protoplanetary disks with a radius of $\lesssim$ 50 au are common around young low-mass stars. We report high resolution ALMA dust continuum observations toward a compact disk around CW Tau at Band 4 ($λ=2.2$ mm), 6 (1.3 mm), 7 (0.89 mm) and 8 (0.75 mm). The SED shows the spectral slope of $2.0\pm0.24$ between 0.75 and 1.3 mm, while it is $3.7\pm0.29$ between 2.17 and 3.56 mm. The steep sl… ▽ More Compact protoplanetary disks with a radius of $\lesssim$ 50 au are common around young low-mass stars. We report high resolution ALMA dust continuum observations toward a compact disk around CW Tau at Band 4 ($λ=2.2$ mm), 6 (1.3 mm), 7 (0.89 mm) and 8 (0.75 mm). The SED shows the spectral slope of $2.0\pm0.24$ between 0.75 and 1.3 mm, while it is $3.7\pm0.29$ between 2.17 and 3.56 mm. The steep slope between 2.17 and 3.56 mm is consistent with that of optically thin emission from small grains ($\lesssim$ 350 ${\rm μm}$). We perform parametric fitting of the ALMA data to characterize the dust disk. Interestingly, if the dust-to-gas mass ratio is 0.01, the Toomre's Q parameter reaches $\sim$ 1-3, suggesting that the CW Tau disk might be marginally gravitationally unstable. The total dust mass is estimated as $\sim250M_{\oplus}$ for the maximum dust size of 140 ${\rm μm}$ that is inferred from the previous Band 7 polarimetric observation and at least $80M_{\oplus}$ even for larger grain sizes. This result shows that the CW Tau disk is quite massive in spite of its smallness. Furthermore, we clearly identify a gap structure located at $\sim20$ au, which might be induced by a giant planet. In spite of these interesting characteristics, the CW Tau disk has normal disk luminosity, size and spectral index at ALMA Band 6, which could be a clue to the mass budget problem in Class II disks. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: 14 pages, 12 figures, accepted for publication in ApJ

Setting the Stage for Planet Formation: Measurements and Implications of the Fundamental Disk Properties

Authors: A. Miotello, I. Kamp, T. Birnstiel, L. I. Cleeves, A. Kataoka

Abstract: The field of planet formation is in an exciting era, where recent observations of disks around low- to intermediate-mass stars made with state of the art interferometers and high-contrast optical and IR facilities have revealed a diversity of substructures, some possibly planet-related. It is therefore important to understand the physical and chemical nature of the protoplanetary building blocks,… ▽ More The field of planet formation is in an exciting era, where recent observations of disks around low- to intermediate-mass stars made with state of the art interferometers and high-contrast optical and IR facilities have revealed a diversity of substructures, some possibly planet-related. It is therefore important to understand the physical and chemical nature of the protoplanetary building blocks, as well as their spatial distribution, to better understand planet formation. Since PPVI, the field has seen tremendous improvements in observational capabilities, enabling both surveys of large samples of disks and high resolution imaging studies of a few bright disks. Improvements in data quality and sample size have, however, opened up many fundamental questions about properties such as the mass budget of disks, its spatial distribution, and its radial extent. Moreover, the vertical structure of disks has been studied in greater detail with spatially resolved observations, providing new insights on vertical layering and temperature stratification, yet also bringing rise to questions about other properties, such as material transport and viscosity. Each one of these properties - disk mass, surface density distribution, outer radius, vertical extent, temperature structure, and transport - is of fundamental interest as they collectively set the stage for disk evolution and corresponding planet formation theories. In this chapter, we will review our understanding of the fundamental properties of disks including the relevant observational techniques to probe their nature, modelling methods, and the respective caveats. Finally, we discuss the implications for theories of disk evolution and planet formation underlining what new questions have since arisen as our observational facilities have improved. △ Less

Submitted 18 March, 2022; originally announced March 2022.

Comments: Review Chapter for Protostars and Planets VII, Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura. Accepted version after interaction with the referees and before community feedback. 30 pages (34 with references), 10 figures

B-fields in Star-Forming Region Observations (BISTRO): Magnetic Fields in the Filamentary Structures of Serpens Main

Authors: Woojin Kwon, Kate Pattle, Sarah Sadavoy, Charles L. H. Hull, Doug Johnstone, Derek Ward-Thompson, James Di Francesco, Patrick M. Koch, Ray Furuya, Yasuo Doi, Valentin J. M. Le Gouellec, Jihye Hwang, A-Ran Lyo, Archana Soam, Xindi Tang, Thiem Hoang, Florian Kirchschlager, Chakali Eswaraiah, Lapo Fanciullo, Kyoung Hee Kim, Takashi Onaka, Vera Könyves, Ji-hyun Kang, Chang Won Lee, Motohide Tamura , et al. (127 additional authors not shown)

Abstract: We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filament… ▽ More We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~μ$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities. △ Less

Submitted 13 January, 2022; originally announced January 2022.

Comments: 18 pages, accepted for publication in ApJ

The JCMT BISTRO Survey: An 850/450$μ$m Polarization Study of NGC 2071IR in OrionB

Authors: A-Ran Lyo, Jongsoo Kim, Sarah Sadavoy, Doug Johnstone, David Berry, Kate Pattle, Woojin Kwon, Pierre Bastien, Takashi Onaka, James Di Francesco, Ji-Hyun Kang, Ray Furuya, Charles L. H. Hull, Motohide Tamura, Patrick M. Koch, Derek Ward-Thompson, Tetsuo Hasegawa, Thiem Hoang, Doris Arzoumanian, Chang Won Lee, Chin-Fei Lee, Do-Young Byun, Florian Kirchschlager, Yasuo Doi, Kee-Tae Kim , et al. (121 additional authors not shown)

Abstract: We present the results of simultaneous 450 $μ$m and 850 $μ$m polarization observations toward the massive star forming region NGC 2071IR, a target of the BISTRO (B-fields in Star-Forming Region Observations) Survey, using the POL-2 polarimeter and SCUBA-2 camera mounted on the James Clerk Maxwell Telescope. We find a pinched magnetic field morphology in the central dense core region, which could b… ▽ More We present the results of simultaneous 450 $μ$m and 850 $μ$m polarization observations toward the massive star forming region NGC 2071IR, a target of the BISTRO (B-fields in Star-Forming Region Observations) Survey, using the POL-2 polarimeter and SCUBA-2 camera mounted on the James Clerk Maxwell Telescope. We find a pinched magnetic field morphology in the central dense core region, which could be due to a rotating toroidal disk-like structure and a bipolar outflow originating from the central young stellar object, IRS 3. Using the modified Davis-Chandrasekhar-Fermi method, we obtain a plane-of-sky magnetic field strength of 563$\pm$421 $μ$G in the central $\sim$0.12 pc region from 850 $μ$m polarization data. The corresponding magnetic energy density of 2.04$\times$10$^{-8}$ erg cm$^{-3}$ is comparable to the turbulent and gravitational energy densities in the region. We find that the magnetic field direction is very well aligned with the whole of the IRS 3 bipolar outflow structure. We find that the median value of polarization fractions, 3.0 \%, at 450 $μ$m in the central 3 arcminute region, which is larger than the median value of 1.2 \% at 850 $μ$m. The trend could be due to the better alignment of warmer dust in the strong radiation environment. We also find that polarization fractions decrease with intensity at both wavelengths, with slopes, determined by fitting a Rician noise model, of $0.59 \pm 0.03$ at 450 $μ$m and $0.36 \pm 0.04$ at 850 $μ$m, respectively. We think that the shallow slope at 850 $μ$m is due to grain alignment at the center being assisted by strong radiation from the central young stellar objects. △ Less

Submitted 28 September, 2021; originally announced September 2021.

Comments: 17 pages, 11 figures, 4 tables

Journal ref: The Astrophysical Journal, 2021, 918, 85

Impact of Differential Dust Settling on the SED and Polarization: Application to the Inner Region of the HL Tau Disk

Authors: Takahiro Ueda, Akimasa Kataoka, Shangjia Zhang, Zhaohuan Zhu, Carlos Carrasco-González, Anibal Sierra

Abstract: The polarimetric observations on the protoplanetary disk around HL Tau have shown the scattering-induced polarization at ALMA Band 7, which indicates that the maximum dust size is $\sim 100~{\rm μm}$, while the Spectral Energy Distribution (SED) has suggested that the maximum dust size is $\sim$ mm. To solve the contradiction, we investigate the impact of differential settling of dust grains on th… ▽ More The polarimetric observations on the protoplanetary disk around HL Tau have shown the scattering-induced polarization at ALMA Band 7, which indicates that the maximum dust size is $\sim 100~{\rm μm}$, while the Spectral Energy Distribution (SED) has suggested that the maximum dust size is $\sim$ mm. To solve the contradiction, we investigate the impact of differential settling of dust grains on the SED and polarization. If the disk is optically thick, longer observing wavelength traces more interior layer which would be dominated by larger grains. We find that, the SED of the center part of the HL Tau disk can be explained with mm-sized grains for a broad range of turbulence strength, while $160~{\rm μm}$-sized grains can explain barely only if the turbulence strength parameter $α_{\rm t}$ is lower than $10^{-5}$. We also find that the observed polarization fraction can be potentially explained with the maximum dust size of $1~{\rm mm}$ if $α_{\rm t}\lesssim10^{-5}$, although models with $160~{\rm μm}$-sized grains are also acceptable. However, if the maximum dust size is $\sim3~{\rm mm}$, the simulated polarization fraction is too low to explain the observations even if the turbulence strength is extremely small, indicating the maximum dust size of $\lesssim1$ mm. The degeneracy between 100 ${\rm μm}$-sized and mm-sized grains can be solved by improving the ALMA calibration accuracy or polarimetric observations at (sub-)cm wavelengths. △ Less

Submitted 13 April, 2021; originally announced April 2021.

Comments: 13 pages, 11 figures. Accepted for publication in ApJ

The C$^{18}$O core mass function toward Orion A: Single-dish observations

Authors: Hideaki Takemura, Fumitaka Nakamura, Shun Ishii, Yoshito Shimajiri, Patricio Sanhueza, Takashi Tsukagoshi, Ryohei Kawabe, Tomoya Hirota, Akimasa Kataoka

Abstract: We have performed an unbiased dense core survey toward the Orion A Giant Molecular Cloud in the C$^{18}$O ($J$=1--0) emission line taken with the Nobeyama Radio Observatory (NRO) 45-m telescope. The effective angular resolution of the map is 26", which corresponds to $\sim$ 0.05 pc at a distance of 414 pc. By using the Herschel-Planck H$_2$ column density map, we calculate the C$^{18}$O fractional… ▽ More We have performed an unbiased dense core survey toward the Orion A Giant Molecular Cloud in the C$^{18}$O ($J$=1--0) emission line taken with the Nobeyama Radio Observatory (NRO) 45-m telescope. The effective angular resolution of the map is 26", which corresponds to $\sim$ 0.05 pc at a distance of 414 pc. By using the Herschel-Planck H$_2$ column density map, we calculate the C$^{18}$O fractional abundance and find that it is roughly constant over the column density range of $\lesssim$ 5 $\times$ 10$^{22}$ cm$^{-3}$, although a trend of C$^{18}$O depletion is determined toward higher column density. Therefore, C$^{18}$O intensity can follow the cloud structure reasonably well. The mean C$^{18}$O abundance in Orion A is estimated to be 5.7$\times$10$^{-7}$, which is about 3 times larger than the fiducial value. We identified 746 C$^{18}$O cores with astrodendro and classified 709 cores as starless cores. We compute the core masses by decomposing the Herschel-Planck dust column density using the relative proportions of the C$^{18}$O integrated intensities of line-of-sight components. Applying this procedure, we attempt to remove the contribution of the background emission, i.e., the ambient gas outside the cores. Then, we derived mass function for starless cores and found that it resembles the stellar initial mass function (IMF). The CMF for starless cores, $dN/dM$, is fitted with a power-law relation of $M^α$ with a power index of $α= -$2.25$\pm$ 0.16 at the high-mass slope ($\gtrsim$ 0.44 $M_\odot$). We also found that the ratio of each core mass to the total mass integrated along the line of sight is significantly large. Therefore, in the previous studies, the core masses derived from the dust image are likely to be overestimated at least by a factor of a few. Accordingly, such previous studies may underestimate the star formation efficiency of individual cores. △ Less

Submitted 25 February, 2021; originally announced March 2021.

Comments: 66 pages, 15 figures, accepted by PASJ

Revealing the diverse magnetic field morphologies in Taurus dense cores with sensitive sub-millimeter polarimetry

Authors: Chakali Eswaraiah, Di Li, Ray S. Furuya, Tetsuo Hasegawa, Derek Ward-Thompson, Keping Qiu, Nagayoshi Ohashi, Kate Pattle, Sarah Sadavoy, Charles L. H. Hull, David Berry, Yasuo Doi, Tao-Chung Ching, Shih-Ping Lai, Jia-Wei Wang, Patrick M. Koch, Jungmi Kwon, Woojin Kwon, Pierre Bastien, Doris Arzoumanian, Simon Coudé, Archana Soam, Lapo Fanciullo, Hsi-Wei Yen, Junhao Liu , et al. (120 additional authors not shown)

Abstract: We have obtained sensitive dust continuum polarization observations at 850 $μ$m in the B213 region of Taurus using POL-2 on SCUBA-2 at the James Clerk Maxwell Telescope (JCMT), as part of the BISTRO (B-fields in STar-forming Region Observations) survey. These observations allow us to probe magnetic field (B-field) at high spatial resolution ($\sim$2000 au or $\sim$0.01 pc at 140 pc) in two protost… ▽ More We have obtained sensitive dust continuum polarization observations at 850 $μ$m in the B213 region of Taurus using POL-2 on SCUBA-2 at the James Clerk Maxwell Telescope (JCMT), as part of the BISTRO (B-fields in STar-forming Region Observations) survey. These observations allow us to probe magnetic field (B-field) at high spatial resolution ($\sim$2000 au or $\sim$0.01 pc at 140 pc) in two protostellar cores (K04166 and K04169) and one prestellar core (Miz-8b) that lie within the B213 filament. Using the Davis-Chandrasekhar-Fermi method, we estimate the B-field strengths in K04166, K04169, and Miz-8b to be 38$\pm$14 $μ$G, 44$\pm$16 $μ$G, and 12$\pm$5 $μ$G, respectively. These cores show distinct mean B-field orientations. B-field in K04166 is well ordered and aligned parallel to the orientations of the core minor axis, outflows, core rotation axis, and large-scale uniform B-field, in accordance with magnetically regulated star formation via ambipolar diffusion taking place in K04166. B-field in K04169 is found to be ordered but oriented nearly perpendicular to the core minor axis and large-scale B-field, and not well-correlated with other axes. In contrast, Miz-8b exhibits disordered B-field which show no preferred alignment with the core minor axis or large-scale field. We found that only one core, K04166, retains a memory of the large-scale uniform B-field. The other two cores, K04169 and Miz-8b, are decoupled from the large-scale field. Such a complex B-field configuration could be caused by gas inflow onto the filament, even in the presence of a substantial magnetic flux. △ Less

Submitted 3 March, 2021; originally announced March 2021.

Comments: 20 pages (10 are main), 5 figures (3 are main), and 3 tables (2 are main); Accepted for its publication in ApJL

Estimate on dust scale height from ALMA dust continuum image of the HD 163296 protoplanetary disk

Authors: Kiyoaki Doi, Akimasa Kataoka

Abstract: We aim at estimating the dust scale height of protoplanetary disks from millimeter continuum observations. First, we present a general expression of intensity of a ring in a protoplanetary disk, and show that we can constrain the dust scale height by the azimuthal intensity variation. Then, we apply the presented methodology to the two distinct rings at 68 au and at 100 au of the protoplanetary di… ▽ More We aim at estimating the dust scale height of protoplanetary disks from millimeter continuum observations. First, we present a general expression of intensity of a ring in a protoplanetary disk, and show that we can constrain the dust scale height by the azimuthal intensity variation. Then, we apply the presented methodology to the two distinct rings at 68 au and at 100 au of the protoplanetary disk around HD 163296. We constrain the dust scale height by comparing the DSHARP high-resolution millimeter dust continuum image with radiative transfer simulations using RADMC-3D. We find that h_d/h_g > 0.84 at the inner ring and h_d/h_g < 0.11 at the outer ring with the 3 sigma uncertainties, where h_d is the dust scale height and h_g is the gas scale height. This indicates that the dust is flared at the inner ring and settled at the outer ring. We further constrain the ratio of turbulence parameter alpha to gas-to-dust-coupling parameter St from the derived dust scale height; alpha/St > 2.4 at the inner ring, and alpha/St < 1.1*10^{-2} at the outer ring. This result shows that the turbulence is stronger or the dust is smaller at the inner ring than at the outer ring. △ Less

Submitted 5 March, 2021; v1 submitted 11 February, 2021; originally announced February 2021.

Comments: 19 pages, 16 figures, accepted to ApJ, error estimation revised

Journal ref: ApJ 912, 164 (2021)

ALMA observation of the protoplanetary disk around WW Cha: faint double-peaked ring and asymmetric structure

Authors: Kazuhiro D. Kanagawa, Jun Hashimoto, Takayuki Muto, Takashi Tsukagoshi, Sanemichi Z. Takahashi, Yasuhiro Hasegawa, Mihoko Konishi, Hideko Nomura, Hauyu Baobab Liu, Ruobing Dong, Akimasa Kataoka, Munetake Momose, Tomohiro Ono, Michael Sitko, Michihiro Takami, Kengo Tomida

Abstract: We present Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 observations of dust continuum emission of the disk around WW Cha. The dust continuum image shows a smooth disk structure with a faint (low-contrast) dust ring, extending from $\sim 40$ au to $\sim 70$ au, not accompanied by any gap. We constructed the simple model to fit the visibility of the observed data by using MCMC method… ▽ More We present Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 observations of dust continuum emission of the disk around WW Cha. The dust continuum image shows a smooth disk structure with a faint (low-contrast) dust ring, extending from $\sim 40$ au to $\sim 70$ au, not accompanied by any gap. We constructed the simple model to fit the visibility of the observed data by using MCMC method and found that the bump (we call the ring without the gap the bump) has two peaks at $40$ au and $70$ au. The residual map between the model and observation indicates asymmetric structures at the center and the outer region of the disk. These asymmetric structures are also confirmed by model-independent analysis of the imaginary part of the visibility. The asymmetric structure at the outer region is consistent with a spiral observed by SPHERE. To constrain physical quantities of the disk (dust density and temperature), we carried out radiative transfer simulations. We found that the midplane temperature around the outer peak is close to the freezeout temperature of CO on water ice ($\sim 30$ K). The temperature around the inner peak is about $50$ K, which is close to the freezeout temperature of H$_2$S and also close to the sintering temperature of several species. We also discuss the size distribution of the dust grains using the spectral index map obtained within the Band 6 data. △ Less

Submitted 17 March, 2021; v1 submitted 25 January, 2021; originally announced January 2021.

Comments: 24 pages, 25 figures, accepted for publication in The Astrophysical Journal

Rotational Disruption of Porous Dust Aggregates due to Gas Flow in Protoplanetary Disks

Authors: Misako Tatsuuma, Akimasa Kataoka

Abstract: We introduce a possible disruption mechanism of dust grains in planet formation by their spinning motion. This mechanism has been discussed as rotational disruption for the interstellar dust grains. We theoretically calculate whether porous dust aggregates can be disrupted by their spinning motion and if it prohibits dust growth in protoplanetary disks. We assume radiative torque and gas-flow torq… ▽ More We introduce a possible disruption mechanism of dust grains in planet formation by their spinning motion. This mechanism has been discussed as rotational disruption for the interstellar dust grains. We theoretically calculate whether porous dust aggregates can be disrupted by their spinning motion and if it prohibits dust growth in protoplanetary disks. We assume radiative torque and gas-flow torque as driving sources of the spinning motion, assume that dust aggregates reach a steady-state rigid rotation, and compare the obtained tensile stress due to the centrifugal force with their tensile strength. We model the irregularly-shaped dust aggregates by introducing a parameter, $γ_\mathrm{ft}$, that mimics the conversion efficiency from force to torque. As a result, we find that porous dust aggregates are rotationally disrupted by their spinning motion induced by gas flow when their mass is larger than $\sim10^8$ g and their volume filling factor is smaller than $\sim 0.01$ in our fiducial model, while relatively compact dust aggregates with volume filling factor more than 0.01 do not face this problem. If we assume the dust porosity evolution, we find that dust aggregates whose Stokes number is $\sim0.1$ can be rotationally disrupted in their growth and compression process. Our results suggest that the growth of dust aggregates may be halted due to rotational disruption or that other compression mechanisms are needed to avoid it. We also note that dust aggregates are not rotationally disrupted when $γ_\mathrm{ft}\leq0.02$ in our fiducial model and the modeling of irregularly-shaped dust aggregates is essential in future work. △ Less

Submitted 7 April, 2021; v1 submitted 13 January, 2021; originally announced January 2021.

Comments: 16 pages, 16 figures, accepted for publication in ApJ

Dust polarized emission observations of NGC 6334; BISTRO reveals the details of the complex but organized magnetic field structure of the high-mass star-forming hub-filament network

Authors: D. Arzoumanian, R. Furuya, T. Hasegawa, M. Tahani, S. Sadavoy, C. L. H. Hull, D. Johnstone, P. M. Koch, S. -i. Inutsuka, Y. Doi, T. Hoang, T. Onaka, K. Iwasaki, Y. Shimajiri, T. Inoue, N. Peretto, P. André, P. Bastien, D. Berry, H. -R. V. Chen, J. Di Francesco, C. Eswaraiah, L. Fanciullo, L. M. Fissel, J. Hwang , et al. (123 additional authors not shown)

Abstract: [Abridged] Filaments and hubs have received special attention recently thanks to studies showing their role in star formation. While the column density and velocity structures of both filaments and hubs have been studied, their magnetic fields (B-field) are not yet characterized. We aim to understand the role of the B-field in the dynamical evolution of the NGC 6334 hub-filament network. We presen… ▽ More [Abridged] Filaments and hubs have received special attention recently thanks to studies showing their role in star formation. While the column density and velocity structures of both filaments and hubs have been studied, their magnetic fields (B-field) are not yet characterized. We aim to understand the role of the B-field in the dynamical evolution of the NGC 6334 hub-filament network. We present new observations of the dust polarized emission at 850$μ$m towards NGC 6334 obtained with the JCMT/POL-2. We study the distribution and dispersion of the polarized intensity ($PI$), the polarization fraction ($PF$), and the B-field angle ($θ_{B}$). We derive the power spectrum of the intensity and $θ_{B}$ along the ridge crest. Our analyses show a complex B-field structure when observed over the whole region ($\sim10$ pc), however, at smaller scales ($\sim1$ pc), $θ_{B}$ varies coherently along the filaments. The observed power spectrum of $θ_{B}$ can be well represented with a power law function with a slope $-1.33\pm0.23$, which is $\sim20\%$ shallower than that of $I$. This result is compatible with the properties of simulated filaments and may indicate the processes at play in the formation of filaments. $θ_{B}$ rotates from being mostly perpendicular to the filament crests to mostly parallel as they merge with the hubs. This variation of $θ_{B}$ may be tracing local velocity flows of matter in-falling onto the hubs. Our analysis suggests a variation of the energy balance along the crests of these filaments, from magnetically critical/supercritical at their far ends to magnetically subcritical near the hubs. We detect an increase of $PF$ towards the high-column density star cluster-forming hubs that may result from the increase of grain alignment efficiency due to stellar radiation from the newborn stars. △ Less

Submitted 23 December, 2020; originally announced December 2020.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 647, A78 (2021)

Observations of magnetic fields surrounding LkH$α$ 101 taken by the BISTRO survey with JCMT-POL-2

Authors: Nguyen Bich Ngoc, Pham Ngoc Diep, Harriet Parsons, Kate Pattle, Thiem Hoang, Derek Ward-Thompson, Le Ngoc Tram, Charles L. H. Hull, Mehrnoosh Tahani, Ray Furuya, Pierre Bastien, Keping Qiu, Tetsuo Hasegawa, Woojin Kwon, Yasuo Doi, Shih-Ping Lai, Simon Coude, David Berry, Tao-Chung Ching, Jihye Hwang, Archana Soam, Jia-Wei Wang, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun , et al. (124 additional authors not shown)

Abstract: We report the first high spatial resolution measurement of magnetic fields surrounding LkH$α$ 101, a part of the Auriga-California molecular cloud. The observations were taken with the POL-2 polarimeter on the James Clerk Maxwell Telescope within the framework of the B-fields In Star-forming Region Observations (BISTRO) survey. Observed polarization of thermal dust emission at 850 $μ$m is found to… ▽ More We report the first high spatial resolution measurement of magnetic fields surrounding LkH$α$ 101, a part of the Auriga-California molecular cloud. The observations were taken with the POL-2 polarimeter on the James Clerk Maxwell Telescope within the framework of the B-fields In Star-forming Region Observations (BISTRO) survey. Observed polarization of thermal dust emission at 850 $μ$m is found to be mostly associated with the red-shifted gas component of the cloud. The magnetic field displays a relatively complex morphology. Two variants of the Davis-Chandrasekhar-Fermi method, unsharp masking and structure function, are used to calculate the strength of magnetic fields in the plane of the sky, yielding a similar result of $B_{\rm POS}\sim 115$ $\mathrmμ$G. The mass-to-magnetic-flux ratio in critical value units, $λ\sim0.3$, is the smallest among the values obtained for other regions surveyed by POL-2. This implies that the LkH$α$ 101 region is sub-critical and the magnetic field is strong enough to prevent gravitational collapse. The inferred $δB/B_0\sim 0.3$ implies that the large scale component of the magnetic field dominates the turbulent one. The variation of the polarization fraction with total emission intensity can be fitted by a power-law with an index of $α=0.82\pm0.03$, which lies in the range previously reported for molecular clouds. We find that the polarization fraction decreases rapidly with proximity to the only early B star (LkH$α$ 101) in the region. The magnetic field tangling and the joint effect of grain alignment and rotational disruption by radiative torques are potential of explaining such a decreasing trend. △ Less

Submitted 8 December, 2020; originally announced December 2020.

Comments: 25 pages, 19 figures, Accepted for publication in The Astrophysical Journal

Modeling of the ALMA HL Tau Polarization by Mixture of Grain Alignment and Self-scattering

Authors: Tomohiro Mori, Akimasa Kataoka

Abstract: Dust polarization at (sub)millimeter wavelengths has been observed for many protoplanetary disks. Theoretically, multiple origins potentially contribute to the polarized emission but it is still uncertain what mechanism is dominant in disk millimeter polarization. To quantitatively address the origin, we perform radiative transfer calculations of the mixture of alignment and self-scattering induce… ▽ More Dust polarization at (sub)millimeter wavelengths has been observed for many protoplanetary disks. Theoretically, multiple origins potentially contribute to the polarized emission but it is still uncertain what mechanism is dominant in disk millimeter polarization. To quantitatively address the origin, we perform radiative transfer calculations of the mixture of alignment and self-scattering induced polarization to reproduce the 3.1 mm polarization of the HL Tau disk, which shows azimuthal pattern in polarization vectors. We find that a mixture of the grain alignment and self-scattering is essential to reproduce the HL Tau 3.1 mm polarization properties. Our model shows that the polarization of the HL Tau at 3.1 mm can be decomposed to be the combination of the self-scattering parallel to the minor axis and the alignment-induced polarization parallel to the major axis, with the orders of $\sim 0.5\%$ fraction for each component. This slightly eases the tight constraints on the grain size of $\sim 70~{\rm~μm}$ to be $\sim 130 {\rm~μm}$ in the previous studies but further modeling is needed. In addition, the grain alignment model requires effectively prolate grains but the physics to reproduce it in protoplanetary disks is still a mystery. △ Less

Submitted 3 December, 2020; originally announced December 2020.

Comments: 12 pages, 5 figures, accepted for publication in ApJ

Solving grain size inconsistency between ALMA polarization and VLA continuum in the Ophiuchus IRS 48 protoplanetary disk

Authors: Satoshi Ohashi, Akimasa Kataoka, Nienke Van der Marel, Charles L. H. Hull, William R. F. Dent, Adriana Pohl, Paola Pinilla, Ewine F. van Dishoeck, Thomas Henning

Abstract: The protoplanetary disk around Ophiuchus IRS 48 shows an azimuthally asymmetric dust distribution in (sub-)millimeter observations, which is interpreted as a vortex, where millimeter/centimeter-sized particles are trapped at the location of the continuum peak. In this paper, we present 860 $μ$m ALMA observations of polarized dust emission of this disk. The polarized emission was detected toward a… ▽ More The protoplanetary disk around Ophiuchus IRS 48 shows an azimuthally asymmetric dust distribution in (sub-)millimeter observations, which is interpreted as a vortex, where millimeter/centimeter-sized particles are trapped at the location of the continuum peak. In this paper, we present 860 $μ$m ALMA observations of polarized dust emission of this disk. The polarized emission was detected toward a part of the disk. The polarization vectors are parallel to the disk minor axis, and the polarization fraction was derived to be $1-2$\%. These characteristics are consistent with models of self-scattering of submillimeter-wave emission, which indicate a maximum grain size of $\sim100$ $μ$m. However, this is inconsistent with the previous interpretation of millimeter/centimeter dust particles being trapped by a vortex. To explain both, ALMA polarization and previous ALMA and VLA observations, we suggest that the thermal emission at 860 $μ$m wavelength is optically thick ($τ_{\rm abs}\sim7.3$) at the dust trap with the maximum observable grain size of $\sim100$ $μ$m rather than an optically thin case with $\sim$ cm dust grains. We note that we cannot rule out that larger dust grains are accumulated near the midplane if the 860 $μ$m thermal emission is optically thick. △ Less

Submitted 29 July, 2020; originally announced July 2020.

Comments: 22 pages, 17 figures, accepted for publication in ApJ

The JCMT BISTRO Survey: Magnetic Fields Associated with a Network of Filaments in NGC 1333

Authors: Yasuo Doi, Tetsuo Hasegawa, Ray S. Furuya, Simon Coudé, Charles L. H. Hull, Doris Arzoumanian, Pierre Bastien, Michael Chun-Yuan Chen, James di Francesco, Rachel Friesen, Martin Houde, Shu-ichiro Inutsuka, Steve Mairs, Masafumi Matsumura, Takashi Onaka, Sarah Sadavoy, Yoshito Shimajiri, Mehrnoosh Tahani, Kohji Tomisaka, Chakali Eswaraiah, Patrick M. Koch, Kate Pattle, Chang Won Lee, Motohide Tamura, David Berry , et al. (113 additional authors not shown)

Abstract: We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary… ▽ More We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multiscale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~1 pc and remains continuous from the scales of filaments (~0.1 pc) to that of protostellar envelopes (~0.005 pc or ~1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data toward the filament network. △ Less

Submitted 21 July, 2020; v1 submitted 30 June, 2020; originally announced July 2020.

Comments: 33 pages, 21 figures, accepted for publication in ApJ

Super-resolution Imaging of the Protoplanetary Disk HD 142527 Using Sparse Modeling

Authors: Masayuki Yamaguchi, Kazunori Akiyama, Takashi Tsukagoshi, Takayuki Muto, Akimasa Kataoka, Fumie Tazaki, Shiro Ikeda, Misato Fukagawa, Mareki Honma, Ryohei Kawabe

Abstract: With an emphasis on improving the fidelity even in super-resolution regimes, new imaging techniques have been intensively developed over the last several years, which may provide substantial improvements to the interferometric observation of protoplanetary disks. In this study, sparse modeling (SpM) is applied for the first time to observational data sets taken by the Atacama Large Millimeter/subm… ▽ More With an emphasis on improving the fidelity even in super-resolution regimes, new imaging techniques have been intensively developed over the last several years, which may provide substantial improvements to the interferometric observation of protoplanetary disks. In this study, sparse modeling (SpM) is applied for the first time to observational data sets taken by the Atacama Large Millimeter/submillimeter Array (ALMA). The two data sets used in this study were taken independently using different array configurations at Band 7 (330 GHz), targeting the protoplanetary disk around HD 142527; one in the shorter-baseline array configuration (~ 430 m), and the other in the longer-baseline array configuration (~ 1570 m). The image resolutions reconstructed from the two data sets are different by a factor of ~ 3. We confirm that the previously known disk structures appear on the images produced by both SpM and CLEAN at the standard beam size. The image reconstructed from the shorter-baseline data using the SpM matches that obtained with the longer-baseline data using CLEAN, achieving a super-resolution image from which a structure finer than the beam size can be reproduced. Our results demonstrate that on-going intensive development in the SpM imaging technique is beneficial to imaging with ALMA. △ Less

Submitted 23 April, 2020; originally announced April 2020.

GW Ori: Interactions Between a Triple-star System and its Circumtriple Disk in Action

Authors: Jiaqing Bi, Nienke van der Marel, Ruobing Dong, Takayuki Muto, Rebecca G. Martin, Jeremy L. Smallwood, Jun Hashimoto, Hauyu Baobab Liu, Hideko Nomura, Yasuhiro Hasegawa, Michihiro Takami, Mihoko Konishi, Munetake Momose, Kazuhiro D. Kanagawa, Akimasa Kataoka, Tomohiro Ono, Michael L. Sitko, Sanemichi Z. Takahashi, Kengo Tomida, Takashi Tsukagoshi

Abstract: GW Ori is a hierarchical triple system which has a rare circumtriple disk. We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and 12CO J=2-1 molecular gas emission of the disk. For the first time, we identify three dust rings in the disk at ~46, 188, and 338 AU, with estimated dust mass of ~70-250 Earth masses, respectively. To our knowledge, the o… ▽ More GW Ori is a hierarchical triple system which has a rare circumtriple disk. We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and 12CO J=2-1 molecular gas emission of the disk. For the first time, we identify three dust rings in the disk at ~46, 188, and 338 AU, with estimated dust mass of ~70-250 Earth masses, respectively. To our knowledge, the outer ring in GW Ori is the largest dust ring ever found in protoplanetary disks. We use visibility modelling of dust continuum to show that the disk has misaligned parts and the innermost dust ring is eccentric. The disk misalignment is also suggested by the CO kinematics modelling. We interpret these substructures as evidence of ongoing dynamical interactions between the triple stars and the circumtriple disk. △ Less

Submitted 29 April, 2020; v1 submitted 7 April, 2020; originally announced April 2020.

Comments: 17 pages, 3+3 figures, 2 tables; accepted by ApJL on 2020-04-29

Scattering-induced intensity reduction: large mass content with small grains in the inner region of the TW Hya disk

Authors: Takahiro Ueda, Akimasa Kataoka, Takashi Tsukagoshi

Abstract: Dust continuum observation is one of the best methods to constrain the properties of protoplanetary disks. Recent theoretical studies have suggested that the dust scattering at the millimeter wavelength potentially reduces the observed intensity, which results in an underestimate in the dust mass. We investigate whether the dust scattering indeed reduces the observed continuum intensity by compari… ▽ More Dust continuum observation is one of the best methods to constrain the properties of protoplanetary disks. Recent theoretical studies have suggested that the dust scattering at the millimeter wavelength potentially reduces the observed intensity, which results in an underestimate in the dust mass. We investigate whether the dust scattering indeed reduces the observed continuum intensity by comparing the ALMA archival data of the TW Hya disk at Band 3, 4, 6, 7 and 9 to models obtained by radiative transfer simulations. We find that the model with scattering by 300 ${\rm μm}$-sized grains well reproduces the observed SED of the central part of the TW Hya disk while the model without scattering is also consistent within the errors of the absolute fluxes. To explain the intensity at Band 3, the dust surface density needs to be $\sim$ 10 ${\rm g\,cm^{-2}}$ at 10 au in the model with scattering, which is 26 times more massive than previously predicted. The model without scattering needs 2.3 times higher dust mass than the model with scattering because it needs lower temperature. At Band 7, scattering reduces the intensity by $\sim$ 35% which makes the disk looks optically thin even though it is optically thick. Our study suggests the TW Hya disk is still capable of forming cores of giant planets at where the current solar system planets exist. △ Less

Submitted 19 March, 2020; originally announced March 2020.

Comments: 17 pages, 15 figures, accepted for publication in ApJ

The detection of dust gap-ring structure in the outer region of the CR Cha protoplanetary disk

Authors: Seongjoong Kim, Sanemichi Takahashi, Hideko Nomura, Takashi Tsukagoshi, Seokho Lee, Takayuki Muto, Ruobing Dong, Yasuhiro Hasegawa, Jun Hashimoto, Kazuhiro Kanagawa, Akimasa Kataoka, Mihoko Konishi, Hauyu Baobab Liu, Munetake Momose, Michael Sitko, Kengo Tomida

Abstract: We observe the dust continuum at 225 GHz and CO isotopologue (12CO, 13CO, and C18O) J=2-1 emission lines toward the CR Cha protoplanetary disk using the Atacama Large Millimeter/Submillimeter Array (ALMA). The dust continuum image shows a dust gap-ring structure in the outer region of the dust disk. A faint dust ring is also detected around 120 au beyond the dust gap. The CO isotopologue lines ind… ▽ More We observe the dust continuum at 225 GHz and CO isotopologue (12CO, 13CO, and C18O) J=2-1 emission lines toward the CR Cha protoplanetary disk using the Atacama Large Millimeter/Submillimeter Array (ALMA). The dust continuum image shows a dust gap-ring structure in the outer region of the dust disk. A faint dust ring is also detected around 120 au beyond the dust gap. The CO isotopologue lines indicate that the gas disk is more extended than the dust disk. The peak brightness temperature of the 13CO line shows a small bump around 130 au while 12CO and C18O lines do not show. We investigate two possible mechanisms for reproducing the observed dust gap-ring structure and a gas temperature bump. First, the observed gap structure can be opened by a Jupiter mass planet using the relation between the planet mass and the gap depth and width. Meanwhile, the radiative transfer calculations based on the observed dust surface density profile show that the observed dust ring could be formed by dust accumulation at the gas temperature bump, that is, the gas pressure bump produced beyond the outer edge of the dust disk. △ Less

Submitted 29 November, 2019; originally announced November 2019.

Comments: 21 pages, figures

Radial variations in grain sizes and dust scale heights in the protoplanetary disk around HD 163296 revealed by ALMA polarization observation

Authors: Satoshi Ohashi, Akimasa Kataoka

Abstract: The HD 163296 disk shows ring and gap substructures with ALMA observations. In addition, this is the only disk where the rings and gaps are spatially resolved in millimeter-wave polarization measurements. In this paper, we conduct radiative transfer modeling that includes self-scattering polarization to constrain the grain size and its distribution. We found that the grain size and dust scale heig… ▽ More The HD 163296 disk shows ring and gap substructures with ALMA observations. In addition, this is the only disk where the rings and gaps are spatially resolved in millimeter-wave polarization measurements. In this paper, we conduct radiative transfer modeling that includes self-scattering polarization to constrain the grain size and its distribution. We found that the grain size and dust scale height are the key parameters for reproducing the radial and azimuthal distributions of the observed polarization signature. Radial variation is mainly determined by grain size. The polarization fraction is high if the particle size is $\sim λ/2π$; it is low if the particle size is larger or smaller than this. In contrast, azimuthal variation in polarization is enhanced if the dust scale height is increased. Based on detailed modeling of the HD 163296 polarization, we found the following radial variations in the grain size and dust scale height. The maximum grain size was 140 microns in the gaps and significantly larger or smaller in the rings. The dust scale height is less than one-third the gas scale height inside the 70 au ring, and two-thirds the gas scale height outside the 70 au ring. Furthermore, we constrained the gas turbulence to be $α\lesssim1.5\times10^{-3}$ in the 50 au gap and $α\sim 0.015-0.3$ in the 90 au gap. The transition of the turbulence strength at the boundary of the 70 au ring indicates the existence of a dead zone. △ Less

Submitted 28 October, 2019; originally announced October 2019.

Comments: 24 pages, 21 figures, accepted for publication in ApJ

Investigating the gas-to-dust ratio in the protoplanetary disk of HD 142527

Authors: Kang-Lou Soon, Munetake Momose, Takayuki Muto, Takashi Tsukagoshi, Akimasa Kataoka, Tomoyuki Hanawa, Misato Fukagawa, Kazuya Saigo, Hiroshi Shibai

Abstract: We present ALMA observations of the $98.5~\mathrm{GHz}$ dust continuum and the $\mathrm{^{13}CO}~J = 1 - 0$ and $\mathrm{C^{18}O}~J = 1 - 0$ line emissions of the protoplanetary disk associated with HD~142527. The $98.5~\mathrm{GHz}$ continuum shows a strong azimuthal-asymmetric distribution similar to that of the previously reported $336~\mathrm{GHz}$ continuum, with a peak emission in dust conce… ▽ More We present ALMA observations of the $98.5~\mathrm{GHz}$ dust continuum and the $\mathrm{^{13}CO}~J = 1 - 0$ and $\mathrm{C^{18}O}~J = 1 - 0$ line emissions of the protoplanetary disk associated with HD~142527. The $98.5~\mathrm{GHz}$ continuum shows a strong azimuthal-asymmetric distribution similar to that of the previously reported $336~\mathrm{GHz}$ continuum, with a peak emission in dust concentrated region in the north. The disk is optically thin in both the $98.5~\mathrm{GHz}$ dust continuum and the $\mathrm{C^{18}O}~J = 1 - 0$ emissions. We derive the distributions of gas and dust surface densities, $Σ_\mathrm{g}$ and $Σ_\mathrm{d}$, and the dust spectral opacity index, $β$, in the disk from ALMA Band 3 and Band 7 data. In the analyses, we assume the local thermodynamic equilibrium and the disk temperature to be equal to the peak brightness temperature of $\mathrm{^{13}CO}~J = 3 - 2$ with a continuum emission. The gas-to-dust ratio, $\mathrm{G/D}$, varies azimuthally with a relation $\mathrm{G/D} \propto Σ_\mathrm{d}^{-0.53}$, and $β$ is derived to be $\approx 1$ and $\approx 1.7$ in the northern and southern regions of the disk, respectively. These results are consistent with the accumulation of larger dust grains in a higher pressure region. In addition, our results show that the peak $Σ_\mathrm{d}$ is located ahead of the peak $Σ_\mathrm{g}$. If the latter corresponds to a vortex of high gas pressure, the results indicate that the dust is trapped ahead of the vortex, as predicted by some theoretical studies. △ Less

Submitted 15 September, 2019; originally announced September 2019.

Comments: 32 pages, 15 figures

An Observational Study for Grain Dynamics in the AS 209 Disk with Submillimeter Polarization

Authors: Tomohiro Mori, Akimasa Kataoka, Satoshi Ohashi, Munetake Momose, Takayuki Muto, Hiroshi Nagai, Takashi Tsukagoshi

Abstract: We present 870 $μ$m ALMA polarization observation toward the Class I\hspace{-.1em}I protoplanetary disk around AS 209, which has concentric, multiple gaps and rings. We successfully detect the polarized emission and find that the polarization orientations and fractions have distinct characteristics between the inner and outer regions. In the inner region, the polarization orientations are parallel… ▽ More We present 870 $μ$m ALMA polarization observation toward the Class I\hspace{-.1em}I protoplanetary disk around AS 209, which has concentric, multiple gaps and rings. We successfully detect the polarized emission and find that the polarization orientations and fractions have distinct characteristics between the inner and outer regions. In the inner region, the polarization orientations are parallel to the minor axis of the disk, which is consistent with the self-scattering model. The mean polarization fraction in the region is $\sim$0.2\%, which is lower than the expected value when the grains have the maximum polarization efficiency, which corresponds to $λ$/2$π$ $\sim$ 140 $μ$m in grain radius. In the outer region, we detect $\sim$1.0\% polarization and find that the polarization orientations are almost in the azimuthal directions. Moreover, the polarization orientations have systematic angular deviations from the azimuthal directions with $Δ$$θ$ $\sim$ 4\fdg5 $\pm$ 1\fdg6. The pattern is consistent with a model that radially drifting dust grains are aligned by the gas flow against the dust grains. We consider possible scenarios of the grain dynamics at the AS 209 ring which can reproduce the polarization pattern. However, the directions of the observed angular deviations are opposite to what is predicted under the fact that the disk rotates clockwise. This poses a question in our understandings of the alignment processes and/or grain dynamics in protoplanetary disks. △ Less

Submitted 24 July, 2019; originally announced July 2019.

Comments: 14 pages, 7 figures

Unveiling Dust Aggregate Structure in Protoplanetary Disks by Millimeter-wave Scattering Polarization

Authors: Ryo Tazaki, Hidekazu Tanaka, Akimasa Kataoka, Satoshi Okuzumi, Takayuki Muto

Abstract: Dust coagulation in a protoplanetary disk is the first step of planetesimal formation. However, a pathway from dust aggregates to planetesimals remains unclear. Both numerical simulations and laboratory experiments have suggested the importance of dust structure in planetesimal formation, but it is not well constrained by observations. We study how dust structure and porosity alters polarimetric i… ▽ More Dust coagulation in a protoplanetary disk is the first step of planetesimal formation. However, a pathway from dust aggregates to planetesimals remains unclear. Both numerical simulations and laboratory experiments have suggested the importance of dust structure in planetesimal formation, but it is not well constrained by observations. We study how dust structure and porosity alters polarimetric images at millimeter wavelength by performing 3D radiative transfer simulations. Aggregates with different porosity and fractal dimension are considered. As a result, we find that dust aggregates with lower porosity and/or higher fractal dimension are favorable to explain observed millimeter-wave scattering polarization of disks. Aggregates with extremely high porosity fail to explain the observations. In addition, we also show that particles with moderate porosity show weak wavelength dependence of scattering polarization, indicating that multi-wavelength polarimetry is useful to constrain dust porosity. Finally, we discuss implications for dust evolution and planetesimal formation in disks. △ Less

Submitted 19 September, 2019; v1 submitted 29 June, 2019; originally announced July 2019.

Comments: 17 pages, 10 Figures, 1 Table, accepted for publication in ApJ

Dust Polarization in Four Protoplanetary Disks at 3 mm: Further Evidence of Multiple Origins

Authors: Rachel E. Harrison, Leslie W. Looney, Ian W. Stephens, Zhi-Yun Li, Haifeng Yang, Akimasa Kataoka, Robert J. Harris, Woojin Kwon, Takayuki Muto, Munetake Momose

Abstract: We present polarimetric observations of four Class II protoplanetary disks (DG Tau, Haro 6-13, RY Tau, and MWC 480) taken with the Atacama Large Millimeter/submillimeter Array (ALMA) at 3 mm. The polarization morphologies observed fall into two distinct categories: azimuthal polarization (DG Tau and Haro 6-13) and polarization parallel to the disk minor axis (RY Tau and MWC 480). The mechanisms re… ▽ More We present polarimetric observations of four Class II protoplanetary disks (DG Tau, Haro 6-13, RY Tau, and MWC 480) taken with the Atacama Large Millimeter/submillimeter Array (ALMA) at 3 mm. The polarization morphologies observed fall into two distinct categories: azimuthal polarization (DG Tau and Haro 6-13) and polarization parallel to the disk minor axis (RY Tau and MWC 480). The mechanisms responsible for disk polarization at millimeter and submillimeter wavelengths are currently under debate. In this Letter, we investigate two mechanisms capable of producing polarized dust emission in disks: self-scattering and grain alignment to the radiation anisotropy. The polarization morphologies of DG Tau and Haro 6-13 are broadly consistent with that expected from radiation alignment (though radiative alignment still does not account for all of the features seen in these disks), while RY Tau and MWC 480 are more consistent with self-scattering. Such a variation in the polarized morphology may provide evidence of dust grain size differences between the sources. △ Less

Submitted 22 June, 2021; v1 submitted 15 May, 2019; originally announced May 2019.

Comments: 10 pages, 5 figures, corrected ALMA project code and corrected formatting of Bacciotti et al. (2018) citation. January 27 2021: Corrected beam sizes for RY Tau and MWC 480. June 22 2021: Corrected percent polarization scale bars in Figures 1, 4, and 5

Journal ref: ApJL 877 (2019) L2

The JCMT BISTRO Survey: The Magnetic Field of the Barnard 1 Star-Forming Region

Authors: Simon Coudé, Pierre Bastien, Martin Houde, Sarah Sadavoy, Rachel Friesen, James Di Francesco, Doug Johnstone, Steve Mairs, Tetsuo Hasegawa, Woojin Kwon, Shih-Ping Lai, Keping Qiu, Derek Ward-Thompson, David Berry, Michael Chun-Yuan Chen, Jason Fiege, Erica Franzmann, Jennifer Hatchell, Kevin Lacaille, Brenda C. Matthews, Gerald H. Moriarty-Schieven, Andy Pon, Philippe André, Doris Arzoumanian, Yusuke Aso , et al. (96 additional authors not shown)

Abstract: We present the POL-2 850 $μ$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We th… ▽ More We present the POL-2 850 $μ$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $μ$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $μ$G and a criticality criterion $λ_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component. △ Less

Submitted 23 April, 2019; v1 submitted 15 April, 2019; originally announced April 2019.

Comments: 22 pages, 7 figures, accepted in the Astrophysical Journal

One Solution to the Mass Budget Problem for Planet Formation: Optically Thick Disks with Dust Scattering

Authors: Zhaohuan Zhu, Shangjia Zhang, Yan-Fei Jiang, Akimasa Kataoka, Tilman Birnstiel, Cornelis P. Dullemond, Sean M. Andrews, Jane Huang, Laura M. Perez, John M. Carpenter, Xue-Ning Bai, David J. Wilner, Luca Ricci

Abstract: ALMA surveys have suggested that the dust in Class II disks may not be enough to explain the averaged solid mass in exoplanets, under the assumption that the mm disk continuum emission is optically thin. This optically thin assumption seems to be supported by recent DSHARP observations where the measured optical depths of spatially resolved disks are mostly less than one. However, we point out tha… ▽ More ALMA surveys have suggested that the dust in Class II disks may not be enough to explain the averaged solid mass in exoplanets, under the assumption that the mm disk continuum emission is optically thin. This optically thin assumption seems to be supported by recent DSHARP observations where the measured optical depths of spatially resolved disks are mostly less than one. However, we point out that dust scattering can considerably reduce the emission from an optically thick region. If that scattering is ignored, the optical depth will be considerably underestimated. An optically thick disk with scattering can be misidentified as an optically thin disk. Dust scattering in more inclined disks can reduce the intensity even further, making the disk look even fainter. The measured optical depth of $\sim$0.6 in several DSHARP disks can be naturally explained by optically thick dust with an albedo of $\sim$0.9 at 1.25 mm. Using the DSHARP opacity, this albedo corresponds to a dust population with the maximum grain size ($s_{max}$) of 0.1-1 mm. For optically thick scattering disks, the measured spectral index $α$ can be either larger or smaller than 2 depending on if the dust albedo increases or decreases with wavelength. Using the DSHARP opacity, $α<2$ corresponds to $s_{max}$ of 0.03-0.3 mm. We describe how this optically thick scattering scenario could explain the observed scaling between submm continuum sizes and luminosities, and might help ease the tension between the dust size constraints from polarization and dust continuum measurements. We suggest that a significant amount of disk mass can be hidden from ALMA observations at short millimeter wavelengths. For compact disks smaller than 30 au, we can easily underestimate the dust mass by more than a factor of 10. Longer wavelength observations (e.g. VLA or SKA) are desired to probe the dust mass in disks. △ Less

Submitted 7 May, 2019; v1 submitted 3 April, 2019; originally announced April 2019.

Comments: 15 pages, 11 figures, Accepted for publication in ApJL

Polarization in Disks

Authors: Ian Stephens, Zhi-Yun Li, Haifeng Yang, Akimasa Kataoka, Leslie W. Looney, Charles L. H. Hull, Manuel Fernández-López, Sarah I. Sadavoy, Woojin Kwon, Satoshi Ohashi, Ryo Tazaki, Dan Li, Thiem Hoang, Gesa H. -M. Bertrang, Carlos Carrasco-González, William R. F. Dent, Satoko Takahashi, Francesca Bacciotti, Felipe O. Alves, Josep M. Girart, Qizhou Zhang, Ramprasad Rao, Adriana Pohl, Marco Padovani, Daniele Galli , et al. (2 additional authors not shown)

Abstract: Polarized dust emission outside of disks reveal the magnetic field morphology of molecular clouds. Within disks, however, polarized dust emission can arise from very different mechanisms (e.g., self-scattering), and each of them are useful for constraining physical properties in the disk. For example, these mechanisms allow us to constrain the disk grain size distributions and grain/disk geometrie… ▽ More Polarized dust emission outside of disks reveal the magnetic field morphology of molecular clouds. Within disks, however, polarized dust emission can arise from very different mechanisms (e.g., self-scattering), and each of them are useful for constraining physical properties in the disk. For example, these mechanisms allow us to constrain the disk grain size distributions and grain/disk geometries, independent from current methods of measuring these parameters. To accurately model these features and disentangle the various polarization mechanisms, multiwavelength observations at very high resolution and sensitivity are required. With significant upgrades to current interferometric facilities, we can understand how grains evolve in disks during the planet formation process. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey

Millimeter-wave polarization due to grain alignment by the gas flow in protoplanetary disks

Authors: Akimasa Kataoka, Satoshi Okuzumi, Ryo Tazaki

Abstract: Dust grains emit intrinsic polarized emission if they are elongated and aligned in the same direction. The direction of the grain alignment is determined by external forces, such as magnetic fields, radiation, and gas flow against the dust grains. In this letter, we apply the concept of the grain alignment by gas flow, which is called mechanical alignment, to the situation of a protoplanetary disk… ▽ More Dust grains emit intrinsic polarized emission if they are elongated and aligned in the same direction. The direction of the grain alignment is determined by external forces, such as magnetic fields, radiation, and gas flow against the dust grains. In this letter, we apply the concept of the grain alignment by gas flow, which is called mechanical alignment, to the situation of a protoplanetary disk. We assume that grains have a certain helicity, which results in the alignment with the minor axis parallel to the grain velocity against the ambient disk gas and discuss the morphology of polarization vectors in a protoplanetary disk. We find that the direction of the polarization vectors depends on the Stokes number, which denotes how well grains are coupled to the gas. If the Stokes number is less than unity, orientation of polarization is in the azimuthal direction since the dust velocity against the gas is in the radial direction. If the Stokes number is as large as unity, the polarization vectors show a leading spiral pattern since the radial and azimuthal components of the gas velocity against the dust grains are comparable. This suggests that if the observed polarization vectors show a leading spiral pattern, it would indicate that Stokes number of dust grains is around unity, which is presumably radially drifting. △ Less

Submitted 8 March, 2019; originally announced March 2019.

Comments: 7 pages, 5 figures, accepted for publication in The Astrophysical Journal Letters

Effect of dust size and structure on scattered light images of protoplanetary discs

Authors: Ryo Tazaki, Hidekazu Tanaka, Takayuki Muto, Akimasa Kataoka, Satoshi Okuzumi

Abstract: We study scattered light properties of protoplanetary discs at near-infrared wavelengths for various dust size and structure by performing radiative transfer simulations. We show that different dust structures might be probed by measuring disk polarisation fraction as long as the dust radius is larger than the wavelength. When the radius is larger than the wavelength, disc scattered light will be… ▽ More We study scattered light properties of protoplanetary discs at near-infrared wavelengths for various dust size and structure by performing radiative transfer simulations. We show that different dust structures might be probed by measuring disk polarisation fraction as long as the dust radius is larger than the wavelength. When the radius is larger than the wavelength, disc scattered light will be highly polarised for highly porous dust aggregates, whereas more compact dust structure tends to show low polarisation fraction. Next, roles of monomer radius and fractal dimension for scattered light colours are studied. We find that, outside the Rayleigh regime, as fractal dimension or monomer radius increases, colours of the effective albedo at near-infrared wavelengths vary from blue to red. Our results imply that discs showing grey or slightly blue colours and high polarisation fraction in near-infrared wavelengths might be explained by the presence of large porous aggregates containing sub-microns sized monomers. △ Less

Submitted 5 March, 2019; originally announced March 2019.

Comments: Accepted for publication in MNRAS, 18 pages, 19 figures

Tensile Strength of Porous Dust Aggregates

Authors: Misako Tatsuuma, Akimasa Kataoka, Hidekazu Tanaka

Abstract: Comets are thought to have information about the formation process of our solar system. Recently, detailed information about comet 67P/Churyumov-Gerasimenko has been found by a spacecraft mission Rosetta. It is remarkable that its tensile strength was estimated. In this paper, we measure and formulate the tensile strength of porous dust aggregates using numerical simulations, motivated by porous d… ▽ More Comets are thought to have information about the formation process of our solar system. Recently, detailed information about comet 67P/Churyumov-Gerasimenko has been found by a spacecraft mission Rosetta. It is remarkable that its tensile strength was estimated. In this paper, we measure and formulate the tensile strength of porous dust aggregates using numerical simulations, motivated by porous dust aggregation model of planetesimal formation. We perform three-dimensional numerical simulations using a monomer interaction model with periodic boundary condition. We stretch out a dust aggregate with a various initial volume filling factor between $10^{-2}$ and 0.5. We find that the tensile stress takes the maximum value at the time when the volume filling factor decreases to about a half of the initial value. The maximum stress is defined to be the tensile strength. We take an average of the results with 10 different initial shapes to smooth out the effects of initial shapes of aggregates. Finally, we numerically obtain the relation between the tensile strength and the initial volume filling factor of dust aggregates. We also use a simple semi-analytical model and successfully reproduce the numerical results, which enables us to apply to a wide parameter range and different materials. The obtained relation is consistent with previous experiments and numerical simulations about silicate dust aggregates. We estimate that the monomer radius of comet 67P has to be about 3.3--220 $\mathrm{μm}$ to reproduce its tensile strength using our model. △ Less

Submitted 21 February, 2019; originally announced February 2019.

Comments: 12 pages, 12 figures, accepted for publication in ApJ

The JCMT BISTRO Survey: The Magnetic Field In The Starless Core $ρ$ Ophiuchus C

Authors: Junhao Liu, Keping Qiu, David Berry, James Di Francesco, Pierre Bastien, Patrick M. Koch, Ray S. Furuya, Kee-Tae Kim, Simon Coudé, Chang Won Lee, Archana Soam, Chakali Eswaraiah, Di Li, Jihye Hwang, A-Ran Lyo, Kate Pattle, Tetsuo Hasegawa, Woojin Kwon, Shih-Ping Lai, Derek Ward-Thompson, Tao-Chung Ching, Zhiwei Chen, Qilao Gu, Dalei Li, Hua-bai Li , et al. (106 additional authors not shown)

Abstract: We report 850~$μ$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $ρ$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic… ▽ More We report 850~$μ$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $ρ$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $μ$G, 136 $\pm$ 69 $μ$G, and 213 $\pm$ 115 $μ$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties. △ Less

Submitted 20 February, 2019; originally announced February 2019.

Comments: 22 pages, 12 figures

JCMT BISTRO survey: Magnetic Fields within the Hub-Filament Structure in IC 5146

Authors: Jia-Wei Wang, Shih-Ping Lai, Chakali Eswaraiah, Kate Pattle, James Di Francesco, Doug Johnstone, Patrick M. Koch, Tie Liu, Motohide Tamura, Ray S. Furuya, Takashi Onaka, Derek Ward-Thompson, Archana Soam, Kee-Tae Kim, Chang Won Lee, Chin-Fei Lee, Steve Mairs, Doris Arzoumanian, Gwanjeong Kim, Thiem Hoang, Jihye Hwang, Sheng-Yuan Liu, David Berry, Pierre Bastien, Tetsuo Hasegawa , et al. (108 additional authors not shown)

Abstract: We present the 850 $μ$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scal… ▽ More We present the 850 $μ$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfvénic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important. △ Less

Submitted 27 March, 2019; v1 submitted 14 December, 2018; originally announced December 2018.

Comments: 24 pages, 15 figures, accepted for publication in ApJ

Does HL Tau Disk Polarization in ALMA Band 3 Come from Radiatively Aligned Grains?

Authors: Haifeng Yang, Zhi-Yun Li, Ian W. Stephens, Akimasa Kataoka, Leslie Looney

Abstract: Disk polarization in (sub)millimeter dust continuum is a rapidly growing field in the ALMA era. It opens up the exciting possibility of detecting and characterizing magnetic fields and grain growth in disks around young stellar objects. However, to use polarization for probing the disk properties, its production mechanism must be ascertained first. To date, the conventional mechanism involving mag… ▽ More Disk polarization in (sub)millimeter dust continuum is a rapidly growing field in the ALMA era. It opens up the exciting possibility of detecting and characterizing magnetic fields and grain growth in disks around young stellar objects. However, to use polarization for probing the disk properties, its production mechanism must be ascertained first. To date, the conventional mechanism involving magnetically aligned grains fails to explain the polarization patterns detected in most disks. This is especially true for the inclined disk of HL Tau in ALMA Band 3 (wavelength $\sim 3$ mm), which has an elliptical polarization pattern. The elliptical pattern was taken as evidence for polarized emission by dust grains aligned with their long axes perpendicular the direction of the radiative flux. We show that the radiatively aligned grains produce a circular, rather than elliptical, polarization pattern even in inclined disks such as HL Tau. An elliptical polarization pattern can be produced if the grains are aligned aerodynamically by the difference in rotation speed between the dust and gas through the Gold mechanism. However, a strong azimuthal variation in polarized intensity is expected for both the radiative and aerodynamic alignment, but not observed in the HL Tau disk in ALMA Band 3. We conclude that neither of these two mechanisms alone can explain the data and the origin of the 3 mm polarization remains a mystery. We speculate that this mystery may be resolved by a combination of both direct emission and scattering by aerodynamically aligned grains. △ Less

Submitted 28 November, 2018; originally announced November 2018.

Comments: 12 pages, 12 figures, accepted by MNRAS

Two Different Grain Size Distributions within the Protoplanetary Disk around HD 142527 Revealed by ALMA Polarization Observation

Authors: Satoshi Ohashi, Akimasa Kataoka, Hiroshi Nagai, Munetake Momose, Takayuki Muto, Tomoyuki Hanawa, Misato Fukagawa, Takashi Tsukagoshi, Kohji Murakawa, Hiroshi Shibai

Abstract: The origin of polarized emission from protoplanetary disks is uncertain. Three mechanisms have been proposed for such polarized emission: grain alignment with magnetic fields, grain alignment with radiation gradients, and self-scattering of thermal dust emission. Aiming to observationally identify the polarization mechanisms, we present ALMA polarization observations of the 0.87 mm dust continuum… ▽ More The origin of polarized emission from protoplanetary disks is uncertain. Three mechanisms have been proposed for such polarized emission: grain alignment with magnetic fields, grain alignment with radiation gradients, and self-scattering of thermal dust emission. Aiming to observationally identify the polarization mechanisms, we present ALMA polarization observations of the 0.87 mm dust continuum emission toward the circumstellar disk around HD 142527 with high spatial resolution. We confirm that the polarization vectors in the northern region are consistent with self-scattering. Furthermore, we show that the polarization vectors in the southern region are consistent with grain alignment by magnetic fields, although self-scattering cannot be ruled out. To understand the differences between the polarization mechanisms, we propose a simple grain size segregation model: small dust grains ($\lesssim$ 100 microns) are dominant and aligned with magnetic fields in the southern region, and middle-sized ($\sim100$ microns) grains in the upper layer emit self-scattered polarized emission in the northern region. The grain size near the middle plane in the northern region cannot be measured because the emission at 0.87 mm is optically thick. However, it can be speculated that larger dust grains ($\gtrsim$ cm) may accumulate near this plane. These results are consistent with those of a previous analysis of the disk, in which large grain accumulation and optically thick emission from the northern region were found. This model is also consistent with theories where smaller dust grains are aligned with magnetic fields. The magnetic fields are toroidal, at least in the southern region. △ Less

Submitted 27 July, 2018; originally announced July 2018.

Comments: 24 pages, 21 figures, 1 table, accepted for publication in ApJ

Magnetic fields towards Ophiuchus-B derived from SCUBA-2 polarization measurements

Authors: Archana Soam, Kate Pattle, Derek Ward-Thompson, Chang Won Lee, Sarah Sadavoy, Patrick M. Koch, Gwanjeong Kim, Jungmi Kwon, Woojin Kwon, Doris Arzoumanian, David Berry, Thiem Hoang, Motohide Tamura, Sang-Sung Lee, Tie Liu, Kee-Tae Kim, Doug Johnstone, Fumitaka Nakamura, A-Ran Lyo, Takashi Onaka, Jongsoo Kim, Ray S. Furuya, Tetsuo Hasegawa, Shih-Ping Lai, Pierre Bastien , et al. (99 additional authors not shown)

Abstract: We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fi… ▽ More We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $μ$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $λ$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical. △ Less

Submitted 7 June, 2018; v1 submitted 16 May, 2018; originally announced May 2018.

Comments: 17 pages, 10 figures, Accepted for publication in ApJ