Showing 1–50 of 101 results for author: Shang, H

Statistical analysis of pulsar flux density distribution

Authors: H. W. Xu, R. S. Zhao, Erbil Gugercinoglu, H. Liu, D. Li, P. Wang, C. H. Niu, C. Miao, X. Zhu, R. W. Tian, W. L. Li, S. D. Wang, Z. F. Tu, Q. J. Zhi, S. J. Dang, L. H. Shang, S. Xiao

Abstract: This study presents a comprehensive analysis of the spectral properties of 886 pulsars across a wide frequency range from 20MHz to 343.5GHz, including a total of 86 millisecond pulsars. The majority of the pulsars exhibit power-law behavior in their spectra, although some exceptions are observed. Five different spectral models, namely simple power-law, broken power-law, low-frequency turn-over, hi… ▽ More This study presents a comprehensive analysis of the spectral properties of 886 pulsars across a wide frequency range from 20MHz to 343.5GHz, including a total of 86 millisecond pulsars. The majority of the pulsars exhibit power-law behavior in their spectra, although some exceptions are observed. Five different spectral models, namely simple power-law, broken power-law, low-frequency turn-over, high-frequency cut-off, and double turn-over, were employed to explore the spectral behaviors. The average spectral index for pulsars modeled with a simple power-law is found to be -1.64 +/-0.80, consistent with previous studies. Additionally, significant correlations between the spectral index and characteristic parameters are observed particularly in millisecond pulsars, while no strong correlation is observed in normal pulsars. Different models show variations in the most influential characteristic parameters associated with the spectral index, indicating diverse dominant radiation mechanisms in millisecond pulsars.Finally, this study identifies 22 pulsars of the Gigahertz-peaked Spectra (GPS) type for the first time based on the Akaike information criterion. △ Less

Submitted 16 April, 2024; v1 submitted 16 April, 2024; originally announced April 2024.

Comments: 39 papers,17figures

Formation of limb-brightened radio jets by angle-dependent energy extraction from rapidly rotating black holes

Authors: Kouichi Hirotani, Hsien Shang, Ruben Krasnopolsky, Kenichi Nishikawa

Abstract: By general relativistic magnetohydrodynamic simulations, it is suggested that the rotational energy of a rapidly rotating black hole (BH) is preferentially extracted along the magnetic field lines threading the event horizon in the middle and lower latitudes. Applying this angle-dependent Poynting flux to the jet downstream, we demonstrate that the jets exhibit limb-brightened structures at variou… ▽ More By general relativistic magnetohydrodynamic simulations, it is suggested that the rotational energy of a rapidly rotating black hole (BH) is preferentially extracted along the magnetic field lines threading the event horizon in the middle and lower latitudes. Applying this angle-dependent Poynting flux to the jet downstream, we demonstrate that the jets exhibit limb-brightened structures at various viewing angles, as observed from Mrk 501, M87, and Cyg A between 5 and 75 degrees, and that the limb-brightening is enhanced when the jet is collimated strongly. It is also found that the jet width perpendicular to the propagation direction shrinks at the projected distance of the altitude where the jet collimates from a conical shape (near the BH) to a parabolic one (in the jet). Comparing with the VLBI observations, we show this collimation takes place within the de-projected altitude of 100 Schwarzschild radii from the BH in the case of the M87 jet. △ Less

Submitted 6 March, 2024; originally announced March 2024.

Comments: 19 pages, 8 figures. The Astrophysical Journal in press

A Unified Model for Bipolar Outflows from Young Stars: Kinematic and Mixing Structures in HH 30

Authors: Tsung-Han Ai, Chun-Fan Liu, Hsien Shang, Doug Johnstone, Ruben Krasnopolsky

Abstract: The young stellar source HH 30 is a textbook example of an ionic optical jet originating from a disk in an edge-on system shown by the HST. It has a remnant envelope in $^{12}$CO observed by ALMA. The optical jet is characterized by its narrow appearance, large line width at the base, and high temperature inferred from line diagnostics. Three featured structures can be identified, most evident in… ▽ More The young stellar source HH 30 is a textbook example of an ionic optical jet originating from a disk in an edge-on system shown by the HST. It has a remnant envelope in $^{12}$CO observed by ALMA. The optical jet is characterized by its narrow appearance, large line width at the base, and high temperature inferred from line diagnostics. Three featured structures can be identified, most evident in the transverse position--velocity diagrams: an extremely--high-velocity (EHV) wide-angle wind component with large spectral widths in the optical, a very--low-velocity (VLV) ambient surrounding medium seen in $^{12}$CO, and a low-velocity (LV) region traced by $^{12}$CO nested both in velocity and location between the primary wind and ambient environment. A layered cavity with multiple shells forms nested morphological and kinematic structures around the optical jet. The atomic gas originating from the innermost region of the disk attains a sufficient temperature and ionization to emit brightly in forbidden lines as an optical jet. The wide-angle portion expands, forming a low-density cavity. The filamentary $^{12}$CO encompassing the wind cavity is mixed and advected inward through the action of the magnetic interplay of the wide-angle wind with the molecular ambient medium. The magnetic interplay results in the layered shells penetrating deeply into the vast cavity of tenuous atomic wind material. The HH 30 system is an ideal manifestation of the unified wind model of \citet{Shang_2020,Shang_2023}, with clearly distinguishable atomic and molecular species mixed through the atomic lightly ionized magnetized wind and the surrounding cold molecular ambient material. △ Less

Submitted 4 February, 2024; originally announced February 2024.

Comments: 19 pages, 9 figures, ApJ in press

An Efficient Algorithm for Astrochemical Systems Using Stoichiometry Matrices

Authors: Kazutaka Motoyama, Ruben Krasnopolsky, Hsien Shang, Kento Aida, Eisaku Sakane

Abstract: Astrochemical simulations are a powerful tool for revealing chemical evolution in the interstellar medium. Astrochemical calculations require efficient processing of large matrices for the chemical networks. The large chemical reaction networks often present bottlenecks for computation because of time derivatives of chemical abundances. We propose an efficient algorithm using a stoichiometry matri… ▽ More Astrochemical simulations are a powerful tool for revealing chemical evolution in the interstellar medium. Astrochemical calculations require efficient processing of large matrices for the chemical networks. The large chemical reaction networks often present bottlenecks for computation because of time derivatives of chemical abundances. We propose an efficient algorithm using a stoichiometry matrix approach in which this time-consuming part is expressed as a loop, unlike the algorithm used in previous studies. Since stoichiometry matrices are sparse in general, the performances of simulations with our algorithm depend on which sparse-matrix storage format is used. We conducted a performance comparison experiment using the common storage formats, including the coordinate (COO) format, the compressed column storage (CCS) format, the compressed row storage (CRS) format, and the Sliced ELLPACK (SELL) format. Experimental results showed that the simulations with the CRS format are the most suitable for astrochemical simulations and about three times faster than those with the algorithm used in previous studies. In addition, our algorithm significantly reduces not only the computation time but also the compilation time. We also explore the beneficial effects of parallelization and sparse-matrix reordering in these algorithms. △ Less

Submitted 8 December, 2023; originally announced December 2023.

Comments: 10 pages, 7 figures, accepted for publication in ApJS

Discovery of four pulsars in a pilot survey at intermediate Galactic latitudes with FAST

Authors: Q. J. Zhi, J. T. Bai, S. Dai, X. Xu, S. J. Dang, L. H. Shang, R. S. Zhao, D. Li, W. W. Zhu, N. Wang, J. P. Yuan, P. Wang, L. Zhang, Y. Feng, J. B. Wang, S. Q. Wang, Q. D. Wu, A. J. Dong, H. Yang, J. Tian, W. Q. Zhong, X. H. Luo, Miroslav D. Filipovi, G. J. Qiao

Abstract: We present the discovery and timing results of four pulsars discovered in a pilot survey at intermediate Galactic latitudes with the Five-hundred Aperture Spherical Telescope (FAST). Among these pulsars, two belong to the category of millisecond pulsars (MSPs) with spin periods of less than 20 ms. The other two fall under the classification of "mildly recycled" pulsars, with massive white dwarfs a… ▽ More We present the discovery and timing results of four pulsars discovered in a pilot survey at intermediate Galactic latitudes with the Five-hundred Aperture Spherical Telescope (FAST). Among these pulsars, two belong to the category of millisecond pulsars (MSPs) with spin periods of less than 20 ms. The other two fall under the classification of "mildly recycled" pulsars, with massive white dwarfs as companions. Remarkably, this small survey, covering an area of 4.7 $deg^2$ , led to the discovery of four recycled pulsars. Such success underscores the immense potential of future surveys at intermediate Galactic latitudes. In order to assess the potential yield of MSPs, we conducted population simulations and found that both FAST and Parkes new phased array feed surveys, focusing on intermediate Galactic latitudes, have the capacity to uncover several hundred new MSPs. △ Less

Submitted 28 December, 2023; v1 submitted 1 November, 2023; originally announced November 2023.

Comments: 7 pages, 4 figures, 2 tables, accepted to ApJ

Exploring the Formation of Resistive Pseudodisks with the GPU Code Astaroth

Authors: Miikka S. Väisälä, Hsien Shang, Daniele Galli, Susana Lizano, Ruben Krasnopolsky

Abstract: Pseudodisks are dense structures formed perpendicular to the direction of the magnetic field during the gravitational collapse of a molecular cloud core. Numerical simulations of the formation of pseudodisks are usually computationally expensive with conventional CPU codes. To demonstrate the proof-of-concept of a fast computing method for this numerically costly problem, we explore the GPU-powere… ▽ More Pseudodisks are dense structures formed perpendicular to the direction of the magnetic field during the gravitational collapse of a molecular cloud core. Numerical simulations of the formation of pseudodisks are usually computationally expensive with conventional CPU codes. To demonstrate the proof-of-concept of a fast computing method for this numerically costly problem, we explore the GPU-powered MHD code Astaroth, a 6th-order finite difference method with low adjustable finite resistivity implemented with sink particles. The formation of pseudodisks is physically and numerically robust and can be achieved with a simple and clean setup for this newly adopted numerical approach for science verification. The method's potential is illustrated by evidencing the dependence on the initial magnetic field strength of specific physical features accompanying the formation of pseudodisks, e.g. the occurrence of infall shocks and the variable behavior of the mass and magnetic flux accreted on the central object. As a performance test, we measure both weak and strong scaling of our implementation to find most efficient way to use the code on a multi-GPU system. Once suitable physics and problem-specific implementations are realized, the GPU-accelerated code is an efficient option for 3-D magnetized collapse problems. △ Less

Submitted 25 October, 2023; originally announced October 2023.

Comments: 29 pages, 1 table, 15 figures, Accepted for publication in the Astrophysical Journal

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): The Warm-Envelope Origin of Hot Corinos

Authors: Shih-Ying Hsu, Sheng-Yuan Liu, Doug Johnstone, Tie Liu, Leonardo Bronfman, Huei-Ru Vivien Chen, Somnath Dutta, David J. Eden, Neal J. Evans II, Naomi Hirano, Mika Juvela, Yi-Jehng Kuan, Woojin Kwon, Chin-Fei Lee, Chang Won Lee, Jeong-Eun Lee, Shanghuo Li, Chun-Fan Liu, Xunchuan Liu, Qiuyi Luo, Sheng-Li Qin, Mark G. Rawlings, Dipen Sahu, Patricio Sanhueza, Hsien Shang , et al. (2 additional authors not shown)

Abstract: Hot corinos are of great interest due to their richness in interstellar complex organic molecules (COMs) and the consequent potential prebiotic connection to solar-like planetary systems. Recent surveys have reported an increasing number of hot corino detections in Class 0/I protostars; however, the relationships between their physical properties and the hot-corino signatures remain elusive. In th… ▽ More Hot corinos are of great interest due to their richness in interstellar complex organic molecules (COMs) and the consequent potential prebiotic connection to solar-like planetary systems. Recent surveys have reported an increasing number of hot corino detections in Class 0/I protostars; however, the relationships between their physical properties and the hot-corino signatures remain elusive. In this study, our objective is to establish a general picture of the detectability of the hot corinos by identifying the origin of the hot-corino signatures in the sample of young stellar objects (YSOs) obtained from the Atacama Large Millimeter/submillimeter Array Survey of Orion Planck Galactic Cold Clumps (ALMASOP) project. We apply spectral energy distribution (SED) modeling to our sample and identify the physical parameters of the modeled YSOs directly, linking the detection of hot-corino signatures to the envelope properties of the YSOs. Imaging simulations of the methanol emission further support this scenario. We, therefore, posit that the observed COM emission originates from the warm inner envelopes of the sample YSOs, based on both the warm region size and the envelope density profile. The former is governed by the source luminosity and is additionally affected by the disk and cavity properties, while the latter is related to the evolutionary stages. This scenario provides a framework for detecting hot-corino signatures toward luminous Class 0 YSOs, with fewer detections observed toward similarly luminous Class I sources. △ Less

