Showing 1–50 of 189 results for author: Myers, P C

Constraining the Stellar Masses and Origin of the Protostellar VLA 1623 System

Authors: Sarah I Sadavoy, Patrick Sheehan, John J. Tobin, Nadia M. Murillo, Richard Teague, Ian Stephens, Thomas Henning, Philip C. Myers, Edwin A. Bergin

Abstract: We present ALMA Band 7 molecular line observations of the protostars within the VLA 1623 system. We map C$^{17}$O (3 - 2) in the circumbinary disk around VLA 1623A and the outflow cavity walls of the collimated outflow. We further detect red-shifted and blue-shifted velocity gradients in the circumstellar disks around VLA 1623B and VLA 1623W that are consistent with Keplerian rotation. We use the… ▽ More We present ALMA Band 7 molecular line observations of the protostars within the VLA 1623 system. We map C$^{17}$O (3 - 2) in the circumbinary disk around VLA 1623A and the outflow cavity walls of the collimated outflow. We further detect red-shifted and blue-shifted velocity gradients in the circumstellar disks around VLA 1623B and VLA 1623W that are consistent with Keplerian rotation. We use the radiative transfer modeling code, pdspy, and simple flared disk models to measure stellar masses of $0.27 \pm 0.03$ M$_\odot$, $1.9^{+0.3}_{-0.2}$ M$_\odot$, and $0.64 \pm 0.06$ M$_\odot$ for the VLA 1623A binary, VLA 1623B, and VLA 1623W, respectively. These results represent the strongest constraints on stellar mass for both VLA 1623B and VLA 1623W, and the first measurement of mass for all stellar components using the same tracer and methodology. We use these masses to discuss the relationship between the young stellar objects (YSOs) in the VLA 1623 system. We find that VLA 1623W is unlikely to be an ejected YSO, as has been previously proposed. While we cannot rule out that VLA 1623W is a unrelated YSO, we propose that it is a true companion star to the VLA 1623A/B system and that the these stars formed in situ through turbulent fragmentation and have had only some dynamical interactions since their inception. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: Accepted to A&A; 16 pages, 12 figures

On the magnetic field properties of protostellar envelopes in Orion

Authors: Bo Huang, Josep M. Girart, Ian W. Stephens, Manuel Fernandez-Lopez, Hector G. Arce, John M. Carpenter, Paulo Cortes, Erin G. Cox, Rachel Friesen, Valentin J. M. Le Gouellec, Charles L. H. Hull, Nicole Karnath, Woojin Kwon, Zhi-Yun Li, Leslie W. Looney, Tom Megeath, Philip C. Myers, Nadia M. Murillo, Jaime E. Pineda, Sarah Sadavoy, Alvaro Sanchez-Monge, Patricio Sanhueza, John J. Tobin, Qizhou Zhang, James M. Jackson , et al. (1 additional authors not shown)

Abstract: We present 870 um polarimetric observations toward 61 protostars in the Orion molecular clouds, with ~400 au (1") resolution using the Atacama Large Millimeter/submillimeter Array. We successfully detect dust polarization and outflow emission in 56 protostars, in 16 of them the polarization is likely produced by self-scattering. Self-scattering signatures are seen in several Class 0 sources, sugge… ▽ More We present 870 um polarimetric observations toward 61 protostars in the Orion molecular clouds, with ~400 au (1") resolution using the Atacama Large Millimeter/submillimeter Array. We successfully detect dust polarization and outflow emission in 56 protostars, in 16 of them the polarization is likely produced by self-scattering. Self-scattering signatures are seen in several Class 0 sources, suggesting that grain growth appears to be significant in disks at earlier protostellar phases. For the rest of the protostars, the dust polarization traces the magnetic field, whose morphology can be approximately classified into three categories: standard-hourglass, rotated-hourglass (with its axis perpendicular to outflow), and spiral-like morphology. 40.0% (+-3.0%) of the protostars exhibit a mean magnetic field direction approximately perpendicular to the outflow on several 100--1000 au scales. However, in the remaining sample, this relative orientation appears to be random, probably due to the complex set of morphologies observed. Furthermore, we classify the protostars into three types based on the C17O (3--2) velocity envelope's gradient: perpendicular to outflow, non-perpendicular to outflow, and unresolved gradient (<1.0~km/s/arcsec). In protostars with a velocity gradient perpendicular to outflow, the magnetic field lines are preferentially perpendicular to outflow, most of them exhibit a rotated hourglass morphology, suggesting that the magnetic field has been overwhelmed by gravity and angular momentum. Spiral-like magnetic fields are associated with envelopes having large velocity gradients, indicating that the rotation motions are strong enough to twist the field lines. All of the protostars with a standard-hourglass field morphology show no significant velocity gradient due to the strong magnetic braking. △ Less

Submitted 15 May, 2024; v1 submitted 11 February, 2024; originally announced February 2024.

Comments: 12 pages, 5 figures (16 pages and 7 figures for the supplementary materials) Accepted for publication on ApJ Letters

Journal ref: ApJL, 963:L31, 2024

Most-Likely DCF Estimates of Magnetic Field Strength

Authors: Philip C. Myers, Ian W. Stephens, Simon Coudé

Abstract: The Davis-Chandrasekhar-Fermi (DCF) method is widely used to evaluate magnetic fields in star-forming regions. Yet it remains unclear how well DCF equations estimate the mean plane-of-the-sky field strength in a map region. To address this question, five DCF equations are applied to an idealized cloud map. Its polarization angles have a normal distribution with dispersion $σ_θ$,and its density and… ▽ More The Davis-Chandrasekhar-Fermi (DCF) method is widely used to evaluate magnetic fields in star-forming regions. Yet it remains unclear how well DCF equations estimate the mean plane-of-the-sky field strength in a map region. To address this question, five DCF equations are applied to an idealized cloud map. Its polarization angles have a normal distribution with dispersion $σ_θ$,and its density and velocity dispersion have negligible variation. Each DCF equation specifies a global field strength $B_{DCF}$ and a distribution of local DCF estimates. The "most-likely" DCF field strength $B_{ml}$ is the distribution mode (Chen et al. 2022), for which a correction factor $β_{ml}$ = $B_{ml}$/$B_{DCF}$ is calculated analytically. For each equation $β_{ml}$ < 1, indicating that $B_{DCF}$ is a biased estimator of $B_{ml}$. The values of $β_{ml}$ are $β_{ml}\approx$ 0.7 when $B_{DCF} \propto {σ_θ}^{-1}$ due to turbulent excitation of Afvénic MHD waves, and $β_{ml}\approx$ 0.9 when $B_{DCF} \propto {σ_θ}^{-1/2}$ due to non-Alfvénic MHD waves. These statistical correction factors $β_{ml}$ have partial agreement with correction factors $β_{sim}$ obtained from MHD simulations. The relative importance of the statistical correction is estimated by assuming that each simulation correction has both a statistical and a physical component. Then the standard, structure function, and original DCF equations appear most accurate because they require the least physical correction. Their relative physical correction factors are 0.1, 0.3, and 0.4 on a scale from 0 to 1. In contrast the large-angle and parallel-$δB$ equations have physical correction factors 0.6 and 0.7. These results may be useful in selecting DCF equations, within model limitations. △ Less

Submitted 14 December, 2023; originally announced December 2023.

Comments: Accepted for publication in The Astrophysical Journal

TRAO Survey of Nearby Filamentary Molecular clouds, the Universal Nursery of Stars (TRAO-FUNS). III. Filaments and dense cores in the NGC 2068 and NGC 2071 regions of Orion B

Authors: Hyunju Yoo, Chang Won Lee, Eun Jung Chung, Shinyoung Kim, Mario Tafalla, Paola Caselli, Philip C. Myers, Kyoung Hee Kim, Tie Liu, Woojin Kwon, Archana Soam, Jongsoo Kim

Abstract: We present the results of molecular line observations performed toward the NGC 2068 and NGC 2071 regions of the Orion B cloud as the TRAO-FUNS project to study the roles of the filamentary structure in the formation of dense cores and stars in the clouds. Gaussian decomposition for the C$^{18}$O spectra with multiple velocity components and application of a Friends-of-Friends algorithm for the dec… ▽ More We present the results of molecular line observations performed toward the NGC 2068 and NGC 2071 regions of the Orion B cloud as the TRAO-FUNS project to study the roles of the filamentary structure in the formation of dense cores and stars in the clouds. Gaussian decomposition for the C$^{18}$O spectra with multiple velocity components and application of a Friends-of-Friends algorithm for the decomposed components allowed us to identify a few tens of velocity coherent filaments. We also identified 48 dense cores from the observations of N$_2$H$^{+}$ using a core finding tool, FellWalker. We made the virial analysis for these filaments and dense cores, finding that the filaments with N$_2$H$^{+}$ dense core are thermally supercritical, and the filaments with larger ratio between the line mass and the thermal critical line mass tend to have more dense cores. We investigated the contribution of the nonthermal motions in dense cores and filaments, showing the dense cores are mostly in transonic/subsonic motions while their natal filaments are mostly in supersonic motions. This may indicates that gas turbulent motions in the filaments have been dissipated at the core scale to form the dense cores there. The filaments with (dynamically evolved) dense cores in infalling motions or with NH$_2$D bright (or chemically evolved) dense cores are all found to be gravitationally critical. Therefore, the criticality of the filament is thought to provide a key condition for its fragmentation, the formation of dense cores, and their kinematical and chemical evolution. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: 43 pages, 34 figures (more figures are available in the online journal), 4 tables, Accepted for publication in Astrophysical Journal

Alignment of dense molecular core morphology and velocity gradients with ambient magnetic fields

Authors: A. Pandhi, R. K. Friesen, L. Fissel, J. E. Pineda, P. Caselli, M. C-Y. Chen, J. Di Francesco, A. Ginsburg, H. Kirk, P. C. Myers, S. S. R. Offner, A. Punanova, F. Quan, E. Redaelli, E. Rosolowsky, S. Scibelli, Y. M. Seo, Y. Shirley

Abstract: Studies of dense core morphologies and their orientations with respect to gas flows and the local magnetic field have been limited to only a small sample of cores with spectroscopic data. Leveraging the Green Bank Ammonia Survey alongside existing sub-millimeter continuum observations and Planck dust polarization, we produce a cross-matched catalogue of 399 dense cores with estimates of core morph… ▽ More Studies of dense core morphologies and their orientations with respect to gas flows and the local magnetic field have been limited to only a small sample of cores with spectroscopic data. Leveraging the Green Bank Ammonia Survey alongside existing sub-millimeter continuum observations and Planck dust polarization, we produce a cross-matched catalogue of 399 dense cores with estimates of core morphology, size, mass, specific angular momentum, and magnetic field orientation. Of the 399 cores, 329 exhibit 2D $\mathrm{v}_\mathrm{LSR}$ maps that are well fit with a linear gradient, consistent with rotation projected on the sky. We find a best-fit specific angular momentum and core size relationship of $J/M \propto R^{1.82 \pm 0.10}$, suggesting that core velocity gradients originate from a combination of solid body rotation and turbulent motions. Most cores have no preferred orientation between the axis of core elongation, velocity gradient direction, and the ambient magnetic field orientation, favouring a triaxial and weakly magnetized origin. We find, however, strong evidence for a preferred anti-alignment between the core elongation axis and magnetic field for protostellar cores, revealing a change in orientation from starless and prestellar populations that may result from gravitational contraction in a magnetically-regulated (but not dominant) environment. We also find marginal evidence for anti-alignment between the core velocity gradient and magnetic field orientation in the L1228 and L1251 regions of Cepheus, suggesting a preferred orientation with respect to magnetic fields may be more prevalent in regions with locally ordered fields. △ Less

Submitted 24 July, 2023; originally announced July 2023.

Comments: 33 pages, 28 figures, accepted to MNRAS

Can Protostellar Outflows Set Stellar Masses?

Authors: Philip C. Myers, Michael M. Dunham, Ian W. Stephens

Abstract: The opening angles of some protostellar outflows appear too narrow to match the expected core-star mass efficiency SFE = 0.3-0.5 if outflow cavity volume traces outflow mass, with a conical shape and a maximum opening angle near 90 deg. However, outflow cavities with paraboloidal shape and wider angles are more consistent with observed estimates of the SFE. This paper presents a model of infall an… ▽ More The opening angles of some protostellar outflows appear too narrow to match the expected core-star mass efficiency SFE = 0.3-0.5 if outflow cavity volume traces outflow mass, with a conical shape and a maximum opening angle near 90 deg. However, outflow cavities with paraboloidal shape and wider angles are more consistent with observed estimates of the SFE. This paper presents a model of infall and outflow evolution based on these properties. The initial state is a truncated singular isothermal sphere which has mass $\approx$1 $M_\odot$, free fall time $\approx$80 kyr, and small fractions of magnetic, rotational, and turbulent energy. The core collapses pressure-free as its protostar and disk launch a paraboloidal wide-angle wind. The cavity walls expand radially and entrain envelope gas into the outflow. The model matches SFE values when the outflow mass increases faster than the protostar mass by a factor 1 - 2, yielding protostar masses typical of the IMF. It matches observed outflow angles if the outflow mass increases at nearly the same rate as the cavity volume. The predicted outflow angles are then typically $\sim$50 deg as they increase rapidly through the stage 0 duration of $\sim$40 kyr. They increase more slowly up to $\sim$110 deg during their stage I duration of $\sim$70 kyr. With these outflow rates and shapes, model predictions appear consistent with observational estimates of typical stellar masses, SFEs, stage durations, and outflow angles, with no need for external mechanisms of core dispersal. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Comments: Accepted for publication by The Astrophysical Journal; 47 pages, 10 figures

