Source author record

Zhi-Yun Li

Zhi-Yun Li appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.SR astro-ph.GA astro-ph.EP astro-ph astro-ph.IM astro-ph.CO astro-ph.HE physics.comp-ph

Catalog footprint

What is connected

92works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2025arXiv

Early Planet Formation in Embedded Disks (eDisk). XX. Constraining the Chemical Tracers of Young Protostellar Sources

Recent studies indicate that the formation of planets in protoplanetary disks begins early in the embedded Class 0/I phases of protostellar evolution. The physical and chemical makeup of the embedded phase can provide valuable insights into the process of star and planet formation. This study aims to provide a thorough overview of the various morphologies for molecular emissions observed on disk scales toward nearby embedded sources. We present high angular resolution (0.1", 15 au) molecular line emissions for $^{12}$CO, $^{13}$CO, C$^{18}$O, SO, SiO, DCN, CH$_3$OH, H$_2$CO, and c-C$_3$H$_2$ towards 19 nearby protostellar sources in the context of the Atacama Large Millimeter/submillimeter Array (ALMA) Large Program "Early Planet Formation in Embedded Disks (eDisk)". Emissions in $^{12}$CO are seen towards all sources and primarily trace outflowing materials. A few sources also show high-velocity jets in SiO emission and high-velocity channel maps of $^{12}$CO. The $^{13}$CO and C$^{18}$O emissions are well-known tracers of high-density regions and trace the inner envelope and disk regions with clear signs of rotation seen at continuum scales. The large-scale emissions of $^{13}$CO also delineate the outflow cavity walls where the outflowing and infalling materials interact with each other, and exposure to UV radiation leads to the formation of hydrocarbons such as c-C$_3$H$_2$. Both DCN and CH$_3$OH, when detected, show compact emissions from the inner envelope and disk regions that peak at the position of the protostar. The CH$_3$OH emissions are contained within the region of DCN emissions, which suggests that CH$_3$OH traces the hot core regions. Likewise, a few sources also display emissions in CH$_3$OH towards the outflow. Both SO and H$_2$CO show complex morphology among the sources, suggesting that they are formed through multiple processes in protostellar systems.

preprint2025arXiv

The Flying Saucer edge-on disc's Near Infrared silhouette revealed by the JWST JEDIce program

Edge-on discs offer a unique opportunity to probe radial and vertical dust and gas distributions in the protoplanetary phase. This study aims to investigate the distribution of micron-sized dust particles in the Flying Saucer (BKLT J162813-243139) in Rho Ophiuchi, leveraging the unique observational conditions of a bright infrared background that enables the edge-on disc to be seen in both silhouette and scattered light at certain, specific wavelengths. As part of the JWST Edge-on Disc Ice program ('JEDIce'), we use NIRSpec IFU observations of the Flying Saucer, serendipitously observed against a PAH-emitting background, to constrain the dust distribution and grain sizes through radiative transfer modelling. Observation of the Flying Saucer in silhouette at 3.29 microns reveals that the midplane radial extent of small dust grains is ~235 au, larger than the large-grain disc extent previously determined to be 190 au from millimetre data. The scattered light observed in emission probes micron sized icy grains at large vertical distances above the midplane. The vertical extent of the disc silhouette is similar at visible, near-IR, and mid-IR wavelengths, corroborating the conclusion that dust settling is inefficient for grains as large as tens of microns, vertically and radially.

preprint2025arXiv

The JWST-NIRCam View of Sagittarius C. III. The Extinction Curve

Determining the infrared extinction curve towards the Galactic centre is crucial for accurately correcting observed data and deriving the underlying stellar populations. However, extinction curves reported in the literature often show discrepancies. We aim to derive the infrared extinction curve towards the Galactic centre based on JWST-NIRCam data for the first time, using observations of the Sagittarius C region in the 1-5 $μ$m range. We determined extinction ratios using two different methods, both based on measuring the reddening vector using the slope of red clump stars, whose intrinsic properties are well known, in observed colour-magnitude diagrams. The extinction curve derived in this work is in good agreement with previous results in the literature. We obtained the following extinction ratios relative to F162M: $A_\mathrm{F115W} : A_\mathrm{F162M} : A_\mathrm{F182M} : A_\mathrm{F212N} : A_\mathrm{F360M} : A_\mathrm{F405N} : A_\mathrm{F470N} : A_\mathrm{F480M} = 1.84 \pm 0.03 : 1.00 : 0.789 \pm 0.005 : 0.607 \pm 0.014 : 0.306 \pm 0.011 : 0.248 \pm 0.017 : 0.240 \pm 0.019 : 0.21 \pm 0.03$. Besides, we found different values of the extinction index for the short- ($λ\sim 1-2.5\,μ$m, $α\sim 2$) and long-wavelength ($λ\sim 2.5-5\,μ$m, $α\sim 1.4$) regimes, with the extinction curve flattening at longer wavelengths. Comparison with extinction curves derived both inside and outside the Galactic centre suggests that the infrared extinction curve does not significantly vary in the central regions, and shows no significant evidence for variations between different lines of sight beyond the inner Galaxy within the uncertainties.

preprint2022arXiv

A VLA View of the Flared, Asymmetric Disk Around the Class 0 Protostar L1527 IRS

We present high resolution Karl G. Jansky Very Large Array (VLA) observations of the protostar L1527 IRS at 7 mm, 1.3 cm, and 2 cm wavelengths. We detect the edge-on dust disk at all three wavelengths and find that it is asymmetric, with the southern side of the disk brighter than the northern side. We confirm this asymmetry through analytic modeling and also find that the disk is flared at 7 mm. We test the data against models including gap features in the intensity profile, and though we cannot rule such models out, they do not provide a statistically significant improvement in the quality of fit to the data. From these fits, we can however place constraints on allowed properties of any gaps that could be present in the true, underlying intensity profile. The physical nature of the asymmetry is difficult to associate with physical features due to the edge-on nature of the disk, but could be related to spiral arms or asymmetries seen in other imaging of more face-on disks.

preprint2022arXiv

Centrifugal Barrier and Super-Keplerian Rotation in Protostellar Disk Formation

With the advent of ALMA, it is now possible to observationally constrain how disks form around deeply embedded protostars. In particular, the recent ALMA C3H2 line observations of the nearby protostar L1527 have been interpreted as evidence for the so-called "centrifugal barrier," where the protostellar envelope infall is gradually decelerated to a stop by the centrifugal force in a region of super-Keplerian rotation. To test the concept of centrifugal barrier, which was originally based on angular momentum conserving-collapse of a rotating test particle around a fixed point mass, we carry out simple axisymmetric hydrodynamic simulations of protostellar disk formation including a minimum set of ingredients: self-gravity, rotation, and a prescribed viscosity that enables the disk to accrete. We find that a super-Keplerian region can indeed exist when the viscosity is relatively large but, unlike the classic picture of centrifugal barrier, the infalling envelope material is not decelerated solely by the centrifugal force. The region has more specific angular momentum than its surrounding envelope material, which points to an origin in outward angular momentum transport in the disk (subject to the constraint of disk expansion by the infalling envelope), rather than the spin-up of the envelope material envisioned in the classic picture as it falls closer to the center in order to conserve angular momentum. For smaller viscosities, the super-Keplerian rotation is weaker or non-existing. We conclude that, despite the existence of super-Keplerian rotation in some parameter regime, the classic picture of centrifugal barrier is not supported by our simulations.

preprint2022arXiv

Disks and Outflows in the Intermediate-mass Star Forming Region NGC 2071 IR

We present ALMA band 6/7 (1.3 mm/0.87 mm) and VLA Ka band (9 mm) observations toward NGC 2071 IR, an intermediate-mass star forming region. We characterize the continuum and associated molecular line emission towards the most luminous protostars, i.e., IRS1 and IRS3, on ~100 au (0. 2") scales. IRS1 is partly resolved in millimeter and centimeter continuum, which shows a potential disk. IRS3 has a well resolved disk appearance in millimeter continuum and is further resolved into a close binary system separated by ~40 au at 9 mm. Both sources exhibit clear velocity gradients across their disk major axes in multiple spectral lines including C18O, H2CO, SO, SO2, and complex organic molecules like CH3OH, 13CH3OH and CH3OCHO. We use an analytic method to fit the Keplerian rotation of the disks, and give constraints on physical parameters with a MCMC routine. The IRS3 binary system is estimated to have a total mass of 1.4-1.5$M_\odot$. IRS1 has a central mass of 3-5$M_\odot$ based on both kinematic modeling and its spectral energy distribution, assuming that it is dominated by a single protostar. For both IRS1 and IRS3, the inferred ejection directions from different tracers, including radio jet, water maser, molecular outflow, and H2 emission, are not always consistent, and for IRS1, these can be misaligned by ~50$^{\circ}$. IRS3 is better explained by a single precessing jet. A similar mechanism may be present in IRS1 as well but an unresolved multiple system in IRS1 is also possible.

preprint2022arXiv

Formation of Dust Rings and Gaps in Non-ideal MHD Disks Through Meridional Gas Flows

Rings and gaps are commonly observed in the dust continuum emission of young stellar disks. Previous studies have shown that substructures naturally develop in the weakly ionized gas of magnetized, non-ideal MHD disks. The gas rings are expected to trap large mm/cm-sized grains through pressure gradient-induced radial dust-gas drift. Using 2D (axisymmetric) MHD simulations that include ambipolar diffusion and dust grains of three representative sizes (1~mm, 3.3~mm, and 1~cm), we show that the grains indeed tend to drift radially relative to the gas towards the centers of the gas rings, at speeds much higher than in a smooth disk because of steeper pressure gradients. However, their spatial distribution is primarily controlled by meridional gas motions, which are typically much faster than the dust-gas drift. In particular, the grains that have settled near the midplane are carried rapidly inwards by a fast accretion stream to the inner edges of the gas rings, where they are lifted up by the gas flows diverted away from the midplane by a strong poloidal magnetic field. The flow pattern in our simulation provides an attractive explanation for the meridional flows recently inferred in HD 163296 and other disks, including both "collapsing" regions where the gas near the disk surface converges towards the midplane and a disk wind. Our study highlights the prevalence of the potentially observable meridional flows associated with the gas substructure formation in non-ideal MHD disks and their crucial role in generating rings and gaps in dust.

preprint2022arXiv

Grain Growth During Protostellar Disk Formation

Recent observations indicate that mm/cm-sized grains may exist in the embedded protostellar disks. How such large grains grow from the micron size (or less) in the earliest phase of star formation remains relatively unexplored. In this study we take a first step to model the grain growth in the protostellar environment, using two-dimensional (2D axisymmetric) radiation hydrodynamic and grain growth simulations. We show that the grain growth calculations can be greatly simplified by the "terminal velocity approximation", where the dust drift velocity relative to the gas is proportional to its stopping time, which is proportional to the grain size. We find that the grain-grain collision from size-dependent terminal velocity alone is too slow to convert a significant fraction of the initially micron-sized grains into mm/cm sizes during the deeply embedded Class 0 phase. Substantial grain growth is achieved when the grain-grain collision speed is enhanced by a factor of 4. The dust growth above and below the disk midplane enables the grains to settle faster towards the midplane, which increases the local dust-to-gas ratio, which, in turn, speeds up further growth there. How this needed enhancement can be achieved is unclear, although turbulence is a strong possibility that deserves further exploration.

preprint2022arXiv

Ice Age : Chemo-dynamical modeling of Cha-MMS1 to predict new solid-phase species for detection with JWST

Chemical models and experiments indicate that interstellar dust grains and their ice mantles play an important role in the production of complex organic molecules (COMs). To date, the most complex solid-phase molecule detected with certainty in the ISM is methanol, but the James Webb Space Telescope (JWST) may be able to identify still larger organic species. In this study, we use a coupled chemo-dynamical model to predict new candidate species for JWST detection toward the young star-forming core Cha-MMS1, combining the gas-grain chemical kinetic code MAGICKAL with a 1-D radiative hydrodynamics simulation using Athena++. With this model, the relative abundances of the main ice constituents with respect to water toward the core center match well with typical observational values, providing a firm basis to explore the ice chemistry. Six oxygen-bearing COMs (ethanol, dimethyl ether, acetaldehyde, methyl formate, methoxy methanol, and acetic acid), as well as formic acid, show abundances as high as, or exceeding, 0.01% with respect to water ice. Based on the modeled ice composition, the infrared spectrum is synthesized to diagnose the detectability of the new ice species. The contribution of COMs to IR absorption bands is minor compared to the main ice constituents, and the identification of COM ice toward the core center of Cha-MMS1 with the JWST NIRCAM/Wide Field Slitless Spectroscopy (2.4-5.0 micron) may be unlikely. However, MIRI observations (5-28 micron) toward COM-rich environments where solid-phase COM abundances exceed 1% with respect to the water ice column density might reveal the distinctive ice features of COMs.

preprint2022arXiv

Magnetic Spirals in Accretion Flows Originated from Misaligned Magnetic Field

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.

preprint2022arXiv

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

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.

preprint2022arXiv

Polarization from Aligned Dust Grains in the $β$ Pic Debris Disk

We present 870 $μ$m ALMA polarization observations of thermal dust emission from the iconic, edge-on debris disk $β$ Pic. While the spatially resolved map does not exhibit detectable polarized dust emission, we detect polarization at the $\sim$3$σ$ level when averaging the emission across the entire disk. The corresponding polarization fraction is $P_\textrm{frac}$ = $0.51 \pm 0.19$%. The polarization position angle $χ$ is aligned with the minor axis of the disk, as expected from models of dust grains aligned via radiative alignment torques (RAT) with respect to a toroidal magnetic field ($B$-RAT) or with respect to the anisotropy in the radiation field ($k$-RAT). When averaging the polarized emission across the outer versus inner thirds of the disk, we find that the polarization arises primarily from the SW third. We perform synthetic observations assuming grain alignment via both $k$-RAT and $B$-RAT. Both models produce polarization fractions close to our observed value when the emission is averaged across the entire disk. When we average the models in the inner versus outer thirds of the disk, we find that $k$-RAT is the likely mechanism producing the polarized emission in $β$ Pic. A comparison of timescales relevant to grain alignment also yields the same conclusion. For dust grains with realistic aspect ratios (i.e., $s > 1.1$), our models imply low grain-alignment efficiencies.

preprint2022arXiv

Probing Magnetic Fields in Protoplanetary Disk Atmospheres through Polarized Optical/IR Light Scattered by Aligned Grains

Magnetic fields play essential roles in protoplanetary disks. Magnetic fields in the disk atmosphere are of particular interest, as they are connected to the wind-launching mechanism. In this work, we study the polarization of the light scattered off of magnetically aligned grains in the disk atmosphere, focusing on the deviation of the polarization orientation from the canonical azimuthal direction, which may be detectable in near-IR polarimetry with instruments such as VLT/SPHERE. We show with a simple disk model that the polarization can even be oriented along the radial (rather than azimuthal) direction, especially in highly inclined disks with toroidally dominated magnetic fields. This polarization reversal is caused by the anisotropy in the polarizibility of aligned grains and is thus a telltale sign of such grains. We show that the near-IR light is scattered mostly by $μ$m-sized grains or smaller at the $τ=1$ surface and such grains can be magnetically aligned if they contain superparamagnetic inclusions. For comparison with observations, we generate synthetic maps of the ratios of $U_ϕ/I$ and $Q_ϕ/I$, which can be used to infer the existence of (magnetically) aligned grains through a negative $Q_ϕ$ (polarization reversal) and/or a significant level of $U_ϕ/I$. We show that two features observed in the existing data, an asymmetric distribution of $U_ϕ$ with respect to the disk minor axis and a spatial distribution of $U_ϕ$ that is predominantly positive or negative, are incompatible with scattering by spherical grains in an axisymmetric disk. They provide indirect evidences for scattering by aligned non-spherical grains.

preprint2022arXiv

The Central 1000 au of a Pre-stellar Core Revealed with ALMA. II. Almost Complete Freeze-out

Pre-stellar cores represent the initial conditions in the process of star and planet formation. Their low temperatures ($<$10 K) allow the formation of thick icy dust mantles, which will be partially preserved in the future protoplanetary disks, ultimately affecting the chemical composition of planetary systems. Previous observations have shown that carbon- and oxygen-bearing species, in particular CO, are heavily depleted in pre-stellar cores due to the efficient molecular freeze-out onto the surface of cold dust grains. However, N-bearing species such as NH$_3$ and, in particular, its deuterated isotopologues, appear to maintain high abundances where CO molecules are mainly in solid phase. Thanks to ALMA, we present here the first clear observational evidence of NH$_2$D freeze-out toward the L1544 pre-stellar core, suggestive of the presence of a"complete-depletion zone" within a $\simeq$1800 au radius, in agreement with astrochemical pre-stellar core model predictions. Our state-of-the-art chemical model coupled with a non-LTE radiative transfer code demonstrates that NH$_2$D becomes mainly incorporated in icy mantles in the central 2000 au and starts freezing-out already at $\simeq$7000 au. Radiative transfer effects within the pre-stellar core cause the NH$_2$D(1$_{11}$-1$_{01}$) emission to appear centrally concentrated, with a flattened distribution within the central $\simeq$3000 au, unlike the 1.3 mm dust continuum emission which shows a clear peak within the central $\simeq$1800 au. This prevented NH$_2$D freeze-out to be detected in previous observations, where the central 1000 au cannot be spatially resolved.

preprint2022arXiv

The Davis-Chandrasekhar-Fermi Method Revisited

Despite the rich observational results on interstellar magnetic fields in star-forming regions, it is still unclear how dynamically significant the magnetic fields are at varying physical scales, because direct measurement of the field strength is observationally difficult. The Davis-Chandrasekhar-Fermi (DCF) method has been the most commonly used method to estimate the magnetic field strength from polarization data. It is based on the assumption that gas turbulent motion is the driving source of field distortion via linear Alfvén waves. In this work, using MHD simulations of star-forming clouds, we test the validity of the assumption underlying the DCF method by examining its accuracy in the real 3D space. Our results suggest that the DCF relation between turbulent kinetic energy and magnetic energy fluctuation should be treated as a statistical result instead of a local property. We then develop and investigate several modifications to the original DCF method using synthetic observations, and propose new recipes to improve the accuracy of DCF-derived magnetic field strength. We further note that the biggest uncertainty in the DCF analysis may come from the linewidth measurement instead of the polarization observation, especially since the line-of-sight gas velocity can be used to estimate the gas volume density, another critical parameter in the DCF method.

preprint2022arXiv

The Magnetic Field in the Milky Way Filamentary Bone G47

Star formation primarily occurs in filaments where magnetic fields are expected to be dynamically important. The largest and densest filaments trace spiral structure within galaxies. Over a dozen of these dense ($\sim$10$^4$\,cm$^{-3}$) and long ($>$10\,pc) filaments have been found within the Milky Way, and they are often referred to as "bones." Until now, none of these bones have had their magnetic field resolved and mapped in their entirety. We introduce the SOFIA legacy project FIELDMAPS which has begun mapping $\sim$10 of these Milky Way bones using the HAWC+ instrument at 214\,$μ$m and 18$\farcs$2 resolution. Here we present a first result from this survey on the $\sim$60\,pc long bone G47. Contrary to some studies of dense filaments in the Galactic plane, we find that the magnetic field is often not perpendicular to the spine (i.e., the center-line of the bone). Fields tend to be perpendicular in the densest areas of active star formation and more parallel or random in other areas. The average field is neither parallel or perpendicular to the Galactic plane nor the bone. The magnetic field strengths along the spine typically vary from $\sim$20 to $\sim$100\,$μ$G. Magnetic fields tend to be strong enough to suppress collapse along much of the bone, but for areas that are most active in star formation, the fields are notably less able to resist gravitational collapse.

preprint2022arXiv

The Twisted Magnetic Field of the Protobinary L483

We present H-band (1.65 $μ$m) and SOFIA HAWC+ 154 $μ$m polarization observations of the low-mass core L483. Our H-band observations reveal a magnetic field that is overwhelmingly in the E-W direction, which is approximately parallel to the bipolar outflow that is observed in scattered IR light and in single-dish $^{12}$CO observations. From our 154 $μ$m data, we infer a $\sim$ 45$^{\circ}$ twist in the magnetic field within the inner 5" (1000 au) of L483. We compare these new observations with published single-dish 350 $μ$m polarimetry and find that the 10,000 au scale H-band data match the smaller scale 350 $μ$m data, indicating that the collapse of L483 is magnetically regulated on these larger scales. We also present high-resolution 1.3 mm ALMA data of L483 which reveals it is a close binary star with a separation of 34 au. The plane of the binary of L483 is observed to be approximately parallel to the twisted field in the inner 1000 au. Comparing this result to the $\sim$ 1000 au protostellar envelope, we find that the envelope is roughly perpendicular to the 1000 au HAWC+ field. Using the data presented, we speculate that L483 initially formed as a wide binary and the companion star migrated to its current position, causing an extreme shift in angular momentum thereby producing the twisted magnetic field morphology observed. More observations are needed to further test this scenario.

preprint2022arXiv

Thermal Emission and Scattering by Aligned Grains: Plane-Parallel Model and Application to Multiwavelength Polarization of the HL Tau Disk

Telescopes are now able to resolve dust polarization across circumstellar disks at multiple wavelengths, allowing the study of the polarization spectrum. Most disks show clear evidence of dust scattering through their unidirectional polarization pattern typically at the shorter wavelength of $\sim 870 μ$m. However, certain disks show an elliptical pattern at $\sim 3$mm, which is likely due to aligned grains. With HL Tau, its polarization pattern at $\sim 1.3$mm shows a transition between the two patterns making it the first example to reveal such transition. We use the T-matrix method to model elongated dust grains and properly treat scattering of aligned non-spherical grains with a plane-parallel slab model. We demonstrate that a change in optical depth can naturally explain the polarization transition of HL Tau. At low optical depths, the thermal polarization dominates, while at high optical depths, dichroic extinction effectively takes out the thermal polarization and scattering polarization dominates. Motivated by results from the plane-parallel slab, we develop a simple technique to disentangle thermal polarization of the aligned grains $T_{0}$ and polarization due to scattering $S$ using the azimuthal variation of the polarization fraction. We find that, with increasing wavelength, the fractional polarization spectrum of the scattering component $S$ decreases, while the thermal component $T_{0}$ increases, which is expected since the optical depth decreases. We find several other sources similar to HL Tau that can be explained by azimuthally aligned scattering prolate grains when including optical depth effects. In addition, we explore how spirally aligned grains with scattering can appear in polarization images.

preprint2021arXiv

Accretion Flows or Outflow Cavities? Uncovering the Gas Dynamics around Lupus 3-MMS

Understanding how material accretes onto the rotationally supported disk from the surrounding envelope of gas and dust in the youngest protostellar systems is important for describing how disks are formed. Magnetohydrodynamic simulations of magnetized, turbulent disk formation usually show spiral-like streams of material (accretion flows) connecting the envelope to the disk. However, accretion flows in these early stages of protostellar formation still remain poorly characterized due to their low intensity and possibly some extended structures are disregarded as being part of the outflow cavity. We use ALMA archival data of a young Class 0 protostar, Lupus 3-MMS, to uncover four extended accretion flow-like structures in C$^{18}$O that follow the edges of the outflows. We make various types of position-velocity cuts to compare with the outflows and find the extended structures are not consistent with the outflow emission, but rather more consistent with a simple infall model. We then use a dendrogram algorithm to isolate five sub-structures in position-position-velocity space. Four out of the five sub-structures fit well ($>$95%) with our simple infall model, with specific angular momenta between $2.7-6.9\times10^{-4}\,$km$\,$s$^{-1}\,$pc and mass-infall rates of $0.5-1.1\times10^{-6}\,M_{\odot}\,$yr$^{-1}$. Better characterization of the physical structure in the supposed "outflow-cavities" is important to disentangle the true outflow cavities and accretion flows.

preprint2021arXiv

The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars V. A Characterization of Protostellar Multiplicity

We characterize protostellar multiplicity in the Orion molecular clouds using ALMA 0.87~mm and VLA 9~mm continuum surveys toward 328 protostars. These observations are sensitive to projected spatial separations as small as $\sim$20~au, and we consider source separations up to 10$^4$~au as potential companions. The overall multiplicity fraction (MF) and companion fraction (CF) for the Orion protostars are 0.30$\pm$0.03 and 0.44$\pm$0.03, respectively, considering separations from 20 to 10$^4$~au. The MFs and CFs are corrected for potential contamination by unassociated young stars using a probabilistic scheme based on the surface density of young stars around each protostar. The companion separation distribution as a whole is double peaked and inconsistent with the separation distribution of solar-type field stars, while the separation distribution of Flat Spectrum protostars is consistent solar-type field stars. The multiplicity statistics and companion separation distributions of the Perseus star-forming region are consistent with those of Orion. Based on the observed peaks in the Class 0 separations at $\sim$100~au and $\sim$10$^3$~au, we argue that multiples with separations $<$500~au are likely produced by both disk fragmentation and turbulent fragmentation with migration, and those at $\ga$10$^3$~au result primarily from turbulent fragmentation. We also find that MFs/CFs may rise from Class 0 to Flat Spectrum protostars between 100 and 10$^3$~au in regions of high YSO density. This finding may be evidence for migration of companions from $>$10$^3$~au to $<$10$^3$~au, and that some companions between 10$^3$ and 10$^4$~au must be (or become) unbound.

preprint2020arXiv

Carbon Chain Chemistry in Hot-Core Regions around Three Massive Young Stellar Objects Associated with 6.7 GHz Methanol Masers

We have carried out observations of CCH ($N=1-0$), CH$_{3}$CN ($J=5-4$), and three $^{13}$C isotopologues of HC$_{3}$N ($J=10-9$) toward three massive young stellar objects (MYSOs), G12.89+0.49, G16.86--2.16, and G28.28--0.36, with the Nobeyama 45-m radio telescope. Combined with previous results on HC$_{5}$N, the column density ratios of $N$(CCH)/$N$(HC$_{5}$N), hereafter the CCH/HC$_{5}$N ratios, in the MYSOs are derived to be $\sim 15$. This value is lower than that in a low-mass warm carbon chain chemistry (WCCC) source by more than one order of magnitude. We compare the observed CCH/HC$_{5}$N ratios with hot-core model calculations (Taniguchi et al. 2019). The observed ratios in the MYSOs can be best reproduced by models when the gas temperature is $\sim 85$ K, which is higher than in L1527, a low-mass WCCC source ($\sim 35$ K). These results suggest that carbon-chain molecules detected around the MYSOs exist at least partially in higher temperature regions than those in low-mass WCCC sources. There is no significant difference in column density among the three $^{13}$C isotopologues of HC$_{3}$N in G12.89+0.49 and G16.86-2.16, while HCC$^{13}$CN is more abundant than the others in G28.28--0.36. We discuss carbon-chain chemistry around the three MYSOs based on the CCH/HC$_{5}$N ratio and the $^{13}$C isotopic fractionation of HC$_{3}$N.

preprint2020arXiv

Detection of Irregular, Sub-mm Opaque Structures in the Orion Molecular Clouds: Protostars within 10000 years of formation?

We report ALMA and VLA continuum observations that potentially identify the four youngest protostars in the Orion Molecular Clouds taken as part of the Orion VANDAM program. These are distinguished by bright, extended, irregular emission at 0.87 mm and 8 mm and are optically thick at 0.87 mm. These structures are distinct from the disk or point-like morphologies seen toward the other Orion protostars. The 0.87 mm emission implies temperatures of 41-170 K, requiring internal heating. The bright 8 mm emission implies masses of 0.5 to 1.2 M_sun assuming standard dust opacity models. One source has a Class 0 companion, while another exhibits substructure indicating a companion-candidate. Three compact outflows are detected, two of which may be driven by companions, with dynamical times of ~300 to ~400 years. The slowest outflow may be driven by a first hydrostatic core. These protostars appear to trace an early phase when the centers of collapsing fragments become optically thick to their own radiation and compression raises the gas temperature. This phase is thought to accompany the formation of hydrostatic cores. A key question is whether these structures are evolving on free fall times of ~100 years, or whether they are evolving on Kelvin-Helmholtz times of several thousand years. The number of these sources imply a lifetime of ~6000 years, in closer agreement with the Kelvin-Helmholtz time. In this case, rotational and/or magnetic support could be slowing the collapse.

preprint2020arXiv

Determining the physical conditions of extremely young Class 0 circumbinary disk around VLA1623A

We present detailed analysis of high-resolution C18O (2-1), SO (88-77), CO (3-2) and DCO+ (3-2) data obtained by the Atacama Large Millimeter/sub-millimeter Array (ALMA) towards a Class 0 Keplerian circumbinary disk around VLA1623A, which represents one of the most complete analysis towards a Class 0 source. From the dendrogram analysis, we identified several accretion flows feeding the circumbinary disk in a highly anisotropic manner. Stream-like SO emission around the circumbinary disk reveals the complicated shocks caused by the interactions between the disk, accretion flows and outflows. A wall-like structure is discovered south of VLA1623B. The discovery of two outflow cavity walls at the same position traveling at different velocities suggests the two outflows from both VLA1623A and VLA1623B overlays on top of each other in the plane of sky. Our detailed flat and flared disk modeling shows that Cycle 2 C18O J = 2-1 data is inconsistent with the combined binary mass of 0.2 Msun as suggested by early Cycle 0 studies. The combined binary mass for VLA1623A should be modified to 0.3 ~ 0.5 Msun.

preprint2020arXiv

Hall Effect in Protostellar Disc Formation and Evolution

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.

preprint2020arXiv

Kinematic Analysis of a Protostellar Multiple System: Measuring the Protostar Masses and Assessing Gravitational Instability in the Disks of L1448 IRS3B and L1448 IRS3A

We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations towards a compact (230~au separation) triple protostar system, L1448 IRS3B, at 879~\micron with \contbeam~resolution. Spiral arm structure within the circum-multiple disk is well resolved in dust continuum toward IRS3B, and we detect the known wide (2300~au) companion, IRS3A, also resolving possible spiral substructure. Using dense gas tracers, C17O, H13CO$+$, and H13CN, we resolve the Keplerian rotation for both the circum-triple disk in IRS3B and the disk around IRS3A. Furthermore, we use the molecular line kinematic data and radiative transfer modeling of the molecular line emission to confirm that the disks are in Keplerian rotation with fitted masses of $1.19^{+0.13}_{-0.07}$ for IRS3B-ab, $1.51^{+0.06}_{-0.07}$~Msun for IRS3A, and place an upper limit on the central protostar mass for the tertiary IRS3B-c of 0.2~Msun. We measure the mass of the fragmenting disk of IRS3B to be 0.29~Msun from the dust continuum emission of the circum-multiple disk and estimate the mass of the clump surrounding IRS3B-c to be 0.07~Msun. We also find that the disk around IRS3A has a mass of 0.04~Msun. By analyzing the Toomre~Q parameter, we find the IRS3A circumstellar disk is gravitationally stable (Q$>$5), while the IRS3B disk is consistent with a gravitationally unstable disk (Q$<$1) between the radii 200-500~au. This coincides with the location of the spiral arms and the tertiary companion IRS3B-c, supporting the hypothesis that IRS3B-c was formed in situ via fragmentation of a gravitationally unstable disk.

preprint2020arXiv

Probing the Temperature Structure of Optically Thick Disks Using Polarized Emission of Aligned Grains

Polarized continuum emission from aligned grains in disks around young stellar objects can be used to probe the magnetic field, radiation anisotropy, or drift between dust and gas, depending on whether the non-spherical grains are aligned magnetically, radiatively or mechanically. We show that it can also be used to probe another key disk property -- the temperature gradient -- along sight lines that are optically thick, independent of the grain alignment mechanism. We first illustrate the technique analytically using a simple 1D slab model, which yields an approximate formula that relates the polarization fraction to the temperature gradient with respect to the optical depth tau at the tau=1 surface. The formula is then validated using models of stellar irradiated disks with and without accretion heating. The promises and challenges of the technique are illustrated with a number of Class 0 and I disks with ALMA dust polarization data, including NGC 1333 IRAS4A1, IRAS 16293B, BHB 07-11, L1527, HH 212 and HH 111. We find, in particular, that the sight lines passing through the near-side of a highly inclined disk trace different temperature gradient directions than those through the far-side, which can lead to a polarization orientation on the near-side that is orthogonal to that on the far-side, and that the HH 111 disk may be such a case. Our technique for probing the disk temperature gradient through dust polarization can complement other methods, particularly those using molecular lines.

preprint2020arXiv

Relative Alignment between Dense Molecular Cores and Ambient Magnetic Field: The Synergy of Numerical Models and Observations

The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular clouds. We find that the results on relative alignment between cores and background magnetic field in synthetic observations slightly disagree with those measured in fully 3D simulation data, which is partly because cores identified in projected 2D maps tend to coexist within filamentary structures, while 3D cores are generally more rounded. In addition, we examine the progression of magnetic field from pc- to core-scale in the simulation, which is consistent with the anisotropic core formation model that gas preferably flow along the magnetic field toward dense cores. When comparing the observed cores identified from the GBT Ammonia Survey (GAS) and Planck polarization-inferred magnetic field orientations, we find that the relative core-field alignment has a regional dependence among different clouds. More specifically, we find that dense cores in the Taurus molecular cloud tend to align perpendicular to the background magnetic field, while those in Perseus and Ophiuchus tend to have random (Perseus) or slightly parallel (Ophiuchus) orientations with respect to the field. We argue that this feature of relative core-field orientation could be used to probe the relative significance of the magnetic field within the cloud.

preprint2020arXiv

The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars IV. Unveiling the Embedded Intermediate-Mass Protostar and Disk within OMC2-FIR3/HOPS-370

We present ALMA (0.87~mm and 1.3~mm) and VLA (9~mm) observations toward the candidate intermediate-mass protostar OMC2-FIR3 (HOPS-370; L$_{bol}$~314~L$_{\odot}$) at $\sim$0.1" (40~au) resolution for the continuum emission and ~0.25" (100 au) resolution of nine molecular lines. The dust continuum observed with ALMA at 0.87~mm and 1.3~mm resolve a near edge-on disk toward HOPS-370 with an apparent radius of ~100 au. The VLA observations detect both the disk in dust continuum and free-free emission extended along the jet direction. The ALMA observations of molecular lines (H$_2$CO, SO, CH$_3$OH, $^{13}$CO, C$^{18}$O, NS, and H$^{13}$CN) reveal rotation of the apparent disk surrounding HOPS-370 orthogonal to the jet/outflow direction. We fit radiative transfer models to both the dust continuum structure of the disk and molecular line kinematics of the inner envelope and disk for the H$_2$CO, CH$_3$OH, NS, and SO lines. The central protostar mass is determined to be $\sim$2.5 M_sun with a disk radius of $\sim$94~au, when fit using combinations of the H$_2$CO, CH$_3$OH, NS, and SO lines, consistent with an intermediate-mass protostar. Modeling of the dust continuum and spectral energy distribution (SED) yields a disk mass of 0.035~M$_{\odot}$ (inferred dust+gas) and a dust disk radius of 62~au, thus the dust disk may have a smaller radius than the gas disk, similar to Class II disks. In order to explain the observed luminosity with the measured protostar mass, HOPS-370 must be accreting at a rate between 1.7 and 3.2$\times$10$^{-5}$~M$_{\odot}$~yr$^{-1}$.

preprint2020arXiv

The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars. A Statistical Characterization of Class 0 and I Protostellar Disks

We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and Flat Spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and Flat Spectrum protostars are 44.9$^{+5.8}_{-3.4}$, 37.0$^{+4.9}_{-3.0}$, and 28.5$^{+3.7}_{-2.3}$ au, respectively, and the mean protostellar dust disk masses are 25.9$^{+7.7}_{-4.0}$, 14.9$^{+3.8}_{-2.2}$, 11.6$^{+3.5}_{-1.9}$ Earth masses, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% out of 379 detected 0.87 mm continuum sources plus 42 non-detections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than that of Class II disks by a factor of $>$4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase.

preprint2020arXiv

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

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.

preprint2020arXiv

Validating Scattering-Induced (Sub)millimeter Disk Polarization through the Spectral Index, Wavelength-Dependent Polarization Pattern, and Polarization Spectrum: The Case of HD 163296

An increasing number of young circumstellar disks show strikingly ordered (sub)millimeter polarization orientations along the minor axis, which is strong evidence for polarization due to scattering by ~0.1 mm sized grains. To test this mechanism further, we model the dust continuum and polarization data of HD 163296, one of the best observed disks with prominent rings and gaps, using the RADMC-3D radiative transfer code. We find that scattering by grains with a maximum size of 90$μ$m can simultaneously reproduce the polarization observed at 0.87 mm (ALMA Band 7) and the unusually low spectral index of $α$ ~ 1.5 between 0.87 and 1.25 mm (ALMA Band 6) in the optically thick inner disk as a result of more efficient scattering at a shorter wavelength. The relatively low spectral index of $α$ ~ 2.5 inferred for the optically thin gaps is also reproduced by the same (relatively small) grains, as a result of telescope beam averaging of the gaps (with an intrinsic $α$ ~ 4) and their adjacent optically thick rings (where $α$ << 2). In this case, the long-standing tension between the grain sizes inferred from polarization and spectral index disappears because the relatively low $α$ values are illusory and do not require large mm-sized grains. In addition, the polarization fraction has a unique pattern of azimuthal variation: higher along the major axis than the minor axis in the gaps but higher along the minor axis than the major axis in the rings. We find a rapidly declining polarization spectrum (with the fraction $p \propto λ^{-3}$ approximately) in the gaps, which becomes flattened or even inverted towards short wavelengths in the optically thick rings. These contrasting behaviors in the rings and gaps provide further tests of scattering-induced polarization that can be tested via multi-wavelength observations that resolve the disk substructure.

preprint2019arXiv

Characterizing magnetic field morphologies in three Serpens protostellar cores with ALMA

With the aim of characterizing the dynamical processes involved in the formation of young protostars, we present high angular resolution ALMA dust polarization observations of the Class 0 protostellar cores Serpens SMM1, Emb 8(N), and Emb 8. With spatial resolutions ranging from 150 to 40 au at 870 μm, we find unexpectedly high values of the polarization fraction along the outflow cavity walls in Serpens Emb8(N). We use 3 mm and 1 mm molecular tracers to investigate outflow and dense gas properties and their correlation with the polarization. These observations allow us to investigate the physical processes involved in the Radiative Alignment Torques (RATs) acting on dust grains along the outflow cavity walls, which experience irradiation from accretion processes and outflow shocks. The inner core of SMM1-a presents a polarization pattern with a poloidal magnetic field at the bases of the two lobes of the bipolar outflow. To the south of SMM1-a we see two polarized filaments, one of which seems to trace the redshifted outflow cavity wall. The other may be an accretion streamer of material infalling onto the central protostar. We propose that the polarized emission we see at millimeter wavelengths along the irradiated cavity walls can be reconciled with the expectations of RAT theory if the aligned grains present at '<' 500 au scales in Class 0 envelopes have grown larger than the 0.1 μm size of ISM dust grains. Our observations allow us to constrain the star-forming sources magnetic field morphologies within the central cores, along the outflow cavity walls, and in possible accretion streamers.

preprint2019arXiv

The Infall Motion in the Low-Mass Protostellar Binary NGC 1333 IRAS 4A1/4A2

We report ALMA observations of NGC 1333 IRAS 4A, a young low-mass protostellar binary, referred as 4A1 and 4A2. With multiple H$_2$CO transitions and HNC (4$-$3) observed at a resolution of 0.25" ($\sim$70 au), we investigate the gas kinematics of 4A1 and 4A2. Our results show that, on the large angular scale ($\sim$10"), 4A1 and 4A2 each display a well-collimated outflow along the N-S direction, and an S-shaped morphology is discerned in the outflow powered by 4A2. On the small scale ($\sim$0.3"), 4A1 and 4A2 exhibit distinct spectral features toward the continuum centroid, with 4A1 showing simple symmetric profiles predominantly in absorption and 4A2 demonstrating rather complicated profiles in emission as well as in absorption. Based on radiative transfer modeling exercises, we find that the physical parameters inferred from earlier low-resolution observations cannot be directly extrapolated down to the 4A1 inner region. Possible reasons for the discrepancies between the observed and modelled profiles are discussed. We constrain the mass infall rate in 4A1 to be at most around 3$\times$10$^{-5}$ M$_{\odot}$ year$^{-1}$ at the layer of 75 au. For the kinematics of the 4A2 inner envelope, the absorbing dips in the H$_2$CO spectra are skewed toward the redshifted side and likely signatures of inward motion. These absorbing dips are relatively narrow. This is, like the case for 4A1, significantly slower than the anticipated inflow speed. We estimate a mass infall rate of 3.1$-$6.2 $\times$ 10$^{-5}$ M$_\odot$ year$^{-1}$ at the layer of 100 au in 4A2.

preprint2019arXiv

The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars I. Identifying and Characterizing the Protostellar Content of the OMC2-FIR4 and OMC2-FIR3 Regions

We present ALMA (0.87~mm) and VLA (9~mm) observations toward OMC2-FIR4 and OMC2-FIR3 within the Orion integral-shaped filament that are thought to be the nearest regions of intermediate mass star formation. We characterize the continuum sources within these regions on $\sim$40~AU (0\farcs1) scales and associated molecular line emission at a factor of $\sim$30 better resolution than previous observations at similar wavelengths. We identify six compact continuum sources within OMC2-FIR4, four in OMC2-FIR3, and one additional source just outside OMC2-FIR4. This continuum emission is tracing the inner envelope and/or disk emission on less than 100~AU scales. HOPS-108 is the only protostar in OMC2-FIR4 that exhibits emission from high-excitation transitions of complex organic molecules (e.g., methanol and other lines) coincident with the continuum emission. HOPS-370 in OMC2-FIR3 with L~$\sim$~360~\lsun, also exhibits emission from high-excitation methanol and other lines. The methanol emission toward these two protostars is indicative of temperatures high enough to thermally evaporate methanol from icy dust grains; overall these protostars have characteristics similar to hot corinos. We do not identify a clear outflow from HOPS-108 in \twco, but find evidence of interaction between the outflow/jet from HOPS-370 and the OMC2-FIR4 region. The multitude of observational constraints indicate that HOPS-108 is likely a low to intermediate-mass protostar in its main mass accretion phase and it is the most luminous protostar in OMC2-FIR4. The high resolution data presented here are essential for disentangling the embedded protostars from their surrounding dusty environments and characterizing them.

preprint2016arXiv

A Triple Protostar System Formed via Fragmentation of a Gravitationally Unstable Disk

Binary and multiple star systems are a frequent outcome of the star formation process, and as a result, almost half of all sun-like stars have at least one companion star. Theoretical studies indicate that there are two main pathways that can operate concurrently to form binary/multiple star systems: large scale fragmentation of turbulent gas cores and filaments or smaller scale fragmentation of a massive protostellar disk due to gravitational instability. Observational evidence for turbulent fragmentation on scales of $>$1000~AU has recently emerged. Previous evidence for disk fragmentation was limited to inferences based on the separations of more-evolved pre-main sequence and protostellar multiple systems. The triple protostar system L1448 IRS3B is an ideal candidate to search for evidence of disk fragmentation. L1448 IRS3B is in an early phase of the star formation process, likely less than 150,000 years in age, and all protostars in the system are separated by $<$200~AU. Here we report observations of dust and molecular gas emission that reveal a disk with spiral structure surrounding the three protostars. Two protostars near the center of the disk are separated by 61 AU, and a tertiary protostar is coincident with a spiral arm in the outer disk at a 183 AU separation. The inferred mass of the central pair of protostellar objects is $\sim$1 M$_{sun}$, while the disk surrounding the three protostars has a total mass of $\sim$0.30 M$_{\sun}$. The tertiary protostar itself has a minimum mass of $\sim$0.085 M$_{sun}$. We demonstrate that the disk around L1448 IRS3B appears susceptible to disk fragmentation at radii between 150~AU and 320~AU, overlapping with the location of the tertiary protostar. This is consistent with models for a protostellar disk that has recently undergone gravitational instability, spawning one or two companion stars.

preprint2016arXiv

An extremely high velocity molecular jet surrounded by an ionized cavity in the protostellar source Serpens SMM1

We report ALMA observations of a one-sided, high-velocity ($\sim$80 km s$^{-1}$) CO($J = 2 \rightarrow 1$) jet powered by the intermediate-mass protostellar source Serpens SMM1-a. The highly collimated molecular jet is flanked at the base by a wide-angle cavity; the walls of the cavity can be seen in both 4 cm free-free emission detected by the VLA and 1.3 mm thermal dust emission detected by ALMA. This is the first time that ionization of an outflow cavity has been directly detected via free-free emission in a very young, embedded Class 0 protostellar source that is still powering a molecular jet. The cavity walls are ionized either by UV photons escaping from the accreting protostellar source, or by the precessing molecular jet impacting the walls. These observations suggest that ionized outflow cavities may be common in Class 0 protostellar sources, shedding further light on the radiation, outflow, and jet environments in the youngest, most embedded forming stars.

preprint2016arXiv

Angular momentum loss in the envelope-disk transition region of HH 111 protostellar system: evidence for magnetic braking?

HH 111 is a Class I protostellar system at a distance of ~ 400 pc, with the central source VLA 1 associated with a rotating disk deeply embedded in a flattened envelope. Here we present the observations of this system at ~ 0.6" (240 AU) resolution in C18O (J=2-1) and 230 GHz continuum obtained with Atacama Large Millimeter/Submillimeter Array, and in SO obtained with Submillimeter Array. The observations show for the first time how a Keplerian rotating disk can be formed inside a flattened envelope. The flattened envelope is detected in C18O, extending out to >~ 2400 AU from the VLA 1 source. It has a differential rotation, with the outer part (>~ 2000 AU) better described by a rotation that has constant specific angular momentum and the innermost part (<~ 160 AU) by a Keplerian rotation. The rotationally supported disk is therefore relatively compact in this system, which is consistent with the dust continuum observations. Most interestingly, if the flow is in steady state, there is a substantial drop in specific angular momentum in the envelope-disk transition region from 2000 AU to 160 AU, by a factor of ~ 3. Such a decrease is not expected outside a disk formed from simple hydrodynamic core collapse, but can happen naturally if the core is significantly magnetized, because magnetic fields can be trapped in the transition region outside the disk by the ram pressure of the protostellar accretion flow, which can lead to efficient magnetic braking. In addition, SO shock emission is detected around the outer radius of the disk and could trace an accretion shock around the disk.

preprint2016arXiv

Axisymmetric Simulations of Hot Jupiter-Stellar Wind Hydrodynamic Interaction

Gas giant exoplanets orbiting at close distances to the parent star are subjected to large radiation and stellar wind fluxes. In this paper, hydrodynamic simulations of the planetary upper atmosphere and its interaction with the stellar wind are carried out to understand the possible flow regimes and how they affect the Lyman-alpha transmission spectrum. Following Tremblin and Chiang, charge exchange reactions are included to explore the role of energetic atoms as compared to thermal particles. In order to understand the role of the tail as compared to the leading edge of the planetary gas, the simulations were carried out under axisymmetry, and photoionization and stellar wind electron impact ionization reactions were included to limit the extent of the neutrals away from the planet. By varying the planetary gas temperature, two regimes are found. At high temperature, a supersonic planetary wind is found, which is turned around by the stellar wind and forms a tail behind the planet. At lower temperatures, the planetary wind is shut off when the stellar wind penetrates inside where the sonic point would have been. In this regime mass is lost by viscous interaction at the boundary between planetary and stellar wind gases. Absorption by cold hydrogen atoms is large near the planetary surface, and decreases away from the planet as expected. The hot hydrogen absorption is in an annulus and typically dominated by the tail, at large impact parameter, rather than by the thin leading edge of the mixing layer near the substellar point.

preprint2016arXiv

Balloon-Borne Submillimeter Polarimetry of the Vela C Molecular Cloud: Systematic Dependence of Polarization Fraction on Column Density and Local Polarization-Angle Dispersion

We present results for Vela C obtained during the 2012 flight of the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol). We mapped polarized intensity across almost the entire extent of this giant molecular cloud, in bands centered at 250, 350, and 500 μm. In this initial paper, we show our 500 μm data smoothed to a resolution of 2.5 arcminutes (approximately 0.5 pc). We show that the mean level of the fractional polarization p and most of its spatial variations can be accounted for using an empirical three-parameter power-law fit, p = p_0 N^(-0.4) S^(-0.6), where N is the hydrogen column density and S is the polarization-angle dispersion on 0.5 pc scales. The decrease of p with increasing S is expected because changes in the magnetic field direction within the cloud volume sampled by each measurement will lead to cancellation of polarization signals. The decrease of p with increasing N might be caused by the same effect, if magnetic field disorder increases for high column density sightlines. Alternatively, the intrinsic polarization efficiency of the dust grain population might be lower for material along higher density sightlines. We find no significant correlation between N and S. Comparison of observed submillimeter polarization maps with synthetic polarization maps derived from numerical simulations provides a promising method for testing star formation theories. Realistic simulations should allow for the possibility of variable intrinsic polarization efficiency. The measured levels of correlation among p, N, and S provide points of comparison between observations and simulations.

preprint2016arXiv

Disk Polarization From Both Emission and Scattering of Magnetically Aligned Grains: The Case of NGC 1333 IRAS4A1

Dust polarization in millimeter (and centimeter) has been mapped in disks around an increasing number of young stellar objects. It is usually thought to come from emission by magnetically aligned (non-spherical) grains, but can also be produced by dust scattering. We present a semi-analytic theory of disk polarization that includes both the direction emission and scattering, with an emphasis on their relative importance and how they are affected by the disk inclination. For face-on disks, both emission and scattering tend to produce polarization in the radial direction, making them difficult to distinguish, although the scattering-induced polarization can switch to the azimuthal direction if the incident radiation is beamed strongly enough in the radial direction in the disk plane. Disk inclination affects the polarizations from emission and scattering differently, especially on the major axis where, in the edge-on limit, the former vanishes while the latter reaches a polarization fraction as large as 1/3. The polarizations from the two competing mechanisms tend to cancel each other on the major axis, producing two low polarization "holes" (one on each side of the center) under certain conditions. We find tantalizing evidence for at least one such "hole" in NGC1333 IRAS4A1, whose polarization observed at 8 mm on the 100 AU scale is indicative of a pattern dominated by scattering close to the center and by direction emission in the outer region. If true, it would imply not only that a magnetic field exists on the disk scale, but that it is strong enough to align large, possibly mm-sized, grains.

preprint2016arXiv

Instrumental performance and results from testing of the BLAST-TNG receiver, submillimeter optics, and MKID arrays

Polarized thermal emission from interstellar dust grains can be used to map magnetic fields in star forming molecular clouds and the diffuse interstellar medium (ISM). The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) flew from Antarctica in 2010 and 2012 and produced degree-scale polarization maps of several nearby molecular clouds with arcminute resolution. The success of BLASTPol has motivated a next-generation instrument, BLAST-TNG, which will use more than 3000 linear polarization sensitive microwave kinetic inductance detectors (MKIDs) combined with a 2.5m diameter carbon fiber primary mirror to make diffraction-limited observations at 250, 350, and 500 $μ$m. With 16 times the mapping speed of BLASTPol, sub-arcminute resolution, and a longer flight time, BLAST-TNG will be able to examine nearby molecular clouds and the diffuse galactic dust polarization spectrum in unprecedented detail. The 250 $μ$m detector array has been integrated into the new cryogenic receiver, and is undergoing testing to establish the optical and polarization characteristics of the instrument. BLAST-TNG will demonstrate the effectiveness of kilo-pixel MKID arrays for applications in submillimeter astronomy. BLAST-TNG is scheduled to fly from Antarctica in December 2017 for 28 days and will be the first balloon-borne telescope to offer a quarter of the flight for "shared risk" observing by the community.