Submitted 11 August, 2023; v1 submitted 10 August, 2023; originally announced August 2023.

Comments: 28 pages, 11 figures

TOI-4860 b, a short-period giant planet transiting an M3.5 dwarf

Authors: J. M. Almenara, X. Bonfils, E. M. Bryant, A. Jordán, G. Hébrard, E. Martioli, A. C. M. Correia, N. Astudillo-Defru, C. Cadieux, L. Arnold, É. Artigau, G. Á. Bakos, S. C. C. Barros, D. Bayliss, F. Bouchy, G. Boué, R. Brahm, A. Carmona, D. Charbonneau, D. R. Ciardi, R. Cloutier, M. Cointepas, N. J. Cook, N. B. Cowan, X. Delfosse , et al. (25 additional authors not shown)

Abstract: We report the discovery and characterisation of a giant transiting planet orbiting a nearby M3.5V dwarf (d = 80.4 pc, $G$ = 15.1 mag, $K$=11.2 mag, R$_\star$ = 0.358 $\pm$ 0.015 R$_\odot$, M$_\star$ = 0.340 $\pm$ 0.009 M$_\odot$). Using the photometric time series from TESS sectors 10, 36, 46, and 63 and near-infrared spectrophotometry from ExTrA, we measured a planetary radius of 0.77 $\pm$ 0.03… ▽ More We report the discovery and characterisation of a giant transiting planet orbiting a nearby M3.5V dwarf (d = 80.4 pc, $G$ = 15.1 mag, $K$=11.2 mag, R$_\star$ = 0.358 $\pm$ 0.015 R$_\odot$, M$_\star$ = 0.340 $\pm$ 0.009 M$_\odot$). Using the photometric time series from TESS sectors 10, 36, 46, and 63 and near-infrared spectrophotometry from ExTrA, we measured a planetary radius of 0.77 $\pm$ 0.03 R$_J$ and an orbital period of 1.52 days. With high-resolution spectroscopy taken by the CFHT/SPIRou and ESO/ESPRESSO spectrographs, we refined the host star parameters ([Fe/H] = 0.27 $\pm$ 0.12) and measured the mass of the planet (0.273 $\pm$ 0.006 M$_J$). Based on these measurements, TOI-4860 b joins the small set of massive planets ($>$80 M$_E$) found around mid to late M dwarfs ($<$0.4 R$_\odot$), providing both an interesting challenge to planet formation theory and a favourable target for further atmospheric studies with transmission spectroscopy. We identified an additional signal in the radial velocity data that we attribute to an eccentric planet candidate ($e=0.66\pm0.09$) with an orbital period of $427\pm7$~days and a minimum mass of $1.66\pm 0.26$ M$_J$, but additional data would be needed to confirm this. △ Less

Submitted 12 January, 2024; v1 submitted 2 August, 2023; originally announced August 2023.

Comments: 16 pages, 14 figures, accepted for publication in A&A

A high-resolution radio study of the L1551 IRS 5 and L1551 NE jets

Authors: A. Feeney-Johansson, S. J. D. Purser, T. P. Ray, C. Carrasco-González, A. Rodríguez-Kamenetzky, J. Eislöffel, J. Lim, R. Galván-Madrid, S. Lizano, L. F. Rodríguez, H. Shang, P. Ho, M. Hoare

Abstract: Using observations with e-MERLIN and the VLA, together with archival data from ALMA, we obtain high-resolution radio images of two binary YSOs: L1551 IRS 5 and L1551 NE, covering a wide range of frequencies from 5 - 336 GHz, and resolving emission from the radio jet on scales of only ~15 au. By comparing these observations to those from a previous epoch, it is shown that there is a high degree of… ▽ More Using observations with e-MERLIN and the VLA, together with archival data from ALMA, we obtain high-resolution radio images of two binary YSOs: L1551 IRS 5 and L1551 NE, covering a wide range of frequencies from 5 - 336 GHz, and resolving emission from the radio jet on scales of only ~15 au. By comparing these observations to those from a previous epoch, it is shown that there is a high degree of variability in the free-free emission from the jets of these sources. In particular, the northern component of L1551 IRS 5 shows a remarkable decline in flux density of a factor of ~5, suggesting that the free-free emission of this source has almost disappeared. By fitting the spectra of the sources, the ionised mass-loss rates of the jets are derived and it is shown that there is significant variability of up to a factor of ~6 on timescales of ~20 years. Using radiative transfer modelling, we also obtained a model image for the jet of the southern component of L1551 IRS 5 to help study the inner region of the ionised high-density jet. The findings favour the X-wind model launched from a very small innermost region. △ Less

Submitted 24 July, 2023; v1 submitted 24 July, 2023; originally announced July 2023.