Models of Rotating Infall for the B335 Protostar

Authors: Neal J. Evans II, Yao-Lun Yang, Joel D. Green, Bo Zhao, James Di Francesco, Jeong-Eun Lee, Jes K. Jørgensen, Minho Choi, Philip C. Myers, Diego Mardones

Abstract: Models of the protostellar source, B335, are developed using axisymmetric three-dimensional models to resolve conflicts found in one-dimensional models. The models are constrained by a large number of observations, including ALMA, Herschel, and Spitzer data. Observations of the protostellar source B335 with ALMA show red-shifted absorption against a central continuum source indicative of infall in… ▽ More Models of the protostellar source, B335, are developed using axisymmetric three-dimensional models to resolve conflicts found in one-dimensional models. The models are constrained by a large number of observations, including ALMA, Herschel, and Spitzer data. Observations of the protostellar source B335 with ALMA show red-shifted absorption against a central continuum source indicative of infall in the HCO$^+$ and HCN $J = 4\rightarrow 3$ transitions. The data are combined with a new estimate of the distance to provide strong constraints to three-dimensional radiative transfer models including a rotating, infalling envelope, outflow cavities, and a very small disk. The models favor ages since the initiation of collapse between $3 \times 10^4$ and $4 \times 10^4$ yr for both the continuum and the lines, resolving a conflict found in one-dimensional models. The models under-predict the continuum emission seen by ALMA, suggesting an additional component such as a pseudo-disk. The best-fitting model is used to convert variations in the 4.5 $μm$ flux in recent years into a model for a variation of a factor of 5-7 in luminosity over the last 8 years. △ Less

Submitted 7 December, 2022; originally announced December 2022.

Comments: 22 pages, 20 figures, in press ApJ

Role of Filamentary Structures in the Formation of Two Dense Cores, L1544 and L694-2

Authors: Shinyoung Kim, Chang Won Lee, Mario Tafalla, Maheswar Gophinathan, Paola Caselli, Philip C. Myers, Eun Jung Chung, Shanghuo Li

