Researcher profile

Hauyu Baobab Liu

Hauyu Baobab Liu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
25works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

25 published item(s)

preprint2024arXiv

Shadowing in the protoplanetary disk of ZZ Tau IRS with HST

An inner component misaligned from an outer component in a protoplanetary disk can result in the former casting shadows on the latter. We present a new instance of shadowing on the outer disk around a very low mass star, ZZ~Tau~IRS. Through the analysis of near-infrared (NIR) archival data at $λ=1.6$~$μ$m acquired with the Wide Field Camera 3 on the Hubble Space Telescope, we identified brightness asymmetries in the top and bottom halves of the highly inclined outer disk, separated by a dark lane. The brighter sides in the top and bottom halves are on the opposite sides, which we attributed to shadows cast by a misaligned inner disk. Radiative transfer modeling of the system with a misaligned angle of 15~deg between the inner and outer disks well reproduced the observations. Additionally, we found an elevated brightness temperature of $^{12}$CO~(3-2) at $r\sim30$~au on the brighter side in NIR wavelengths in the top half disk, which can be explained by the shadowing effect too. While the origin of the misaligned inner disk remains unclear, future monitoring observations to search for temporal variations in brightness asymmetries will likely provide useful clues.

preprint2022arXiv

Clustered Formation of Massive Stars within an Ionized Rotating Disk

We present ALMA observations with a 800 au resolution and radiative-transfer modelling of the inner part ($r\approx6000$ au) of the ionized accretion flow around a compact star cluster in formation at the center of the luminous ultra-compact (UC) HII region G10.6-0.4. We modeled the flow with an ionized Keplerian disk with and without radial motions in its outer part, or with an external Ulrich envelope. The MCMC fits to the data give total stellar masses $M_\star$ from 120 to $200~M_\odot$, with much smaller ionized-gas masses $M_\mathrm{ion-gas} = 0.2$ to $0.25~M_\odot$. The stellar mass is distributed within the gravitational radius $R_g\approx 1000$ to 1500 au, where the ionized gas is bound. The viewing inclination angle from the face-on orientation is $i = 49$ to $56~°$. Radial motions at radii $r > R_g$ converge to $v_{r,0} \approx 8.7$ km/s, or about the speed of sound of ionized gas, indicating that this gas is marginally unbound at most. From additional constraints on the ionizing-photon rate and far-IR luminosity of the region, we conclude that the stellar cluster consists of a few massive stars with $M_\mathrm{star} = 32$ to $60~M_\odot$, or one star in this range of masses accompanied by a population of lower-mass stars. Any active accretion of ionized gas onto the massive (proto)stars is residual. The inferred cluster density is very large, comparable to that reported at similar scales in the Galactic Center. Stellar interactions are likely to occur within the next Myr.

preprint2022arXiv

Formation of dust clumps with sub-Jupiter mass and cold shadowed region in gravitationally unstable disk around Class 0/I protostar in L1527 IRS

We have investigated the protostellar disk around a Class 0/I protostar, L1527 IRS, using multi-wavelength observations of the dust continuum emission at $λ=0.87$, 2.1, 3.3, and 6.8 mm obtained by the Atacama Large Millimeter/submillimeter Array (ALMA) and the Jansky Very Large Array (VLA). Our observations achieved a spatial resolution of $3-13$ au and revealed an edge-on disk structure with a size of $\sim80-100$ au. The emission at 0.87 and 2.1 mm is found to be optically thick within a projected disk radius of $ r_{\rm proj}\lesssim50$ au. The emission at 3.3 and 6.8 mm shows that the power-law index of the dust opacity ($β$) is $β\sim1.7$ around $ r_{\rm proj}\sim 50$ au, suggesting that grain growth has not yet begun. The dust temperature ($T_{\rm dust}$) shows a steep decrease with $T_{\rm dust}\propto r_{\rm proj}^{-2}$ outside of the VLA clumps previously identified at $r_{\rm proj}\sim20$ au. Furthermore, the disk is gravitationally unstable at $r_{\rm proj}\sim20$ au, as indicated by a Toomre {\it Q} parameter value of $Q\lesssim1.0$. These results suggest that the VLA clumps are formed via gravitational instability, which creates a shadow on the outside of the substructure, resulting in the sudden drop in temperature. The derived dust masses for the VLA clumps are $\gtrsim0.1$ $M_{\rm J}$. Thus, we suggest that Class 0/I disks can be massive enough to be gravitationally unstable, which might be the origin of gas-giant planets in a 20 au radius. Furthermore, the protostellar disks can be cold due to shadowing.