preprint2016arXiv

Interferometric mapping of magnetic fields: The ALMA view of the massive star forming clump W43-MM1

Here we present the first results from ALMA observations of 1 mm polarized dust emission towards the W43-MM1 high mass star forming clump. We have detected a highly fragmented filament with source masses ranging from 14Msun to 312Msun, where the largest fragment, source A, is believed to be one of the most massive in our Galaxy. We found a smooth, ordered, and detailed polarization pattern throughout the filament which we used to derived magnetic field morphologies and strengths for 12 out of the 15 fragments detected ranging from 0.2 to 9 mG. The dynamical equilibrium of each fragment was evaluated finding that all the fragments are in a super-critical state which is consistent with previously detected infalling motions towards W43-MM1. Moreover, there are indications suggesting that the field is being dragged by gravity as the whole filament is collapsing.

preprint2016arXiv

Protostellar Disk Formation Enabled by Removal of Small Dust Grains

It has been shown that a realistic level of magnetization of dense molecular cloud cores can suppress the formation of a rotationally supported disk (RSD) through catastrophic magnetic braking in the axisymmetric ideal MHD limit. In this study, we present conditions for the formation of RSDs through non-ideal MHD effects computed self-consistently from an equilibrium chemical network. We find that removing from the standard MRN distribution the large population of very small grains (VSGs) of ~10 $Å$ to few 100 $Å$ that dominate the coupling of the bulk neutral matter to the magnetic field increases the ambipolar diffusivity by ~1--2 orders of magnitude at densities below 10$^{10}$ cm$^{-3}$. The enhanced ambipolar diffusion (AD) in the envelope reduces the amount of magnetic flux dragged by the collapse into the circumstellar disk-forming region. Therefore, magnetic braking is weakened and more angular momentum can be retained. With continuous high angular momentum inflow, RSDs of tens of AU are able to form, survive, and even grow in size, depending on other parameters including cosmic-ray ionization rate, magnetic field strength, and rotation speed. Some disks become self-gravitating and evolve into rings in our 2D (axisymmetric) simulations, which have the potential to fragment into (close) multiple systems in 3D. We conclude that disk formation in magnetized cores is highly sensitive to chemistry, especially to grain sizes. A moderate grain coagulation/growth to remove the large population of VSGs, either in the prestellar phase or during free-fall collapse, can greatly promote AD and help formation of tens of AU RSDs.