Abstract: We present mapping results of two prestellar cores, L1544 and L694-2, embedded in filamentary clouds in C$^{18}$O (3-2), $^{13}$CO (3-2), $^{12}$CO (3-2), HCO$^+$ (4-3), and H$^{13}$CO$^+$ (4-3) with the JCMT telescope to examine the role of the filamentary structures in the formation of dense cores in the clouds, with new distance estimates for L1544 ($175_{-3}^{+4}$ pc) and L694-2 (… ▽ More We present mapping results of two prestellar cores, L1544 and L694-2, embedded in filamentary clouds in C$^{18}$O (3-2), $^{13}$CO (3-2), $^{12}$CO (3-2), HCO$^+$ (4-3), and H$^{13}$CO$^+$ (4-3) with the JCMT telescope to examine the role of the filamentary structures in the formation of dense cores in the clouds, with new distance estimates for L1544 ($175_{-3}^{+4}$ pc) and L694-2 ($203_{-7}^{+6}$ pc). From these observations, we found that the non-thermal velocity dispersion of two prestellar cores and their surrounding clouds are smaller than or comparable to the sound speed. This may indicate that the turbulence has already been dissipated for both filaments and cores during their formation time. We also found a $λ/4$ shift between the periodic oscillations in the velocity and the column density distributions implying the possible presence of gravitational core-forming flow motion along the axis of the filament. The mass accretion rates due to these flow motions are estimated to be 2-3 M$_\odot$ Myr$^{-1}$, being comparable to that for Serpens cloud but much smaller than those for the Hub filaments, cluster, or high mass forming filaments by 1 or 2 order of magnitudes. From this study, we suggest that the filaments in our targets might be formed from the shock-compression of colliding clouds, and then the cores are formed by gravitational fragmentation of the filaments to evolve to the prestellar stage. We conclude that the filamentary structures in the clouds play an important role in the entire process of formation of dense cores and their evolution. △ Less

Submitted 29 September, 2022; originally announced September 2022.

Comments: 30 pages, 22 figures, 6 tables, accepted for publication in ApJ

The Twisted Magnetic Field of the Protobinary L483

Authors: Erin G. Cox, Giles Novak, Sarah Sadavoy, Leslie W. Looney, Dennis Lee, Marc Berthoud, Tyler L. Bourke, Simon Coudé, Frankie Encalada, Laura M. Fissel, Rachel Harrison, Martin Houde, Zhi-Yun Li, Philip C. Myers, Kate Pattle, Fabio P. Santos, Ian W. Stephens, Hailin Wang, Sebastian Wolf

Abstract: 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 i… ▽ More 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. △ Less

Submitted 1 June, 2022; originally announced June 2022.

Comments: Accepted for publication in ApJ

Journal ref: 2022 ApJ 932 34

Virial Clumps in Central Molecular Zone Clouds

Authors: Philip C. Myers, H Perry Hatchfield, Cara Battersby

Abstract: CMZoom survey observations with the Submillimeter Array are analyzed to describe the virial equilibrium (VE) and star-forming potential of 755 clumps in 22 clouds in the Central Molecular Zone (CMZ) of the Milky Way. In each cloud, nearly all clumps follow the column-density-mass trend N~M^s, where s = 0.38 +- 0.03 is near the pressure-bounded limit s=1/3. This trend is expected when gravitational… ▽ More CMZoom survey observations with the Submillimeter Array are analyzed to describe the virial equilibrium (VE) and star-forming potential of 755 clumps in 22 clouds in the Central Molecular Zone (CMZ) of the Milky Way. In each cloud, nearly all clumps follow the column-density-mass trend N~M^s, where s = 0.38 +- 0.03 is near the pressure-bounded limit s=1/3. This trend is expected when gravitationally unbound clumps in VE have similar velocity dispersion and external pressure. Nine of these clouds also harbor one or two distinctly more massive clumps. These properties allow a VE model of bound and unbound clumps in each cloud, where the most massive clump has the VE critical mass. These models indicate that 213 clumps have velocity dispersion 1-2 km s^(-1), mean external pressure 0.5-4 x 10^8 cm^(-3) K, bound clump fraction 0.06, and typical virial parameter alpha=4-15. These mostly unbound clumps may be in VE with their turbulent cloud pressure, possibly driven by inflow from the galactic bar. In contrast, most Sgr B2 clumps are bound according to their associated sources and N-M trends. When the CMZ clumps are combined into mass distributions, their typical power-law slope is analyzed with a stopped accretion model. It also indicates that most clumps are unbound and cannot grow significantly, due to their similar time scales of accretion and dispersal, ~0.2 Myr. Thus, virial and dynamical analyses of the most extensive clump census available indicate that star formation in the CMZ may be suppressed by a significant deficit of gravitationally bound clumps. △ Less

Submitted 28 February, 2022; originally announced February 2022.

Comments: Accepted for publication, The Astrophysical Journal

The Magnetic Field in the Milky Way Filamentary Bone G47

Authors: Ian W. Stephens, Philip C. Myers, Catherine Zucker, James M. Jackson, B-G Andersson, Rowan Smith, Archana Soam, Cara Battersby, Patricio Sanhueza, Taylor Hogge, Howard A. Smith, Giles Novak, Sarah Sadavoy, Thushara Pillai, Zhi-Yun Li, Leslie W. Looney, Koji Sugitani, Simon Coude, Andres Guzman, Alyssa Goodman, Takayoshi Kusune, Fabio P. Santos, Leah Zuckerman, Frankie Encalada

Abstract: 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 magne… ▽ More 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. △ Less

Submitted 8 February, 2022; v1 submitted 27 January, 2022; originally announced January 2022.

Comments: Accepted to ApJL; typo in author list metadata corrected

Evolution and Kinematics of Protostellar Envelopes in the Perseus Molecular Cloud

Authors: Daniel J. Heimsoth, Ian W. Stephens, Hector G. Arce, Tyler L. Bourke, Philip C. Myers, Michael M. Dunham

Abstract: We present a comprehensive analysis on the evolution of envelopes surrounding protostellar systems in the Perseus molecular cloud using data from the MASSES survey. We focus our attention to the C$^{18}$O(2--1) spectral line, and we characterize the shape, size, and orientation of 54 envelopes and measure their fluxes, velocity gradients, and line widths. To look for evolutionary trends, we compar… ▽ More We present a comprehensive analysis on the evolution of envelopes surrounding protostellar systems in the Perseus molecular cloud using data from the MASSES survey. We focus our attention to the C$^{18}$O(2--1) spectral line, and we characterize the shape, size, and orientation of 54 envelopes and measure their fluxes, velocity gradients, and line widths. To look for evolutionary trends, we compare these parameters to the bolometric temperature Tbol, a tracer of protostellar age. We find evidence that the angular difference between the elongation angle of the C$^{18}$O envelope and the outflow axis direction generally becomes increasingly perpendicular with increasing Tbol, suggesting the envelope evolution is directly affected by the outflow evolution. We show that this angular difference changes at Tbol = $53 \pm 20$ K, which includes the conventional delineation between Class 0 and I protostars of 70K. We compare the C$^{18}$O envelopes with larger gaseous structures in other molecular clouds and show that the velocity gradient increases with decreasing radius ($|\mathcal{G}| \sim R^{-0.72 \pm 0.06}$). From the velocity gradients we show that the specific angular momentum follows a power law fit $J/M \propto R^{1.83 \pm 0.05}$ for scales from 1pc down to $\sim$500 au, and we cannot rule out a possible flattening out at radii smaller than $\sim$1000 au. △ Less

Submitted 18 December, 2021; originally announced December 2021.

Comments: ApJ Accepted, 19 pages, 12 figures

Are massive dense clumps truly sub-virial? A new analysis using Gould Belt ammonia data

Authors: Ayushi Singh, Christopher D. Matzner, Rachel K. Friesen, Peter G. Martin, Jaime E. Pineda, Erik W. Rosolowsky, Felipe Alves, Ana Chacón-Tanarro, Hope How-Huan Chen, Michael Chun-Yuan Chen, Spandan Choudhury, James Di Francesco, Jared Keown, Helen Kirk, Anna Punanova, Youngmin Seo, Yancy Shirley, Adam Ginsburg, Stella S. R. Offner, Héctor G. Arce, Paola Caselli, Alyssa A. Goodman, Philip C. Myers, Elena Redaelli

Abstract: Dynamical studies of dense structures within molecular clouds often conclude that the most massive clumps contain too little kinetic energy for virial equilibrium, unless they are magnetized to an unexpected degree. This raises questions about how such a state might arise, and how it might persist long enough to represent the population of massive clumps. In an effort to re-examine the origins of… ▽ More Dynamical studies of dense structures within molecular clouds often conclude that the most massive clumps contain too little kinetic energy for virial equilibrium, unless they are magnetized to an unexpected degree. This raises questions about how such a state might arise, and how it might persist long enough to represent the population of massive clumps. In an effort to re-examine the origins of this conclusion, we use ammonia line data from the Green Bank Ammonia Survey and Planck-calibrated dust emission data from Herschel to estimate the masses and kinetic and gravitational energies for dense clumps in the Gould Belt clouds. We show that several types of systematic error can enhance the appearance of low kinetic-to-gravitational energy ratios: insufficient removal of foreground and background material; ignoring the kinetic energy associated with velocity differences across a resolved cloud; and over-correcting for stratification when evaluating the gravitational energy. Using an analysis designed to avoid these errors, we find that the most massive Gould Belt clumps harbor virial motions, rather than sub-virial ones. As a byproduct, we present a catalog of masses, energies, and virial energy ratios for 85 Gould Belt clumps. △ Less

Submitted 11 August, 2021; originally announced August 2021.

Comments: Submitted to ApJ

TRAO Survey of the nearby filamentary molecular clouds, the universal nursery of stars (TRAO FUNS). II. Filaments and Dense cores in IC 5146

Authors: Eun Jung Chung, Chang Won Lee, Shinyoung Kim, Maheswar Gopinathan, Mario Tafalla, Paola Caselli, Philip C. Myers, Tie Liu, Hyunju Yoo, Kyoung Hee Kim, Mi-Ryang Kim, Archana Soam, Jungyeon Cho, Woojin Kwon, Changhoon Lee, Hyunwoo Kang

Abstract: We present the results on the physical properties of filaments and dense cores in IC 5146, as a part of the TRAO FUNS project. We carried out On-The-Fly mapping observations using the Taeduk Radio Astronomy Observatory (TRAO) 14m telescope covering about 1 square degree of the area of IC 5146 using various molecular lines. We identified 14 filaments (24 in total, including sub-filaments) from the… ▽ More We present the results on the physical properties of filaments and dense cores in IC 5146, as a part of the TRAO FUNS project. We carried out On-The-Fly mapping observations using the Taeduk Radio Astronomy Observatory (TRAO) 14m telescope covering about 1 square degree of the area of IC 5146 using various molecular lines. We identified 14 filaments (24 in total, including sub-filaments) from the C$^{18}$O (1-0) data cube and 22 dense cores from the $\rm N_{2}H^{+}$ (1-0) data. We examined the filaments' gravitational criticality, turbulence properties, accretion rate from filaments to dense cores, and relative evolutionary stages of cores. Most filaments in IC 5146 are gravitationally supercritical within the uncertainty, and most dense cores are formed in them. We found that dense cores in the hubs show a systemic velocity shift of ~0.3 km/s between the $\rm N_{2}H^{+}$ and C$^{18}$O gas. Besides, these cores are subsonic or transonic, while the surrounding filament gas is transonic or supersonic, indicating that the cores in the hubs are likely formed by the turbulence dissipation in the colliding turbulent filaments and the merging is still ongoing. We estimated the mass accretion rate of $15 - 35~M_{\odot}~\rm Myr^{-1}$ from the filaments to the dense cores, and the required time scales to collect the current core mass are consistent with the lifetime of the dense cores. The structures of filaments and dense cores in the hub can form by a collision of turbulent converging flows, and mass flow along the filaments to the dense cores may play an important role in forming dense cores. △ Less

Submitted 16 September, 2021; v1 submitted 7 June, 2021; originally announced June 2021.

Comments: 41 pages, 23 figures, 6 tables, Accepted for publication in ApJ

The Single-Cloud Star Formation Relation

Authors: Riwaj Pokhrel, Robert A. Gutermuth, Mark R. Krumholz, Christoph Federrath, Mark Heyer, Shivan Khullar, S. Thomas Megeath, Philip C. Myers, Stella S. R. Offner, Judith L. Pipher, William J. Fischer, Thomas Henning, Joseph L. Hora

Abstract: One of the most important and well-established empirical results in astronomy is the Kennicutt-Schmidt (KS) relation between the density of interstellar gas and the rate at which that gas forms stars. A tight correlation between these quantities has long been measured at galactic scales. More recently, using surveys of YSOs, a KS relationship has been found within molecular clouds relating the sur… ▽ More One of the most important and well-established empirical results in astronomy is the Kennicutt-Schmidt (KS) relation between the density of interstellar gas and the rate at which that gas forms stars. A tight correlation between these quantities has long been measured at galactic scales. More recently, using surveys of YSOs, a KS relationship has been found within molecular clouds relating the surface density of star formation to the surface density of gas; however, the scaling of these laws varies significantly from cloud to cloud. In this Letter, we use a recently developed, high-accuracy catalog of young stellar objects from $\textit{Spitzer}$ combined with high-dynamic-range gas column density maps of twelve nearby ($<$1.5 kpc) molecular clouds from $\textit{Herschel}$ to re-examine the KS relation within individual molecular clouds. We find a tight, linear correlation between clouds' star formation rate per unit area and their gas surface density normalized by the gas free-fall time. The measured intracloud KS relation, which relates star formation rate to the volume density, extends over more than two orders of magnitude within each cloud and is nearly identical in each of the twelve clouds, implying a constant star formation efficiency per free-fall time $ε_{\rm ff}\approx 0.026$. The finding of a universal correlation within individual molecular clouds, including clouds that contain no massive stars or massive stellar feedback, favors models in which star formation is regulated by local processes such as turbulence or stellar feedback such as protostellar outflows, and disfavors models in which star formation is regulated only by galaxy properties or supernova feedback on galactic scales. △ Less

Submitted 9 April, 2021; originally announced April 2021.

Comments: Accepted in ApJL. 16 pages, 4 figures

Magnetic Properties of Star-Forming Dense Cores

Authors: Philip C. Myers, Shantanu Basu

Abstract: Magnetic and energetic properties are presented for 17 dense cores within a few hundred pc of the Sun. Their plane-of-sky field strengths are estimated from the dispersion of polarization directions, following Davis, Chandrasekhar and Fermi (DCF). Their ratio of mass to magnetic critical mass is 0.5-3, indicating nearly critical field strengths. The field strength B_pos is correlated with column d… ▽ More Magnetic and energetic properties are presented for 17 dense cores within a few hundred pc of the Sun. Their plane-of-sky field strengths are estimated from the dispersion of polarization directions, following Davis, Chandrasekhar and Fermi (DCF). Their ratio of mass to magnetic critical mass is 0.5-3, indicating nearly critical field strengths. The field strength B_pos is correlated with column density N as B_pos~N^p, where p=1.05+-0.08, and with density n as B_pos~n^q, where q=0.66+-0.05. These magnetic properties are consistent with those derived from Zeeman studies (Crutcher et al. 2010), with less scatter. Relations between virial mass M_V, magnetic critical mass M_B, and Alfven amplitude sigma_B/B match the observed range of M/M_B for cores observed to be nearly virial, with M/M_V=0.5-2, with moderate Alfven amplitudes, and with sigma_B/B=0.1-0.4. The B-N and B-n correlations in the DCF and Zeeman samples can be explained when such bound, Alfvenic, and nearly-critical cores have central concentration and spheroidal shape. For these properties, B~N because M/M_B is nearly constant compared to the range of N, and B~n^(2/3) because M^(1/3) is nearly constant compared to the range of n^(2/3). The observed core fields which follow B~n^(2/3) need not be much weaker than gravity, in contrast to core fields which follow B~n^(2/3) due to spherical contraction at constant mass (Mestel 1966). Instead, the nearly critical values of M/M_B suggest that the observed core fields are nearly as strong as possible, among values which allow gravitational contraction. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Comments: Accepted for publication, The Astrophysical Journal

Gas infalling motions in the envelopes of Very Low Luminosity Objects

Authors: Mi-Ryang Kim, Chang Won Lee, Maheswar, G, Philip C. Myers, Gwanjeong Kim

Abstract: We present the results of a single dish survey toward 95 VeLLOs in optically thick (HCN 1-0) and thin ($\rm N_2H^+$ 1-0) lines performed for the purpose of understanding the physical processes of inward motions in the envelopes of the VeLLOs and characterizing their true nature. The normalized velocity differences ($δV_{HCN}$) between the peak velocities of the two lines were derived for 41 VeLLOs… ▽ More We present the results of a single dish survey toward 95 VeLLOs in optically thick (HCN 1-0) and thin ($\rm N_2H^+$ 1-0) lines performed for the purpose of understanding the physical processes of inward motions in the envelopes of the VeLLOs and characterizing their true nature. The normalized velocity differences ($δV_{HCN}$) between the peak velocities of the two lines were derived for 41 VeLLOs detected in both lines. The $δV$ distribution of these VeLLOs is found to be significantly skewed to the blue, indicating the dominance of infalling motions in their envelopes. The infall speeds were derived for 15 infall candidates by using the HILL5 radiative transfer model. The speeds were in the range of 0.03 $\rm km~s^{-1}$ to 0.3 $\rm km~s^{-1}$, with a median value of 0.16 $\rm km~s^{-1}$, being consistent with the gravitational free-fall speeds from pressure-free envelopes. The mass infall rates calculated from the infall speeds are mostly of the order of $10^{-6} M_{\odot}~yr^{-1}$ with a median value of $\rm 3.4 \pm 1.5 \times 10^{-6} M_{\odot}~yr^{-1}$. These are found to be also consistent with the values predicted with the inside-out collapse model and show a fairly good correlation with the internal luminosities of the VeLLOs. This again indicates that the infall motions observed toward the VeLLOs are likely to be due to the gravitational infall motions in their envelopes. Our study suggests that most of the VeLLOs are potentially faint protostars while two of the VeLLOs could possibly be proto-brown dwarf candidates. △ Less

Submitted 22 February, 2021; originally announced February 2021.

Magnetized filamentary gas flows feeding the young embedded cluster in Serpens South

Authors: Thushara G. S. Pillai, Dan P. Clemens, Stefan Reissl, Philip C. Myers, Jens Kauffmann, Enrique Lopez-Rodriguez, Felipe O. Alves, Gabriel P. Franco, Jonathan D. Henshaw, Karl M. Menten, Fumitaka Nakamura, Daniel Seifried, Koji Sugitani, Helmut Wiesemeyer

Abstract: Observations indicate that molecular clouds are strongly magnetized, and that magnetic fields influence the formation of stars. A key observation supporting the conclusion that molecular clouds are significantly magnetized is that the orientation of their internal structure is closely related to that of the magnetic field. At low column densities the structure aligns parallel with the field, where… ▽ More Observations indicate that molecular clouds are strongly magnetized, and that magnetic fields influence the formation of stars. A key observation supporting the conclusion that molecular clouds are significantly magnetized is that the orientation of their internal structure is closely related to that of the magnetic field. At low column densities the structure aligns parallel with the field, whereas at higher column densities, the gas structure is typically oriented perpendicular to magnetic fields, with a transition at visual extinctions $A_V\gtrsim{}3~\rm{}mag$. Here we use far-infrared polarimetric observations from the HAWC+ polarimeter on SOFIA to report the discovery of a further transition in relative orientation, i.e., a return to parallel alignment at $A_V\gtrsim{}21~\rm{}mag$ in parts of the Serpens South cloud. This transition appears to be caused by gas flow and indicates that magnetic supercriticality sets in near $A_V\gtrsim{}21~\rm{}mag$, allowing gravitational collapse and star cluster formation to occur even in the presence of relatively strong magnetic fields. △ Less

Submitted 29 September, 2020; originally announced September 2020.

Comments: 16 pages, 4 figures, Published in Nature Astronomy (August 2020). This is the authors' version before final edits. Link to the NA publication: https://www.nature.com/articles/s41550-020-1172-6

Journal ref: Nature Astronomy, August 2020

Ubiquitous $\rm NH_3$ supersonic component in L1688 coherent cores

Authors: Spandan Choudhury, Jaime E. Pineda, Paola Caselli, Adam Ginsburg, Stella S. R. Offner, Erik Rosolowsky, Rachel K. Friesen, Felipe O. Alves, Ana Chacón-Tanarro, Anna Punanova, Elena Redaelli, Helen Kirk, Philip C. Myers, Peter G. Martin, Yancy Shirley, Michael Chun-Yuan Chen, Alyssa A. Goodman, James Di Francesco

Abstract: Context : Star formation takes place in cold dense cores in molecular clouds. Earlier observations have found that dense cores exhibit subsonic non-thermal velocity dispersions. In contrast, CO observations show that the ambient large-scale cloud is warmer and has supersonic velocity dispersions. Aims : We aim to study the ammonia ($\rm NH_3$) molecular line profiles with exquisite sensitivity tow… ▽ More Context : Star formation takes place in cold dense cores in molecular clouds. Earlier observations have found that dense cores exhibit subsonic non-thermal velocity dispersions. In contrast, CO observations show that the ambient large-scale cloud is warmer and has supersonic velocity dispersions. Aims : We aim to study the ammonia ($\rm NH_3$) molecular line profiles with exquisite sensitivity towards the coherent cores in L1688 in order to study their kinematical properties in unprecedented detail. Methods : We used $\rm NH_3$ (1,1) and (2,2) data from the first data release (DR1) in the Green Bank Ammonia Survey (GAS). We first smoothed the data to a larger beam of 1' to obtain substantially more extended maps of velocity dispersion and kinetic temperature, compared to the DR1 maps. We then identified the coherent cores in the cloud and analysed the averaged line profiles towards the cores. Results : For the first time, we detected a faint (mean $\rm NH_3$(1,1) peak brightness $<$0.25 K in $T_{MB}$), supersonic component towards all the coherent cores in L1688. We fitted two components, one broad and one narrow, and derived the kinetic temperature and velocity dispersion of each component. The broad components towards all cores have supersonic linewidths ($\mathcal{M}_S \ge 1$). This component biases the estimate of the narrow dense core component's velocity dispersion by $\approx$28% and the kinetic temperature by $\approx$10%, on average, as compared to the results from single-component fits. Conclusions : Neglecting this ubiquitous presence of a broad component towards all coherent cores causes the typical single-component fit to overestimate the temperature and velocity dispersion. This affects the derived detailed physical structure and stability of the cores estimated from $\rm NH_3$ observations. △ Less

Submitted 20 July, 2020; v1 submitted 14 July, 2020; originally announced July 2020.

Comments: Accepted for publication in Astronomy & Astrophysics on 06/07/2020. 15 pages, 16 figures, 1 table. Language edits from previous version

Journal ref: A&A 640, L6 (2020)

Star-Gas Surface Density Correlations in Twelve Nearby Molecular Clouds I: Data Collection and Star-Sampled Analysis

Authors: Riwaj Pokhrel, Robert A. Gutermuth, Sarah K. Betti, Stella S. R. Offner, Philip C. Myers, S. Thomas Megeath, Alyssa D. Sokol, Babar Ali, Lori Allen, Tom S. Allen, Michael M. Dunham, William J. Fischer, Thomas Henning, Mark Heyer, Joseph L. Hora, Judith L. Pipher, John J. Tobin, Scott J. Wolk

Abstract: We explore the relation between the stellar mass surface density and the mass surface density of molecular hydrogen gas in twelve nearby molecular clouds that are located at $<$1.5 kpc distance. The sample clouds span an order of magnitude range in mass, size, and star formation rates. We use thermal dust emission from $Herschel$ maps to probe the gas surface density and the young stellar objects… ▽ More We explore the relation between the stellar mass surface density and the mass surface density of molecular hydrogen gas in twelve nearby molecular clouds that are located at $<$1.5 kpc distance. The sample clouds span an order of magnitude range in mass, size, and star formation rates. We use thermal dust emission from $Herschel$ maps to probe the gas surface density and the young stellar objects from the most recent $Spitzer$ Extended Solar Neighborhood Archive (SESNA) catalog to probe the stellar surface density. Using a star-sampled nearest neighbor technique to probe the star-gas surface density correlations at the scale of a few parsecs, we find that the stellar mass surface density varies as a power-law of the gas mass surface density, with a power-law index of $\sim$2 in all the clouds. The consistent power-law index implies that star formation efficiency is directly correlated with gas column density, and no gas column density threshold for star formation is observed. We compare the observed correlations with the predictions from an analytical model of thermal fragmentation, and with the synthetic observations of a recent hydrodynamic simulation of a turbulent star-forming molecular cloud. We find that the observed correlations are consistent for some clouds with the thermal fragmentation model and can be reproduced using the hydrodynamic simulations. △ Less

Submitted 11 May, 2020; originally announced May 2020.

Comments: 39 pages, 10 figures, 3 tables. Accepted in ApJ

Magnetic Field Structure in Spheroidal Star-Forming Clouds. II. Estimating Field Structure from Observed Maps

Authors: Philip C. Myers, Ian W. Stephens, Sayantan Auddy, Shantanu Basu, Tyler L. Bourke, Charles L. H. Hull

Abstract: This paper presents models to estimate the structure of density and magnetic field strength in spheroidal condensations, from maps of their column density and their polarization of magnetically aligned dust grains. The density model is obtained by fitting a column density map with an embedded p = 2 Plummer spheroid of any aspect ratio and inclination. The magnetic properties are based on the densi… ▽ More This paper presents models to estimate the structure of density and magnetic field strength in spheroidal condensations, from maps of their column density and their polarization of magnetically aligned dust grains. The density model is obtained by fitting a column density map with an embedded p = 2 Plummer spheroid of any aspect ratio and inclination. The magnetic properties are based on the density model, on the Davis-Chandrasekhar-Fermi (DCF) model of Alfvénic fluctuations, and on the Spheroid Flux Freezing (SFF) model of mass and flux conservation in Paper I. The field strength model has the resolution of the column density map, which is finer than the resolution of the DCF estimate of field strength. The models are applied to ALMA observations of the envelope of the protostar BHR71 IRS1. Column density fits give the density model, from (2.0 +- 0.4) x 10^5 cm^-3 to (7 +- 1) x 10^7 cm^-3 . The density model predicts the field directions map, which fits the polarization map best within 1100 au, with standard deviation of angle differences 17°. In this region the DCF mean field strength is 0.7 +- 0.2 mG and the envelope mass is supercritical, with ratio of mass to magnetic critical mass 1.5 +- 0.4. The SFF field strength profile scales with the DCF field strength, from 60 x 10μG to 4+-1 mG. The spatial resolution of the SFF field strength estimate is finer than the DCF resolution by a factor ~7, and the peak SFF field strength exceeds the DCF field strength by a factor ~5. △ Less

Submitted 8 May, 2020; originally announced May 2020.

Comments: Accepted for publication in The Astrophysical Journal

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

Authors: Che-Yu Chen, Erica A. Behrens, Jasmin E. Washington, Laura M. Fissel, Rachel K. Friesen, Zhi-Yun Li, Jaime E. Pineda, Adam Ginsburg, Helen Kirk, Samantha Scibelli, Felipe Alves, Elena Redaelli, Paola Caselli, Anna Punanova, James Di Francesco, Erik Rosolowsky, Stella S. R. Offner, Peter G. Martin, Ana Chacón-Tanarro, Hope H. -H. Chen, Michael C. -Y. Chen, Jared Keown, Youngmin Seo, Yancy Shirley, Hector G. Arce , et al. (4 additional authors not shown)

Abstract: 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 cl… ▽ More 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. △ Less

Submitted 24 March, 2020; originally announced March 2020.

Comments: 18 pages, 11 figures, accepted for publication in MNRAS

CS Depletion in Prestellar Cores

Authors: Shinyoung Kim, Chang Won Lee, Maheswar Gopinathan, Mario Tafalla, Jungjoo Sohn, Gwanjeong Kim, Mi-Ryang Kim, Archana Soam, Philip C. Myers

Abstract: The CS molecule is known to be absorbed onto dust in the cold and dense conditions, causing it to get significantly depleted in the central region of cores. This study is aimed to investigate the depletion of the CS molecule using the optically thin C$^{34}$S molecular line observations. We mapped five prestellar cores, L1544, L1552, L1689B, L694-2, and L1197 using two molecular lines, C$^{34}$S… ▽ More The CS molecule is known to be absorbed onto dust in the cold and dense conditions, causing it to get significantly depleted in the central region of cores. This study is aimed to investigate the depletion of the CS molecule using the optically thin C$^{34}$S molecular line observations. We mapped five prestellar cores, L1544, L1552, L1689B, L694-2, and L1197 using two molecular lines, C$^{34}$S $(J=2-1)$ and N$_2$H$^+$ $(J=1-0)$ with the NRO 45-m telescope, doubling the number of cores where the CS depletion was probed using C$^{34}$S. In most of our targets, the distribution of C$^{34}$S emission shows features that suggest that the CS molecule is generally depleted in the center of the prestellar cores. The radial profile of the CS abundance with respect to H$_2$ directly measured from the CS emission and the Herschel dust emission indicates that the CS molecule is depleted by a factor of $\sim$3 toward the central regions of the cores with respect to their outer regions. The degree of the depletion is found to be even more enhanced by an order of magnitude when the contaminating effect introduced by the presence of CS molecules in the surrounding envelope that lie along the line-of-sight is removed. Except for L1197 which is classified as relatively the least evolved core in our targets based on its observed physical parameters, we found that the remaining four prestellar cores are suffering from significant CS depletion at their central region regardless of the relative difference in their evolutionary status. △ Less

Submitted 17 February, 2020; originally announced February 2020.

Comments: 20 pages, 9 figures, Accepted for publication in ApJ

Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) -- Full Data Release