preprint2022arXiv

Gap Opening and Inner Disk Structure in the Strongly Accreting Transition Disk of DM Tau

Large inner dust gaps in transition disks are frequently posited as evidence of giant planets sculpting gas and dust in the disk, or the opening of a gap by photoevaporative winds. Although the former hypothesis is strongly supported by the observations of planets and deep depletions in gas within the gap some disks, many T Tauri stars hosting transition disks accrete at rates typical for an undepleted disk, raising the question of how gap opening occurs in these objects. We thus present an analysis of the structure of the transition disk around the T Tauri star DM Tau, which is strongly accreting ($\sim 10^{-8.3}~\mathrm{M}_\odot~ \mathrm{yr}^{-1}$) and turbulent ($α=0.078 \pm 0.02$). Using the DALI thermochemical code, we fit disk models to simultaneously reproduce the accretion rate, high level of turbulence, the gas traced by ALMA band 6 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O J=2--1 lines, and the observed dust emission from the mm continuum and spectral energy distribution. We find a shallow depletion in gas surface density of $\sim 10$ relative to the outer disk and a gas rich inner disk is consistent with the observations. The planet mass of $<1$ M$_\mathrm{Jup}$ implied by the gap depth is in tension with predictions for dust trapping in a highly viscous disk, which requires a more massive planet of of $\sim10$M$_\mathrm{Jup}$. Photoevaporative models including a dead zone can qualitatively reproduce some features of the DM Tau disk, but still struggle to explain the high accretion rates and the observed mm continuum flux.

preprint2022arXiv

Tracking the ice mantle history in the Solar-type Protostars of NGC 1333 IRAS 4