Comments: 13 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 677, A97 (2023)

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): A forming quadruple system with continuum `ribbons' and intricate outflows

Authors: Qiu-yi Luo, Tie Liu, Aaron T. Lee, Stella S. R. Offner, James di Francesco, Doug Johnstone, Mika Juvela, Paul F. Goldsmith, Sheng-Li Qin, Xiaofeng Mai, Xun-chuan Liu, Patricio Sanhueza, Feng-Wei Xu, Ken'ichi Tatematsu, Somnath Dutta, Huei-Ru Vivien Chen, Shanghuo Li, Aiyuan Yang, Sheng-Yuan Liu, Chin-Fei Lee, Naomi Hirano, Chang Won Lee, Dipen Sahu, Hsien Shang, Shih-Ying Hsu , et al. (9 additional authors not shown)

Abstract: One of the most poorly understood aspects of low-mass star formation is how multiple-star systems are formed. Here we present the results of Atacama Large Millimeter/submillimeter Array (ALMA) Band-6 observations towards a forming quadruple protostellar system, G206.93-16.61E2, in the Orion B molecular cloud. ALMA 1.3 mm continuum emission reveals four compact objects, of which two are Class I you… ▽ More One of the most poorly understood aspects of low-mass star formation is how multiple-star systems are formed. Here we present the results of Atacama Large Millimeter/submillimeter Array (ALMA) Band-6 observations towards a forming quadruple protostellar system, G206.93-16.61E2, in the Orion B molecular cloud. ALMA 1.3 mm continuum emission reveals four compact objects, of which two are Class I young stellar objects (YSOs), and the other two are likely in prestellar phase. The 1.3 mm continuum emission also shows three asymmetric ribbon-like structures that are connected to the four objects, with lengths ranging from $\sim$500 au to $\sim$2200 au. By comparing our data with magneto-hydrodynamic (MHD) simulations, we suggest that these ribbons trace accretion flows and also function as gas bridges connecting the member protostars. Additionally, ALMA CO J=2-1 line emission reveals a complicated molecular outflow associated with G206.93-16.61E2 with arc-like structures suggestive of an outflow cavity viewed pole-on. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Comments: The paper was accepted by APJL

Density Structure of Centrally Concentrated Prestellar Cores from Multi-scale Observations

Authors: Dipen Sahu, Sheng-Yuan Liu, Doug Johnstone, Tie Liu, Neal J. Evans II, Naomi Hirano, Kenichi Tatematsu, James Di Francesco, Chin-Fei Lee, Kee-Tae Kim, Somnath Dutta, Shih-Ying Hsu, Shanghuo Li, Qiu-Yi Luo, Patricio Sanhueza, Hsien Shang, Alessio Traficante, Mika Juvela, Chang Won Lee, David J. Eden, Paul F. Goldsmith, Leonardo Bronfman, Woojin Kwon, Jeong-Eun Lee, Yi-Jehng Kuan , et al. (1 additional authors not shown)

Abstract: Starless cores represent the initial stage of evolution toward (proto)star formation, and a subset of them, known as prestellar cores, with high density (~ 10^6 cm^-3 or higher) and being centrally concentrated are expected to be embryos of (proto)stars. Determining the density profile of prestellar cores, therefore provides an important opportunity to gauge the initial conditions of star formatio… ▽ More Starless cores represent the initial stage of evolution toward (proto)star formation, and a subset of them, known as prestellar cores, with high density (~ 10^6 cm^-3 or higher) and being centrally concentrated are expected to be embryos of (proto)stars. Determining the density profile of prestellar cores, therefore provides an important opportunity to gauge the initial conditions of star formation. In this work, we perform rigorous modeling to estimate the density profiles of three nearly spherical prestellar cores among a sample of five highly dense cores detected by our recent observations. We employed multi-scale observational data of the (sub)millimeter dust continuum emission including those obtained by SCUBA-2 on the JCMT with a resolution of ~5600 au and by multiple ALMA observations with a resolution as high as ~480 au. We are able to consistently reproduce the observed multi-scale dust continuum images of the cores with a simple prescribed density profile, which bears an inner region of flat density and a r^-2 profile toward the outer region. By utilizing the peak density and the size of the inner flat region as a proxy for the dynamical stage of the cores, we find that the three modeled cores are most likely unstable and prone to collapse. The sizes of the inner flat regions, as compact as ~500 au, signify them being the highly evolved prestellar cores rarely found to date. △ Less

Submitted 14 February, 2023; originally announced February 2023.

Comments: Accepted for publication in ApJ

A Unified Model for Bipolar Outflows from Young Stars: Apparent Magnetic Jet Acceleration

Authors: Hsien Shang, Ruben Krasnopolsky, Chun-Fan Liu

Abstract: We explore a new, efficient mechanism that can power toroidally magnetized jets up to two to three times their original terminal velocity after they enter a self-similar phase of magnetic acceleration. Underneath the elongated outflow lobe formed by a magnetized bubble, a wide-angle free wind, through the interplay with its ambient toroid, is compressed and accelerated around its axial jet. The ex… ▽ More We explore a new, efficient mechanism that can power toroidally magnetized jets up to two to three times their original terminal velocity after they enter a self-similar phase of magnetic acceleration. Underneath the elongated outflow lobe formed by a magnetized bubble, a wide-angle free wind, through the interplay with its ambient toroid, is compressed and accelerated around its axial jet. The extremely magnetic bubble can inflate over its original size, depending on the initial Alfvén Mach number $M_A$ of the launched flow. The shape-independent slope $\partial{}v_r/\partial{}r=2/3t$ is a salient feature of the self-similarity in the acceleration phase. Peculiar kinematic signatures are observable in the position--velocity (PV) diagrams and can combine with other morphological signatures as probes for the density-collimated jets arising in toroidally dominated magnetized winds. The apparent second acceleration is powered by the decrease of the toroidal magnetic field but operates far beyond the scales of the primary magnetocentrifugal launch region and the free asymptotic terminal state. Rich implications may connect the jets arising from the youngest protostellar outflows such as HH 211 and HH 212 and similar systems with parsec-scale jets across the mass and evolutionary spectra. △ Less

Submitted 20 January, 2023; originally announced January 2023.

Comments: 21 pages, 8 figures, 1 table, to appear in Astrophysical Journal Letters (2023)

Journal ref: The Astrophysical Journal Letters, 945:L1 (21pp), 2023 March 1

A Unified Model for Bipolar Outflows from Young Stars: Kinematic Signatures of Jets, Winds, and Their Magnetic Interplay with the Ambient Toroids

Authors: Hsien Shang, Chun-Fan Liu, Ruben Krasnopolsky, Liang-Yao Wang

Abstract: Kinematic signatures of the jet, winds, multicavities, and episodic shells arising in the unified model of bipolar outflows developed in Shang et al.\ (2020), in which an outflow forms by radially directed, wide-angle toroidally magnetized winds interacting with magnetized isothermal toroids, are extracted in the form of position--velocity diagrams. Elongated outflow lobes, driven by magnetized wi… ▽ More Kinematic signatures of the jet, winds, multicavities, and episodic shells arising in the unified model of bipolar outflows developed in Shang et al.\ (2020), in which an outflow forms by radially directed, wide-angle toroidally magnetized winds interacting with magnetized isothermal toroids, are extracted in the form of position--velocity diagrams. Elongated outflow lobes, driven by magnetized winds and their interplay with the environment, are dominated by extended bubble structures with mixing layers beyond the conventional thin-shell models. The axial cylindrically stratified density jet carries a broad profile near the base, across the projected velocity of the wide-angle wind, and narrows down along the axis with the collimated flow. The reverse shock encloses the magnetized free wind, forms an innermost cavity, and deflects the flow pattern. Shear, Kelvin--Helmholtz instabilities, and pseudopulses add fine and distinctive features between the jet--shell components, and the fluctuating jet velocities. The broad webbed velocity features connect the extremely high and the low velocities across the multicavities, mimicking nested outflowing slower-wind components. Rings and ovals in the perpendicular cuts trace multicavities at different heights, and the compressed ambient gap regions enrich the low-velocity features with protruding spikes. Our kinematic signatures capture the observed systematics of the high-, intermediate-, and low-velocity components from Class 0 to II jet--outflow systems in molecular and atomic lines. The nested shells observed in HH 212, HH 30, and DG Tau B are naturally explained. Outflows as bubbles are ubiquitous and form an inevitable integrative outcome of the interaction between wind and ambient media. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: 66 pages, 26 figures, 3 tables, to appear in the Astrophysical Journal (2023)

Journal ref: The Astrophysical Journal, 944:230 (55pp), 2023 February 20

A sub-Neptune planet around TOI-1695 discovered and characterized with SPIRou and TESS

Authors: F. Kiefer, G. Hébrard, E. Martioli, E. Artigau, R. Doyon, J. -F. Donati, C. Cadieux, A. Carmona, D. R. Ciardi, P. I. Cristofari, L. de Almeida, P. Figueira, E. Gaidos, E. Gonzales, A. Lecavelier, K. G. Stassun, L. Arnold, B. Benneke, I. Boisse, X. Bonfils, N. J. Cook, P. Cortés-Zuleta, X. Delfosse, J. Dias do Nascimento, M. Fausnaugh , et al. (17 additional authors not shown)

Abstract: TOI-1695 is a V-mag=13 M-dwarf star from the northern hemisphere at 45$\,$pc from the Sun, around which a 3.134-day periodic transit signal from a super-Earth candidate was identified in TESS photometry. With a transit depth of 1.3$\,$mmag, the radius of candidate TOI-1695.01 was estimated by the TESS pipeline to be 1.82$\,$R$_\oplus$ with an equilibrium temperature of $\sim 620\,$K. We successful… ▽ More TOI-1695 is a V-mag=13 M-dwarf star from the northern hemisphere at 45$\,$pc from the Sun, around which a 3.134-day periodic transit signal from a super-Earth candidate was identified in TESS photometry. With a transit depth of 1.3$\,$mmag, the radius of candidate TOI-1695.01 was estimated by the TESS pipeline to be 1.82$\,$R$_\oplus$ with an equilibrium temperature of $\sim 620\,$K. We successfully detect a reflex motion of the star and establish it is due to a planetary companion at an orbital period consistent with the photometric transit period thanks to a year-long radial-velocity monitoring of TOI-1695 by the SPIRou infrared spectropolarimeter. We use and compare different methods to reduce and analyse those data. We report a 5.5-$σ$ detection of the planetary signal, giving a mass of $5.5 \pm 1.0\,$M$_\oplus$ and a radius of $2.03 \pm 0.18\,$R$_\oplus$. We derive a mean equilibrium planet temperature of $590 \pm 90\,$K. The mean density of this small planet of $3.6 \pm 1.1\,$g$\,$cm$^{-3}$ is similar (1.7-$σ$ lower) than that of the Earth. It leads to a non-negligible fraction of volatiles in its atmosphere with $f_{H,He}=0.28^{+0.46}_{-0.23}$% or $f_\text{water}=23 \pm 12$%. TOI-1695$\,$b is a new sub-Neptune planet at the border of the M-dwarf radius valley that can help test formation scenarios for super-Earth/sub-Neptune-like planets. △ Less

Submitted 14 November, 2022; v1 submitted 11 November, 2022; originally announced November 2022.

Comments: 15 pages (+9 appendix pages), 19 figures, 8 tables, accepted for publication in A&A

Journal ref: A&A 670, A136 (2023)

Time-Variable Jet Ejections from RW Aur A, RY Tau and DG Tau

Authors: Michihiro Takami, Hans Moritz Guenther, P. Christian Schneider, Tracy L. Beck, Jennifer L. Karr, Youichi Ohyama, Roberto Galvan-Madrid, Taichi Uyama, Marc White, Konstantin Grankin, Deirdre Coffey, Chun-Fan Liu, Misato Fukagawa, Nadine Manset, Wen-Ping Chen, Tae-Soo Pyo, Hsien Shang, Thomas P. Ray, Masaaki Otsuka, Mei-Yin Chou

Abstract: We present Gemini-NIFS, VLT-SINFONI and Keck-OSIRIS observations of near-infrared [Fe II] emission associated with the well-studied jets from three active T Tauri stars; RW Aur A, RY Tau and DG Tau taken from 2012-2021. We primarily covered the redshifted jet from RW Aur A, and the blueshifted jets from RY Tau and DG Tau, to investigate long-term time variabilities potentially related to the activ… ▽ More We present Gemini-NIFS, VLT-SINFONI and Keck-OSIRIS observations of near-infrared [Fe II] emission associated with the well-studied jets from three active T Tauri stars; RW Aur A, RY Tau and DG Tau taken from 2012-2021. We primarily covered the redshifted jet from RW Aur A, and the blueshifted jets from RY Tau and DG Tau, to investigate long-term time variabilities potentially related to the activities of mass accretion and/or the stellar magnetic fields. All of these jets consist of several moving knots with tangential velocities of 70-240 km s-1, ejected from the star with different velocities and at irregular time intervals. Via comparison with literature, we identify significant differences in tangential velocities for the DG Tau jet between 1985-2008 and 2008-2021. The sizes of the individual knots appear to increase with time, and in turn, their peak brightnesses in the 1.644-micron emission decreased up to a factor of ~30 during the epochs of our observations. A variety of the decay timescales measured in the [Fe II] 1.644 micron emission can be attributed to different pre-shock conditions if the moving knots are unresolved shocks. However, our data do not exclude the possibility that these knots are due to non-uniform density/temperature distributions with another heating mechanism, or in some cases due to stationary shocks without proper motions. Spatially resolved observations of these knots with significantly higher angular resolutions are necessary to better understand their physical nature. △ Less

Submitted 19 October, 2022; originally announced October 2022.

Comments: 24 pages, 11 figures, 8 tables, accepted for ApJS

Two-dimensional Particle-in-Cell simulations of axisymmetric black hole magnetospheres: angular dependence of the Blandford-Znajek flux

Authors: Kouichi Hirotani, Hsien Shang, Ruben Krasnopolsky, Ken-ichi Nishikawa

Abstract: We examine the temporary evolution of axisymmetric magnetospheres around rapidly rotating black holes (BHs), by applying our two-dimensional particle-in-cell simulation code. Assuming a stellar-mass BH, we find that the created pairs fail to screen the electric field along the magnetic field, provided that the mass accretion rate is much small compared to the Eddington limit. Magnetic islands are… ▽ More We examine the temporary evolution of axisymmetric magnetospheres around rapidly rotating black holes (BHs), by applying our two-dimensional particle-in-cell simulation code. Assuming a stellar-mass BH, we find that the created pairs fail to screen the electric field along the magnetic field, provided that the mass accretion rate is much small compared to the Eddington limit. Magnetic islands are created by reconnection near the equator and migrate toward the event horizon, expelling magnetic flux tubes from the BH vicinity during a large fraction of time. When the magnetic islands stick to the horizon due to redshift and virtually vanish, a strong magnetic field penetrates the horizon, enabling efficient extraction of energy from the BH. During this flaring phase, a BH gap appears around the inner light surface with a strong meridional return current toward the equator within the ergosphere. If the mass accretion rate is 0.025 percent of the Eddington limit, the BH's spin-down luminosity becomes 16-19 times greater than its analytical estimate during the flares, although its long-term average is only 6 percent of it. We demonstrate that the extracted energy flux concentrates along the magnetic field lines threading the horizon in the middle latitudes. It is implied that this meridional concentration of the Poynting flux may result in the formation of limb-brightened jets from low-accreting BH systems. △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: Accepted for publication in ApJ; 25 pages, 17 figures, 1 movie (fig. 8, 18MB, filename=Aph_Npe_515.gif)

Journal ref: The Astrophysical Journal, 943:164 (19pp), 2023 February 1

A Q-band line survey towards Orion KL using the Tianma radio telescope

Authors: Xunchuan Liu, Tie Liu, Zhiqiang Shen, Sheng-Li Qin, Qiuyi Luo, Yu Cheng, Qilao Gu, Tianwei Zhang, Fengyao Zhu, Sheng-Yuan Liu, Xing Lu, Rongbing Zhao, Weiye Zhong, Yajun Wu, Juan Li, Zhang Zhao, Jinqing Wang, Qinghui Liu, Bo Xia, Bin Li, Li Fu, Zhen Yan, Chao Zhang, Lingling Wang, Qian Ye , et al. (7 additional authors not shown)

Abstract: We have conducted a line survey towards Orion KL using the Q-band receiver of Tianma 65 m radio telescope (TMRT), covering 34.8--50 GHz with a velocity resolution between 0.79 km s$^{-1}$ and 0.55 km s$^{-1}$ respectively. The observations reach a sensitivity on the level of 1-8 mK, proving that the TMRT is sensitive for conducting deep line surveys. In total, 597 Gaussian features are extracted.… ▽ More We have conducted a line survey towards Orion KL using the Q-band receiver of Tianma 65 m radio telescope (TMRT), covering 34.8--50 GHz with a velocity resolution between 0.79 km s$^{-1}$ and 0.55 km s$^{-1}$ respectively. The observations reach a sensitivity on the level of 1-8 mK, proving that the TMRT is sensitive for conducting deep line surveys. In total, 597 Gaussian features are extracted. Among them, 177 radio recombination lines (RRLs) are identified, including 126, 40 and 11 RRLs of hydrogen, helium and carbon, with a maximum $Δn$ of 16, 7, and 3, respectively. The carbon RRLs are confirmed to originate from photodissociation regions with a $V_{\rm LSR}\sim$9 km s$^{-1}$. In addition, 371 molecular transitions of 53 molecular species are identified. Twenty-one molecular species of this survey were not firmly detected in the Q band by Rizzo et al. (2017), including species such as H$_2$CS, HCOOH, C$_2$H$_5$OH, H$_2^{13}$CO, H$_2$CCO, CH$_3$CHO, CH$_2$OCH$_2$, HCN $v_2=1$, and CH$_3$OCHO $v_t=1$. In particular, the vibrationally excited states of ethyl cyanide (C$_2$H$_5$CN $v$13/$v$21) are for the first time firmly detected in the Q band. NH$_3$ (15,15) and (16,16) are identified, and they are so far the highest transitions of the NH$_3$ inversion lines detected towards Orion KL. All the identified lines can be reproduced by a radiative transfer model. △ Less

Submitted 7 September, 2022; originally announced September 2022.

Comments: 51 pages, 18 figures, accepted by ApJS

The JCMT Transient Survey: Single Epoch Transients and Variability of Faint Sources

Authors: Doug Johnstone, Bhavana Lalchand, Steve Mairs, Hsien Shang, Wen Ping Chen, Geoffrey C. Bower, Gregory J. Herczeg, Jeong-Eun Lee, Jan Forbrich, Bo-Yan Chen, Carlos Contreras-Pena, Yong-Hee Lee, Wooseok Park, Colton Broughton, Spencer Plovie, The JCMT Transient Team

Abstract: Short-duration flares at millimeter wavelengths provide unique insights into the strongest magnetic reconnection events in stellar coronae, and combine with longer-term variability to introduce complications to next-generation cosmology surveys. We analyze 5.5 years of JCMT Transient Survey 850 micron submillimeter monitoring observations toward eight Gould Belt star-forming regions to search for… ▽ More Short-duration flares at millimeter wavelengths provide unique insights into the strongest magnetic reconnection events in stellar coronae, and combine with longer-term variability to introduce complications to next-generation cosmology surveys. We analyze 5.5 years of JCMT Transient Survey 850 micron submillimeter monitoring observations toward eight Gould Belt star-forming regions to search for evidence of transient events or long-duration variability from faint sources. The eight regions (30 arcmin diameter fields), including ~1200 infrared-selected YSOs, have been observed on average 47 times with integrations of approximately half an hour, or one day total spread over 5.5 years. Within this large data set, only two robust faint source detections are recovered: JW 566 in OMC 2/3 and MGM12 2864 in NGC 2023. JW 566, a Class II TTauri binary system previously identified as an extraordinary submillimeter flare, remains unique, the only clear single-epoch transient detection in this sample with a flare eight times bright than our ~4.5 sigma detection threshold of 55 mJy/beam. The lack of additional recovered flares intermediate between JW 566 and our detection limit is puzzling, if smaller events are more common than larger events. In contrast, the other submillimeter variable identified in our analysis, Source 2864, is highly variable on all observed timescales. Although Source 2864 is occasionally classified as a YSO, the source is most likely a blazar. The degree of variability across the electromagnetic spectrum may be used to aid source classification. △ Less

Submitted 16 August, 2022; originally announced August 2022.

Comments: Accepted for publication in the Astrophysicall Journal. 17 pages, 6 figures, 2 tables

Formation of episodic jets and associated flares from black hole accretion systems

Authors: Miljenko Cemeljic, Hai Yang, Feng Yuan, Hsien Shang

Abstract: Episodic ejections of blobs (episodic jets) are widely observed in black hole sources and usually associated with flares. In this paper, by performing and analyzing three dimensional general relativity magnetohydrodynamical numerical simulations of accretion flows, we investigate their physical mechanisms. We find that magnetic reconnection occurs in the accretion flow, likely due to the turbulent… ▽ More Episodic ejections of blobs (episodic jets) are widely observed in black hole sources and usually associated with flares. In this paper, by performing and analyzing three dimensional general relativity magnetohydrodynamical numerical simulations of accretion flows, we investigate their physical mechanisms. We find that magnetic reconnection occurs in the accretion flow, likely due to the turbulent motion and differential rotation of the accretion flow, resulting in flares and formation of flux ropes. Flux ropes formed inside of 10-15 gravitational radii are found to mainly stay within the accretion flow, while flux ropes formed beyond this radius are ejected outward by magnetic forces and form the episodic jets. These results confirm the basic scenario proposed in Yuan et al.(2009). Moreover, our simulations find that the predicted velocity of the ejected blobs is in good consistency with observations of Sgr A*, M81, and M87. The whole processes are found to occur quasi-periodically, with the period being the orbital time at the radius where the flux rope is formed. The predicted period of flares and ejections is consistent with those found from the light curves or image of Sgr A*, M87, and PKS 1510-089. The possible applications to protostellar accretion systems are discussed. △ Less

Submitted 21 May, 2022; originally announced May 2022.

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

Resolving the collimation zone of an intermediate-mass protostar

Authors: Adriana R. Rodríguez-Kamenetzky, Carlos Carrasco-González, Luis Felipe Rodríguez Jorge, Tom P. Ray, Alberto Sanna, Luca Moscadelli, Melvin Hoare, Roberto Galván-Madrid, Hsien Shang, Susana Lizano, Jochen Eislöffel, Jeremy Lim, José M. Torrelles, Paul Ho, Anton Feeney-Johansson

Abstract: We report new VLA and e-MERLIN high resolution and sensitivity images of the Triple Radio continuum Source in the Serpens star forming region. These observations allowed us to perform a deep multi-frequency, multi-epoch study by exploring the innermost regions (<~100 au) of an intermediate-mass YSO for the first time, with a physical resolution of ~15 au. The kinematic analysis of knots recently e… ▽ More We report new VLA and e-MERLIN high resolution and sensitivity images of the Triple Radio continuum Source in the Serpens star forming region. These observations allowed us to perform a deep multi-frequency, multi-epoch study by exploring the innermost regions (<~100 au) of an intermediate-mass YSO for the first time, with a physical resolution of ~15 au. The kinematic analysis of knots recently ejected by the protostar indicates that the jet is undergoing episodic variations in velocity. In addition, our multi-frequency images reveal striking characteristics, e.g., a highly collimated ionized stream that would be launched at a radial distance of ~0.4 au from the protostar, and a narrow (~28 au wide) ionized cavity that would be excited by the interaction of a wide-angle component with the surrounding toroid of infalling material. In light of these results, we propose the scenario in which both a highly-collimated jet and a wide-angle wind coexist to be the most plausible to explain our observations, either launched by the X-wind or X- plus Disk-wind mechanism. △ Less

Submitted 20 May, 2022; originally announced May 2022.

Comments: 13 pages, 3 figures. Accepted for publication in APJL

MSC Class: 85-02 ACM Class: J.2

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): Deriving Inclination Angle and Velocity of the Protostellar Jets from their SiO Knots

Authors: Kai-Syun Jhan, Chin-Fei Lee, Doug Johnstone, Tie Liu, Sheng-Yuan Liu, Naomi Hirano, Kenichi Tatematsu, Somnath Dutta, Anthony Moraghan, Hsien Shang, Jeong-Eun Lee, Shanghuo Li, Chun-Fan Liu, Shih-Ying Hsu, Woojin Kwon, Dipen Sahu, Xun-Chuan Liu, Kee-Tae Kim, Qiuyi Luo, Sheng-Li Qin, Patricio Sanhueza, Leonardo Bronfman, Zhang Qizhou, David Eden, Alessio Traficante , et al. (1 additional authors not shown)

Abstract: We have selected six sources (G209.55-19.68S2, G205.46-14.56S1$_{-}$A, G203.21-11.20W2, G191.90-11.21S, G205.46-14.56S3, and G206.93-16.61W2) from the Atacama Large Millimeter/submillimeter Array Survey of Orion Planck Galactic Cold Clumps (ALMASOP), in which these sources have been mapped in the CO (J=2-1), SiO (J=5-4), and C$^{18}$O (J=2-1) lines. These sources have high-velocity SiO jets surrou… ▽ More We have selected six sources (G209.55-19.68S2, G205.46-14.56S1$_{-}$A, G203.21-11.20W2, G191.90-11.21S, G205.46-14.56S3, and G206.93-16.61W2) from the Atacama Large Millimeter/submillimeter Array Survey of Orion Planck Galactic Cold Clumps (ALMASOP), in which these sources have been mapped in the CO (J=2-1), SiO (J=5-4), and C$^{18}$O (J=2-1) lines. These sources have high-velocity SiO jets surrounded by low-velocity CO outflows. The SiO jets consist of a chain of knots. These knots have been thought to be produced by semi-periodical variations in jet velocity. Therefore, we adopt a shock-forming model, which uses such variations to estimate the inclination angle and velocity of the jets. We also derive the inclination angle of the CO outflows using the wide-angle wind-driven shell model, and find it to be broadly consistent with that of the associated SiO jets. In addition, we apply this shock-forming model to another three protostellar sources with SiO jets in the literature -- HH 211, HH 212, and L1448C(N) -- and find that their inclination angle and jet velocity are consistent with those previously estimated from proper motion and radial velocity studies. △ Less

Submitted 27 April, 2022; originally announced April 2022.

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): Evidence for a Molecular Jet Launched at an Unprecedented Early Phase of Protostellar evolution

Authors: Somnath Dutta, Chin-Fei Lee, Naomi Hirano, Tie Liu, Doug Johnstone, Sheng-Yuan Liu, Kenichi Tatematsu, Paul F. Goldsmith, Dipen Sahu, Neal J. Evans, Patricio Sanhueza, Woojin Kwon, Sheng-Li Qin, Manash Ranjan Samal, Qizhou Zhang, Kee-Tae Kim, Hsien Shang, Chang Won Lee, Anthony Moraghan, Kai-Syun Jhan, Shanghuo Li, Jeong-Eun Lee, Alessio Traficante, Mika Juvela, Leonardo Bronfman , et al. (9 additional authors not shown)

Abstract: Protostellar outflows and jets play a vital role in star formation as they carry away excess angular momentum from the inner disk surface, allowing the material to be transferred toward the central protostar. Theoretically, low velocity and poorly collimated outflows appear from the beginning of the collapse, at the first hydrostatic core (FHSC) stage. With growing protostellar core mass, high-den… ▽ More Protostellar outflows and jets play a vital role in star formation as they carry away excess angular momentum from the inner disk surface, allowing the material to be transferred toward the central protostar. Theoretically, low velocity and poorly collimated outflows appear from the beginning of the collapse, at the first hydrostatic core (FHSC) stage. With growing protostellar core mass, high-density jets are launched which entrain an outflow from the infalling envelope. Until now, molecular jets have been observed at high velocity ($\gtrsim$ 100 km/s) in early Class\,0 protostars. We, for the first time, detect a dense molecular jet in SiO emission with small-velocity ($\sim$ 4.2 km\,s$^{-1}$, deprojected $\sim$ 24 km\,s$^{-1}$) from source G208.89-20.04Walma (hereafter, G208Walma) using ALMA Band\,6 observations. This object has some characteristics of FHSCs, such as a small outflow/jet velocity, extended 1.3\,mm continuum emission, and N$_2$D$^+$ line emission. Additional characteristics, however, are typical of early protostars: collimated outflow and SiO jet. The full extent of the outflow corresponds to a dynamical time scale of $\sim$ 930$^{+200}_{-100}$ years. The spectral energy distribution also suggests a very young source having an upper limit of T$_{bol}$ $\sim$ 31 K and L$_{bol}$ $\sim$ 0.8 L$_\sun$. We conclude that G208Walma is likely in the transition phase from FHSC to protostar, and the molecular jet has been launched within a few hundred years of initial collapse. Therefore, G208Walma may be the earliest object discovered in the protostellar phase with a molecular jet. △ Less

Submitted 13 April, 2022; originally announced April 2022.

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

Detection of strong scattering close to the eclipse region of PSR B1957+20

Authors: J. T. Bai, S. Dai, Q. J. Zhi, W. A. Coles, D. Li, W. W. Zhu, G. Hobbs, G. J. Qiao, N. Wang, J. P. Yuan, M. D. Filipovic, J. B. Wang, Z. C. Pan, L. H. Shang, S. J. Dang, S. Q. Wang, C. C. Miao

Abstract: We present the first measurement of pulse scattering close to the eclipse region of PSR B1957+20, which is in a compact binary system with a low-mass star. We measured pulse scattering time-scales up to 0.2 ms close to the eclipse and showed that it scales with the dispersion measure (DM) excess roughly as $τ\proptoΔ{\rm DM}^{2}$. Our observations provide the first evidence of strong scattering du… ▽ More We present the first measurement of pulse scattering close to the eclipse region of PSR B1957+20, which is in a compact binary system with a low-mass star. We measured pulse scattering time-scales up to 0.2 ms close to the eclipse and showed that it scales with the dispersion measure (DM) excess roughly as $τ\proptoΔ{\rm DM}^{2}$. Our observations provide the first evidence of strong scattering due to multi-path propagation effects in the eclipsing material. We show that Kolmogorov turbulence in the eclipsing material with an inner scale of $\sim100$ m and an outer scale of the size of the eclipse region can naturally explain the observation. Our results show that the eclipsing material in such systems can be highly turbulent and suggest that scattering is one of the main eclipsing mechanisms at around 1.4 GHz. △ Less

Submitted 28 March, 2022; originally announced March 2022.

Comments: 8 pages, 4 figures, MNRAS accepted

Magneto-Centrifugal Origin for Protostellar Jets Validated through Detection of Radial Flow at the Jet Base

Authors: Chin-Fei Lee, Zhi-Yun Li, Hsien Shang, Naomi Hirano

Abstract: Jets can facilitate the mass accretion onto the protostars in star formation. They are believed to be launched from accretion disks around the protostars by magneto-centrifugal force, as supported by the detections of rotation and magnetic field in some of them. Here we report a radial flow of the textbook-case protostellar jet HH 212 at the base to further support this jet-launching scenario. Thi… ▽ More Jets can facilitate the mass accretion onto the protostars in star formation. They are believed to be launched from accretion disks around the protostars by magneto-centrifugal force, as supported by the detections of rotation and magnetic field in some of them. Here we report a radial flow of the textbook-case protostellar jet HH 212 at the base to further support this jet-launching scenario. This radial flow validates a central prediction of the magneto-centrifugal theory of jet formation and collimation, namely, the jet is the densest part of a wide-angle wind that flows radially outward at distances far from the (small, sub-au) launching region. Additional evidence for the radially flowing wide-angle component comes from its ability to reproduce the structure and kinematics of the shells detected around the HH 212 jet. This component, which can transport material from inner to outer disk, could account for the chondrules and Ca-Al-rich inclusions detected in the Solar System at large distances. △ Less

Submitted 2 March, 2022; originally announced March 2022.

Comments: 18 pages, 6 figures

Magnetic Spirals in Accretion Flows Originated from Misaligned Magnetic Field

Authors: Weixiao Wang, Miikka S. Väisälä, Hsien Shang, Ruben Krasnopolsky, Zhi-Yun Li, Ka Ho Lam, Feng Yuan

Abstract: Misalignment between rotation and magnetic field has been suggested to be one type of physical mechanisms which can easen the effects of magnetic braking during collapse of cloud cores leading to formation of protostellar disks. However, its essential factors are poorly understood. Therefore, we perform a more detailed analysis of the physics involved. We analyze existing simulation data to measur… ▽ More Misalignment between rotation and magnetic field has been suggested to be one type of physical mechanisms which can easen the effects of magnetic braking during collapse of cloud cores leading to formation of protostellar disks. However, its essential factors are poorly understood. Therefore, we perform a more detailed analysis of the physics involved. We analyze existing simulation data to measure the system torques, mass accretion rates and Toomre Q parameters. We also examine the presence of shocks in the system. While advective torques are generally the strongest, we find that magnetic and gravitational torques can play substantial roles in how angular momentum is transferred during the disk formation process. Magnetic torques can shape the accretion flows, creating two-armed magnetized inflow spirals aligned with the magnetic field. We find evidence of an accretion shock that is aligned according to the spiral structure of the system. Inclusion of ambipolar diffusion as explored in this work has shown a slight influence in the small scale structures but not in the main morphology. We discuss potential candidate systems where some of these phenomena could be present. △ Less

Submitted 22 January, 2022; originally announced January 2022.

Comments: 26 pages, 26 figures, 2 tables, Accepted for publication in the Astrophysical Journal

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): A Hot Corino Survey toward Protostellar Cores in the Orion Cloud

Authors: Shih-Ying Hsu, Sheng-Yuan Liu, Tie Liu, Dipen Sahu, Chin-Fei Lee, Kenichi Tatematsu, Kee-Tae Kim, Naomi Hirano, Yao-Lun Yang, Doug Johnstone, Hongli Liu, Mika Juvela, Leonardo Bronfman, Huei-Ru Vivien Chen, Somnath Dutta, David J. Eden, Kai-Syun Jhan, Yi-Jehng Kuan, Chang Won Lee, Jeong-Eun Lee, Shanghuo Li, Chun-Fan Liu, Sheng-Li Qin, Patricio Sanhueza, Hsien Shang , et al. (3 additional authors not shown)

Abstract: The presence of complex organic molecules (COMs) in the interstellar medium (ISM) is of great interest since it may link to the origin and prevalence of life in the universe. Aiming to investigate the occurrence of COMs and their possible origins, we conducted a chemical census toward a sample of protostellar cores as part of the ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP) project.… ▽ More The presence of complex organic molecules (COMs) in the interstellar medium (ISM) is of great interest since it may link to the origin and prevalence of life in the universe. Aiming to investigate the occurrence of COMs and their possible origins, we conducted a chemical census toward a sample of protostellar cores as part of the ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP) project. We report the detection of 11 hot corino sources, which exhibit compact emissions from warm and abundant COMs, among 56 Class 0/I protostellar cores. All the hot corino sources discovered are likely Class 0 and their sizes of the warm region ($>$ 100 K) are comparable to 100 au. The luminosity of the hot corino sources exhibits positive correlations with the total number of methanol and the extent of its emissions. Such correlations are consistent with the thermal desorption picture for the presence of hot corino and suggest that the lower luminosity (Class 0) sources likely have a smaller region with COMs emissions. With the same sample selection method and detection criteria being applied, the detection rates of the warm methanol in the Orion cloud (15/37) and the Perseus cloud (28/50) are statistically similar when the cloud distances and the limited sample size are considered. Observing the same set of COM transitions will bring a more informative comparison between the cloud properties. △ Less

Submitted 27 March, 2022; v1 submitted 7 January, 2022; originally announced January 2022.

Comments: 34 pages. 14 figures. 1 figure set. 1 machine-readable table. Accepted for publication in ApJ

Detection of a dense SiO jet in the evolved protostellar phase

Authors: Somnath Dutta, Chin-Fei Lee, Doug Johnstone, Tie Liu, Naomi Hirano, Sheng-Yuan Liu, Jeong-Eun Lee, Hsien Shang, Ken'ichi Tatematsu, Kee-Tae Kim, Dipen Sahu, Patricio Sanhueza, James Di Francesco, Kai-Syun Jhan, Chang Won Lee, Woojin Kwon, Shanghuo Li, Leonardo Bronfman, Hong-li Liu, Alessio Traficante, Yi-Jehng Kuan, Shih-Ying Hsu, Anthony Moraghan, Chun-Fan Liu, David Eden , et al. (3 additional authors not shown)

Abstract: Jets and outflows trace the accretion history of protostars. High-velocity molecular jets have been observed from several protostars in the early Class\,0 phase of star formation, detected with the high-density tracer SiO. Until now, no clear jet has been detected with SiO emission from isolated evolved Class\,I protostellar systems. We report a prominent dense SiO jet from a Class\,I source G205S… ▽ More Jets and outflows trace the accretion history of protostars. High-velocity molecular jets have been observed from several protostars in the early Class\,0 phase of star formation, detected with the high-density tracer SiO. Until now, no clear jet has been detected with SiO emission from isolated evolved Class\,I protostellar systems. We report a prominent dense SiO jet from a Class\,I source G205S3 (HOPS\,315: T$_{bol}$ $\sim$ 180 K, spectral index $\sim$ 0.417), with a moderately high mass-loss rate ($\sim$ 0.59 $\times$ 10$^{-6}$ M$_\odot$ yr$^{-1}$) estimated from CO emission. Together, these features suggest that G205S3 is still in a high accretion phase, similar to that expected of Class\,0 objects. We compare G205S3 to a representative Class\,0 system G206W2 (HOPS\,399) and literature Class\,0/I sources to explore the possible explanations behind the SiO emission seen at the later phase. We estimate a high inclination angle ($\sim$ 40$^\circ$) for G205S3 from CO emission, which may expose the infrared emission from the central core and mislead the spectral classification. However, the compact 1.3\,mm continuum, C$^{18}$O emission, location in the bolometric luminosity to sub-millimeter fluxes diagram, outflow force ($\sim$ 3.26 $\times$ 10$^{-5}$ M$_\odot$km s$^{-1}$/yr) are also analogous to that of Class\,I systems. We thus consider G205S3 to be at the very early phase of Class\,I, and in the late phase of ``high-accretion". The episodic ejection could be due to the presence of an unknown binary, a planetary companion, or dense clumps, where the required mass for such high accretion could be supplied by a massive circumbinary disk. △ Less

Submitted 26 October, 2021; originally announced October 2021.

Comments: 13 pages, 5 figures, Accepted for publication in the Astrophysical Journal

The JCMT Transient Survey: Four Year Summary of Monitoring the Submillimeter Variability of Protostars

Authors: Yong-Hee Lee, Doug Johnstone, Jeong-Eun Lee, Gregory Herczeg, Steve Mairs, Carlos Contreras-Peña, Jennifer Hatchell, Tim Naylor, Graham S. Bell, Tyler L. Bourke, Colton Broughton, Logan Francis, Aashish Gupta, Daniel Harsono, Sheng-Yuan Liu, Geumsook Park, Spencer Plovie, Gerald H. Moriarty-Schieven, Aleks Scholz, Tanvi Sharma, Paula Stella Teixeira, Yao-Te Wang, Yuri Aikawa, Geoffrey C. Bower, Huei-Ru Vivien Chen , et al. (27 additional authors not shown)

Abstract: We present the four-year survey results of monthly submillimeter monitoring of eight nearby ($< 500 $pc) star-forming regions by the JCMT Transient Survey. We apply the Lomb-Scargle Periodogram technique to search for and characterize variability on 295 submillimeter peaks brighter than 0.14 Jy beam$^{-1}$, including 22 disk sources (Class II), 83 protostars (Class 0/I), and 190 starless sources.… ▽ More We present the four-year survey results of monthly submillimeter monitoring of eight nearby ($< 500 $pc) star-forming regions by the JCMT Transient Survey. We apply the Lomb-Scargle Periodogram technique to search for and characterize variability on 295 submillimeter peaks brighter than 0.14 Jy beam$^{-1}$, including 22 disk sources (Class II), 83 protostars (Class 0/I), and 190 starless sources. We uncover 18 secular variables, all of them protostars. No single-epoch burst or drop events and no inherently stochastic sources are observed. We classify the secular variables by their timescales into three groups: Periodic, Curved, and Linear. For the Curved and Periodic cases, the detectable fractional amplitude, with respect to mean peak brightness, is $\sim4$ % for sources brighter than $\sim$ 0.5 Jy beam$^{-1}$. Limiting our sample to only these bright sources, the observed variable fraction is 37 % (16 out of 43). Considering source evolution, we find a similar fraction of bright variables for both Class 0 and Class I. Using an empirically motivated conversion from submillimeter variability to variation in mass accretion rate, six sources (7 % of our full sample) are predicted to have years-long accretion events during which the excess mass accreted reaches more than 40 % above the total quiescently accreted mass: two previously known eruptive Class I sources, V1647 Ori and EC 53 (V371 Ser), and four Class 0 sources, HOPS 356, HOPS 373, HOPS 383, and West 40. Considering the full protostellar ensemble, the importance of episodic accretion on few years timescale is negligible, only a few percent of the assembled mass. However, given that this accretion is dominated by events of order the observing time-window, it remains uncertain as to whether the importance of episodic events will continue to rise with decades-long monitoring. △ Less

Submitted 22 July, 2021; originally announced July 2021.

Comments: Accepted for publication in the Astrophysical Journal

Characterization of dense Planck clumps observed with Herschel and SCUBA-2

Authors: E. Mannfors, M. Juvela, L. Bronfman, D. J. Eden, Jinhua He, Gwanjeon Kim, Kee-Tae Kim, H. Kirppu, T. Liu, J. Montillaud, H. Parsons, Patricio Sanhueza, Hsien Shang, A. Soam, K. Tatematsu, A. Traficante, M. S. Väisälä, Chang Won Lee

Abstract: We aim to characterize a diverse selection of dense, potentially star-forming cores, clumps, and clouds within the Milky Way in terms of their dust emission and SF activity. We studied 53 fields that have been observed in the JCMT SCUBA-2 continuum survey SCOPE and have been mapped with Herschel. We estimated dust properties by fitting Herschel observations with modified blackbody functions, studi… ▽ More We aim to characterize a diverse selection of dense, potentially star-forming cores, clumps, and clouds within the Milky Way in terms of their dust emission and SF activity. We studied 53 fields that have been observed in the JCMT SCUBA-2 continuum survey SCOPE and have been mapped with Herschel. We estimated dust properties by fitting Herschel observations with modified blackbody functions, studied the relationship between dust temperature and dust opacity spectral index $β$, and estimated column densities. We extracted clumps from the SCUBA-2 850 $μ$m maps with the FellWalker algorithm and examined their masses and sizes. Clumps are associated with young stellar objects found in several catalogs. We estimated the gravitational stability of the clumps with virial analysis. The clumps are categorized as unbound starless, prestellar, or protostellar. We find 529 dense clumps, typically with high column densities from (0.3-4.8)$\times 10^{22}$ cm$^{-2}$, with a mean of (1.5$\pm$0.04)$\times10^{22}$ cm$^{-2}$, low temperatures ($T\sim $10-20 K), and estimated submillimeter $β$ =1.7$\pm$0.1. We detect a slight increase in opacity spectral index toward millimeter wavelengths. Masses of the sources range from 0.04 $M_\odot$ to 4259 $M_\odot$. Mass, linear size, and temperature are correlated with distance. Furthermore, the estimated gravitational stability is dependent on distance, and more distant clumps appear more virially bound. Finally, we present a catalog of properties of the clumps.Our sources present a large array of SF regions, from high-latitude, nearby diffuse clouds to large SF complexes near the Galactic center. Analysis of these regions will continue with the addition of molecular line data, which will allow us to study the densest regions of the clumps in more detail. △ Less

Submitted 18 June, 2021; originally announced June 2021.

Comments: 60 pages (article 18 pages), 73 figures To be published in A&A

Journal ref: A&A 654, A123 (2021)

Revealing Ionization Conditions of Sz 102 with Spatially Resolved [Ne III] Microjets

Authors: Chun-Fan Liu, Hsien Shang, Gregory J. Herczeg, Frederick M. Walter

Abstract: Forbidden neon emission lines from small-scale microjets can probe high-energy processes in low-mass young stellar systems. We obtained spatially resolved [Ne III] spectra of the microjets from the classical T Tauri Star Sz 102 using the Hubble Space Telescope Imaging Spectrograph (HST/STIS) at a spatial resolution of ~0".1. The blueshifted and redshifted [Ne III] emission both peak in intensity w… ▽ More Forbidden neon emission lines from small-scale microjets can probe high-energy processes in low-mass young stellar systems. We obtained spatially resolved [Ne III] spectra of the microjets from the classical T Tauri Star Sz 102 using the Hubble Space Telescope Imaging Spectrograph (HST/STIS) at a spatial resolution of ~0".1. The blueshifted and redshifted [Ne III] emission both peak in intensity within ~0".1 of the star and gradually decay along the flow outward to ~0".24. The spatial distribution and extent of the [Ne III] microjet is consistent with a jet that is ionized close to the base and subsequently recombines on a longer timescale than the flow time. Ca II H and K lines are also detected from the redshifted microjet with a line full-width at half-maximum of ~170 km/s, consistent with those of other forbidden emission lines, atop a 300-km/s wide stellar component. The launching radius of the Sz 102 jet, inferred from the observed line centroids and the range of inclination angles and stellar masses from the literature, is on the order of ~0.03 au. The possible proximity of the launching region to the star allows immediate ionization without distance dilution from the circumstellar ionization sources, most likely keV X-ray flares generated by magnetic reconnection events in the star-disk system, to sustain the observed [Ne III] flux. △ Less

Submitted 30 December, 2020; originally announced December 2020.

Comments: 16 pages, 10 figures, 2 tables; accepted by ApJ

Interaction of large- and small-scale dynamos in isotropic turbulent flows from GPU-accelerated simulations

Authors: Miikka S. Väisälä, Johannes Pekkilä, Maarit J. Käpylä, Matthias Rheinhardt, Hsien Shang, Ruben Krasnopolsky

Abstract: Magnetohydrodynamical (MHD) dynamos emerge in many different astrophysical situations where turbulence is present, but the interaction between large-scale (LSD) and small-scale dynamos (SSD) is not fully understood. We performed a systematic study of turbulent dynamos driven by isotropic forcing in isothermal MHD with magnetic Prandtl number of unity, focusing on the exponential growth stage. Both… ▽ More Magnetohydrodynamical (MHD) dynamos emerge in many different astrophysical situations where turbulence is present, but the interaction between large-scale (LSD) and small-scale dynamos (SSD) is not fully understood. We performed a systematic study of turbulent dynamos driven by isotropic forcing in isothermal MHD with magnetic Prandtl number of unity, focusing on the exponential growth stage. Both helical and non-helical forcing was employed to separate the effects of LSD and SSD in a periodic domain. Reynolds numbers (Rm) up to $\approx 250$ were examined and multiple resolutions used for convergence checks. We ran our simulations with the Astaroth code, designed to accelerate 3D stencil computations on graphics processing units (GPUs) and to employ multiple GPUs with peer-to-peer communication. We observed a speedup of $\approx 35$ in single-node performance compared to the widely used multi-CPU MHD solver Pencil Code. We estimated the growth rates both from the averaged magnetic fields and their power spectra. At low Rm, LSD growth dominates, but at high Rm SSD appears to dominate in both helically and non-helically forced cases. Pure SSD growth rates follow a logarithmic scaling as a function of Rm. Probability density functions of the magnetic field from the growth stage exhibit SSD behaviour in helically forced cases even at intermediate Rm. We estimated mean-field turbulence transport coefficients using closures like the second-order correlation approximation (SOCA). They yield growth rates similar to the directly measured ones and provide evidence of $α$ quenching. Our results are consistent with the SSD inhibiting the growth of the LSD at moderate Rm, while the dynamo growth is enhanced at higher Rm. △ Less

Submitted 16 December, 2020; originally announced December 2020.

Comments: 22 pages, 23 figures, 2 tables, Accepted for publication in the Astrophysical Journal

Report number: NORDITA-2020-067

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): Detection of extremely high density compact structure of prestellar cores and multiple substructures within

Authors: Dipen Sahu, Sheng-Yuan Liu, Tie Liu, Neal J. Evans II, Naomi Hirano, Ken'ichi Tatematsu, Chin-Fei Lee, Kee-Tae Kim, Somnath Dutta, Dana Alina, Leonardo Bronfman, Maria Cunningham, David J. Eden, Guido Garay, Paul F. Goldsmith, Jinhua He, Shih-Ying Hsu, Kai-Syun Jhan, Doug Johnstone, Mika Juvela, Gwanjeong Kim, Yi-Jehng Kuan, Woojin Kwon, Chang Won Lee, Jeong-Eun Lee , et al. (16 additional authors not shown)

Abstract: Prestellar cores are self-gravitating dense and cold structures within molecular clouds where future stars are born. They are expected, at the stage of transitioning to the protostellar phase, to harbor centrally concentrated dense (sub)structures that will seed the formation of a new star or the binary/multiple stellar systems. Characterizing this critical stage of evolution is key to our underst… ▽ More Prestellar cores are self-gravitating dense and cold structures within molecular clouds where future stars are born. They are expected, at the stage of transitioning to the protostellar phase, to harbor centrally concentrated dense (sub)structures that will seed the formation of a new star or the binary/multiple stellar systems. Characterizing this critical stage of evolution is key to our understanding of star formation. In this work, we report the detection of high density (sub)structures on the thousand-au scale in a sample of dense prestellar cores. Through our recent ALMA observations towards the Orion molecular cloud, we have found five extremely dense prestellar cores, which have centrally concentrated regions $\sim$ 2000 au in size, and several $10^7$ $cm^{-3}$ in average density. Masses of these centrally dense regions are in the range of 0.30 to 6.89 M$_\odot$. {\it For the first time}, our higher resolution observations (0.8$'' \sim $ 320 au) further reveal that one of the cores shows clear signatures of fragmentation; such individual substructures/fragments have sizes of 800 -1700 au, masses of 0.08 to 0.84 M$_\odot$, densities of $2 - 8\times 10^7$ $cm^{-3}$ and separations of $\sim 1200$ au. The substructures are massive enough ($\gtrsim 0.1~M_\odot$) to form young stellar objects and are likely examples of the earliest stage of stellar embryos which can lead to widely ($\sim$ 1200 au) separated multiple systems. △ Less

Submitted 12 February, 2023; v1 submitted 15 December, 2020; originally announced December 2020.

Comments: Published in The Astrophysical Journal Letters (ApJL)

Two-dimensional Particle-in-Cell simulations of axisymmetric black hole magnetospheres

Authors: Kouichi Hirotani, Ruben Krasnopolsky, Hsien Shang, Ken-ichi Nishikawa, Michael Watson

Abstract: We investigate the temporal evolution of an axisymmetric magnetosphere around a rapidly rotating, stellar-mass black hole, applying a two-dimensional particle-in-cell simulation scheme. Adopting a homogeneous pair production, and assuming that the mass accretion rate is much less than the Eddington limit, we find that the black hole's rotational energy is preferentially extracted from the middle l… ▽ More We investigate the temporal evolution of an axisymmetric magnetosphere around a rapidly rotating, stellar-mass black hole, applying a two-dimensional particle-in-cell simulation scheme. Adopting a homogeneous pair production, and assuming that the mass accretion rate is much less than the Eddington limit, we find that the black hole's rotational energy is preferentially extracted from the middle latitudes, and that this outward energy flux exhibits an enhancement that lasts approximately 160 dynamical time scales. It is demonstrated that the magnetohydrodynamic approximations cannot be justified in such a magnetically-dominated magnetosphere, because the Ohm's law completely breaks down, and because the charge-separated electron-positron plasmas are highly non-neutral. An implication is given regarding the collimation of relativistic jets. △ Less

Submitted 13 December, 2020; originally announced December 2020.

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

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP) II. Survey overview: a first look at 1.3 mm continuum maps and molecular outflows

Authors: Somnath Dutta, Chin-Fei Lee, Tie Liu, Naomi Hirano, Sheng-Yuan Liu, Kenichi Tatematsu, Kee-Tae Kim, Hsien Shang, Dipen Sahu, Gwanjeong Kim, Anthony Moraghan, Kai-Syun Jhan, Shih-Ying Hsu, Neal J. Evans, Doug Johnstone, Derek Ward-Thompson, Yi-Jehng Kuan, Chang Won Lee, Jeong-Eun Lee, Alessio Traficante, Mika Juvela, Charlotte Vastel, Qizhou Zhang, Patricio Sanhueza, Archana Soam , et al. (10 additional authors not shown)

Abstract: Planck Galactic Cold Clumps (PGCCs) are contemplated to be the ideal targets to probe the early phases of star formation. We have conducted a survey of 72 young dense cores inside PGCCs in the Orion complex with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.3\,mm (band 6) using three different configurations (resolutions $\sim$ 0$\farcs$35, 1$\farcs$0, and 7$\farcs$0) to statistical… ▽ More Planck Galactic Cold Clumps (PGCCs) are contemplated to be the ideal targets to probe the early phases of star formation. We have conducted a survey of 72 young dense cores inside PGCCs in the Orion complex with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.3\,mm (band 6) using three different configurations (resolutions $\sim$ 0$\farcs$35, 1$\farcs$0, and 7$\farcs$0) to statistically investigate their evolutionary stages and sub-structures. We have obtained images of the 1.3\,mm continuum and molecular line emission ($^{12}$CO, and SiO) at an angular resolution of $\sim$ 0$\farcs$35 ($\sim$ 140\,au) with the combined arrays. We find 70 substructures within 48 detected dense cores with median dust-mass $\sim$ 0.093\,M$_{\sun}$ and deconvolved size $\sim$ 0$\farcs$27. Dense substructures are clearly detected within the central 1000\,au of four candidate prestellar cores. The sizes and masses of the substructures in continuum emission are found to be significantly reduced with protostellar evolution from Class\,0 to Class\,I. We also study the evolutionary change in the outflow characteristics through the course of protostellar mass accretion. A total of 37 sources exhibit CO outflows, and 20 ($>$50\%) show high-velocity jets in SiO. The CO velocity-extents ($Δ$Vs) span from 4 to 110 km/s with outflow cavity opening angle width at 400\,au ranging from $[Θ_{obs}]_{400}$ $\sim$ 0$\farcs$6 to 3$\farcs$9, which corresponds to 33$\fdg$4$-$125$\fdg$7. For the majority of the outflow sources, the $Δ$Vs show a positive correlation with $[Θ_{obs}]_{400}$, suggesting that as protostars undergo gravitational collapse, the cavity opening of a protostellar outflow widens and the protostars possibly generate more energetic outflows. △ Less

Submitted 27 October, 2020; originally announced October 2020.

Comments: 43 pages, 25 figures (14 in the main text, 9 in the appendix), accepted for publication in the Astrophysical Journal Supplement Series (ApJS)

The Interplay between Ambipolar Diffusion and Hall Effect on Magnetic Field Decoupling and Protostellar Disc Formation

Authors: Bo Zhao, Paola Caselli, Zhi-Yun Li, Ruben Krasnopolsky, Hsien Shang, Ka Ho Lam

Abstract: Non-ideal MHD effects have been shown recently as a robust mechanism of averting the magnetic braking "catastrophe" and promoting protostellar disc formation. However, the magnetic diffusivities that determine the efficiency of non-ideal MHD effects are highly sensitive to microphysics. We carry out non-ideal MHD simulations to explore the role of microphysics on disc formation and the interplay b… ▽ More Non-ideal MHD effects have been shown recently as a robust mechanism of averting the magnetic braking "catastrophe" and promoting protostellar disc formation. However, the magnetic diffusivities that determine the efficiency of non-ideal MHD effects are highly sensitive to microphysics. We carry out non-ideal MHD simulations to explore the role of microphysics on disc formation and the interplay between ambipolar diffusion (AD) and Hall effect during the protostellar collapse. We find that removing the smallest grain population ($\lesssim$10 nm) from the standard MRN size distribution is sufficient for enabling disc formation. Further varying the grain sizes can result in either a Hall-dominated or an AD-dominated collapse; both form discs of tens of AU in size regardless of the magnetic field polarity. The direction of disc rotation is bimodal in the Hall dominated collapse but unimodal in the AD-dominated collapse. We also find that AD and Hall effect can operate either with or against each other in both radial and azimuthal directions, yet the combined effect of AD and Hall is to move the magnetic field radially outward relative to the infalling envelope matter. In addition, microphysics and magnetic field polarity can leave profound imprints both on observables (e.g., outflow morphology, disc to stellar mass ratio) and on the magnetic field characteristics of protoplanetary discs. Including Hall effect relaxes the requirements on microphysics for disc formation, so that prestellar cores with cosmic-ray ionization rate of $\lesssim$2--3$\times10^{-16}$ s$^{-1}$ can still form small discs of $\lesssim$10 AU radius. We conclude that disc formation should be relatively common for typical prestellar core conditions, and that microphysics in the protostellar envelope is essential to not only disc formation, but also protoplanetary disc evolution. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Comments: 21 pages, 16 figures, submitted to MNRAS

Hall Effect in Protostellar Disc Formation and Evolution

Authors: Bo Zhao, Paola Caselli, Zhi-Yun Li, Ruben Krasnopolsky, Hsien Shang, Ka Ho Lam

Abstract: The Hall effect is recently shown to be efficient in magnetized dense molecular cores, and could lead to a bimodal formation of rotationally supported discs (RSDs) in the first core phase. However, how such Hall dominated systems evolve in the protostellar accretion phase remains unclear. We carry out 2D axisymmetric simulations including Hall effect and Ohmic dissipation, with realistic magnetic… ▽ More The Hall effect is recently shown to be efficient in magnetized dense molecular cores, and could lead to a bimodal formation of rotationally supported discs (RSDs) in the first core phase. However, how such Hall dominated systems evolve in the protostellar accretion phase remains unclear. We carry out 2D axisymmetric simulations including Hall effect and Ohmic dissipation, with realistic magnetic diffusivities computed from our equilibrium chemical network. We find that Hall effect only becomes efficient when the large population of very small grains (VSGs: $\lesssim$10 nm) is removed from the standard MRN size distribution. With such an enhanced Hall effect, however, the bimodality of disc formation does not continue into the main accretion phase. The outer part of the initial $\sim$40 AU disc formed in the anti-aligned configuration (${\bf Ω\cdot B}<0$) flattens into a thin rotationally supported Hall current sheet as Hall effect moves the poloidal magnetic field radially inward relative to matter, leaving only the inner $\lesssim$10--20 AU RSD. In the aligned configuration (${\bf Ω\cdot B}>0$), disc formation is suppressed initially but a counter-rotating disc forms subsequently due to efficient azimuthal Hall drift. The counter-rotating disc first grows to $\sim$30 AU as Hall effect moves the magnetic field radially outward, but only the inner $\lesssim$10 AU RSD is long-lived like in the anti-aligned case. Besides removing VSGs, cosmic ray ionization rate should be below a few 10$^{-16}$ s$^{-1}$ for Hall effect to be efficient in disc formation. We conclude that Hall effect produces small $\lesssim$10--20 AU discs regardless of the polarity of the magnetic field, and that radially outward diffusion of magnetic fields remains crucial for disc formation and growth. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Comments: 22 pages, 18 figures

Journal ref: MNRAS 492, 3375Z (2020)

Possible Time Correlation Between Jet Ejection and Mass Accretion for RW Aur A

Authors: Michihiro Takami, Tracy L. Beck, P. Christian Schneider, Hans Moritz Guenther, Marc White, Konstantin Grankin, Jennifer L. Karr, Youichi Ohyama, Deirdre Coffey, Hauyu Baobab Liu, Roberto Galvan-Madrid, Chun-Fan Liu, Misato Fukagawa, Nadine Manset, Wen-Ping Chen, Tae-Soo Pyo, Hsien Shang, Thomas P. Ray, Masaaki Otsuka, Mei-Yin Chou

Abstract: For the active T-Taur star RW Aur A we have performed long-term (~10 yr) monitoring observations of (1) jet imaging in the [Fe II] 1.644-micron emission line using Gemini-NIFS and VLT-SINFONI; (2) optical high-resolution spectroscopy using CFHT-ESPaDOnS; and (3) V-band photometry using the CrAO 1.25-m telescope and AAVSO. The latter two observations confirm the correlation of time variabilities be… ▽ More For the active T-Taur star RW Aur A we have performed long-term (~10 yr) monitoring observations of (1) jet imaging in the [Fe II] 1.644-micron emission line using Gemini-NIFS and VLT-SINFONI; (2) optical high-resolution spectroscopy using CFHT-ESPaDOnS; and (3) V-band photometry using the CrAO 1.25-m telescope and AAVSO. The latter two observations confirm the correlation of time variabilities between (A) the Ca II 8542 A and O I 7772 A line profiles associated with magnetospheric accretion, and (B) optical continuum fluxes. The jet images and their proper motions show that four knot ejections occurred at the star over the past ~15 years with an irregular interval of 2-6 years. The time scale and irregularity of these intervals are similar to those of the dimming events seen in the optical photometry data. Our observations show a possible link between remarkable (Delta_V < -1 mag.) photometric rises and jet knot ejections. Observations over another few years may confirm or reject this trend. If confirmed, this would imply that the location of the jet launching region is very close to the star (r <<0.1 au) as predicted by some jet launching models. Such a conclusion would be crucial for understanding disk evolution within a few au of the star, and therefore possible ongoing planet formation at these radii. △ Less

Submitted 31 July, 2020; originally announced July 2020.

Comments: 14 pages, 5 figures, Accepted by Astrophysical Journal

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): I. Detection of New Hot Corinos with ACA

Authors: Shih-Ying Hsu, Sheng-Yuan Liu, Tie Liu, Dipen Sahu, Naomi Hirano, Chin-Fei Lee, Kenichi Tatematsu, Gwanjeong Kim, Mika Juvela, Patricio Sanhueza, Jinhua He, Doug Johnstone, Sheng-Li Qin, Leonardo J. Bronfman, Huei-Ru Chen, Somnath Dutta, David Eden, Kai-Syun Jhan, Kee-Tae Kim, Yi-Jehng Kuan, Woojin Kwon, Chang Won Lee, Jeong-Eun Lee, Anthony Moraghan, Mark Rawlings , et al. (7 additional authors not shown)

Abstract: We report the detection of four new hot corino sources, G211.47-19.27S, G208.68-19.20N1, G210.49-19.79W and G192.12-11.10 from a survey study of Planck Galactic Cold Clumps in the Orion Molecular Cloud Complex with the Atacama Compact Array (ACA). Three sources had been identified as low mass Class 0 protostars in the Herschel Orion Protostar Survey (HOPS). One source in the lambda Orionis region… ▽ More We report the detection of four new hot corino sources, G211.47-19.27S, G208.68-19.20N1, G210.49-19.79W and G192.12-11.10 from a survey study of Planck Galactic Cold Clumps in the Orion Molecular Cloud Complex with the Atacama Compact Array (ACA). Three sources had been identified as low mass Class 0 protostars in the Herschel Orion Protostar Survey (HOPS). One source in the lambda Orionis region is firstly reported as a protostellar core. We have observed abundant complex organic molecules (COMs), primarily methanol but also other oxygen-bearing COMs (in G211.47-19.27S and G208.68-19.20N1) and the molecule of prebiotic interest NH2CHO (in G211.47-19.27S), signifying the presence of hot corinos. While our spatial resolution is not sufficient for resolving most of the molecular emission structure, the large linewidth and high rotational temperature of COMs suggest that they likely reside in the hotter and innermost region immediately surrounding the protostar. In G211.47-19.27S, the D/H ratio of methanol ([CH2DOH]/[CH3OH]) and the 12C/13C ratio of methanol ([CH3OH]/[13CH3OH]) are comparable to those of other hot corinos. Hydrocarbons and long carbon-chain molecules such as c-C3H2 and HCCCN are also detected in the four sources, likely tracing the outer and cooler molecular envelopes. △ Less

Submitted 29 June, 2020; originally announced June 2020.

Comments: 37 pages, 51 figures, to be published in ApJ

Transition from ordered pinched to warped magnetic field on a 100 au scale in the Class 0 protostar B335

Authors: Hsi-Wei Yen, Bo Zhao, Patrick Koch, Ruben Krasnopolsky, Zhi-Yun Li, Nagayoshi Ohashi, Hsien Shang, Shigehisa Takakuwa, Ya-Wen Tang

Abstract: We present our observational results of the 0.87 mm polarized dust emission in the Class 0 protostar B335 obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) at a 0.2" (20 au) resolution. We compared our data at 0.87 mm with those at 1.3 mm from the ALMA archive. The observed polarization orientations at the two wavelengths are consistent within the uncertainty, and the polarizat… ▽ More We present our observational results of the 0.87 mm polarized dust emission in the Class 0 protostar B335 obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) at a 0.2" (20 au) resolution. We compared our data at 0.87 mm with those at 1.3 mm from the ALMA archive. The observed polarization orientations at the two wavelengths are consistent within the uncertainty, and the polarization percentages are systematically higher at 1.3 mm than 0.87 mm by a factor of ~1.7, suggesting that the polarized emission originates from magnetically aligned dust grains. We inferred the magnetic field orientations from the observed polarization orientations. We found that the magnetic field changes from ordered and highly pinched to more complicated and asymmetric structures within the inner 100 au scale of B335, and the magnetic field connects to the center along the equatorial plane as well as along the directions which are ~40-60 degrees from the equatorial plane. We performed non-ideal MHD simulations of collapsing dense cores. We found that similar magnetic field structures appear in our simulations of dense cores with the magnetic field and rotational axis slightly misaligned by 15 degrees but not in those with the aligned magnetic field and rotational axis. Our results suggest that the midplane of the inner envelope within the inner 100 au scale of B335 could be warped because of the misaligned magnetic field and rotational axis, and the magnetic field could be dragged by the warped accretion flows. △ Less

Submitted 9 March, 2020; originally announced March 2020.

Comments: 11 pages, 6 figures, Accepted by ApJ

Direct Calculation of Self-Gravitational Force for Infinitesimally Thin Gaseous Disks Using Adaptive Mesh Refinement

Authors: Yao-Huan Tseng, Hsien Shang, Chien-Chang Yen

Abstract: Yen et al. (2012) advanced a direct approach for the calculation of self-gravitational force to second order accuracy based on uniform grid discretization. This method improves the accuracy of N-body calculation by using exact integration of kernel functions and employing the Fast Fourier Transform (FFT) to reduce complexity of computation to nearly linear. This direct approach is free of artifici… ▽ More Yen et al. (2012) advanced a direct approach for the calculation of self-gravitational force to second order accuracy based on uniform grid discretization. This method improves the accuracy of N-body calculation by using exact integration of kernel functions and employing the Fast Fourier Transform (FFT) to reduce complexity of computation to nearly linear. This direct approach is free of artificial boundary conditions, however, the applicability is limited by the uniform discretization of grids. We report here an advancement in the direct method with the implementation of adaptive mesh refinement (AMR) and maintaining second-order accuracy, which breaks the barrier set by uniform grid discretization. The adoption of graphic process units (GPUs) can significantly speed up the computation and make application of this method possible to astrophysical systems of gaseous disk galaxies and protoplanetary disks. △ Less

Submitted 8 January, 2020; originally announced January 2020.

Comments: 16 pages

Journal ref: The Astrophysical Journal Supplement Series, 244:26 (13pp), 2019 October

Planet formation: The case for large efforts on the computational side

Authors: Wladimir Lyra, Thomas Haworth, Bertram Bitsch, Simon Casassus, Nicolás Cuello, Thayne Currie, Andras Gáspár, Hannah Jang-Condell, Hubert Klahr, Nathan Leigh, Giuseppe Lodato, Mordecai-Mark Mac Low, Sarah Maddison, George Mamatsashvili, Colin McNally, Andrea Isella, Sebastián Pérez, Luca Ricci, Debanjan Sengupta, Dimitris Stamatellos, Judit Szulágyi, Richard Teague, Neal Turner, Orkan Umurhan, Jacob White , et al. (32 additional authors not shown)

Abstract: Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been a… ▽ More Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: White paper submitted to the Astro2020 decadal survey

SCOPE: SCUBA-2 Continuum Observations of Pre-protostellar Evolution - Survey Description and Compact Source Catalogue

Authors: D. J. Eden, Tie Liu, Kee-Tae Kim, S. -Y. Liu, K. Tatematsu, J. Di Francesco, K. Wang, Y. Wu, M. A. Thompson, G. A. Fuller, Di Li, I. Ristorcelli, Sung-ju Kang, N. Hirano, D. Johnstone, Y. Lin, J. H. He, P. M. Koch, Patricio Sanhueza, S. -L. Qin, Q. Zhang, P. F. Goldsmith, N. J. Evans II, J. Yuan, C. -P. Zhang , et al. (136 additional authors not shown)

Abstract: We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14… ▽ More We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other. △ Less

Submitted 26 February, 2019; originally announced February 2019.

Comments: 16 pages, 13 figures, 1 table. Accepted for publication in MNRAS

Ejection History of the IRAS 04166+2706 Molecular Jet

Authors: Liang-Yao Wang, Hsien Shang, Tzu-Yang Chiang

Abstract: The high-velocity molecular jet driven by Class 0 protostar IRAS 04166+2706 exhibits a unique saw-tooth velocity pattern. It consists of a series of well-aligned symmetric knots with similar averaged speeds, whose speeds at peaks of emission decreases roughly linearly away from the origin. Recent ALMA observations of knots R6 and B6 reveal kinematic behavior with expansion velocity increasing line… ▽ More The high-velocity molecular jet driven by Class 0 protostar IRAS 04166+2706 exhibits a unique saw-tooth velocity pattern. It consists of a series of well-aligned symmetric knots with similar averaged speeds, whose speeds at peaks of emission decreases roughly linearly away from the origin. Recent ALMA observations of knots R6 and B6 reveal kinematic behavior with expansion velocity increasing linearly from the axis to the edge. This pattern can be formed by a spherically expanding wind with axial density concentration. In this picture, the diverging velocity profile naturally possesses an increasing expansion velocity away from the axis, resulting in a tooth-like feature on the position-velocity diagram through projection. Such geometric picture predicts a correspondence between the slopes of the teeth and the outflow inclination angles, and the same inclination angle of 52$^\circ$ of the IRAS 04166+2706 can generally explain the whole pattern. Aided by numerical simulations in the framework of unified wind model by Shang et al. (2006), the observed velocity pattern can indeed be generated. A proper geometrical distribution of the jet and wind material is essential to the reconstruction the ejection history of the system. △ Less

Submitted 15 February, 2019; originally announced February 2019.

Comments: 15 pages, 12 figures, accepted by ApJ

Time Evolution of 3D Disk Formation with Misaligned Magnetic Field and Rotation Axes

Authors: Miikka S. Väisälä, Hsien Shang, Ruben Krasnopolsky, Sheng-Yuan Liu, Ka Ho Lam, Zhi-Yun Li

Abstract: Distinguishing diagnostic observational signatures produced by MHD models is essential in understanding the physics for the formation of protostellar disks in the ALMA era. Developing suitable tools along with time evolution will facilitate better identification of diagnostic features. With a ray-tracing based radiative transfer code Perspective, we explore time evolution of MHD models carried out… ▽ More Distinguishing diagnostic observational signatures produced by MHD models is essential in understanding the physics for the formation of protostellar disks in the ALMA era. Developing suitable tools along with time evolution will facilitate better identification of diagnostic features. With a ray-tracing based radiative transfer code Perspective, we explore time evolution of MHD models carried out in Li, Krasnopolsky & Shang (2013) - most of which have $90^\circ$ misalignment between the rotational axis and the magnetic field. Four visible object types can be characterized, origins of which are dependent on the initial conditions. Our results show complex spiraling density, velocity and polarization structures. The systems are under constant change, but many of those distinctive features are present already early on, and they grow more visible in time, but most could not be identified from the data without examining their change in time. The results suggest that spiraling pseudodisk structures could function as an effective observation signature of the formation process, and we witness accretion in the disk with eccentric orbits which appear as spiral-like perturbation from simple circular Keplerian orbits. Magnetically aligned polarization appears purely azimuthal in the disk and magnetic field can lead to precession of the disk. △ Less

Submitted 30 December, 2018; originally announced December 2018.

Comments: 28 pages, 19 Figures and 3 Tables. Accepted for Publication in ApJ

Journal ref: ApJ, 873:114 (22pp), 2019 March 10

The formation of rings and gaps in wind-launching non-ideal MHD disks: three-dimensional simulations

Authors: Scott S. Suriano, Zhi-Yun Li, Ruben Krasnopolsky, Takeru K. Suzuki, Hsien Shang

Abstract: Previous axisymmetric investigations in two dimensions (2D) have shown that rings and gaps develop naturally in non-ideal magnetohydrodynamic (MHD) disk-wind systems, especially in the presence of ambipolar diffusion (AD). Here we extend these 2D simulations to three dimensions (3D) and find that rings and gaps are still formed in the presence of a moderately strong ambipolar diffusion. The rings… ▽ More Previous axisymmetric investigations in two dimensions (2D) have shown that rings and gaps develop naturally in non-ideal magnetohydrodynamic (MHD) disk-wind systems, especially in the presence of ambipolar diffusion (AD). Here we extend these 2D simulations to three dimensions (3D) and find that rings and gaps are still formed in the presence of a moderately strong ambipolar diffusion. The rings and gaps form from the same basic mechanism that was identified in the 2D simulations, namely, the redistribution of the poloidal magnetic flux relative to the disk material as a result of the reconnection of a sharply pinched poloidal magnetic field lines. Thus, the less dense gaps are more strongly magnetized with a large poloidal magnetic field compared to the less magnetized (dense) rings. The rings and gaps start out rather smoothly in 3D simulations that have axisymmetric initial conditions. Non-axisymmetric variations arise spontaneously at later times, but they do not grow to such an extent as to disrupt the rings and gaps. These disk substructures persist to the end of the simulations, lasting up to 3000 orbital periods at the inner edge of the simulated disk. The longevity of the perturbed yet still coherent rings make them attractive sites for trapping large grains that would otherwise be lost to rapid radial migration due to gas drag. As the ambipolar diffusivity decreases, both the disk and the wind become increasingly turbulent, driven by the magnetorotational instability, with tightly-wound spiral arms becoming more prominent in the disk. △ Less

Submitted 16 December, 2018; v1 submitted 4 October, 2018; originally announced October 2018.

Comments: 19 pages, 13 figures, accepted by MNRAS

ALMA observations of the very young Class 0 protostellar system HH 211-mms: a 30-au dusty disk with a disk-wind traced by SO?

Authors: Chin-Fei Lee, Zhi-Yun Li, Naomi Hirano, Hsien Shang, Paul T. P. Ho, Qizhou Zhang

Abstract: HH 211-mms is one of the youngest Class 0 protostellar systems in Perseus at ~ 235 pc away. We have mapped its central region at up to ~ 7 AU (0.03") resolution. A dusty disk is seen deeply embedded in a flattened envelope, with an intensity jump in dust continuum at ~ 350 GHz. It is nearly edge-on and is almost exactly perpendicular to the jet axis. It has a size of ~ 30 au along the major axis.… ▽ More HH 211-mms is one of the youngest Class 0 protostellar systems in Perseus at ~ 235 pc away. We have mapped its central region at up to ~ 7 AU (0.03") resolution. A dusty disk is seen deeply embedded in a flattened envelope, with an intensity jump in dust continuum at ~ 350 GHz. It is nearly edge-on and is almost exactly perpendicular to the jet axis. It has a size of ~ 30 au along the major axis. It is geometrically thick, indicating that the (sub)millimeter light emitting grains have yet to settle to the midplane. Its inner part is expected to have transformed into a Keplerian rotating disk with a radius of ~ 10 au. A rotating disk atmosphere and a compact rotating bipolar outflow are detected in SO. The outflow fans out from the inner disk surfaces and is rotating in the same direction as the flattened envelope, and hence could trace a disk wind carrying away angular momentum from the inner disk. From the rotation of the disk atmosphere, the protostellar mass is estimated to be <~ 50 M_Jup. Together with results from the literature, our result favors a model where the disk radius grows linearly with the protostellar mass, as predicted by models of pre-stellar dense core evolution that asymptotes to an $r^{-1}$ radial profile for both the column density and angular velocity. △ Less

Submitted 14 July, 2018; originally announced July 2018.

Comments: 20 pags, 6 figures

A 100 au-Wide Bipolar Rotating Shell Emanating From The HH 212Protostellar Disk: A Disk Wind?

Authors: Chin-Fei Lee, Zhi-Yun Li, Claudio Codella, Paul T. Ho, Linda Podio, Naomi Hirano, Hsien Shang, Neal J. Turner, Qizhou Zhang

Abstract: HH 212 is a Class 0 protostellar system found to host a "hamburger"-shaped dusty disk with a rotating disk atmosphere and a collimated SiO jet at a distance of ~ 400 pc. Recently, a compact rotating outflow has been detected in SO and SO2 toward the center along the jet axis at ~ 52 au (0.13") resolution. Here we resolve the compact outflow into a small-scale wide-opening rotating outflow shell an… ▽ More HH 212 is a Class 0 protostellar system found to host a "hamburger"-shaped dusty disk with a rotating disk atmosphere and a collimated SiO jet at a distance of ~ 400 pc. Recently, a compact rotating outflow has been detected in SO and SO2 toward the center along the jet axis at ~ 52 au (0.13") resolution. Here we resolve the compact outflow into a small-scale wide-opening rotating outflow shell and a collimated jet, with the observations in the same S-bearing molecules at ~ 16 au (0.04") resolution. The collimated jet is aligned with the SiO jet, tracing the shock interactions in the jet. The wide-opening outflow shell is seen extending out from the inner disk around the SiO jet and has a width of ~ 100 au. It is not only expanding away from the center, but also rotating around the jet axis. The specific angular momentum of the outflow shell is ~ 40 au km/s. Simple modeling of the observed kinematics suggests that the rotating outflow shell can trace either a disk wind or disk material pushed away by an unseen wind from the inner disk or protostar. We also resolve the disk atmosphere in the same S-bearing molecules, confirming the Keplerian rotation there. △ Less

Submitted 10 February, 2018; originally announced February 2018.

Comments: 29 pages, 11 figures

The formation of rings and gaps in magnetically coupled disk-wind systems: ambipolar diffusion and reconnection

Authors: Scott S. Suriano, Zhi-Yun Li, Ruben Krasnopolsky, Hsien Shang

Abstract: Radial substructures in circumstellar disks are now routinely observed by ALMA. There is also growing evidence that disk winds drive accretion in such disks. We show through 2D (axisymmetric) simulations that rings and gaps develop naturally in magnetically coupled disk-wind systems on the scale of tens of au, where ambipolar diffusion (AD) is the dominant non-ideal MHD effect. In simulations wher… ▽ More Radial substructures in circumstellar disks are now routinely observed by ALMA. There is also growing evidence that disk winds drive accretion in such disks. We show through 2D (axisymmetric) simulations that rings and gaps develop naturally in magnetically coupled disk-wind systems on the scale of tens of au, where ambipolar diffusion (AD) is the dominant non-ideal MHD effect. In simulations where the magnetic field and matter are moderately coupled, the disk remains relatively laminar with the radial electric current steepened by AD into a thin layer near the midplane. The toroidal magnetic field sharply reverses polarity in this layer, generating a large magnetic torque that drives fast accretion, which drags the poloidal field into a highly pinched radial configuration. The reconnection of this pinched field creates magnetic loops where the net poloidal magnetic flux (and thus the accretion rate) is reduced, yielding dense rings. Neighbouring regions with stronger poloidal magnetic fields accrete faster, forming gaps. In better magnetically coupled simulations, the so-called `avalanche accretion streams' develop continuously near the disk surface, rendering the disk-wind system more chaotic. Nevertheless, prominent rings and gaps are still produced, at least in part, by reconnection, which again enables the segregation of the poloidal field and the disk material similar to the more diffusive disks. However, the reconnection is now driven by the non-linear growth of MRI channel flows. The formation of rings and gaps in rapidly accreting yet laminar disks has interesting implications for dust settling and trapping, grain growth, and planet formation. △ Less