Authors: Ian W. Stephens, Tyler L. Bourke, Michael M. Dunham, Philip C. Myers, Riwaj Pokhrel, John J. Tobin, Héctor G. Arce, Sarah I. Sadavoy, Eduard I. Vorobyov, Jaime E. Pineda, Stella S. R. Offner, Katherine I. Lee, Lars E. Kristensen, Jes K. Jørgensen, Mark A. Gurwell, Alyssa A. Goodman

Abstract: We present and release the full dataset for the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. This survey used the Submillimeter Array (SMA) to image the 74 known protostars within the Perseus molecular cloud. The SMA was used in two array configurations to capture outflows for scales $>$30$^{\prime\prime}$ ($>$9000 au) and to probe scales down to $\sim$1… ▽ More We present and release the full dataset for the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. This survey used the Submillimeter Array (SMA) to image the 74 known protostars within the Perseus molecular cloud. The SMA was used in two array configurations to capture outflows for scales $>$30$^{\prime\prime}$ ($>$9000 au) and to probe scales down to $\sim$1$^{\prime\prime}$ ($\sim$300 au). The protostars were observed with the 1.3 mm and 850 $��$m receivers simultaneously to detect continuum at both wavelengths and molecular line emission from CO(2-1), $^{13}$CO(2-1), C$^{18}$O(2-1), N$_2$D$^+$(3-2), CO(3-2), HCO$^+$(4-3), and H$^{13}$CO$^+$(4-3). Some of the observations also used the SMA's recently upgraded correlator, SWARM, whose broader bandwidth allowed for several more spectral lines to be observed (e.g., SO, H$_2$CO, DCO$^+$, DCN, CS, CN). Of the main continuum and spectral tracers observed, 84% of the images and cubes had emission detected. The median C$^{18}$O(2-1) linewidth is $\sim$1.0 km s$^{-1}$, which is slightly higher than those measured with single-dish telescopes at scales of 3000-20000 au. Of the 74 targets, six are suggested to be first hydrostatic core candidates, and we suggest that L1451-mm is the best candidate. We question a previous continuum detection toward L1448 IRS2E. In the SVS13 system, SVS13A certainly appears to be the most evolved source, while SVS13C appears to be hotter and more evolved than SVS13B. The MASSES survey is the largest publicly available interferometric continuum and spectral line protostellar survey to date, and is largely unbiased as it only targets protostars in Perseus. All visibility ($uv$) data and imaged data are publicly available at https://dataverse.harvard.edu/dataverse/full_MASSES/. △ Less

Submitted 19 November, 2019; originally announced November 2019.

Comments: Accepted to ApJS

Dust Polarization Toward Embedded Protostars in Ophiuchus with ALMA. III. Survey Overview

Authors: Sarah I. Sadavoy, Ian W. Stephens, Philip C. Myers, Leslie Looney, John Tobin, Woojin Kwon, Benoit Commercon, Dominique Segura-Cox, Thomas Henning, Patrick Hennebelle

Abstract: We present 0.25 arcsec resolution (35 au) ALMA 1.3 mm dust polarization observations for 37 young stellar objects (YSOs) in the Ophiuchus cloud. These data encompass all the embedded protostars in the cloud and several Flat and Class II objects to produce the largest, homogeneous study of dust polarization on disk scales to date. The goal of this study is to study dust polarization down to disk sc… ▽ More We present 0.25 arcsec resolution (35 au) ALMA 1.3 mm dust polarization observations for 37 young stellar objects (YSOs) in the Ophiuchus cloud. These data encompass all the embedded protostars in the cloud and several Flat and Class II objects to produce the largest, homogeneous study of dust polarization on disk scales to date. The goal of this study is to study dust polarization down to disk scales. We find that 14/37 (38%) of the YSOs are detected in polarization. Nine of these sources have uniform polarization angles and four sources have azimuthal polarization. The sources with uniform polarization tend to have steeper inclinations (> 60 degree) than those with azimuthal polarization (< 60 degree). The majority (9/14) of the detected sources have polarization morphologies and disk properties consistent with dust self-scattering in optically thick disks. The remaining sources may be instead tracing magnetic fields. Their inferred field directions from rotating the polarization vectors by 90 degree are mainly poloidal or hourglass shaped. We find no evidence of a strong toroidal field component toward any of our disks. For the 23 YSOs that are undetected in polarization, roughly half of them have 3-sigma upper limits of < 2%. These sources also tend to have inclinations < 60 degree and they are generally compact. Since lower inclination sources tend to have azimuthal polarization, these YSOs may be undetected in polarization due to unresolved polarization structure within our beam. We propose that disks with inclinations > 60 degree are the best candidates for future polarization studies of dust self-scattering as these systems will generally show uniform polarization vectors that do not require very high resolution to resolve. We release the continuum and polarization images for all the sources with this publication. Data from the entire survey can be obtained from Dataverse. △ Less

Submitted 5 September, 2019; originally announced September 2019.

Comments: Accepted to ApJS, 52 pages, 52 figures, full data release on Dataverse

KFPA Examinations of Young STellar Object Natal Environments (KEYSTONE): Hierarchical Ammonia Structures in Galactic Giant Molecular Clouds