To understand the origin of the diversity observed in exoplanetary systems, it is crucial to characterize the early stages of their formation, represented by Solar-type protostars. Likely, the gaseous chemical content of these objects directly depends on the composition of the dust grain mantles formed before the collapse. Directly retrieving the ice mantle composition is challenging, but it can be done indirectly by observing the major components, such as NH3 and CH3OH at cm wavelengths, once they are released into the gas-phase during the warm protostellar stage. We observed several CH3OH and NH3 lines toward three Class 0 protostars in NGC1333 (IRAS 4A1, IRAS 4A2, and IRAS 4B), at high angular resolution (1&#34;; ~300 au) with the VLA interferometer at 24-26 GHz. Using a non-LTE LVG analysis, we derived a similar NH3/CH3OH abundance ratio in the three protostars (<0.5, 0.015-0.5, and 0.003-0.3 for IRAS 4A1, 4A2, and 4B, respectively). Hence, we infer they were born from pre-collapse material with similar physical conditions. Comparing the observed abundance ratios with astrochemical model predictions, we constrained the dust temperature at the time of the mantle formation to be ~17 K, which coincides with the average temperature of the southern NGC 1333 diffuse cloud. We suggest that a brutal event started the collapse that eventually formed IRAS 4A1, 4A2 and 4B, which,therefore, did not experience the usual pre-stellar core phase. This event could be the clash of a bubble with NGC 1333 south, that has previously been evoked in the literature.

preprint2021arXiv

ALMA observations of NGC 6334S. II. Subsonic and Transonic Narrow Filaments in a High-mass Star Formation Cloud

We present a study of narrow filaments toward a massive infrared dark cloud, NGC 6334S, using the Atacama Large Millimeter/submillimeter Array (ALMA). Thirteen gas filaments are identified using the H$^{13}$CO$^{+}$ line, while a single continuum filament is revealed by the continuum emission. The filaments present a compact radial distribution with a median filament width of $\sim$0.04 pc narrower than the previously proposed `quasi-universal&#39; 0.1~pc filament width. The higher spatial resolution observations and higher-density gas tracer tend to identify even narrower and lower mass filaments. The filament widths are roughly twice the size of embedded cores. The gas filaments are largely supported by thermal motions. The nonthermal motions are predominantly subsonic and transonic in both identified gas filaments and embedded cores, which may imply that stars are likely born in environments of low turbulence. A fraction of embedded objects show a narrower velocity dispersion compared with their corresponding natal filaments, which may indicate that the turbulent dissipation is taking place in these embedded cores. The physical properties (mass, mass per unit length, gas kinematics, and width) of gas filaments are analogous to those of narrow filaments found in low- to high-mass star-forming regions. The more evolved sources are found to be farther away from the filaments, a situation that may have resulted from the relative motions between the YSOs and their natal filaments.

preprint2021arXiv

Distinguishing between different mechanisms of FU-Orionis-type luminosity outbursts

Aims. Accretion and luminosity bursts triggered by three distinct mechanisms: the magnetorotational instability in the inner disk regions, clump infall in gravitationally fragmented disks and close encounters with an intruder star, were studied to determine the disk kinematic characteristics that can help to distinguish between these burst mechanisms. Methods. Numerical hydrodynamics simulations in the thin-disk limit were employed to model the bursts in disk environments that are expected for each burst mechanism. Results. We found that the circumstellar disks featuring accretion bursts can bear kinematic features that are distinct for different burst mechanisms, which can be useful when identifying the burst origin. The disks in the stellar encounter and clump-infall models are characterized by tens of percent deviations from the Keplerian rotation, whie the disks in the MRI models are characterized only a few percent deviation, which is mostly caused by the gravitational instability that fuels the MRI bursts. Velocity channel maps also show distinct kinks and wiggles, which are caused by gas disk flows that are peculiar to each considered burst mechanism. The deviations of velocity channels in the burst-hosting disks from a symmetric pattern typical of Keplerian disks are strongest for the clump-infall and collision models, and carry individual features that may be useful for the identification of the corresponding burst mechanism. The considered burst mechanisms produce a variety of light curves with the burst amplitudes varying in the Δm=2.5-3.7 limits, except for the clump-infall model where Δm can reach 5.4, although the derived numbers may be affected by a small sample and boundary conditions. Conclusions. Burst triggering mechanisms are associated with distinct kinematic features in the burst-hosting disks that may be used for their identification. Abridged.

preprint2021arXiv

Evolution of dust in protoplanetary disks of eruptive stars

Luminosity bursts in young FU Orionis-type stars warm up the surrounding disks of gas and dust, thus inflicting changes on their morphological and chemical composition. In this work, we aim at studying the effects that such bursts may have on the spatial distribution of dust grain sizes and the corresponding spectral index in protoplanetary disks. We use the numerical hydrodynamics code FEOSAD, which simulates the co-evolution of gas, dust, and volatiles in a protoplanetary disk, taking dust growth and back reaction on gas into account. The dependence of the maximum dust size on the water ice mantles is explicitly considered. The burst is initialized by increasing the luminosity of the central star to 100-300 L_sun for a time period of 100 yr. The water snowline shifts during the burst to a larger distance, resulting in the drop of the maximum dust size interior to the snowline position because of more efficient fragmentation of bare grains. After the burst, the water snowline shifts quickly back to its preburst location followed by renewed dust growth. The timescale of dust regrowth after the burst depends on the radial distance so that the dust grains at smaller distances reach the preburst values faster than the dust grains at larger distances. As a result, a broad peak in the radial distribution of the spectral index in the millimeter dust emission develops at \approx 10 au, which shifts further out as the disk evolves and dust grains regrow to preburst values at progressively larger distances. This feature is most pronounced in evolved axisymmetric disks rather than in young gravitationally unstable counterparts, although young disks may still be good candidates if gravitational instability is suppressed. Abridged.

preprint2021arXiv

Millimeter-sized Dust Grains Appear Surviving the Water-sublimating Temperature in the Inner 10 au of the FU Ori Disk

Previous observations have shown that the $\lesssim$10 au, $\gtrsim$400 K hot inner disk of the archetypal accretion outburst young stellar object, FU Ori, is dominated by viscous heating. To constrain dust properties in this region, we have performed radio observations toward this disk using the Karl G. Jansky Very Large Array (JVLA) in 2020 June-July, September, and November. We also performed complementary optical photometric monitoring observations. We found that the dust thermal emission from the hot inner disk mid-plane of FU Ori has been approximately stationary and the maximum dust grain size is $\gtrsim$1.6 mm in this region. If the hot inner disk of FU Ori which is inward of the 150-170 K water snowline is turbulent (e.g., corresponding to a Sunyaev & Shakura viscous $α_{t}\gtrsim$0.1), or if the actual maximum grain size is still larger than the lower limit we presently constrain, then as suggested by the recent analytical calculations and the laboratory measurements, water-ice free dust grains may be stickier than water-ice coated dust grains in protoplanetary disks. Additionally, we find that the free-free emission and the Johnson B and V bands magnitudes of these binary stars are brightening in 2016-2020. The optical and radio variability might be related to the dynamically evolving protostellar or disk accretion activities. Our results highlight that hot inner disks of outbursting objects are important laboratories for testing models of dust grain growth. Given the active nature of such systems, to robustly diagnose the maximum dust grain sizes, it is important to carry out coordinated multi-wavelength radio observations.

preprint2021arXiv

The evolution of temperature and density structures of OB cluster-forming molecular clumps

OB star clusters originate from parsec-scale massive molecular clumps. We aim to understand the evolution of temperature and density structures on the intermediate-scale ($\lesssim$0.1-1 pc) extended gas of massive clumps. We performed $\sim$0.1 pc resolution observations (SMA+APEX) of multiple molecular line tracers (e.g., CH$_{3}$CCH, H$_{2}$CS, CH$_{3}$CN, CH$_{3}$OH) which cover a wide range of excitation conditions, towards a sample of eight massive clumps. Based on various radiative transfer models, we constrain the gas temperature and density structures and establish an evolutionary picture, aided by a spatially-dependent virial analysis and abundance ratios of multiple species. We determine temperature radial profiles varying between 30-200 K over a continuous scale, from the center of the clumps out to 0.3-0.4 pc radii. The clumps&#39; radial gas density profiles, described by radial power-laws with slopes between -0.6 and $\sim$-1.5, are steeper for more evolved sources, as suggested by results based on both dust continuum, representing the bulk of the gas ($\sim$10$^{4}$ cm$^{-3}$), and CH$_{3}$OH lines probing the dense gas ($\gtrsim$10$^{6}$-10$^{8}$ cm$^{-3}$) regime. The density contrast between the dense gas and the bulk gas increases with evolution, and may be indicative of spatially and temporally varying star formation efficiencies. The radial profiles of the virial parameter show a global variation towards a sub-virial state as the clump evolves. The line-widths decline with increasing radius around the central core region and increase in the outer envelope, with a slope shallower than the case of the supersonic turbulence ($\,\propto\,$$r^{0.5}$) and the subsonic Kolmogorov scaling ($\,\propto\,$$r^{0.33}$). In the context of clump evolution, we also find that the abundance ratios of [CCH]/[CH$_{3}$OH] and [CH$_{3}$CN]/[CH$_{3}$OH] show correlations with clump $L/M$.

preprint2020arXiv

ALMA observations of NGC 6334S $-$ I: Forming massive stars and cluster in subsonic and transonic filamentary clouds

We present Atacama Large Millimeter/submillimeter Array (ALMA) and Karl G. Jansky Very Large Array (JVLA) observations of the massive infrared dark cloud NGC 6334S (also known as IRDC G350.56+0.44), located at the southwestern end of the NGC 6334 molecular cloud complex. The H$^{13}$CO$^{+}$ and the NH$_{2}$D lines covered by the ALMA observations at a $\sim$3$^{\prime\prime}$ angular resolution ($\sim$0.02 pc) reveal that the spatially unresolved non-thermal motions are predominantly subsonic and transonic, a condition analogous to that found in low-mass star-forming molecular clouds. The observed supersonic non-thermal velocity dispersions in massive star forming regions, often reported in the literature, might be significantly biased by poor spatial resolutions that broaden the observed line widths due to unresolved motions within the telescope beam. Our 3~mm continuum image resolves 49 dense cores, whose masses range from 0.17 to 14 $M_{\odot}$. The majority of them are resolved with multiple velocity components. Our analyses of these gas velocity components find an anti-correlation between the gas mass and the virial parameter. This implies that the more massive structures tend to be more gravitationally unstable. Finally, we find that the external pressure in the NGC 6334S cloud is important in confining these dense structures, and may play a role in the formation of dense cores, and subsequently, the embedded young stars.

preprint2020arXiv

Formation of the Hub-Filament System G33.92+0.11: Local Interplay between Gravity, Velocity, and Magnetic Field

The formation of filaments in molecular clouds is an important process in star formation. Hub-filament systems (HFSs) are a transition stage connecting parsec-scale filaments and proto-clusters. Understanding the origin of HFSs is crucial to reveal how star formation proceeds from clouds to cores. Here, we report JCMT POL-2 850 $μ$m polarization and IRAM 30-m C$^{18}$O (2-1) line observations toward the massive HFS G33.92+0.11. The 850 $μ$m continuum map reveals four major filaments converging to the center of G33.92+0.11 with numerous short filaments connecting to the major filaments at local intensity peaks. We estimate the local orientations of filaments, magnetic field, gravity, and velocity gradients from observations, and we examine their correlations based on their local properties. In the high-density areas, our analysis shows that the filaments tend to align with the magnetic field and local gravity. In the low-density areas, we find that the local velocity gradients tend to be perpendicular to both the magnetic field and local gravity, although the filaments still tend to align with local gravity. A global virial analysis suggests that the gravitational energy overall dominates the magnetic and kinematic energy. Combining local and global aspects, we conclude that the formation of G33.92+0.11 is predominantly driven by gravity, dragging and aligning the major filaments and magnetic field on the way to the inner dense center. Traced by local velocity gradients in the outer diffuse areas, ambient gas might be accreted onto the major filaments directly or via the short filaments.

preprint2020arXiv

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

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

preprint2020arXiv

Hot Corinos Chemical Diversity: Myth or Reality?

After almost 20 years of hunting, only about a dozen hot corinos, hot regions enriched in interstellar complex organic molecules (iCOMs), are known. Of them, many are binary systems with the two components showing drastically different molecular spectra. Two obvious questions arise. Why are hot corinos so difficult to find and why do their binary components seem chemically different? The answer to both questions could be a high dust opacity that would hide the molecular lines. To test this hypothesis, we observed methanol lines at centimeter wavelengths, where dust opacity is negligible, using the Very Large Array interferometer. We targeted the NGC 1333 IRAS 4A binary system, for which one of the two components, 4A1, has a spectrum deprived of iCOMs lines when observed at millimeter wavelengths, while the other component, 4A2, is very rich in iCOMs. We found that centimeter methanol lines are similarly bright toward 4A1 and 4A2. Their non-LTE analysis indicates gas density and temperature ($\geq2\times10^6$ cm$^{-3}$ and 100--190 K), methanol column density ($\sim10^{19}$ cm$^{-2}$) and extent ($\sim$35 au in radius) similar in 4A1 and 4A2, proving that both are hot corinos. Furthermore, the comparison with previous methanol line millimeter observations allows us to estimate the optical depth of the dust in front of 4A1 and 4A2, respectively. The obtained values explain the absence of iCOMs line emission toward 4A1 at millimeter wavelengths and indicate that the abundances toward 4A2 are underestimated by $\sim$30\%. Therefore, centimeter observations are crucial for the correct study of hot corinos, their census, and their molecular abundances.

preprint2020arXiv

Ionized Gas in the NGC 5253 Supernebula:High Spatial and Spectral Resolution Observations with the JVLA and TEXES

The youngest, closest and most compact embedded massive star cluster known excites the supernebula in the nearby dwarf galaxy NGC 5253. It is a crucial target and test case for studying the birth and evolution of the most massive star clusters. We present observations of the ionized gas in this source with high spatial and spectral resolution. The data includes continuum images of free-free emission with ~0.15&#39;&#39; resolution made with the JVLA at 15, 22 and 33 GHz, and a full data cube of the [SIV]10.5 micron fine-structure emission line with ~4.5 km/s velocity resolution and 0.3&#39;&#39; beam, obtained with TEXES on Gemini North. We find that 1) the ionized gas extends out from the cluster in arms or jets, and 2) the ionized gas comprises two components offset both spatially and in velocity. We discuss mechanisms that may have created the observed velocity field; possibilities include large-scale jets or a subcluster falling onto the main source.

preprint2020arXiv

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

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.

preprint2020arXiv

Ring formation by coagulation of dust aggregates in early phase of disk evolution around a protostar

Ring structures are observed by (sub-)millimeter dust continuum emission in various circumstellar disks from early stages of Class 0 and I to late stage of Class II young stellar objects (YSOs). In this paper, we study one of the possible scenarios of such ring formation in early stage, which is coagulation of dust aggregates. The dust grains grow in an inside-out manner because the growth timescale is roughly proportional to the orbital period. The boundary of the dust evolution can be regarded as the growth front, where the growth time is comparable to the disk age. With radiative transfer calculations based on the dust coagulation model, we find that the growth front can be observed as a ring structure because dust surface density is sharply changed at this position. Furthermore, we confirm that the observed ring positions in the YSOs with an age of $\lesssim1$ Myr are consistent with the growth front. The growth front could be important to create the ring structure in particular for early stage of the disk evolution such as Class 0 and I sources.

preprint2020arXiv

Substructure Formation in a Protostellar Disk of L1527 IRS

We analyze multi-frequency, high-resolution continuum data obtained by ALMA and JVLA to study detailed structure of the dust distribution in the infant disk of a Class~0/I source, L1527 IRS. We find three clumps aligning in the north-south direction in the $7 {\rm \,mm}$ radio continuum image. The three clumps remain even after subtracting free-free contamination, which is estimated from the $1.3{\rm \,cm}$ continuum observations. The northern and southern clumps are located at a distance of $\sim 15{\rm \,au}$ from the central clump and are likely optically thick at $7{\rm \,mm}$ wavelength. The clumps have similar integrated intensities. The symmetric physical properties could be realized when a dust ring or spiral arms around the central protostar is projected to the plane of the sky. We demonstrates for the first time that such substructure may form even in the disk-forming stage, where the surrounding materials actively accrete toward a disk-protostar system.

preprint2020arXiv

The MUSTANG-2 Galactic Plane Survey (MGPS90) pilot

We report the results of a pilot program for a Green Bank Telescope (GBT) MUSTANG Galactic Plane survey at 3 mm (90 GHz), MGPS90. The survey achieves a typical $1σ$ depth of $1-2$ mJy beam$^{-1}$ with a 9&#34; beam. We describe the survey parameters, quality assessment process, cataloging, and comparison with other data sets. We have identified 709 sources over seven observed fields selecting some of the most prominent millimeter-bright regions between $0°< \ell < 50°$ (total area $\approx 7.5 °^2$). The majority of these sources have counterparts at other wavelengths. By applying flux selection criteria to these sources, we successfully recovered several known hypercompact HII (HCHII) regions, but did not confirm any new ones. We identify 126 sources that have mm-wavelength counterparts but do not have cm-wavelength counterparts and are therefore candidate HCHII regions; of these, 10 are morphologically compact and are strong candidates for new HCHII regions. Given the limited number of candidates in the extended area in this survey compared to the relatively large numbers seen in protoclusters W51 and W49, it appears that most HCHII regions exist within dense protoclusters. Comparing the counts of HCHII to ultracompact HII (UCHII) regions, we infer the HCHII region lifetime is 16-46% that of the UCHII region lifetime. We additionally separated the 3 mm emission into dust and free-free emission by comparing with archival 870 $μ$m and 20 cm data. In the selected pilot fields, most ($\gtrsim80$%) of the 3 mm emission comes from plasma, either through free-free or synchrotron emission.

preprint2019arXiv

Physical Properties of the star-forming clusters in NGC 6334

We have studied the high-mass star-forming complex NGC 6334 with ALMA in the continuum emission at a frequency of 87.6 GHz, achieving a spatial resolution of 1300 au. Detecting 142 compact sources distributed over the whole observed area, we then used machine learning algorithms to group the compact cores in different clusters. A total of four main clusters were identified: NGC 6334-E, NGC 6334-I, NGC 6334-I(N) and NGC 6334-I(NW). The typical separations between cluster members (4000-12000 au) together with the core masses (0.2-100 $M_{\odot}$) are in agreement with turbulent fragmentation at scales of 0.1 pc. We find that the CMFs (core mass functions) show an excess of massive cores compared to the IMF. Typical uncertainties in temperature or unresolved multiplicity may mimic the observed excess of massive cores. Evidence of mass segregation in NGC 6334-I and NGC 6334-I(N), with the most massive cores located closer to the center, was found. Correlations between the physical properties of the clusters and their evolutionary stage were searched for, and found to show a larger separation between cores in the more evolved clusters, favouring the role of gas expulsion and stellar ejection over time. From these results, we suggest that NGC 6334-I(N) was primordially segregated, while NGC 6334-I may have become mass-segregated over time due to dynamical effects. Finally, the lack of massive cores in the most evolved cluster suggests that the gas reservoir may be already exhausted, while the less evolved clusters still have access to a larger gas mass reservoir. In general, the fragmentation in NGC 6334 at large scales (about 1 pc) is governed by turbulent pressure, while at smaller scales (a few hundred au), thermal pressure regulates the fragmentation process.

preprint2019arXiv

Resolving Linear Polarization due to Emission and Extinction of Aligned Dust Grains on NGC1333 IRAS4A with JVLA and ALMA

We report high angular resolution observations of linearly polarized dust emission towards the Class 0 young stellar object (YSO) NGC1333 IRAS4A (hereafter, IRAS4A) using the Karl G. Jansky Very Large Array (JVLA) at K (11.5-16.7 mm), Ka (8.1-10.3 mm), and Q bands (6.3-7.9 mm), and using the Atacama Large Millimeter Array (ALMA) at Band 6 (1.2 mm) and Band 7 (0.85-0.89 mm). On 100-1000 AU scales, all of these observations consistently trace the hourglass shaped magnetic field topology as shown in the previous studies. In the innermost 100 AU region of IRAS4A1, the polarization position angles (E-field) detected at 6.3-16.7 mm are consistent, however, are nearly 90 degrees offset from those detected at 1.2 mm and 0.85-0.89 mm. Such a 90 degree offset may be explained by that the inner $\sim$100 AU area is optically thick at wavelengths shorter than $\sim$1.5 mm, whereby the observations probe the absorption of aligned dust against the weakly or unpolarized warm dust emission from the innermost region. This can also consistently explain why the highest angular resolution ALMA images at Band 7 show that the polarization percentage increases with dust brightness temperature in the inner $\sim$100 AU region of IRAS4A1. Following this interpretation and assuming that the dust grains are aligned with the magnetic fields, the inferred magnetic field position angle based on the 90$^{\circ}$ rotated at 6.3-7.9 mm in the central peak of IRAS4A1 is $\sim-22^{\circ}$, which is approximately consistent with the outflow direction $\sim-9^{\circ}$.

preprint2019arXiv

Resolving the FU Ori System with ALMA: Interacting Twin Disks?

FU Orionis objects are low-mass pre-main sequence stars characterized by dramatic outbursts of several magnitudes in brightness. These outbursts are linked to episodic accretion events in which stars gain a significant portion of their mass. The physical processes behind these accretion events are not yet well understood. The archetypical FU Ori system, FU Orionis, is composed of two young stars with detected gas and dust emission. The continuum emitting regions have not been resolved until now. Here, we present 1.3 mm observations of the FU Ori binary system with ALMA. The disks are resolved at 40 mas resolution. Radiative transfer modeling shows that the emission from FU Ori north (primary) is consistent with a dust disk with a characteristic radius of $\sim$11 au. The ratio between major and minor axes shows that the inclination of the disk is $\sim$37 deg. FU Ori south is consistent with a dust disk of similar inclination and size. Assuming the binary orbit shares the same inclination angle as the disks, the deprojected distance between north and south components is 0.6&#39;&#39;, i.e. $\sim$250 au. Maps of $^{12}$CO emission show a complex kinematic environment with signatures disk rotation at the location of the northern component, and also (to a lesser extent) for FU Ori south. The revised disk geometry allows us to update FU Ori accretion models (Zhu et al.), yielding a stellar mass and mass accretion rate of FU Ori north of 0.6 M$_{\odot}$ and 3.8$\times10^{-5}$ M$_{\odot}$ yr$^{-1}$, respectively.

preprint2019arXiv

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

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

preprint2018arXiv

Investigating fragmentation of gas structures in OB cluster-forming molecular clump G33.92+0.11 with 1000 AU resolution observations of ALMA

We report new, $\sim$1000 AU spatial resolution observations of 225 GHz dust continuum emission towards the OB cluster-forming molecular clump G33.92+0.11. On parsec scales, this molecular clump presents a morphology with several arm-like dense gas structures surrounding the two central massive ($\gtrsim$100 $M_{\odot}$) cores. From the new, higher resolution observations, we identified 28 localized, spatially compact dust continuum emission sources, which may be candidates of young stellar objects. Only one of them is not embedded within known arm-like (or elongated) dense gas structures. The spatial separations of these compact sources can be very well explained by Jeans lengths. We found that G33.92+0.11 may be consistently described by a marginally centrifugally supported, Toomre unstable accretion flow which is approximately in a face-on projection. The arm-like overdensities are natural consequence of the Toomre instability, which can fragment to form young stellar objects in shorter time scales than the timescale of the global clump contraction. On our resolved spatial scales, there is not yet evidence that the fragmentation is halted by turbulence, magnetic field, or stellar feedback.

preprint2013arXiv

Probability Friends-of-Friends (PFOF) Group Finder: Performance Study and Observational Data Applications on Photometric Surveys

(Abridged) In tandem with observational datasets, we utilize realistic mock catalogs, based on a semi-analytic galaxy formation model, constructed specifically for Pan-STARRS1 Medium Deep Surveys in order to assess the performance of the Probability Friends-of-Friends (PFOF, Liu et al.) group finder, and aim to develop a grouping optimization method applicable to surveys like Pan-STARRS1. Producing mock PFOF group catalogs under a variety of photometric redshift accuracies (σΔz/(1+zs)), we find that catalog purities and completenesses from ``good&#39; σΔz/(1+zs)) ~ 0.01) to ``poor&#39; σΔz/(1+zs)) ~ 0.07) photo-zs gradually degrade respectively from 77% and 70% to 52% and 47%. To avoid model dependency of the mock for use on observational data we apply a ``subset optimization&#39; approach, using spectroscopic-redshift group data from the target field to train the group finder for application to that field, as an alternative method for the grouping optimization. We demonstrate this approach using these spectroscopically identified groups as the training set, i.e. zCOSMOS groups for PFOF searches within PS1 Medium Deep Field04 (PS1MD04) and DEEP2 EGS groups for searches in PS1MD07. We ultimately apply PFOF to four datasets spanning the photo-z uncertainty range from 0.01 to 0.06 in order to quantify the dependence of group recovery performance on photo-z accuracy. We find purities and completenesses calculated from observational datasets broadly agree with their mock analogues. Further tests of the PFOF algorithm are performed via matches to X-ray clusters identified within the PS1MD04 and COSMOS footprints. Across over a decade in group mass, we find PFOF groups match ~85% of X-ray clusters in COSMOS and PS1MD04, but at a lower statistical significance in the latter.