Submitted 16 March, 2018; v1 submitted 17 December, 2017; originally announced December 2017.

Comments: 21 pages, 16 figures, accepted by MNRAS

A Rotating Protostellar Jet Launched from the Innermost Disk of HH 212

Authors: Chin-Fei Lee, Paul. T. P Ho, Zhi-Yun Li, Naomi Hirano, Qizhou Zhang, Hsien Shang

Abstract: The central problem in forming a star is the angular momentum in the circumstellar disk which prevents material from falling into the central stellar core. An attractive solution to the "angular momentum problem" appears to be the ubiquitous (low-velocity and poorly-collimated) molecular outflows and (high-velocity and highly-collimated) protostellar jets accompanying the earliest phase of star fo… ▽ More The central problem in forming a star is the angular momentum in the circumstellar disk which prevents material from falling into the central stellar core. An attractive solution to the "angular momentum problem" appears to be the ubiquitous (low-velocity and poorly-collimated) molecular outflows and (high-velocity and highly-collimated) protostellar jets accompanying the earliest phase of star formation that remove angular momentum at a range of disk radii. Previous observations suggested that outflowing material carries away the excess angular momentum via magneto-centrifugally driven winds from the surfaces of circumstellar disks down to ~ 10 AU scales, allowing the material in the outer disk to transport to the inner disk. Here we show that highly collimated protostellar jets remove the residual angular momenta at the ~ 0.05 AU scale, enabling the material in the innermost region of the disk to accrete toward the central protostar. This is supported by the rotation of the jet measured down to ~ 10 AU from the protostar in the HH 212 protostellar system. The measurement implies a jet launching radius of ~ 0.05_{-0.02}^{+0.05} AU on the disk, based on the magneto-centrifugal theory of jet production, which connects the properties of the jet measured at large distances to those at its base through energy and angular momentum conservation. △ Less