Authors: Jared Keown, James Di Francesco, Erik Rosolowsky, Ayushi Singh, Charles Figura, Helen Kirk, L. D. Anderson, Michael Chun-Yuan Chen, Davide Elia, Rachel Friesen, Adam Ginsburg, A. Marston, Stefano Pezzuto, Eugenio Schisano, Sylvain Bontemps, Paola Caselli, Hong-Li Liu, Steven Longmore, Frederique Motte, Philip C. Myers, Stella S. R. Offner, Patricio Sanhueza, Nicola Schneider, Ian Stephens, James Urquhart , et al. (1 additional authors not shown)

Abstract: We present initial results from the K-band focal plane array Examinations of Young STellar Object Natal Environments (KEYSTONE) survey, a large project on the 100-m Green Bank Telescope mapping ammonia emission across eleven giant molecular clouds at distances of $0.9-3.0$ kpc (Cygnus X North, Cygnus X South, M16, M17, MonR1, MonR2, NGC2264, NGC7538, Rosette, W3, and W48). This data release includ… ▽ More We present initial results from the K-band focal plane array Examinations of Young STellar Object Natal Environments (KEYSTONE) survey, a large project on the 100-m Green Bank Telescope mapping ammonia emission across eleven giant molecular clouds at distances of $0.9-3.0$ kpc (Cygnus X North, Cygnus X South, M16, M17, MonR1, MonR2, NGC2264, NGC7538, Rosette, W3, and W48). This data release includes the NH$_3$ (1,1) and (2,2) maps for each cloud, which are modeled to produce maps of kinetic temperature, centroid velocity, velocity dispersion, and ammonia column density. Median cloud kinetic temperatures range from $11.4\pm2.2$ K in the coldest cloud (MonR1) to $23.0\pm6.5$ K in the warmest cloud (M17). Using dendrograms on the NH$_3$ (1,1) integrated intensity maps, we identify 856 dense gas clumps across the eleven clouds. Depending on the cloud observed, $40-100\%$ of the clumps are aligned spatially with filaments identified in H$_2$ column density maps derived from SED-fitting of dust continuum emission. A virial analysis reveals that 523 of the 835 clumps ($\sim63\%$) with mass estimates are bound by gravity alone. We find no significant difference between the virial parameter distributions for clumps aligned with the dust-continuum filaments and those unaligned with filaments. In some clouds, however, hubs or ridges of dense gas with unusually high mass and low virial parameters are located within a single filament or at the intersection of multiple filaments. These hubs and ridges tend to host water maser emission, multiple 70$μ$m-detected protostars, and have masses and radii above an empirical threshold for forming massive stars. △ Less

Submitted 29 August, 2019; v1 submitted 27 August, 2019; originally announced August 2019.

Comments: Accepted for publication in ApJ

The specific angular momentum radial profile in dense cores: improved initial conditions for disk formation

Authors: Jaime E. Pineda, Bo Zhao, Anika Schmiedeke, Dominique M. Segura-Cox, Paola Caselli, Philip C. Myers, John Tobin, Michael Dunham

Abstract: The determination of the specific angular momentum radial profile, $j(r)$, in the early stages of star formation is crucial to constrain star and circumstellar disk formation theories. The specific angular momentum is directly related to the largest Keplerian disk possible, and it could constrain the angular momentum removal mechanism. We determine $j(r)$ towards two Class 0 objects and a first hy… ▽ More The determination of the specific angular momentum radial profile, $j(r)$, in the early stages of star formation is crucial to constrain star and circumstellar disk formation theories. The specific angular momentum is directly related to the largest Keplerian disk possible, and it could constrain the angular momentum removal mechanism. We determine $j(r)$ towards two Class 0 objects and a first hydrostatic core candidate in the Perseus cloud, which is consistent across all three sources and well fit with a single power-law relation between 800 and 10,000\,au: $j_{fit}(r)=10^{-3.60\pm0.15}\left(r/\textrm{1,000au}\right)^{1.80\pm0.04}$ km s$^{-1}$ pc. This power-law relation is in between solid body rotation ($\propto r^2$) and pure turbulence ($\propto r^{1.5}$). This strongly suggests that even at 1,000\,au, the influence of the dense core's initial level of turbulence or the connection between core and the molecular cloud is still present. The specific angular momentum at 10,000\,au is $\approx3\times$ higher than previously estimated, while at 1,000\,au it is lower by $2\times$. We do not find a region of conserved specific angular momentum, although it could still be present at a smaller radius. We estimate an upper limit to the largest Keplerian disk radius of 60\,au, which is small but consistent with published upper limits. Finally, these results suggest that more realistic initial conditions for numerical simulations of disk formation are needed. Some possible solutions include: a) use a larger simulation box to include some level of driven turbulence or connection to the parental cloud, or b) incorporate the observed $j(r)$ to setup the dense core kinematics initial conditions. △ Less

Submitted 25 June, 2019; v1 submitted 13 June, 2019; originally announced June 2019.

Comments: Accepted to ApJ, 13 pages, 7 figures

TRAO Survey of Nearby Filamentary Molecular clouds, the Universal Nursery of Stars (TRAO FUNS) I. Dynamics and Chemistry of L1478 in the California Molecular Cloud

Authors: Eun Jung Chung, Chang Won Lee, Shinyoung Kim, Gwanjeong Kim, Paola Caselli, Mario Tafalla, Philip C. Myers, Archana Soam, Tie Liu, Maheswar Gopinathan, Miryang Kim, Kyoung Hee Kim, Woojin Kwon, Hyunwoo Kang, Changhoon Lee

Abstract: "TRAO FUNS" is a project to survey Gould Belt's clouds in molecular lines. This paper presents its first results on the central region of the California molecular cloud, L1478. We performed On-The-Fly mapping observations using the Taedeok Radio Astronomy Observatory (TRAO) 14m single dish telescope equipped with a 16 multi-beam array covering $\sim$1.0 square degree area of this region using C… ▽ More "TRAO FUNS" is a project to survey Gould Belt's clouds in molecular lines. This paper presents its first results on the central region of the California molecular cloud, L1478. We performed On-The-Fly mapping observations using the Taedeok Radio Astronomy Observatory (TRAO) 14m single dish telescope equipped with a 16 multi-beam array covering $\sim$1.0 square degree area of this region using C$^{18}$O (1-0) mainly tracing low density cloud and about 460 square arcminute area using N$_{2}$H$^{+}$ (1-0) mainly tracing dense cores. CS (2-1) and SO $(3_{2}-2_{1})$ were also used simultaneously to map $\sim$440 square arcminute area of this region. We identified 10 filaments by applying the dendrogram technique to the C$^{18}$O data-cube and 8 dense N$_{2}$H$^{+}$ cores by using {\sc FellWalker}. Basic physical properties of filaments such as mass, length, width, velocity field, and velocity dispersion are derived. It is found that L1478 consists of several filaments with slightly different velocities. Especially the filaments which are supercritical are found to contain dense cores detected in N$_{2}$H$^{+}$. Comparison of non-thermal velocity dispersions derived from C$^{18}$O and N$_{2}$H$^{+}$ for the filaments and dense cores indicates that some of dense cores share similar kinematics with those of the surrounding filaments while several dense cores have different kinematics with those of their filaments. This suggests that the formation mechanism of dense cores and filaments can be different in individual filaments depending on their morphologies and environments. △ Less

Submitted 23 April, 2019; originally announced April 2019.

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

Efficient methanol production on the dark side of a prestellar core

Authors: Jorma Harju, Jaime E. Pineda, Anton I. Vasyunin, Paola Caselli, Stella S. R. Offner, Alyssa A. Goodman, Mika Juvela, Olli Sipilae, Alexandre Faure, Romane Le Gal, Pierre Hily-Blant, Joao Alves, Luca Bizzocchi, Andreas Burkert, Hope Chen, Rachel K. Friesen, Rolf Guesten, Philip C. Myers, Anna Punanova, Claire Rist, Erik Rosolowsky, Stephan Schlemmer, Yancy Shirley, Silvia Spezzano, Charlotte Vastel , et al. (1 additional authors not shown)

Abstract: We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particul… ▽ More We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particularly efficient. Methanol and sulphur monoxide are most abundant in a narrow zone that follows the eastern side of the core. This side is sheltered from the stronger external radiation field coming from the west. We show that photodissociation on the illuminated side can give rise to an asymmetric methanol distribution, but that the stark contrast observed in H-MM1 is hard to explain without assuming enhanced desorption on the shaded side. The region of the brightest emission has a wavy structure that rolls up at one end. This is the signature of Kelvin-Helmholtz instability occurring in sheared flows. We suggest that in this zone, methanol and sulphur are released as a result of grain-grain collisions induced by shear vorticity. △ Less

Submitted 4 May, 2020; v1 submitted 27 March, 2019; originally announced March 2019.

Comments: Accepted for publication in the Astrophysical Journal

The Mass Evolution of Protostellar Disks and Envelopes in the Perseus Molecular Cloud

Authors: Bridget C. Andersen, Ian W. Stephens, Michael M. Dunham, Riwaj Pokhrel, Jes K. Jørgensen, Søren Frimann, Dominique Segura-Cox, Philip C. Myers, Tyler L. Bourke, John J. Tobin, Łukasz Tychoniec

Abstract: In the standard picture for low-mass star formation, a dense molecular cloud undergoes gravitational collapse to form a protostellar system consisting of a new central star, a circumstellar disk, and a surrounding envelope of remaining material. The mass distribution of the system evolves as matter accretes from the large-scale envelope through the disk and onto the protostar. While this general p… ▽ More In the standard picture for low-mass star formation, a dense molecular cloud undergoes gravitational collapse to form a protostellar system consisting of a new central star, a circumstellar disk, and a surrounding envelope of remaining material. The mass distribution of the system evolves as matter accretes from the large-scale envelope through the disk and onto the protostar. While this general picture is supported by simulations and indirect observational measurements, the specific timescales related to disk growth and envelope dissipation remain poorly constrained. In this paper we conduct a rigorous test of a method introduced by Jørgensen et al. (2009) to obtain mass measurements of disks and envelopes around embedded protostars with observations that do not resolve the disk (resolution of $\sim$1000\,AU). Using unresolved data from the recent Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey, we derive disk and envelope mass estimates for $59$ protostellar systems in the Perseus molecular cloud. We compare our results to independent disk mass measurements from the VLA Nascent Disk and Multiplicity (VANDAM) survey and find a strong linear correlation, suggesting that accurate disk masses can be measured from unresolved observations. Then, leveraging the size of the MASSES sample, we find no significant trend in protostellar mass distribution as a function of age, as approximated from bolometric temperatures. These results may indicate that the disk mass of a protostar is set near the onset of the Class 0 protostellar stage and remains roughly constant throughout the Class I protostellar stage. △ Less

Submitted 15 February, 2019; originally announced February 2019.

Comments: Accepted to ApJ

Magnetic Field Structure of Dense Cores using Spectroscopic Methods

Authors: Sayantan Auddy, Philip C. Myers, Shantanu Basu, Jorma Harju, Jaime E. Pineda, Rachel K. Friesen

Abstract: We develop a new ''core field structure'' (CFS) model to predict the magnetic field strength and magnetic field fluctuation profile of dense cores using gas kinematics. We use spatially resolved observations of the nonthermal velocity dispersion from the Green Bank Ammonia survey along with column density maps from SCUBA-2 to estimate the magnetic field strength across seven dense cores located in… ▽ More We develop a new ''core field structure'' (CFS) model to predict the magnetic field strength and magnetic field fluctuation profile of dense cores using gas kinematics. We use spatially resolved observations of the nonthermal velocity dispersion from the Green Bank Ammonia survey along with column density maps from SCUBA-2 to estimate the magnetic field strength across seven dense cores located in the L1688 region of Ophiuchus. The CFS model predicts the profile of the relative field fluctuation, which is related to the observable dispersion in direction of the polarization vectors. Within the context of our model we find that all the cores have a transcritical mass-to-flux ratio. △ Less

Submitted 28 January, 2019; originally announced January 2019.

Comments: 18 pages, 13 figures, to appear in ApJ

Catalogue of High Protostellar Surface Density Regions in Nearby Embedded Clusters

Authors: Juan Li, Philip C. Myers, Helen Kirk, Robert A. Gutermuth, Michael M. Dunham, Riwaj Pokhrel

Abstract: We analyze high-quality stellar catalogs for 24 young and nearby (within 1 kpc) embedded clusters and present a catalogue of 32 groups which have a high concentration of protostars. The median effective radius of these groups is 0.17 pc. The median protostellar and pre-main sequence star surface densities are 46 M_{\odot} pc^{-2} and 11 M_{\odot} pc^{-2}, respectively. We estimate the age of these… ▽ More We analyze high-quality stellar catalogs for 24 young and nearby (within 1 kpc) embedded clusters and present a catalogue of 32 groups which have a high concentration of protostars. The median effective radius of these groups is 0.17 pc. The median protostellar and pre-main sequence star surface densities are 46 M_{\odot} pc^{-2} and 11 M_{\odot} pc^{-2}, respectively. We estimate the age of these groups using a model of constant birthrate and random accretion stopping and find a median value of 0.25 Myr. Some groups in Aquila and Serpens, Corona Australia and Ophichus L1688 show high protostellar surface density and high molecular gas surface density, which seem to be undergoing vigorous star formation. These groups provide an excellent opportunity to study initial conditions of clustered star formation. Comparison of protostellar and pre-main-sequence stellar surface densities reveal continuous low-mass star formation of these groups over several Myr in some clouds. For groups with typical protostellar separations of less than 0.4 pc, we find that these separations agree well with the thermal Jeans fragmentation scale. On the other hand, for groups with typical protostellar separations larger than 0.4 pc, these separations are always larger than the associated Jeans length. △ Less