preprint2016arXiv

Revolutionizing our View of Protostellar Multiplicity and Disks: The VLA Nascent Disk and Multiplicity (VANDAM) Survey of the Perseus Molecular Cloud

There is substantial evidence for disk formation taking place during the early stages of star formation and for most stars being born in multiple systems; however, protostellar multiplicity and disk searches have been hampered by low resolution, sample bias, and variable sensitivity. We have conducted an unbiased, high-sensitivity Karl G. Jansky Very Large Array (VLA) survey toward all known protostars (n = 94) in the Perseus molecular cloud (d~230 pc), with a resolution of ~15 AU (0.06") at 8 mm. We have detected candidate protostellar disks toward 17 sources (with 12 of those in the Class 0 stage) and we have found substructure on < 50AU scales for three Class 0 disk candidates, possibly evidence for disk fragmentation. We have discovered 16 new multiple systems (or new components) in this survey; the new systems have separations < 500 AU and 3 by < 30 AU. We also found a bi-modal distribution of separations, with peaks at ~75 AU and ~3000 AU, suggestive of formation through two distinct mechanisms: disk and turbulent fragmentation. The results from this survey demonstrate the necessity and utility of uniform, unbiased surveys of protostellar systems at millimeter and centimeter wavelengths.

preprint2016arXiv

Submillimeter Polarization Spectrum in the Vela C Molecular Cloud

Polarization maps of the Vela C molecular cloud were obtained at 250, 350, and 500um during the 2012 flight of the balloon-borne telescope BLASTPol. These measurements are used in conjunction with 850um data from Planck to study the submillimeter spectrum of the polarization fraction for this cloud. The spectrum is relatively flat and does not exhibit a pronounced minimum at λ~350um as suggested by previous measurements of other molecular clouds. The shape of the spectrum does not depend strongly on the radiative environment of the dust, as quantified by the column density or the dust temperature obtained from Herschel data. The polarization ratios observed in Vela C are consistent with a model of a porous clumpy molecular cloud being uniformly heated by the interstellar radiation field.

preprint2016arXiv

The VLA Nascent Disk and Multiplicity Survey of Perseus Protostars (VANDAM). II. Multiplicity of Protostars in the Perseus Molecular Cloud

We present a multiplicity study of all known protostars (94) in the Perseus molecular cloud from a Karl G. Jansky Very Large Array (VLA) survey at Ka-band (8 mm and 1 cm) and C-band (4 cm and 6.6 cm). The observed sample has a bolometric luminosity range between 0.1 L$_{\odot}$ and $\sim$33 L$_{\odot}$, with a median of 0.7 L$_{\odot}$. This multiplicity study is based on the Ka-band data, having a best resolution of $\sim$0.065" (15 AU) and separations out to $\sim$43" (10000 AU) can be probed. The overall multiplicity fraction (MF) is found to be of 0.40$\pm$0.06 and the companion star fraction (CSF) is 0.71$\pm$0.06. The MF and CSF of the Class 0 protostars are 0.57$\pm$0.09 and 1.2$\pm$0.2, and the MF and CSF of Class I protostars are both 0.23$\pm$0.08. The distribution of companion separations appears bi-modal, with a peak at $\sim$75 AU and another peak at $\sim$3000 AU. Turbulent fragmentation is likely the dominant mechanism on $>$1000 AU scales and disk fragmentation is likely to be the dominant mechanism on $<$200 AU scales. Toward three Class 0 sources we find companions separated by $<$30 AU. These systems have the smallest separations of currently known Class 0 protostellar binary systems. Moreover, these close systems are embedded within larger (50 AU to 400 AU) structures and may be candidates for ongoing disk fragmentation.

preprint2016arXiv

The VLA Nascent Disk and Multiplicity Survey: First Look at Resolved Candidate Disks around Class 0 and I Protostars in the Perseus Molecular Cloud

We present the first dust emission results toward a sample of seven protostellar disk candidates around Class 0 and I sources in the Perseus molecular cloud from the VLA Nascent Disk and Multiplicity (VANDAM) survey with ~0.05'' or 12 AU resolution. To examine the surface brightness profiles of these sources, we fit the Ka-band 8 mm dust-continuum data in the u,v-plane to a simple, parametrized model based on the Shakura-Sunyaev disk model. The candidate disks are well-fit by a model with a disk-shaped profile and have masses consistent with known Class 0 and I disks. The inner-disk surface densities of the VANDAM candidate disks have shallower density profiles compared to disks around more evolved Class II systems. The best-fit model radii of the seven early-result candidate disks are R_c > 10 AU; at 8 mm, the radii reflect lower limits on the disk size since dust continuum emission is tied to grain size and large grains radially drift inwards. These relatively large disks, if confirmed kinematically, are inconsistent with theoretical models where the disk size is limited by strong magnetic braking to < 10 AU at early times.

preprint2015arXiv

High Resolution 8 mm and 1 cm Polarization of IRAS 4A from the VLA Nascent Disk and Multiplicity (VANDAM) Survey

Magnetic fields can regulate disk formation, accretion and jet launching. Until recently, it has been difficult to obtain high resolution observations of the magnetic fields of the youngest protostars in the critical region near the protostar. The VANDAM survey is observing all known protostars in the Perseus Molecular Cloud. Here we present the polarization data of IRAS 4A. We find that with ~ 0.2'' (50 AU) resolution at λ = 8.1 and 10.3 mm, the inferred magnetic field is consistent with a circular morphology, in marked contrast with the hourglass morphology seen on larger scales. This morphology is consistent with frozen-in field lines that were dragged in by rotating material entering the infall region. The field morphology is reminiscent of rotating circumstellar material near the protostar. This is the first polarization detection of a protostar at these wavelengths. We conclude from our observations that the dust emission is optically thin with β ~ 1.3, suggesting that mm/cm-sized grains have grown and survived in the short lifetime of the protostar.

preprint2015arXiv

Local Magnetic Field Role in Star Formation

We highlight distinct and systematic observational features of magnetic field morphologies in polarized submm dust continuum. We illustrate this with specific examples and show statistical trends from a sample of 50 star-forming regions.

preprint2015arXiv

Mass Assembly of Stellar Systems and Their Evolution with the SMA (MASSES). Multiplicity and the Physical Environment in L1448N

We present continuum and molecular line observations at 230 GHz and 345 GHz from the Sub-millimeter Array (SMA) toward three protostars in the Perseus L1448N region. The data are from the large project "Mass Assembly of Stellar Systems and their Evolution with the SMA" (MASSES). Three dust continuum sources, Source B, Source NW, and Source A, are detected at both frequencies. These sources have corresponding emission peaks in C18O (J=2-1), 13CO (J=2-1), and HCO+ (J=4-3), and have offsets with N2D+ (J=3-2) peaks. High angular resolution data from a complimentary continuum survey with the Karl G. Jansky Very Large Array show that Source B is associated with three 8 mm continuum objects, Source NW with two, and Source A remains single. These results suggest that multiplicity in L1448N exists at different spatial scales from a few thousand AU to < 100 AU. Velocity gradients in each source obtained from two-dimensional fits to the SMA C18O emission are found to be perpendicular to within 20 degrees of the outflow directions as revealed by 12CO (J=2-1). We have observed that Sources B and NW with multiplicity have higher densities than Source A without multiplicity. This suggests that thermal Jeans fragmentation can be relevant in the fragmentation process. However, we have not observed a difference in the ratio between rotational and gravitational energy between sources with and without multiplicity. We also have not observed a trend between non-thermal velocity dispersions and the level of fragmentation. Our study has provided the first direct and comprehensive comparison between multiplicity and core properties in low-mass protostars, although based on small number statistics.

preprint2015arXiv

Testing protostellar disk formation models with ALMA observations

Abridged: Recent simulations have explored different ways to form accretion disks around low-mass stars. We aim to present observables to differentiate a rotationally supported disk from an infalling rotating envelope toward deeply embedded young stellar objects and infer their masses and sizes. Two 3D magnetohydrodynamics (MHD) formation simulations and 2D semi-analytical model are studied. The dust temperature structure is determined through continuum radiative transfer RADMC3D modelling. A simple temperature dependent CO abundance structure is adopted and synthetic spectrally resolved submm rotational molecular lines up to $J_{\rm u} = 10$ are simulated. All models predict similar compact components in continuum if observed at the spatial resolutions of 0.5-1$"$ (70-140 AU) typical of the observations to date. A spatial resolution of $\sim$14 AU and high dynamic range ($> 1000$) are required to differentiate between RSD and pseudo-disk in the continuum. The peak-position velocity diagrams indicate that the pseudo-disk shows a flatter velocity profile with radius than an RSD. On larger-scales, the CO isotopolog single-dish line profiles are similar and are narrower than the observed line widths of low-$J$ lines, indicating significant turbulence in the large-scale envelopes. However a forming RSD can provide the observed line widths of high-$J$ lines. Thus, either RSDs are common or a higher level of turbulence ($b \sim 0.8 \ {\rm km \ s^{-1}}$ ) is required in the inner envelope compared with the outer part. Multiple spatially and spectrally resolved molecular line observations are needed. The continuum data give a better estimate on disk masses whereas the disk sizes can be estimated from the spatially resolved molecular lines observations. The general observable trends are similar between the 2D semi-analytical models and 3D MHD RSD simulations.

Zhi-Yun Li

What is connected

Connect this record

See the researcher in context

Building this map preview

92 published item(s)

Early Planet Formation in Embedded Disks (eDisk). XX. Constraining the Chemical Tracers of Young Protostellar Sources

The Flying Saucer edge-on disc's Near Infrared silhouette revealed by the JWST JEDIce program

The JWST-NIRCam View of Sagittarius C. III. The Extinction Curve

A VLA View of the Flared, Asymmetric Disk Around the Class 0 Protostar L1527 IRS

Centrifugal Barrier and Super-Keplerian Rotation in Protostellar Disk Formation

Disks and Outflows in the Intermediate-mass Star Forming Region NGC 2071 IR

Formation of Dust Rings and Gaps in Non-ideal MHD Disks Through Meridional Gas Flows

Grain Growth During Protostellar Disk Formation

Ice Age : Chemo-dynamical modeling of Cha-MMS1 to predict new solid-phase species for detection with JWST

Magnetic Spirals in Accretion Flows Originated from Misaligned Magnetic Field

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

Polarization from Aligned Dust Grains in the $β$ Pic Debris Disk

Probing Magnetic Fields in Protoplanetary Disk Atmospheres through Polarized Optical/IR Light Scattered by Aligned Grains

The Central 1000 au of a Pre-stellar Core Revealed with ALMA. II. Almost Complete Freeze-out

The Davis-Chandrasekhar-Fermi Method Revisited

The Magnetic Field in the Milky Way Filamentary Bone G47

The Twisted Magnetic Field of the Protobinary L483

Thermal Emission and Scattering by Aligned Grains: Plane-Parallel Model and Application to Multiwavelength Polarization of the HL Tau Disk

Accretion Flows or Outflow Cavities? Uncovering the Gas Dynamics around Lupus 3-MMS

The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars V. A Characterization of Protostellar Multiplicity

Carbon Chain Chemistry in Hot-Core Regions around Three Massive Young Stellar Objects Associated with 6.7 GHz Methanol Masers

Detection of Irregular, Sub-mm Opaque Structures in the Orion Molecular Clouds: Protostars within 10000 years of formation?

Determining the physical conditions of extremely young Class 0 circumbinary disk around VLA1623A

Hall Effect in Protostellar Disc Formation and Evolution

Kinematic Analysis of a Protostellar Multiple System: Measuring the Protostar Masses and Assessing Gravitational Instability in the Disks of L1448 IRS3B and L1448 IRS3A

Probing the Temperature Structure of Optically Thick Disks Using Polarized Emission of Aligned Grains

Relative Alignment between Dense Molecular Cores and Ambient Magnetic Field: The Synergy of Numerical Models and Observations

The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars IV. Unveiling the Embedded Intermediate-Mass Protostar and Disk within OMC2-FIR3/HOPS-370

The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars. A Statistical Characterization of Class 0 and I Protostellar Disks

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

Validating Scattering-Induced (Sub)millimeter Disk Polarization through the Spectral Index, Wavelength-Dependent Polarization Pattern, and Polarization Spectrum: The Case of HD 163296

Characterizing magnetic field morphologies in three Serpens protostellar cores with ALMA

The Infall Motion in the Low-Mass Protostellar Binary NGC 1333 IRAS 4A1/4A2

The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars I. Identifying and Characterizing the Protostellar Content of the OMC2-FIR4 and OMC2-FIR3 Regions

A Triple Protostar System Formed via Fragmentation of a Gravitationally Unstable Disk

An extremely high velocity molecular jet surrounded by an ionized cavity in the protostellar source Serpens SMM1

Angular momentum loss in the envelope-disk transition region of HH 111 protostellar system: evidence for magnetic braking?

Axisymmetric Simulations of Hot Jupiter-Stellar Wind Hydrodynamic Interaction

Balloon-Borne Submillimeter Polarimetry of the Vela C Molecular Cloud: Systematic Dependence of Polarization Fraction on Column Density and Local Polarization-Angle Dispersion

Disk Polarization From Both Emission and Scattering of Magnetically Aligned Grains: The Case of NGC 1333 IRAS4A1

Instrumental performance and results from testing of the BLAST-TNG receiver, submillimeter optics, and MKID arrays

Interferometric mapping of magnetic fields: The ALMA view of the massive star forming clump W43-MM1

Protostellar Disk Formation Enabled by Removal of Small Dust Grains

Revolutionizing our View of Protostellar Multiplicity and Disks: The VLA Nascent Disk and Multiplicity (VANDAM) Survey of the Perseus Molecular Cloud

Submillimeter Polarization Spectrum in the Vela C Molecular Cloud

The VLA Nascent Disk and Multiplicity Survey of Perseus Protostars (VANDAM). II. Multiplicity of Protostars in the Perseus Molecular Cloud

The VLA Nascent Disk and Multiplicity Survey: First Look at Resolved Candidate Disks around Class 0 and I Protostars in the Perseus Molecular Cloud

High Resolution 8 mm and 1 cm Polarization of IRAS 4A from the VLA Nascent Disk and Multiplicity (VANDAM) Survey

Local Magnetic Field Role in Star Formation

Mass Assembly of Stellar Systems and Their Evolution with the SMA (MASSES). Multiplicity and the Physical Environment in L1448N

Testing protostellar disk formation models with ALMA observations

Confronting Outflow-Regulated Cluster Formation Model with Observations

Dense molecular cloud cores as a source of micrometer-sized grains in galaxies

Discovery of Two Rare Rigidly-Rotating Magnetosphere Stars in the APOGEE Survey

Magnetic Fields and Massive Star Formation

Magnetohydrodynamic Simulations of Hot Jupiter Upper Atmospheres

On the Role of Pseudodisk Warping and Reconnection in Protostellar Disk Formation in Turbulent Magnetized Cores

Recollimation boundary layers as X-ray sources in young stellar jets

Spatially Resolved Magnetic Field Structure in the Disk of a T Tauri Star

Star Cluster Formation and Feedback

The Earliest Stages of Star and Planet Formation: Core Collapse, and the Formation of Disks and Outflows

The Importance of the Magnetic Field from an SMA-CSO-Combined Sample of Star-Forming Regions

The Link between Magnetic Fields and Cloud/Star Formation

The Magnetic Field in the Class 0 Protostellar Disk of L1527

The Next Generation BLAST Experiment

The next-generation BLASTPol experiment

The VLA Nascent Disk And Multiplicity (VANDAM) Survey of Perseus Protostars. Resolving the Sub-Arcsecond Binary System in NGC 1333 IRAS2A

Alignment Between Flattened Protostellar Infall Envelopes and Ambient Magnetic Fields

Condition for the formation of micron-sized dust grains in dense molecular cloud cores

Does Magnetic Field-Rotation Misalignment Solve the Magnetic Braking Catastrophe in Protostellar Disk Formation?

Earliest Stages of Protocluster Formation: Substructure and Kinematics of Starless Cores in Orion

Effect of Magnetic Misalignment on Protobinary Evolution

Hα Absorption in Transiting Exoplanet Atmospheres

Fragmentation and OB Star Formation in High-Mass Molecular Hub-Filament System