Submitted 20 June, 2017; originally announced June 2017.

Comments: 23 pages, 4 figures

Journal ref: Nature Astronomy, Volume 1, id. 0152 (2017)

Formation and Atmosphere of Complex Organic Molecules of the HH 212 Protostellar Disk

Authors: Chin-Fei Lee, Zhi-Yun Li, Paul T. P. Ho, Naomi Hirano, Qizhou Zhang, Hsien Shang

Abstract: HH 212 is a nearby (400 pc) Class 0 protostellar system recently found to host a "hamburger"-shaped dusty disk with a radius of ~ 60 AU, deeply embedded in an infalling-rotating flattened envelope. We have spatially resolved this envelope-disk system with the Atacama Large Millimeter/submillimeter Array at up to ~ 16 AU (0.04") resolution. The envelope is detected in HCO+ J=4-3 down to the dusty d… ▽ More HH 212 is a nearby (400 pc) Class 0 protostellar system recently found to host a "hamburger"-shaped dusty disk with a radius of ~ 60 AU, deeply embedded in an infalling-rotating flattened envelope. We have spatially resolved this envelope-disk system with the Atacama Large Millimeter/submillimeter Array at up to ~ 16 AU (0.04") resolution. The envelope is detected in HCO+ J=4-3 down to the dusty disk. Complex organic molecules (COMs) and doubly deuterated formaldehyde (D2CO) are detected above and below the dusty disk within ~ 40 AU of the central protostar. The COMs are methanol (CH3OH), deuterated methanol (CH2DOH), methyl mercaptan (CH3SH), and formamide (NH2CHO, a prebiotic precursor). We have modeled the gas kinematics in HCO+ and COMs, and found a centrifugal barrier at a radius of ~ 44 AU, within which a Keplerian rotating disk is formed. This indicates that HCO+ traces the infalling-rotating envelope down to centrifugal barrier and COMs trace the atmosphere of a Keplerian rotating disk within the centrifugal barrier. The COMs are spatially resolved for the first time, both radially and vertically, in the atmosphere of a disk in the earliest, Class 0 phase of star formation. Our spatially resolved observations of COMs favor their formation in the disk rather than a rapidly infalling (warm) inner envelope. The abundances and spatial distributions of the COMs provide strong constraints on models of their formation and transport in low-mass star formation. △ Less

Submitted 19 June, 2017; originally announced June 2017.

Comments: 32 pages, 9 figures