Submitted 3 December, 2018; originally announced December 2018.

Comments: 29 pages, 24 figures, to be published in ApJ

Dust Polarization Toward Embedded Protostars in Ophiuchus with ALMA. II. IRAS 16293-2422

Authors: Sarah I. Sadavoy, Philip C. Myers, Ian W. Stephens, John Tobin, Woojin Kwon, Dominique Segura-Cox, Thomas Henning, Benoit Commercon, Leslie Looney

Abstract: We present high-resolution (~ 35 au) ALMA Band 6 1.3 mm dust polarization observations of IRAS 16293. These observations spatially resolve the dust polarization across the two protostellar sources and toward the filamentary structures between them. The dust polarization and inferred magnetic field have complicated structures throughout the region. In particular, we find that the magnetic field is… ▽ More We present high-resolution (~ 35 au) ALMA Band 6 1.3 mm dust polarization observations of IRAS 16293. These observations spatially resolve the dust polarization across the two protostellar sources and toward the filamentary structures between them. The dust polarization and inferred magnetic field have complicated structures throughout the region. In particular, we find that the magnetic field is aligned parallel to three filamentary structures. We characterize the physical properties of the filamentary structure that bridges IRAS 16293A and IRAS 16293B and estimate a magnetic field strength of 23-78 mG using the Davis-Chandrasekhar-Fermi method. We construct a toy model for the bridge material assuming that the young stars dominate the mass and gravitational potential of the system. We find that the expected gas flow to each star is of comparable order to the Alfven speed, which suggests that the field may be regulating the gas flow. We also find that the bridging material should be depleted in ~ 1000 yr. If the bridge is part of the natal filament that formed the stars, then it must have accreted new material. Alternatively, the bridge could be a transient structure. Finally, we show that the 1.3 mm polarization morphology of the optically thick IRAS 16293B system is qualitatively similar to dust self-scattering. Based on similar polarization measurements at 6.9 mm, we propose that IRAS 16293B has produced a substantial population of large dust grains with sizes between 200 and 2000 um. △ Less

Submitted 13 December, 2018; v1 submitted 7 November, 2018; originally announced November 2018.

Comments: Accepted to ApJ, 15 pages, 9 figures; minor corrections for typos

Magnetic Field Structure in Spheroidal Star-Forming Clouds

Authors: Philip C. Myers, Shantanu Basu, Sayantan Auddy

Abstract: A model of magnetic field structure is presented to help test the prevalence of flux freezing in star-forming clouds of various shapes, orientations, and degrees of central concentration, and to estimate their magnetic field strength. The model is based on weak-field flux freezing in centrally condensed Plummer spheres and spheroids of oblate and prolate shape. For a spheroid of given density cont… ▽ More A model of magnetic field structure is presented to help test the prevalence of flux freezing in star-forming clouds of various shapes, orientations, and degrees of central concentration, and to estimate their magnetic field strength. The model is based on weak-field flux freezing in centrally condensed Plummer spheres and spheroids of oblate and prolate shape. For a spheroid of given density contrast, aspect ratio, and inclination, the model estimates the local field strength and direction, and the global field pattern of hourglass shape. Comparisons with a polarization simulation indicate typical angle agreement within 1 - 10 degrees. Scalable analytic expressions are given to match observed polarization patterns, and to provide inputs to radiative transfer codes for more accurate predictions. The model may apply to polarization observations of dense cores, elongated filamentary clouds, and magnetized circumstellar disks. △ Less

Submitted 9 October, 2018; originally announced October 2018.

Comments: accepted, Astrophysical Journal

High-resolution ALMA Study of the Proto-Brown-Dwarf Candidate L328-IRS

Authors: Chang Won Lee, Gwanjeong Kim, Philip C. Myers, Masao Saito, Shinyoung Kim, Woojin Kwon, A-Ran Lyo, Archana Soam, Mi-Ryang Kim

Abstract: This paper presents our observational attempts to precisely measure the central mass of a proto-brown dwarf candidate, L328-IRS, in order to investigate whether L328-IRS is in the substellar mass regime. Observations were made for the central region of L328-IRS with the dust continuum and CO isotopologue line emission at ALMA band 6, discovering the detailed outflow activities and a deconvolved di… ▽ More This paper presents our observational attempts to precisely measure the central mass of a proto-brown dwarf candidate, L328-IRS, in order to investigate whether L328-IRS is in the substellar mass regime. Observations were made for the central region of L328-IRS with the dust continuum and CO isotopologue line emission at ALMA band 6, discovering the detailed outflow activities and a deconvolved disk structure of a size of $\sim$87 AU $\times$ 37 AU. We investigated the rotational velocities as a function of the disk radius, finding that its motions between 130 AU and 60 AU are partially fitted with a Keplerian orbit by a stellar object of $\sim0.30 ~\rm\,M_\odot$, while the motions within 60 AU do not follow any Keplerian orbit at all. This makes it difficult to lead a reliable estimation of the mass of L328-IRS. Nonetheless, our ALMA observations were useful enough to well constrain the inclination angle of the outflow cavity of L328-IRS as $\sim 66^{\circ}$, enabling us to better determine the mass accretion rate of $\sim 8.9\times 10^{-7} \rm\,M_\odot~yr^{-1}$. From assumptions that the internal luminosity of L328-IRS is mostly due to this mass accretion process in the disk, or that L328-IRS has mostly accumulated the mass through this constant accretion rate during its outflow activity, its mass was estimated to be $\sim 0.012 - 0.023\rm\,M_\odot$, suggesting L328-IRS to be a substellar object. However, we leave our identification of L328-IRS as a proto-brown dwarf to be tentative because of various uncertainties especially regarding the mass accretion rate. △ Less

Submitted 27 September, 2018; originally announced September 2018.

Comments: Accepted in ApJ, 24 pages, 12 figures

Droplets I: Pressure-Dominated Sub-0.1 pc Coherent Structures in L1688 and B18

Authors: Hope How-Huan Chen, Jaime E. Pineda, Alyssa A. Goodman, Andreas Burkert, Stella S. R. Offner, Rachel K. Friesen, Philip C. Myers, Felipe Alves, Hector G. Arce, Paola Caselli, Ana Chacon-Tanarro, Michael Chun-Yuan Chen, James Di Francesco, Adam Ginsburg, Jared Keown, Helen Kirk, Peter G. Martin, Christopher Matzner, Anna Punanova, Elena Redaelli, Erik Rosolowsky, Samantha Scibelli, Young Min Seo, Yancy Shirley, Ayushi Singh

Abstract: We present the observation and analysis of newly discovered coherent structures in the L1688 region of Ophiuchus and the B18 region of Taurus. Using data from the Green Bank Ammonia Survey (GAS), we identify regions of high density and near-constant, almost-thermal, velocity dispersion. Eighteen coherent structures are revealed, twelve in L1688 and six in B18, each of which shows a sharp "transiti… ▽ More We present the observation and analysis of newly discovered coherent structures in the L1688 region of Ophiuchus and the B18 region of Taurus. Using data from the Green Bank Ammonia Survey (GAS), we identify regions of high density and near-constant, almost-thermal, velocity dispersion. Eighteen coherent structures are revealed, twelve in L1688 and six in B18, each of which shows a sharp "transition to coherence" in velocity dispersion around its periphery. The identification of these structures provides a chance to study the coherent structures in molecular clouds statistically. The identified coherent structures have a typical radius of 0.04 pc and a typical mass of 0.4 Msun, generally smaller than previously known coherent cores identified by Goodman et al. (1998), Caselli et al. (2002), and Pineda et al. (2010). We call these structures "droplets." We find that unlike previously known coherent cores, these structures are not virially bound by self-gravity and are instead predominantly confined by ambient pressure. The droplets have density profiles shallower than a critical Bonnor-Ebert sphere, and they have a velocity (VLSR) distribution consistent with the dense gas motions traced by NH3 emission. These results point to a potential formation mechanism through pressure compression and turbulent processes in the dense gas. We present a comparison with a magnetohydrodynamic simulation of a star-forming region, and we speculate on the relationship of droplets with larger, gravitationally bound coherent cores, as well as on the role that droplets and other coherent structures play in the star formation process. △ Less

Submitted 15 May, 2019; v1 submitted 26 September, 2018; originally announced September 2018.

Comments: Accepted by ApJ in April, 2019

Journal ref: 2019ApJ...877...93C

Mass Assembly of Stellar Systems and their Evolution with the SMA - 1.3 mm Subcompact Data Release

Authors: Ian W. Stephens, Michael M. Dunham, Philip C. Myers, Riwaj Pokhrel, Tyler L. Bourke, Eduard I. Vorobyov, John J. Tobin, Sarah I. Sadavoy, Jaime E. Pineda, Stella S. R. Offner, Katherine I. Lee, Lars E. Kristensen, Jes K. Jørgensen, Alyssa A. Goodman, Héctor G. Arce, Mark Gurwell

Abstract: We present the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey, which uses the Submillimeter Array (SMA) interferometer to map the continuum and molecular lines for all 74 known Class 0/I protostellar systems in the Perseus molecular cloud. The primary goal of the survey is to observe an unbiased sample of young protostars in a single molecular cloud so that we ca… ▽ More We present the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey, which uses the Submillimeter Array (SMA) interferometer to map the continuum and molecular lines for all 74 known Class 0/I protostellar systems in the Perseus molecular cloud. The primary goal of the survey is to observe an unbiased sample of young protostars in a single molecular cloud so that we can characterize the evolution of protostars. This paper releases the MASSES 1.3 mm data from the subcompact configuration ($\sim$4$^{\prime\prime}$ or $\sim$1000 au resolution), which is the SMA's most compact array configuration. We release both $uv$ visibility data and imaged data for the spectral lines CO(2-1),$^{13}$CO(2-1), C$^{18}$O(2-1), and N$_2$D$^+$(3-2), as well as for the 1.3 mm continuum. We identify the tracers that are detected toward each source. We also show example images of continuum and CO(2-1) outflows, analyze C$^{18}$O(2-1) spectra, and present data from the SVS 13 star-forming region. The calculated envelope masses from the continuum show a decreasing trend with bolometric temperature (a proxy for age). Typical C$^{18}$O(2-1) linewidths are 1.45 km s$^{-1}$, which is higher than the C$^{18}$O linewidths detected toward Perseus filaments and cores. We find that N$_2$D$^+$(3-2) is significantly more likely to be detected toward younger protostars. We show that the protostars in SVS 13 are contained within filamentary structures as traced by C$^{18}$O(2-1) and N$_2$D$^+$(3-2). We also present the locations of SVS 13A's high velocity (absolute line-of-sight velocities $>$150 km s$^{-1}$) red and blue outflow components. Data can be downloaded from https://dataverse.harvard.edu/dataverse/MASSES . △ Less

Submitted 19 June, 2018; originally announced June 2018.

Comments: Accepted ApJS

Dust Polarization Toward Embedded Protostars in Ophiuchus with ALMA. I. VLA 1623

Authors: Sarah I. Sadavoy, Philip C. Myers, Ian W. Stephens, John Tobin, Benoit Commercon, Thomas Henning, Leslie Looney, Woojin Kwon, Dominique Segura-Cox, Robert Harris

Abstract: We present high resolution (~ 30 au) ALMA Band 6 dust polarization observations of VLA 1623. The VLA 1623 data resolve compact ~ 40 au inner disks around the two protobinary sources, VLA 1623-A and VLA 1623-B, and also an extended ~ 180 au ring of dust around VLA 1623-A. This dust ring was previously identified as a large disk in lower-resolution observations. We detect highly structured dust pola… ▽ More We present high resolution (~ 30 au) ALMA Band 6 dust polarization observations of VLA 1623. The VLA 1623 data resolve compact ~ 40 au inner disks around the two protobinary sources, VLA 1623-A and VLA 1623-B, and also an extended ~ 180 au ring of dust around VLA 1623-A. This dust ring was previously identified as a large disk in lower-resolution observations. We detect highly structured dust polarization toward the inner disks and the extended ring with typical polarization fractions ~ 1.7% and ~ 2.4%, respectively. The two components also show distinct polarization morphologies. The inner disks have uniform polarization angles aligned with their minor axes. This morphology is consistent with expectations from dust scattering. By contrast, the extended dust ring has an azimuthal polarization morphology not previously seen in lower-resolution observations. We find that our observations are well fit by a static, oblate spheroid model with a flux-frozen, poloidal magnetic field. We propose that the polarization traces magnetic grain alignment likely from flux freezing on large scales and magnetic diffusion on small scales. Alternatively, the azimuthal polarization may be attributed to grain alignment by the anisotropic radiation field. If the grains are radiatively aligned, then our observations indicate that large (~ 100 um) dust grains grow quickly at large angular extents. Finally, we identify significant proper motion of VLA 1623 using our observations and those in the literature. This result indicates that the proper motion of nearby systems must be corrected for when combining ALMA data from different epochs. △ Less

Submitted 5 June, 2018; v1 submitted 16 April, 2018; originally announced April 2018.

Comments: Accepted in ApJ; 15 pages, 13 figures

Filamentary Fragmentation and Accretion in High-Mass Star-Forming Molecular Clouds

Authors: Xing Lu, Qizhou Zhang, Hauyu Baobab Liu, Patricio Sanhueza, Ken'ichi Tatematsu, Siyi Feng, Howard A. Smith, Philip C. Myers, T. K. Sridharan, Qiusheng Gu

Abstract: Filamentary structures are ubiquitous in high-mass star-forming molecular clouds. Their relation with high-mass star formation is still to be understood. Here we report interferometric observations toward 8 filamentary high-mass star-forming clouds. A total of 50 dense cores are identified in these clouds, most of which present signatures of high-mass star formation. Five of them are not associate… ▽ More Filamentary structures are ubiquitous in high-mass star-forming molecular clouds. Their relation with high-mass star formation is still to be understood. Here we report interferometric observations toward 8 filamentary high-mass star-forming clouds. A total of 50 dense cores are identified in these clouds, most of which present signatures of high-mass star formation. Five of them are not associated with any star formation indicators, hence are prestellar core candidates. Evolutionary phases of these cores and their linewidths, temperatures, NH$_3$ abundances, and virial parameters are found to be correlated. In a sub-sample of 4 morphologically well-defined filaments, we find that their fragmentation can not be solely explained by thermal or turbulence pressure support. We also investigate distributions of gas temperatures and non-thermal motions along the filaments, and find a spatial correlation between non-thermal linewidths and star formation activities. We find evidence of gas flows along these filaments, and derive an accretion rate along filaments of $\sim$10$^{-4}$ M$_\odot$ yr$^{-1}$. These results suggest a strong relationship between massive filaments and high-mass star formation, through i) filamentary fragmentation in very early evolutionary phases to form dense cores, ii) accretion flows along filaments that are important for the growth of dense cores and protostars, and iii) enhancement of non-thermal motion in the filaments by the feedback or accretion during star formation. △ Less

Submitted 18 January, 2018; originally announced January 2018.

Comments: 31 pages, 11 figures, 8 tables, ApJ Accepted. 3-D interactive figures to visualize outflows (see Figure 7) will be published online, and are available at https://xinglunju.github.io/outflows.html

Hierarchical fragmentation in the Perseus molecular cloud: From the cloud scale to protostellar objects

Authors: Riwaj Pokhrel, Philip C. Myers, Michael M. Dunham, Ian W. Stephens, Sarah I. Sadavoy, Qizhou Zhang, Tyler L. Bourke, John J. Tobin, Katherine I. Lee, Robert A. Gutermuth, Stella S. R. Offner

Abstract: We present a study of hierarchical structure in the Perseus molecular cloud, from the scale of the entire cloud ($\gtrsim$10 pc) to smaller clumps ($\sim$1 pc), cores ($\sim$0.05-0.1 pc), envelopes ($\sim$300-3000 AU) and protostellar objects ($\sim$15 AU). We use new observations from the Submillimeter Array (SMA) large project "Mass Assembly of Stellar Systems and their Evolution with the SMA (M… ▽ More We present a study of hierarchical structure in the Perseus molecular cloud, from the scale of the entire cloud ($\gtrsim$10 pc) to smaller clumps ($\sim$1 pc), cores ($\sim$0.05-0.1 pc), envelopes ($\sim$300-3000 AU) and protostellar objects ($\sim$15 AU). We use new observations from the Submillimeter Array (SMA) large project "Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES)" to probe the envelopes, and recent single-dish and interferometric observations from the literature for the remaining scales. This is the first study to analyze hierarchical structure over five scales in the same cloud complex. We compare the number of fragments with the number of Jeans masses in each scale to calculate the Jeans efficiency, or the ratio of observed to expected number of fragments. The velocity dispersion is assumed to arise either from purely thermal motions, or from combined thermal and non-thermal motions inferred from observed spectral line widths. For each scale, thermal Jeans fragmentation predicts more fragments than observed, corresponding to inefficient thermal Jeans fragmentation. For the smallest scale, thermal plus non-thermal Jeans fragmentation also predicts too many protostellar objects. However at each of the larger scales thermal plus non-thermal Jeans fragmentation predicts fewer than one fragment, corresponding to no fragmentation into envelopes, cores, and clumps. Over all scales, the results are inconsistent with complete Jeans fragmentation based on either thermal or thermal plus non-thermal motions. They are more nearly consistent with inefficient thermal Jeans fragmentation, where the thermal Jeans efficiency increases from the largest to the smallest scale. △ Less

Submitted 13 December, 2017; originally announced December 2017.

Comments: Accepted in ApJ. 21 pages, 15 figures, 4 tables

Intensity-Corrected Herschel Observations of Nearby Isolated Low-Mass Clouds

Authors: Sarah I. Sadavoy, Eric Keto, Tyler L. Bourke, Michael M. Dunham, Philip C. Myers, Ian W. Stephens, James Di Francesco, Kristi Webb, Amelia Stutz, Ralf Launhardt, John Tobin

Abstract: We present intensity-corrected Herschel maps at 100 um, 160 um, 250 um, 350 um, and 500 um for 56 isolated low-mass clouds. We determine the zero-point corrections for Herschel PACS and SPIRE maps from the Herschel Science Archive (HSA) using Planck data. Since these HSA maps are small, we cannot correct them using typical methods. Here, we introduce a technique to measure the zero-point correctio… ▽ More We present intensity-corrected Herschel maps at 100 um, 160 um, 250 um, 350 um, and 500 um for 56 isolated low-mass clouds. We determine the zero-point corrections for Herschel PACS and SPIRE maps from the Herschel Science Archive (HSA) using Planck data. Since these HSA maps are small, we cannot correct them using typical methods. Here, we introduce a technique to measure the zero-point corrections for small Herschel maps. We use radial profiles to identify offsets between the observed HSA intensities and the expected intensities from Planck. Most clouds have reliable offset measurements with this technique. In addition, we find that roughly half of the clouds have underestimated HSA-SPIRE intensities in their outer envelopes relative to Planck, even though the HSA-SPIRE maps were previously zero-point corrected. Using our technique, we produce corrected Herschel intensity maps for all 56 clouds and determine their line-of-sight average dust temperatures and optical depths from modified black body fits. The clouds have typical temperatures of ~ 14-20 K and optical depths of ~ 1e-5 - 1e-3. Across the whole sample, we find an anti-correlation between temperature and optical depth. We also find lower temperatures than what was measured in previous Herschel studies, which subtracted out a background level from their intensity maps to circumvent the zero-point correction. Accurate Herschel observations of clouds are key to obtain accurate density and temperature profiles. To make such future analyses possible, intensity-corrected maps for all 56 clouds are publicly available in the electronic version. △ Less

Submitted 30 November, 2017; originally announced December 2017.

Comments: Accepted to ApJ, 35 pages, 15 figures, 4 Tables. Intensity-corrected FITS images and temperature/optical depths maps are available with the electronic version

The Green Bank Ammonia Survey: Unveiling the Dynamics of the Barnard 59 star-forming Clump

Authors: E. Redaelli, F. O. Alves, P. Caselli, J. E. Pineda, R. K. Friesen, A. Chacón-Tanarro, C. D. Matzner, A. Ginsburg, E. Rosolowsky, J. Keown, S. S. R. Offner, J. Di Francesco, H. Kirk, P. C. Myers, A. Hacar, A. Cimatti, H. H. Chen, M. C. Chen, Y. M. Seo, K. I. Lee

Abstract: Understanding the early stages of star formation is a research field of ongoing development, both theoretically and observationally. In this context, molecular data have been continuously providing observational constraints on the gas dynamics at different excitation conditions and depths in the sources. We have investigated the Barnard 59 core, the only active site of star formation in the Pipe N… ▽ More Understanding the early stages of star formation is a research field of ongoing development, both theoretically and observationally. In this context, molecular data have been continuously providing observational constraints on the gas dynamics at different excitation conditions and depths in the sources. We have investigated the Barnard 59 core, the only active site of star formation in the Pipe Nebula, to achieve a comprehensive view of the kinematic properties of the source. These information were derived by simultaneously fitting ammonia inversion transition lines (1,1) and (2,2). Our analysis unveils the imprint of protostellar feedback, such as increasing line widths, temperature and turbulent motions in our molecular data. Combined with complementary observations of dust thermal emission, we estimate that the core is gravitationally bound following a virial analysis. If the core is not contracting, another source of internal pressure, most likely the magnetic field, is supporting it against gravitational collapse and limits its star formation efficiency. △ Less

Submitted 6 November, 2017; originally announced November 2017.

Comments: 18 pages, 18 figures

The Green Bank Ammonia Survey: Observations of Hierarchical Dense Gas Structures in Cepheus-L1251

Authors: Jared Keown, James Di Francesco, Helen Kirk, Rachel K. Friesen, Jaime E. Pineda, Erik Rosolowsky, Adam Ginsburg, Stella S. R. Offner, Paola Caselli, Felipe Alves, Ana Chacón-Tanarro, Anna Punanova, Elena Redaelli, Young Min Seo, Christopher D. Matzner, Michael Chun-Yuan Chen, Alyssa A. Goodman, How-Huan Chen, Yancy Shirley, Ayushi Singh, Hector G. Arce, Peter Martin, Philip C. Myers

Abstract: We use Green Bank Ammonia Survey observations of NH$_3$ (1,1) and (2,2) emission with 32'' FWHM resolution from a ~ 10 pc$^{2}$ portion of the Cepheus-L1251 molecular cloud to identify hierarchical dense gas structures. Our dendrogram analysis of the NH$_3$ data results in 22 top-level structures, which reside within 13 lower-level, parent structures. The structures are compact (0.01 pc… ▽ More We use Green Bank Ammonia Survey observations of NH$_3$ (1,1) and (2,2) emission with 32'' FWHM resolution from a ~ 10 pc$^{2}$ portion of the Cepheus-L1251 molecular cloud to identify hierarchical dense gas structures. Our dendrogram analysis of the NH$_3$ data results in 22 top-level structures, which reside within 13 lower-level, parent structures. The structures are compact (0.01 pc $\lesssim R_{eff} \lesssim$ 0.1 pc) and are spatially correlated with the highest H$_2$ column density portions of the cloud. We also compare the ammonia data to a catalog of dense cores identified by higher-resolution (18.2'' FWHM) Herschel Space Observatory observations of dust continuum emission from Cepheus-L1251. Maps of kinetic gas temperature, velocity dispersion, and NH$_3$ column density, derived from detailed modeling of the NH$_3$ data, are used to investigate the stability and chemistry of the ammonia-identified and Herschel-identified structures. We show that the dust and dense gas in the structures have similar temperatures, with median $T_{dust}$ and $T_K$ measurements of 11.7 $\pm$ 1.1 K and 10.3 $\pm$ 2.0 K, respectively. Based on a virial analysis, we find that the ammonia-identified structures are gravitationally dominated, yet may be in or near a state of virial equilibrium. Meanwhile, the majority of the Herschel-identified dense cores appear to be not bound by their own gravity and instead confined by external pressure. CCS $(2_0-1_0)$ and HC$_5$N $(9-8)$ emission from the region reveal broader line widths and centroid velocity offsets when compared to the NH$_3$ (1,1) emission in some cases, likely due to these carbon-based molecules tracing the turbulent outer layers of the dense cores. △ Less

Submitted 12 October, 2017; originally announced October 2017.

Comments: Accepted for publication in ApJ

The Green Bank Ammonia Survey: Dense Cores Under Pressure in Orion A

Authors: Helen Kirk, Rachel K. Friesen, Jaime E. Pineda, Erik Rosolowsky, Stella S. R. Offner, Christopher D. Matzner, Philip C. Myers, James Di Francesco, Paola Caselli, Felipe O. Alves, Ana Chacón-Tanarro, How-Huan Chen, Michael Chun-Yuan Chen, Jared Keown, Anna Punanova, Young Min Seo, Yancy Shirley, Adam Ginsburg, Christine Hall, Ayushi Singh, Héctor G. Arce, Alyssa A. Goodman, Peter Martin, Elena Redaelli

Abstract: We use gas temperature and velocity dispersion data from the Green Bank Ammonia Survey and core masses and sizes from the James Clerk Maxwell Telescope Gould Belt Survey to estimate the virial states of dense cores within the Orion A molecular cloud. Surprisingly, we find that almost none of the dense cores are sufficiently massive to be bound when considering only the balance between self-gravity… ▽ More We use gas temperature and velocity dispersion data from the Green Bank Ammonia Survey and core masses and sizes from the James Clerk Maxwell Telescope Gould Belt Survey to estimate the virial states of dense cores within the Orion A molecular cloud. Surprisingly, we find that almost none of the dense cores are sufficiently massive to be bound when considering only the balance between self-gravity and the thermal and non-thermal motions present in the dense gas. Including the additional pressure binding imposed by the weight of the ambient molecular cloud material and additional smaller pressure terms, however, suggests that most of the dense cores are pressure confined. △ Less

Submitted 17 August, 2017; originally announced August 2017.

Comments: Accepted for publication in ApJ Aug 10, 2017

Alignment Between Protostellar Outflows and Filamentary Structure

Authors: Ian W. Stephens, Michael M. Dunham, Philip C. Myers, Riwaj Pokhrel, Sarah I. Sadavoy, Eduard I. Vorobyov, John J. Tobin, Jaime E. Pineda, Stella S. R. Offner, Katherine I. Lee, Lars E. Kristensen, Jes K. Jørgensen, Alyssa A. Goodman, Tyler L. Bourke, Héctor G. Arce, Adele L. Plunkett

Abstract: We present new Submillimeter Array (SMA) observations of CO(2-1) outflows toward young, embedded protostars in the Perseus molecular cloud as part of the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. For 57 Perseus protostars, we characterize the orientation of the outflow angles and compare them with the orientation of the local filaments as derived from… ▽ More We present new Submillimeter Array (SMA) observations of CO(2-1) outflows toward young, embedded protostars in the Perseus molecular cloud as part of the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. For 57 Perseus protostars, we characterize the orientation of the outflow angles and compare them with the orientation of the local filaments as derived from $Herschel$ observations. We find that the relative angles between outflows and filaments are inconsistent with purely parallel or purely perpendicular distributions. Instead, the observed distribution of outflow-filament angles are more consistent with either randomly aligned angles or a mix of projected parallel and perpendicular angles. A mix of parallel and perpendicular angles requires perpendicular alignment to be more common by a factor of $\sim$3. Our results show that the observed distributions probably hold regardless of the protostar's multiplicity, age, or the host core's opacity. These observations indicate that the angular momentum axis of a protostar may be independent of the large-scale structure. We discuss the significance of independent protostellar rotation axes in the general picture of filament-based star formation. △ Less

Submitted 25 July, 2017; originally announced July 2017.

Comments: Accepted to ApJ

An Estimation of the Star Formation Rate in the Perseus Complex

Authors: Seyma Mercimek, Philip C. Myers, Katherine I. Lee, Sarah I. Sadavoy

Abstract: We present the results of our investigation of the star-forming potential in the Perseus star-forming complex. We build on previous starless core, protostellar core, and young stellar object (YSO) catalogs from Spitzer, Herschel, and SCUBA observations in the literature. We place the cores and YSOs within seven star-forming clumps based on column densities greater than 5x10^21 cm^-2. We calculate… ▽ More We present the results of our investigation of the star-forming potential in the Perseus star-forming complex. We build on previous starless core, protostellar core, and young stellar object (YSO) catalogs from Spitzer, Herschel, and SCUBA observations in the literature. We place the cores and YSOs within seven star-forming clumps based on column densities greater than 5x10^21 cm^-2. We calculate the mean density and free-fall time for 69 starless cores as 5.55x10^-19 gcm^-3 and 0.1 Myr,respectively, and we estimate the star formation rate for the near future as 150 Msun Myr^-1. According to Bonnor Ebert stability analysis, we find that majority of starless cores in Perseus are unstable. Broadly, these cores can collapse to form the next generation of stars. We found a relation between starless cores and YSOs, where the numbers of young protostars (Class 0 + Class I) are similar to the numbers of starless cores. This similarity, which shows a one-to-one relation, suggests that these starless cores may form the next generation of stars with approximately the same formation rate as the current generation, as identified by the Class 0 and Class I protostars. It follows that if such a relation between starless cores and any YSO stage exists, the SFR values of these two populations must be nearly constant. In brief, we propose that this one-to-one relation is an important factor in better understanding the star formation process within a cloud. △ Less

Submitted 25 April, 2017; originally announced April 2017.

Comments: Published in AJ, 2017, 153, 214; 14 pages, 7 figures

Journal ref: 2017AJ....153..214M

The Green Bank Ammonia Survey (GAS): First Results of NH3 mapping the Gould Belt

Authors: Rachel K. Friesen, Jaime E. Pineda, Erik Rosolowsky, Felipe Alves, Ana Chacón-Tanarro, Hope How-Huan Chen, Michael Chun-Yuan Chen, James Di Francesco, Jared Keown, Helen Kirk, Anna Punanova, Youngmin Seo, Yancy Shirley, Adam Ginsburg, Christine Hall, Stella S. R. Offner, Ayushi Singh, Héctor G. Arce, Paola Caselli, Alyssa A. Goodman, Peter G. Martin, Christopher Matzner, Philip C. Myers, Elena Redaelli

Abstract: We present an overview of the first data release (DR1) and first-look science from the Green Bank Ammonia Survey (GAS). GAS is a Large Program at the Green Bank Telescope to map all Gould Belt star-forming regions with $A_V \gtrsim 7$ mag visible from the northern hemisphere in emission from NH$_3$ and other key molecular tracers. This first release includes the data for four regions in Gould Belt… ▽ More We present an overview of the first data release (DR1) and first-look science from the Green Bank Ammonia Survey (GAS). GAS is a Large Program at the Green Bank Telescope to map all Gould Belt star-forming regions with $A_V \gtrsim 7$ mag visible from the northern hemisphere in emission from NH$_3$ and other key molecular tracers. This first release includes the data for four regions in Gould Belt clouds: B18 in Taurus, NGC 1333 in Perseus, L1688 in Ophiuchus, and Orion A North in Orion. We compare the NH$_3$ emission to dust continuum emission from Herschel, and find that the two tracers correspond closely. NH$_3$ is present in over 60\% of lines-of-sight with $A_V \gtrsim 7$ mag in three of the four DR1 regions, in agreement with expectations from previous observations. The sole exception is B18, where NH$_3$ is detected toward ~ 40\% of lines-of-sight with $A_V \gtrsim 7$ mag. Moreover, we find that the NH$_3$ emission is generally extended beyond the typical 0.1 pc length scales of dense cores. We produce maps of the gas kinematics, temperature, and NH$_3$ column densities through forward modeling of the hyperfine structure of the NH$_3$ (1,1) and (2,2) lines. We show that the NH$_3$ velocity dispersion, $σ_v$, and gas kinetic temperature, $T_K$, vary systematically between the regions included in this release, with an increase in both the mean value and spread of $σ_v$ and $T_K$ with increasing star formation activity. The data presented in this paper are publicly available. △ Less

Submitted 20 April, 2017; originally announced April 2017.

Comments: 33 pages, 27 figures, accepted to ApJS. Datasets are publicly available: https://dataverse.harvard.edu/dataverse/GAS_DR1

The Herschel Orion Protostar Survey: Luminosity and Envelope Evolution

Authors: William J. Fischer, S. Thomas Megeath, Elise Furlan, Babar Ali, Amelia M. Stutz, John J. Tobin, Mayra Osorio, Thomas Stanke, P. Manoj, Charles A. Poteet, Joseph J. Booker, Lee Hartmann, Thomas L. Wilson, Philip C. Myers, Dan M. Watson

Abstract: The Herschel Orion Protostar Survey obtained well-sampled 1.2 - 870 micron spectral energy distributions (SEDs) of over 300 protostars in the Orion molecular clouds, home to most of the young stellar objects (YSOs) in the nearest 500 pc. We plot the bolometric luminosities and temperatures for 330 Orion YSOs, 315 of which have bolometric temperatures characteristic of protostars. The histogram of… ▽ More The Herschel Orion Protostar Survey obtained well-sampled 1.2 - 870 micron spectral energy distributions (SEDs) of over 300 protostars in the Orion molecular clouds, home to most of the young stellar objects (YSOs) in the nearest 500 pc. We plot the bolometric luminosities and temperatures for 330 Orion YSOs, 315 of which have bolometric temperatures characteristic of protostars. The histogram of bolometric temperature is roughly flat; 29% of the protostars are in Class 0. The median luminosity decreases by a factor of four with increasing bolometric temperature; consequently, the Class 0 protostars are systematically brighter than the Class I protostars, with a median luminosity of 2.3 L_sun as opposed to 0.87 L_sun. At a given bolometric temperature, the scatter in luminosities is three orders of magnitude. Using fits to the SEDs, we analyze how the luminosities corrected for inclination and foreground reddening relate to the mass in the inner 2500 AU of the best-fit model envelopes. The histogram of envelope mass is roughly flat, while the median corrected luminosity peaks at 15 L_sun for young envelopes and falls to 1.7 L_sun for late-stage protostars with remnant envelopes. The spread in luminosity at each envelope mass is three orders of magnitude. Envelope masses that decline exponentially with time explain the flat mass histogram and the decrease in luminosity, while the formation of a range of stellar masses explains the dispersion in luminosity. △ Less

Submitted 19 April, 2017; originally announced April 2017.

Comments: Accepted by ApJ, 16 pages, 14 figures

ALMA Observations of Starless Core Substructure in Ophiuchus

Authors: Helen Kirk, Michael M. Dunham, James Di Francesco, Doug Johnstone, Stella S. R. Offner, Sarah I. Sadavoy, John J. Tobin, Hector G. Arce, Tyler L. Bourke, Steve Mairs, Philip C. Myers, Jaime E. Pineda, Scott Schnee, Yancy L. Shirley

Abstract: Compact substructure is expected to arise in a starless core as mass becomes concentrated in the central region likely to form a protostar. Additionally, multiple peaks may form if fragmentation occurs. We present ALMA Cycle 2 observations of 60 starless and protostellar cores in the Ophiuchus molecular cloud. We detect eight compact substructures which are >15 arcsec from the nearest Spitzer YSO.… ▽ More Compact substructure is expected to arise in a starless core as mass becomes concentrated in the central region likely to form a protostar. Additionally, multiple peaks may form if fragmentation occurs. We present ALMA Cycle 2 observations of 60 starless and protostellar cores in the Ophiuchus molecular cloud. We detect eight compact substructures which are >15 arcsec from the nearest Spitzer YSO. Only one of these has strong evidence for being truly starless after considering ancillary data, e.g., from Herschel and X-ray telescopes. An additional extended emission structure has tentative evidence for starlessness. The number of our detections is consistent with estimates from a combination of synthetic observations of numerical simulations and analytical arguments. This result suggests that a similar ALMA study in the Chamaeleon I cloud, which detected no compact substructure in starless cores, may be due to the peculiar evolutionary state of cores in that cloud. △ Less

Submitted 1 March, 2017; originally announced March 2017.

Comments: accepted for publication in ApJ. 34 pages

Star-forming filament models

Authors: Philip C. Myers

Abstract: New models of star-forming filamentary clouds are presented, to quantify their properties and to predict their evolution. These 2D axisymmetric models describe filaments having no core, one low-mass core, and one cluster-forming core. They are based on Plummer-like cylinders and spheroids, bounded by a constant-density surface of finite extent. In contrast to 1D Plummer-like models, they have spec… ▽ More New models of star-forming filamentary clouds are presented, to quantify their properties and to predict their evolution. These 2D axisymmetric models describe filaments having no core, one low-mass core, and one cluster-forming core. They are based on Plummer-like cylinders and spheroids, bounded by a constant-density surface of finite extent. In contrast to 1D Plummer-like models, they have specific values of length and mass, they approximate observed column density maps, and their distributions of column density (N-pdfs) are pole-free. Each model can estimate the star-forming potential of a core-filament system, by identifying the zone of gas dense enough to form low-mass stars, and by counting the number of enclosed thermal Jeans masses. This analysis suggests that the Musca Center filament may be near the start of its star-forming life, with enough dense gas to make its first ~3 protostars, while the Coronet filament is near the midpoint of its star formation, with enough dense gas to add ~ 8 protostars to its ~20 known stars. In contrast L43 appears near the end of its star-forming life, since it lacks enough dense gas to add any new protostars to the 2 YSOs already known. △ Less

Submitted 8 February, 2017; v1 submitted 7 February, 2017; originally announced February 2017.

Comments: Accepted by The Astrophysical Journal