Showing 1–47 of 47 results for author: Friesen, R K

The stability of dense cores near the Serpens South protocluster

Authors: Rachel K. Friesen, Emma Jarvis

Abstract: Most stars form in clusters and groups rather than in isolation. We present $\lesssim 5^{\prime\prime}$ angular resolution ($\sim 2000$ au, or 0.01 pc) Very Large Array NH$_3$ (1,1), (2,2), and (3,3) and 1.3 cm continuum emission observations of the dense gas within the Serpens South protocluster and extended filaments to the north and south. We identify 94 dense cores using a dendrogram analysis… ▽ More Most stars form in clusters and groups rather than in isolation. We present $\lesssim 5^{\prime\prime}$ angular resolution ($\sim 2000$ au, or 0.01 pc) Very Large Array NH$_3$ (1,1), (2,2), and (3,3) and 1.3 cm continuum emission observations of the dense gas within the Serpens South protocluster and extended filaments to the north and south. We identify 94 dense cores using a dendrogram analysis of the NH$_3$ (1,1) integrated intensity. Gas temperatures $T_K$ and non-thermal linewidths $σ_\mathrm{NT}$ both increase towards the centre of the young stellar cluster, in the dense gas generally and in the cores specifically. We find that most cores (54\%) are super-virial, with gravitationally bound cores located primarily in the filaments. Cores in the protocluster have higher virial parameters by a factor $\sim 1.7$, driven primarily by the increased core $σ_\mathrm{NT}$ values. These cores cannot collapse to form stars unless they accrete additional mass or their core internal motions are reduced. The southern filament shows a significant velocity gradient previously interpreted as mass flow toward the cluster. We find more complex kinematics in the northern filament. We find a strong correlation between $σ_\mathrm{NT}$ and $T_K$, and argue that the enhanced temperatures and non-thermal motions are due to mechanical heating and interaction between the protocluster-driven outflows and the dense gas. Filament-led accretion may also contribute to the increased $σ_\mathrm{NT}$ values. Assuming a constant fraction of core mass ends up in the young stars, future star formation in the Serpens South protocluster will shift to higher masses by a factor $\sim 2$. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 27 pages, 12 figures, accepted to ApJ

The Reservoir of the Per-emb-2 Streamer

Authors: Kotomi Taniguchi, Jaime E Pineda, Paola Caselli, Tomomi Shimoikura, Rachel K. Friesen, Dominique M. Segura-Cox, Anika Schmiedeke

Abstract: Streamers bring gas from outer regions to protostellar systems and could change the chemical composition around protostars and protoplanetary disks. We have carried out mapping observations of carbon-chain species (HC$_3$N, HC$_5$N, CCH, and CCS) in the 3mm and 7mm bands toward the streamer flowing to the Class 0 young stellar object (YSO) Per-emb-2 with the Nobeyama 45m radio telescope. A region… ▽ More Streamers bring gas from outer regions to protostellar systems and could change the chemical composition around protostars and protoplanetary disks. We have carried out mapping observations of carbon-chain species (HC$_3$N, HC$_5$N, CCH, and CCS) in the 3mm and 7mm bands toward the streamer flowing to the Class 0 young stellar object (YSO) Per-emb-2 with the Nobeyama 45m radio telescope. A region with a diameter of $\sim0.04$ pc is located north with a distance of $\sim 20,500$ au from the YSO. The streamer connects to this north region which is the origin of the streamer. The reservoir has high density and low temperature ($n_{\rm {H}_2} \approx 1.9 \times 10^4$ cm$^{-3}$, $T_{\rm {kin}} = 10$ K), which are similar to those of early stage starless cores. By comparisons with the observed abundance ratios of CCS/HC$_3$N to the chemical simulations, the reservoir and streamer are found to be chemically young. The total mass available for the streamer is derived to be $24-34$ M$_{\odot}$. If all of the gas in the reservoir will accrete onto the Per-emb-2 protostellar system, the lifetime of the streamer has been estimated at ($1.1 - 3.2$)$\times10^{5}$ yr, suggesting that the mass accretion via the streamer would continue until the end of the Class I stage. △ Less

Submitted 7 March, 2024; v1 submitted 29 February, 2024; originally announced February 2024.

Comments: Accepted for The Astrophysical Journal (ApJ), 18 pages, 10 figures, 6 tables

Deuterium fractionation in cold dense cores in the low-mass star forming region L1688

Authors: I. V. Petrashkevich, A. F. Punanova, P. Caselli, O. Sipilä, J. E. Pineda, R. K. Friesen, M. G. Korotaeva, A. I. Vasyunin

Abstract: In this work, we study deuterium fractionation in four starless cores in the low-mass star-forming region L1688 in the Ophiuchus molecular cloud. We study how the deuterium fraction ($R_D$) changes with environment, compare deuteration of ions and neutrals, core centre and its envelope, and attempt to reproduce the observed results with a gas-grain chemical model. We chose high and low gas density… ▽ More In this work, we study deuterium fractionation in four starless cores in the low-mass star-forming region L1688 in the Ophiuchus molecular cloud. We study how the deuterium fraction ($R_D$) changes with environment, compare deuteration of ions and neutrals, core centre and its envelope, and attempt to reproduce the observed results with a gas-grain chemical model. We chose high and low gas density tracers to study both core centre and the envelope. With the IRAM 30m antenna, we mapped N$_2$H$^+$(1-0), N$_2$D$^+$(1-0), H$^{13}$CO$^+$ (1-0) and (2-1), DCO$^+$(2-1), and $p$-NH$_2$D(1$_{11}$-1$_{01}$) towards the chosen cores. The missing $p$-NH$_3$ and N$_2$H$^+$(1-0) data were taken from the literature. To measure the molecular hydrogen column density, dust and gas temperature within the cores, we used the Herschel/SPIRE dust continuum emission data, the GAS survey data (ammonia), and the COMPLETE survey data to estimate the upper limit on CO depletion. We present the deuterium fraction maps for three species towards four starless cores. Deuterium fraction of the core envelopes traced by DCO$^+$/H$^{13}$CO$^+$ is one order of magnitude lower ($\sim$0.08) than that of the core central parts traced by the nitrogen-bearing species ($\sim$0.5). Deuterium fraction increases with the gas density as indicated by high deuterium fraction of high gas density tracers and low deuterium fraction of lower gas density tracers and by the decrease of $R_D$ with core radii, consistent with the predictions of the chemical model. Our model results show a good agreement with observations for $R_D$(N$_2$D$^+$/N$_2$H$^+$) and R$_D$(DCO$^+$/HCO$^+$) and underestimate the $R_D$(NH$_2$D/NH$_3$). △ Less

Submitted 10 January, 2024; v1 submitted 9 January, 2024; originally announced January 2024.

Comments: 30 pages with 22 figures, accepted for publication in MNRAS

MSC Class: 70F15 (Primary) ACM Class: J.2

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

Turbulence and Accretion: a High-resolution Study of the B5 Filaments

Authors: Michael Chun-Yuan Chen, James Di Francesco, Jaime E. Pineda, Stella S. Offner, Rachel K. Friesen

Abstract: High-resolution observations of the Perseus B5 "core" have previously revealed that this subsonic region actually consists of several filaments that are likely in the process of forming a quadruple stellar system. Since subsonic filaments are thought to be produced at the $\sim 0.1$ pc sonic scale by turbulent compression, a detailed kinematic study is crucial to test such a scenario in the contex… ▽ More High-resolution observations of the Perseus B5 "core" have previously revealed that this subsonic region actually consists of several filaments that are likely in the process of forming a quadruple stellar system. Since subsonic filaments are thought to be produced at the $\sim 0.1$ pc sonic scale by turbulent compression, a detailed kinematic study is crucial to test such a scenario in the context of core and star formation. Here we present a detailed kinematic follow-up study of the B5 filaments at a 0.009 pc resolution using the VLA and GBT combined observations fitted with multi-component spectral models. Using precisely identified filament spines, we find a remarkable resemblance between the averaged width profiles of each filament and Plummer-like functions, with filaments possessing FWHM widths of $\sim 0.03$ pc. The velocity dispersion profiles of the filaments also show decreasing trends towards the filament spines. Moreover, the velocity gradient field in B5 appears to be locally well ordered ($\sim 0.04$ pc) but globally complex, with kinematic behaviors suggestive of inhomogeneous turbulent accretion onto filaments and longitudinal flows towards a local overdensity along one of the filaments. △ Less

Submitted 1 September, 2022; originally announced September 2022.

Comments: 27 pages, 12 figures, published in The Astrophysical Journal

Journal ref: The Astrophysical Journal, Volume 935, Issue 1, id.57, 16 pp (2022)

From Bubbles and Filaments to Cores and Disks: Gas Gathering and Growth of Structure Leading to the Formation of Stellar Systems

Authors: Jaime E. Pineda, Doris Arzoumanian, Philippe André, Rachel K. Friesen, Annie Zavagno, Seamus D. Clarke, Tsuyoshi Inoue, Che-Yu Chen, Yueh-Ning Lee, Juan D. Soler, Michael Kuffmeier

Abstract: The study of the development of structures on multiple scales in the cold interstellar medium has experienced rapid expansion in the past decade, on both the observational and the theoretical front. Spectral line studies at (sub-)millimeter wavelengths over a wide range of physical scales have provided unique probes of the kinematics of dense gas in star-forming regions, and have been complemented… ▽ More The study of the development of structures on multiple scales in the cold interstellar medium has experienced rapid expansion in the past decade, on both the observational and the theoretical front. Spectral line studies at (sub-)millimeter wavelengths over a wide range of physical scales have provided unique probes of the kinematics of dense gas in star-forming regions, and have been complemented by extensive, high dynamic range dust continuum surveys of the column density structure of molecular cloud complexes, while dust polarization maps have highlighted the role of magnetic fields. This has been accompanied by increasingly sophisticated numerical simulations including new physics (e.g., supernova driving, cosmic rays, non-ideal magneto-hydrodynamics, radiation pressure) and new techniques such as zoom-in simulations allowing multi-scale studies. Taken together, these new data have emphasized the anisotropic growth of dense structures on all scales, from giant ISM bubbles driven by stellar feedback on $\sim$50-100 pc scales through parsec-scale molecular filaments down to $<$0.1 pc dense cores and $<$1000 au protostellar disks. Combining observations and theory, we present a coherent picture for the formation and evolution of these structures and synthesize a comprehensive physical scenario for the initial conditions and early stages of star and disk formation. △ Less

Submitted 23 March, 2023; v1 submitted 8 May, 2022; originally announced May 2022.

Comments: 42 pages, 20 figures. To appear in Protostars and Planets VII, Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura. Small changes to match proofs and the updated format

An Interferometric View of H-MM1. I. Direct Observation of NH3 Depletion

Authors: Jaime E. Pineda, Jorma Harju, Paola Caselli, Olli Sipilä, Mika Juvela, Charlotte Vastel, Erik Rosolowsky, Andreas Burkert, Rachel K. Friesen, Yancy Shirley, María José Maureira, Spandan Choudhury, Dominique M. Segura-Cox, Rolf Güsten, Anna Punanova, Luca Bizzocchi, Alyssa A. Goodman

Abstract: Spectral lines of ammonia, NH$_3$, are useful probes of the physical conditions in dense molecular cloud cores. In addition to advantages in spectroscopy, ammonia has also been suggested to be resistant to freezing onto grain surfaces, which should make it a superior tool for studying the interior parts of cold, dense cores. Here we present high-resolution NH$_3$ observations with the Very Large A… ▽ More Spectral lines of ammonia, NH$_3$, are useful probes of the physical conditions in dense molecular cloud cores. In addition to advantages in spectroscopy, ammonia has also been suggested to be resistant to freezing onto grain surfaces, which should make it a superior tool for studying the interior parts of cold, dense cores. Here we present high-resolution NH$_3$ observations with the Very Large Array (VLA) and Green Bank Telescope (GBT) towards a prestellar core. These observations show an outer region with a fractional NH$_3$ abundance of X(NH$_3$) = (1.975$\pm$0.005)$\times 10^{-8}$ ($\pm 10\%$ systematic), but it also reveals that after all, the X(NH$_3$) starts to decrease above a H$_2$ column density of $\approx 2.6 \times 10^{22}$ cm$^{-2}$. We derive a density model for the core and find that the break-point in the fractional abundance occurs at the density n(H$_2$) $\sim 2\times10^5$ cm$^{-3}$, and beyond this point the fractional abundance decreases with increasing density, following the power law $n^{-1.1}$. This power-law behavior is well reproduced by chemical models where adsorption onto grains dominates the removal of ammonia and related species from the gas at high densities. We suggest that the break-point density changes from core to core depending on the temperature and the grain properties, but that the depletion power law is anyway likely to be close to $n^{-1}$ owing to the dominance of accretion in the central parts of starless cores. △ Less

Submitted 2 May, 2022; originally announced May 2022.

Comments: Accepted for publication in AJ. 18 Pages, 16 Figures, 3 Tables

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

The Core Mass Function in the Orion Nebula Cluster Region: What Determines the Final Stellar Masses?

Authors: Hideaki Takemura, Fumitaka Nakamura, Shuo Kong, Héctor G. Arce, John M. Carpenter, Volker Ossenkopf-Okada, Ralf Klessen, Patricio Sanhueza, Yoshito Shimajiri, Takashi Tsukagoshi, Ryohei Kawabe, Shun Ishii, Kazuhito Dobashi, Tomomi Shimoikura, Paul F. Goldsmith, Álvaro Sánchez-Monge, Jens Kauffmann, Thushara Pillai, Paolo Padoan, Adam Ginsberg, Rowan J. Smith, John Bally, Steve Mairs, Jaime E. Pineda, Dariusz C. Lis , et al. (7 additional authors not shown)

Abstract: Applying dendrogram analysis to the CARMA-NRO C$^{18}$O ($J$=1--0) data having an angular resolution of $\sim$ 8", we identified 692 dense cores in the Orion Nebula Cluster (ONC) region. Using this core sample, we compare the core and initial stellar mass functions in the same area to quantify the step from cores to stars. About 22 \% of the identified cores are gravitationally bound. The derived… ▽ More Applying dendrogram analysis to the CARMA-NRO C$^{18}$O ($J$=1--0) data having an angular resolution of $\sim$ 8", we identified 692 dense cores in the Orion Nebula Cluster (ONC) region. Using this core sample, we compare the core and initial stellar mass functions in the same area to quantify the step from cores to stars. About 22 \% of the identified cores are gravitationally bound. The derived core mass function (CMF) for starless cores has a slope similar to Salpeter's stellar initial mass function (IMF) for the mass range above 1 $M_\odot$, consistent with previous studies. Our CMF has a peak at a subsolar mass of $\sim$ 0.1 $M_\odot$, which is comparable to the peak mass of the IMF derived in the same area. We also find that the current star formation rate is consistent with the picture in which stars are born only from self-gravitating starless cores. However, the cores must gain additional gas from the surroundings to reproduce the current IMF (e.g., its slope and peak mass), because the core mass cannot be accreted onto the star with a 100\% efficiency. Thus, the mass accretion from the surroundings may play a crucial role in determining the final stellar masses of stars. △ Less

Submitted 25 February, 2021; originally announced March 2021.

Comments: 13 pages, 5 figures, accepted by ApJL

Transition from Coherent Cores to Surrounding Cloud in L1688

Authors: Spandan Choudhury, Jaime E. Pineda, Paola Caselli, Stella S. R. Offner, Erik Rosolowsky, Rachel K. Friesen, Elena Redaelli, Ana Chacón-Tanarro, Yancy Shirley, Anna Punanova, Helen Kirk

Abstract: Stars form in cold dense cores showing subsonic velocity dispersions. The parental molecular clouds display higher temperatures and supersonic velocity dispersions. The transition from core to cloud has been observed in velocity dispersion, but temperature and abundance variations are unknown. We aim to study the transition from cores to ambient cloud in temperature and velocity dispersion using a… ▽ More Stars form in cold dense cores showing subsonic velocity dispersions. The parental molecular clouds display higher temperatures and supersonic velocity dispersions. The transition from core to cloud has been observed in velocity dispersion, but temperature and abundance variations are unknown. We aim to study the transition from cores to ambient cloud in temperature and velocity dispersion using a single tracer. We use NH3 (1,1) and (2,2) maps in L1688 from the Green Bank Ammonia Survey, smoothed to 1', and determine the physical properties from fits. We identify the coherent cores and study the changes in temperature and velocity dispersion from cores to the surrounding cloud. We obtain a kinetic temperature map tracing the extended cloud, improving from previous maps tracing mostly the cores. The cloud is 4-6 K warmer than the cores, and shows a larger velocity dispersion (diff. = 0.15-0.25 km/s). Comparing to Herschel-based measurements, we find that cores show kinetic temperature $\approx$1.8 K lower than the dust temperature; while the gas temperature is higher than the dust temperature in the cloud. We find an average p-NH3 fractional abundance (with respect to H2) of $(4.2\pm0.2) \times 10^{-9}$ towards the coherent cores, and $(1.4\pm0.1) \times 10^{-9}$ outside the core boundaries. Using stacked spectra, we detect two components, one narrow and one broad, towards cores and their neighbourhoods. We find the turbulence in the narrow component to be correlated to the size of the structure (Pearson-r=0.54). With these unresolved regional measurements, we obtain a turbulence-size relation of $σ_{v,NT}\propto r^{0.5}$, similar to previous findings using multiple tracers. We discover that the subsonic component extends up to 0.15 pc beyond the typical coherent boundaries, unveiling larger extents of the coherent cores and showing gradual transition to coherence over ~0.2 pc. △ Less

Submitted 12 February, 2021; originally announced February 2021.

Comments: 37 pages, 33 figures, 1 table. Accepted for publication in A&A

Journal ref: A&A 648, A114 (2021)

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)

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

Velocity-coherent Filaments in NGC 1333: Evidence for Accretion Flow?

Authors: Michael Chun-Yuan Chen, James Di Francesco, Erik Rosolowsky, Jared Keown, Jaime E. Pineda, Rachel K. Friesen, Paola Caselli, How-Huan Chen, Christopher D. Matzner, Stella S. Offner, Anna Punanova, Elena Redaelli, Samantha Scibelli, Yancy Shirley

Abstract: Recent observations of global velocity gradients across and along molecular filaments have been interpreted as signs of gas accreting onto and along these filaments, potentially feeding star-forming cores and proto-clusters. The behavior of velocity gradients in filaments, however, has not been studied in detail, particularly on small scales (< 0.1 pc). In this paper, we present MUFASA, an efficie… ▽ More Recent observations of global velocity gradients across and along molecular filaments have been interpreted as signs of gas accreting onto and along these filaments, potentially feeding star-forming cores and proto-clusters. The behavior of velocity gradients in filaments, however, has not been studied in detail, particularly on small scales (< 0.1 pc). In this paper, we present MUFASA, an efficient, robust, and automatic method to fit ammonia lines with multiple velocity components, generalizable to other molecular species. We also present CRISPy, a Python package to identify filament spines in 3D images (e.g., position-position-velocity cubes), along with a complementary technique to sort fitted velocity components into velocity-coherent filaments. In NGC 1333, we find a wealth of velocity gradient structures on a beam-resolved scale of ~0.05 pc. Interestingly, these local velocity gradients are not randomly oriented with respect to filament spines and their perpendicular, i.e., radial, component decreases in magnitude towards the spine for many filaments. Together with remarkably constant velocity gradients on larger scales along many filaments, these results suggest a scenario in which gas falling onto filaments is progressively damped and redirected to flow along these filaments. △ Less

Submitted 26 February, 2020; originally announced February 2020.

Comments: Accepted to ApJ

Journal ref: ApJ 891 (2020) 84

Droplets II: Internal Velocity Structures and Potential Rotational Motions in Pressure-dominated Coherent Structures

Authors: Hope How-Huan Chen, Jaime E. Pineda, Stella S. R. Offner, Alyssa A. Goodman, Andreas Burkert, Rachel K. Friesen, Erik Rosolowsky, Samantha Scibelli, Yancy Shirley

Abstract: We present an analysis of the internal velocity structures of the newly identified sub-0.1 pc coherent structures, droplets, in L1688 and B18. By fitting 2D linear velocity fields to the observed maps of velocity centroids, we determine the magnitudes of linear velocity gradients and examine the potential rotational motions that could lead to the observed velocity gradients. The results show that… ▽ More We present an analysis of the internal velocity structures of the newly identified sub-0.1 pc coherent structures, droplets, in L1688 and B18. By fitting 2D linear velocity fields to the observed maps of velocity centroids, we determine the magnitudes of linear velocity gradients and examine the potential rotational motions that could lead to the observed velocity gradients. The results show that the droplets follow the same power-law relation between the velocity gradient and size found for larger-scale dense cores. Assuming that rotational motion giving rise to the observed velocity gradient in each core is a solid-body rotation of a rotating body with a uniform density, we derive the "net rotational motions" of the droplets. We find a ratio between rotational and gravitational energies, $β$, of $\sim 0.046$ for the droplets, and when including both droplets and larger-scale dense cores, we find $β\sim 0.039$. We then examine the alignment between the velocity gradient and the major axis of each droplet, using methods adapted from the histogram of relative orientations (HRO) introduced by Soler et al. (2013). We find no definitive correlation between the directions of velocity gradients and the elongations of the cores. Lastly, we discuss physical processes other than rotation that may give rise to the observed velocity field. △ Less

Submitted 18 October, 2019; v1 submitted 12 August, 2019; originally announced August 2019.

Comments: accepted to ApJ in Oct, 2019; in press

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

SOFIA {\it FORCAST} Photometry of 12 Extended Green Objects in the Milky Way

Authors: Allison P. M. Towner, Crystal L. Brogan, Todd R. Hunter, Claudia J. Cyganowski, Rachel K. Friesen

Abstract: Massive young stellar objects are known to undergo an evolutionary phase in which high mass accretion rates drive strong outflows. A class of objects believed to trace this phase accurately is the GLIMPSE Extended Green Object (EGO) sample, so named for the presence of extended 4.5 $μ$m emission on sizescales of $\sim$0.1 pc in \textit{Spitzer} images. We have been conducting a multi-wavelength ex… ▽ More Massive young stellar objects are known to undergo an evolutionary phase in which high mass accretion rates drive strong outflows. A class of objects believed to trace this phase accurately is the GLIMPSE Extended Green Object (EGO) sample, so named for the presence of extended 4.5 $μ$m emission on sizescales of $\sim$0.1 pc in \textit{Spitzer} images. We have been conducting a multi-wavelength examination of a sample of 12 EGOs with distances of 1 to 5 kpc. In this paper, we present mid-infrared images and photometry of these EGOs obtained with the SOFIA telescope, and subsequently construct SEDs for these sources from the near-IR to sub-millimeter regimes using additional archival data. We compare the results from greybody models and several publicly-available software packages which produce model SEDs in the context of a single massive protostar. The models yield typical \rstar\/$\sim$10 \rsun, \tstar\/$\sim$10$^3$ to 10$^4$ K, and \lstar\/$\sim$1~$-$~40~$\times$~10$^3$ \lsun; the median $L/M$ for our sample is 24.7 \lsun/\msun. Model results rarely converge for \rstar\/ and \tstar, but do for \lstar, which we take to be an indication of the multiplicity and inherently clustered nature of these sources even though, typically, only a single source dominates in the mid-infrared. The median $L/M$ value for the sample suggests that these objects may be in a transitional stage between the commonly described ``IR-quiet'' and ``IR-bright'' stages of MYSO evolution. The median $T_{dust}$ for the sample is less conclusive, but suggests that these objects are either in this transitional stage or occupy the cooler (and presumably younger) part of the IR-bright stage. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: 43 pages, 18 Figures, 10 Tables, 2 Appendices

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

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

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

Submitted 20 February, 2019; originally announced February 2019.

Comments: 22 pages, 12 figures

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

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

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

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

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

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

ALMA detections of the youngest protostars in Ophiuchus

Authors: Rachel K. Friesen, Andy Pon, Tyler L. Bourke, P. Caselli, J. Di Francesco, Jes K. Jørgensen, Jaime E. Pineda

Abstract: We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.1 mm dust continuum and CO 2-1 emission toward six dense cores within the Ophiuchus molecular cloud. We detect compact, sub-arcsecond continuum structures toward three targets, two of which (Oph A N6 and SM1) are located in the Ophiuchus A ridge. Two targets, SM1 and GSS 30, contain two compact sources within the ALMA… ▽ More We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.1 mm dust continuum and CO 2-1 emission toward six dense cores within the Ophiuchus molecular cloud. We detect compact, sub-arcsecond continuum structures toward three targets, two of which (Oph A N6 and SM1) are located in the Ophiuchus A ridge. Two targets, SM1 and GSS 30, contain two compact sources within the ALMA primary beam. We argue that several of the compact structures are small ($R \lesssim 80$ au) accretion disks around young protostars, due to their resolved, elongated structures, coincident radio and x-ray detections, or bipolar outflow detections. While CO line wings extend to $\pm 10-20$ km s$^{-1}$ for the more evolved sources GSS 30 IRS3 and IRS1, CO emission toward other sources, where detected, only extends a few km s$^{-1}$ from the cloud $v_\mathrm{LSR}$. The dust spectral index toward the compact objects suggests that the disks are either optically thick at 1.1 mm, or that significant grain growth has already occurred. We identify, for the first time, a single compact continuum source ($R \sim 100$ au) toward N6 embedded within a larger continuum structure. SM1N is extended in the continuum but is highly centrally concentrated, with a density profile that follows a $r^{-1.3}$ power law within 200 au, and additional structure suggested by the uv-data. Both N6 and SM1N show no clear bipolar outflows with velocities greater than a few km s$^{-1}$ from the cloud velocity. These sources are candidates to be the youngest protostars or first hydrostatic cores in the Ophiuchus molecular cloud. △ Less

Submitted 9 November, 2018; originally announced November 2018.

Comments: 19 pages, 10 figures, accepted to ApJ

Science with an ngVLA: Deuteration in starless and prestellar cores

Authors: Rachel K. Friesen, Maria T. Beltrán, Paola Caselli, Robin T. Garrod

Abstract: In dense starless and protostellar cores, the relative abundance of deuterated species to their non-deuterated counterparts can become orders of magnitude greater than in the local interstellar medium. This enhancement proceeds through multiple pathways in the gas phase and on dust grains, where the chemistry is strongly dependent on the physical conditions. In this Chapter, we discuss how sensiti… ▽ More In dense starless and protostellar cores, the relative abundance of deuterated species to their non-deuterated counterparts can become orders of magnitude greater than in the local interstellar medium. This enhancement proceeds through multiple pathways in the gas phase and on dust grains, where the chemistry is strongly dependent on the physical conditions. In this Chapter, we discuss how sensitive, high resolution observations with the ngVLA of emission from deuterated molecules will trace both the dense gas structure and kinematics on the compact physical scales required to track the gravitational collapse of star-forming cores and the subsequent formation of young protostars and circumstellar accretion regions. Simultaneously, such observations will play a critical role in tracing the chemical history throughout the various phases of star and planet formation. Many low-J transitions of key deuterated species, along with their undeuterated counterparts, lie within the 60-110 GHz frequency window, the lower end of which is largely unavailable with current facilities and instrumentation. The combination of sensitivity and angular resolution provided only by the ngVLA will enable unparalleled detailed studies of the physics and chemistry of the earliest stages of star formation. △ Less

Submitted 16 October, 2018; originally announced October 2018.

Comments: To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA)

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

The extraordinary outburst in the massive protostellar system NGC6334I-MM1: Flaring of the water masers in a north-south bipolar outflow driven by MM1B

Authors: C. L. Brogan, T. R. Hunter, C. J. Cyganowski, J. O. Chibueze, R. K. Friesen, T. Hirota, G. C. MacLeod, B. A. McGuire, A. M. Sobolev

Abstract: We compare multi-epoch sub-arcsecond VLA imaging of the 22 GHz water masers toward the massive protocluster NGC6334I observed before and after the recent outburst of MM1B in (sub)millimeter continuum. Since the outburst, the water maser emission toward MM1 has substantially weakened. Simultaneously, the strong water masers associated with the synchrotron continuum point source CM2 have flared by a… ▽ More We compare multi-epoch sub-arcsecond VLA imaging of the 22 GHz water masers toward the massive protocluster NGC6334I observed before and after the recent outburst of MM1B in (sub)millimeter continuum. Since the outburst, the water maser emission toward MM1 has substantially weakened. Simultaneously, the strong water masers associated with the synchrotron continuum point source CM2 have flared by a mean factor of 6.5 (to 4.2 kJy) with highly-blueshifted features (up to 70 km/s from LSR) becoming more prominent. The strongest flaring water masers reside 3000 au north of MM1B and form a remarkable bow shock pattern whose vertex coincides with CM2 and tail points back to MM1B. Excited OH masers trace a secondary bow shock located ~120 au downstream. ALMA images of CS (6-5) reveal a highly-collimated north-south structure encompassing the flaring masers to the north and the non-flaring masers to the south seen in projection toward the MM3-UCHII region. Proper motions of the southern water masers over 5.3 years indicate a bulk projected motion of 117 km/s southward from MM1B with a dynamical time of 170 yr. We conclude that CM2, the water masers, and many of the excited OH masers trace the interaction of the high velocity bipolar outflow from MM1B with ambient molecular gas. The previously-excavated outflow cavity has apparently allowed the radiative energy of the current outburst to propagate freely until terminating at the northern bow shock where it strengthened the masers. Additionally, water masers have been detected toward MM7 for the first time, and a highly-collimated CS (6-5) outflow has been detected toward MM4. △ Less

Submitted 11 September, 2018; originally announced September 2018.

Comments: 19 pages, 9 figures, accepted for publication by The Astrophysical Journal

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

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

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

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

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

A First Look at BISTRO Observations of The $ρ$ Oph-A core

Authors: Jungmi Kwon, Yasuo Doi, Motohide Tamura, Masafumi Matsumura, Kate Pattle, David Berry, Sarah Sadavoy, Brenda C. Matthews, Derek Ward-Thompson, Tetsuo Hasegawa, Ray S. Furuya, Andy Pon, James Di Francesco, Doris Arzoumanian, Saeko S. Hayashi, Koji S. Kawabata, Takashi Onaka, Minho Choi, Miju Kang, Thiem Hoang, Chang Won Lee, Sang-Sung Lee, Hong-Li Liu, Tie Liu, Shu-Ichiro Inutsuka , et al. (97 additional authors not shown)

Abstract: We present 850 $μ$m imaging polarimetry data of the $ρ$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We iden… ▽ More We present 850 $μ$m imaging polarimetry data of the $ρ$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $ρ$ Oph A core; some of them can be part of a coherent magnetic field structure in the $ρ$ Oph region. The results are consistent with previous observations of the brightest regions of $ρ$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components. △ Less

Submitted 24 April, 2018; originally announced April 2018.

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

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

VLA Survey of Dense Gas in Extended Green Objects: Prevalence of 25 GHz Methanol Masers

Authors: A. P. M. Towner, C. L. Brogan, T. R. Hunter, C. J. Cyganowski, B. A. McGuire, R. Indebetouw, R. K. Friesen, C. J. Chandler

Abstract: We present $\sim1-4"$ resolution Very Large Array (VLA) observations of four CH$_3$OH $J_2-J_1$-$E$ 25~GHz transitions ($J$=3, 5, 8, 10) along with 1.3~cm continuum toward 20 regions of active massive star formation containing Extended Green Objects (EGOs), 14 of which we have previously studied with the VLA in the Class~I 44~GHz and Class~II 6.7~GHz maser lines (Cyganowski et al. 2009). Sixteen r… ▽ More We present $\sim1-4"$ resolution Very Large Array (VLA) observations of four CH$_3$OH $J_2-J_1$-$E$ 25~GHz transitions ($J$=3, 5, 8, 10) along with 1.3~cm continuum toward 20 regions of active massive star formation containing Extended Green Objects (EGOs), 14 of which we have previously studied with the VLA in the Class~I 44~GHz and Class~II 6.7~GHz maser lines (Cyganowski et al. 2009). Sixteen regions are detected in at least one 25~GHz line ($J$=5), with 13 of 16 exhibiting maser emission. In total, we report 34 new sites of CH$_3$OH maser emission and ten new sites of thermal CH$_3$OH emission, significantly increasing the number of 25~GHz Class I CH$_3$OH masers observed at high angular resolution. We identify probable or likely maser counterparts at 44~GHz for all 15 of the 25~GHz masers for which we have complementary data, providing further evidence that these masers trace similar physical conditions despite uncorrelated flux densities. The sites of thermal and maser emission of CH$_3$OH are both predominantly associated with the 4.5 $μ$m emission from the EGO, and the presence of thermal CH$_3$OH emission is accompanied by 1.3~cm continuum emission in 9 out of 10 cases. Of the 19 regions that exhibit 1.3~cm continuum emission, it is associated with the EGO in 16 cases (out of a total of 20 sites), 13 of which are new detections at 1.3~cm. Twelve of the 1.3~cm continuum sources are associated with 6.7~GHz maser emission and likely trace deeply-embedded massive protostars. △ Less

Submitted 23 May, 2017; v1 submitted 18 May, 2017; originally announced May 2017.

First results from BISTRO -- a SCUBA-2 polarimeter survey of the Gould Belt

Authors: Derek Ward-Thompson, Kate Pattle, Pierre Bastien, Ray S. Furuya, Woojin Kwon, Shih-Ping Lai, Keping Qiu, David Berry, Minho Choi, Simon Coudé, James Di Francesco, Thiem Hoang, Erica Franzmann, Per Friberg, Sarah F. Graves, Jane S. Greaves, Martin Houde, Doug Johnstone, Jason M. Kirk, Patrick M. Koch, Jungmi Kwon, Chang Won Lee, Di Li, Brenda C. Matthews, Joseph C. Mottram , et al. (89 additional authors not shown)

Abstract: We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will a… ▽ More We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments. △ Less

Submitted 27 April, 2017; originally announced April 2017.

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

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 fragmentation and stability of hierarchical structure in Serpens South

Authors: R. K. Friesen, T. L. Bourke, J. Di Francesco, R. Gutermuth, P. C. Myers

Abstract: Filamentary structures are ubiquitous in molecular clouds, and have been recently argued to play an important role in regulating the size and mass of embedded clumps through fragmentation and mass accretion. Here, we reveal the dynamical state and fragmentation of filamentary molecular gas associated with the Serpens South protocluster through analysis of wide (~4 x 4 pc) observations of NH3 (1,1)… ▽ More Filamentary structures are ubiquitous in molecular clouds, and have been recently argued to play an important role in regulating the size and mass of embedded clumps through fragmentation and mass accretion. Here, we reveal the dynamical state and fragmentation of filamentary molecular gas associated with the Serpens South protocluster through analysis of wide (~4 x 4 pc) observations of NH3 (1,1) and (2,2) inversion transitions with the Green Bank Telescope. Detailed modeling of the NH3 lines reveals that the kinematics of the cluster and surrounding filaments are complex. We identify hierarchical structure using a dendrogram analysis of the NH3 emission. The distance between neighbour structures that are embedded within the same parent structure is generally greater than expected from a spherical Jeans analysis, and is in better agreement with cylindrical fragmentation models. The NH3 line width-size relation is flat, and average gas motions are sub- or trans-sonic over all physical scales observed. Subsonic regions extend far beyond the typical 0.1 pc scale previously identified in star-forming cores. As a result, we find a strong trend of decreasing virial parameter with increasing structure mass in Serpens South. Extremely low virial parameters on the largest scales probed by our data suggest that the previously observed, ordered magnetic field is insufficient to support the region against collapse, in agreement with large radial infall motions previously measured toward some of the filaments. A more complex magnetic field configuration in the dense gas, however, may be able to support the filaments. △ Less

Submitted 31 October, 2016; originally announced October 2016.

Comments: 20 pages, 12 figures. Accepted to ApJ

Deuteration of ammonia in the starless core Ophiuchus/H-MM1

Authors: Jorma Harju, Fabien Daniel, Olli Sipilä, Paola Caselli, Jaime E. Pineda, Rachel K. Friesen, Anna Punanova, Rolf Güsten, Laurent Wiesenfeld, Philip C. Myers, Alexandre Faure, Pierre Hily-Blant, Claire Rist, Erik Rosolowsky, Stephan Schlemmer, Yancy L. Shirley

Abstract: Ammonia and its deuterated isotopologues probe physical conditions in dense molecular cloud cores. With the aim of testing the current understanding of the spin-state chemistry of these molecules, we observed spectral lines of NH3, NH2D, NHD2, ND3, and N2D+ towards a dense, starless core in Ophiuchus with the APEX, GBT, and IRAM 30-m telescopes. The observations were interpreted using a gas-grain… ▽ More Ammonia and its deuterated isotopologues probe physical conditions in dense molecular cloud cores. With the aim of testing the current understanding of the spin-state chemistry of these molecules, we observed spectral lines of NH3, NH2D, NHD2, ND3, and N2D+ towards a dense, starless core in Ophiuchus with the APEX, GBT, and IRAM 30-m telescopes. The observations were interpreted using a gas-grain chemistry model combined with radiative transfer calculations. The chemistry model distinguishes between the different nuclear spin states of light hydrogen molecules, ammonia, and their deuterated forms. High deuterium fractionation ratios with NH2D/NH3=0.4, NHD2/NH2D=0.2, and ND3/NHD2=0.06 were found in the core. The observed ortho/para ratios of NH2D and NHD2 are close to the corresponding nuclear spin statistical weights. The chemistry model can approximately reproduce the observed abundances, but predicts uniformly too low ortho/para-NH2D, and too large ortho/para-NHD2 ratios. The longevity of N2H+ and NH3 in dense gas, which is prerequisite to their strong deuteration, can be attributed to the chemical inertia of N2 on grain surfaces. The discrepancies between the chemistry model and the observations are likely to be caused by the fact that the model assumes complete scrambling in principal gas-phase deuteration reactions of ammonia, which means that all the nuclei are mixed in reactive collisions. If, instead, these reactions occur through proton hop/hydrogen abstraction processes, statistical spin ratios are to be expected. The present results suggest that while the deuteration of ammonia changes with physical conditions and time, the nuclear spin ratios of ammonia isotopologues do not probe the evolutionary stage of a cloud. △ Less

Submitted 22 December, 2016; v1 submitted 19 April, 2016; originally announced April 2016.

Comments: to appear in Astronomy & Astrophysics

Journal ref: A&A 600, A61 (2017)

Contraction Signatures Toward Dense Cores in the Perseus Molecular Cloud

Authors: J. L. Campbell, R. K. Friesen, P. G. Martin, P. Caselli, J. Kauffmann, J. E. Pineda

Abstract: We report the results of an HCO+ (3-2) and N2D+ (3-2) molecular line survey performed toward 91 dense cores in the Perseus molecular cloud using the James Clerk Maxwell Telescope, to identify the fraction of starless and protostellar cores with systematic radial motions. We quantify the HCO+ asymmetry using a dimensionless asymmetry parameter $δ_v$, and identify 20 cores with significant blue or r… ▽ More We report the results of an HCO+ (3-2) and N2D+ (3-2) molecular line survey performed toward 91 dense cores in the Perseus molecular cloud using the James Clerk Maxwell Telescope, to identify the fraction of starless and protostellar cores with systematic radial motions. We quantify the HCO+ asymmetry using a dimensionless asymmetry parameter $δ_v$, and identify 20 cores with significant blue or red line asymmetries in optically-thick emission indicative of collapsing or expanding motions, respectively. We separately fit the HCO+ profiles with an analytic collapse model and determine contraction (expansion) speeds toward 22 cores. Comparing the $δ_v$ and collapse model results, we find that $δ_v$ is a good tracer of core contraction if the optically-thin emission is aligned with the model-derived systemic velocity. The contraction speeds range from subsonic (0.03 km/s) to supersonic (0.4 km/s), where the supersonic contraction speeds may trace global rather than local core contraction. Most cores have contraction speeds significantly less than their free-fall speeds. Only 7 of 28 starless cores have spectra well-fit by the collapse model, which more than doubles (15 of 28) for protostellar cores. Starless cores with masses greater than the Jeans mass (M/M$_J$ > 1) are somewhat more likely to show contraction motions. We find no trend of optically-thin non-thermal line width with M/M$_J$, suggesting that any undetected contraction motions are small and subsonic. Most starless cores in Perseus are either not in a state of collapse or expansion, or are in a very early stage of collapse. △ Less

Submitted 26 January, 2016; originally announced January 2016.

Comments: Accepted to ApJ, 18 pages, 12 figures, 4 tables

Dust Emissivity in the Star-Forming Filament OMC 2/3

Authors: S. I. Sadavoy, A. M. Stutz, S. Schnee, B. S. Mason, J. Di Francesco, R. K. Friesen

Abstract: We present new measurements of the dust emissivity index, beta, for the high-mass, star-forming OMC 2/3 filament. We combine 160-500 um data from Herschel with long-wavelength observations at 2 mm and fit the spectral energy distributions across a ~ 2 pc long, continuous section of OMC 2/3 at 15000 AU (0.08 pc) resolution. With these data, we measure beta and reconstruct simultaneously the filtere… ▽ More We present new measurements of the dust emissivity index, beta, for the high-mass, star-forming OMC 2/3 filament. We combine 160-500 um data from Herschel with long-wavelength observations at 2 mm and fit the spectral energy distributions across a ~ 2 pc long, continuous section of OMC 2/3 at 15000 AU (0.08 pc) resolution. With these data, we measure beta and reconstruct simultaneously the filtered-out large-scale emission at 2 mm. We implement both variable and fixed values of beta, finding that beta = 1.7 - 1.8 provides the best fit across most of OMC 2/3. These beta values are consistent with a similar analysis carried out with filtered Herschel data. Thus, we show that beta values derived from spatial filtered emission maps agree well with those values from unfiltered data at the same resolution. Our results contradict the very low beta values (~ 0.9) previously measured in OMC 2/3 between 1.2 mm and 3.3 mm data, which we attribute to elevated fluxes in the 3.3 mm observations. Therefore, we find no evidence or rapid, extensive dust grain growth in OMC 2/3. Future studies with Herschel data and complementary ground-based long-wavelength data can apply our technique to obtain robust determinations of beta in nearby cold molecular clouds. △ Less

Submitted 16 February, 2016; v1 submitted 25 January, 2016; originally announced January 2016.

Comments: Accepted to A&A, 12 pages, 12 figures, 1 table. Minor changes to the text and additional information about downloading the 2 mm observations

Physical Conditions of the Earliest Phases of Massive Star Formation: Single-Dish and Interferometric Observations of Ammonia and CCS in Infrared Dark Clouds

Authors: William J. Dirienzo, Crystal Brogan, Remy Indebetouw, Claire J. Chandler, Rachel K. Friesen, Kathryn E. Devine

Abstract: Infrared Dark Clouds (IRDCs) harbor the earliest phases of massive star formation, and many of the compact cores in IRDCs, traced by millimeter continuum or by molecular emission in high critical density lines, host massive young stellar objects (YSOs). We used the Robert C. Byrd Green Bank Telescope (GBT) and the Karl G. Jansky Very Large Array (VLA) to map NH$_{3}$ and CCS in nine IRDCs to revea… ▽ More Infrared Dark Clouds (IRDCs) harbor the earliest phases of massive star formation, and many of the compact cores in IRDCs, traced by millimeter continuum or by molecular emission in high critical density lines, host massive young stellar objects (YSOs). We used the Robert C. Byrd Green Bank Telescope (GBT) and the Karl G. Jansky Very Large Array (VLA) to map NH$_{3}$ and CCS in nine IRDCs to reveal the temperature, density, and velocity structures and explore chemical evolution in the dense ($>10^{22}$ cm$^{-2}$) gas. Ammonia is an excellent molecular tracer for these cold, dense environments. The internal structure and kinematics of the IRDCs include velocity gradients, filaments, and possibly colliding clumps that elucidate the formation process of these structures and their YSOs. We find a wide variety of substructure including filaments and globules at distinct velocities, sometimes overlapping at sites of ongoing star formation. It appears that these IRDCs are still being assembled from molecular gas clumps even as star formation has already begun, and at least three of them appear consistent with the morphology of ``hub-filament structures'' discussed in the literature. Furthermore, we find that these clumps are typically near equipartition between gravitational and kinetic energies, so these structures may survive for multiple free-fall times.Keywords: molecular data -- ISM: clouds -- (ISM:) dust, extinction -- ISM: molecules -- Stars: formation -- radio lines: ISM △ Less

Submitted 5 November, 2015; v1 submitted 7 August, 2015; originally announced August 2015.

Comments: 83 pages, 23 figures, published in The Astronomical Journal. Version 2 contains minor revisions to match the published version of the article. References were updated. Figures and tables were moved from the end to where they are first referenced in the article

Revealing H2D+ depletion and compact structure in starless and protostellar cores with ALMA

Authors: R. K. Friesen, J. Di Francesco, T. L. Bourke, P. Caselli, J. K. Jørgensen, J. E. Pineda, M. Wong

Abstract: We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the submillimeter dust continuum and H2D+ 1_{10}-1_{11} emission toward two evolved, potentially protostellar cores within the Ophiuchus molecular cloud, Oph A SM1 and SM1N. The data reveal small-scale condensations within both cores, with mass upper limits of M <~ 0.02M_Sun (~ 20 M_Jup). The SM1 condensation is consist… ▽ More We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the submillimeter dust continuum and H2D+ 1_{10}-1_{11} emission toward two evolved, potentially protostellar cores within the Ophiuchus molecular cloud, Oph A SM1 and SM1N. The data reveal small-scale condensations within both cores, with mass upper limits of M <~ 0.02M_Sun (~ 20 M_Jup). The SM1 condensation is consistent with a nearly-symmetric Gaussian source with a width of only 37 AU. The SM1N condensation is elongated, and extends 500 AU along its major axis. No evidence for substructure is seen in either source. A Jeans analysis indicates these sources are unlikely to fragment, suggesting that both will form single stars. H2D+ is only detected toward SM1N, offset from the continuum peak by ~150-200 AU. This offset may be due to either heating from an undetected, young, low luminosity protostellar source or first hydrostatic core, or HD (and consequently H2D+) depletion in the cold centre of the condensation. We propose that SM1 is protostellar, and that the condensation detected by ALMA is a warm (T ~ 30-50 K) accretion disk. The less concentrated emission of the SM1N condensation suggests that it is still starless, but we cannot rule out the presence of a low-luminosity source, perhaps surrounded by a pseudodisk. These data reveal observationally the earliest stages of the formation of circumstellar accretion regions, and agree with theoretical predictions that disk formation can occur very early in the star formation process, coeval with or just after the formation of a first hydrostatic core or protostar. △ Less

Submitted 14 October, 2014; originally announced October 2014.

Comments: 13 pages, 4 figures, accepted by ApJ

Abundant cyanopolyynes as a probe of infall in the Serpens South cluster-forming region

Authors: R. K. Friesen, L. Medeiros, S. Schnee, T. L. Bourke, J. Di Francesco, R. Gutermuth, P. C. Myers

Abstract: We have detected bright HC7N J = 21-20 emission toward multiple locations in the Serpens South cluster-forming region using the K-Band Focal Plane Array at the Robert C. Byrd Green Bank Telescope. HC7N is seen primarily toward cold filamentary structures that have yet to form stars, largely avoiding the dense gas associated with small protostellar groups and the main central cluster of Serpens Sou… ▽ More We have detected bright HC7N J = 21-20 emission toward multiple locations in the Serpens South cluster-forming region using the K-Band Focal Plane Array at the Robert C. Byrd Green Bank Telescope. HC7N is seen primarily toward cold filamentary structures that have yet to form stars, largely avoiding the dense gas associated with small protostellar groups and the main central cluster of Serpens South. Where detected, the HC7N abundances are similar to those found in other nearby star forming regions. Toward some HC7N `clumps', we find consistent variations in the line centroids relative to NH3 (1,1) emission, as well as systematic increases in the HC7N non-thermal line widths, which we argue reveal infall motions onto dense filaments within Serpens South with minimum mass accretion rates of M ~ 2-5 M_sun Myr^-1. The relative abundance of NH3 to HC7N suggests that the HC7N is tracing gas that has been at densities n ~ 10^4 cm^-3, for timescales t < 1-2 x 10^5 yr. Since HC7N emission peaks are rarely co-located with those of either NH3 or continuum, it is likely that Serpens South is not particularly remarkable in its abundance of HC7N, but instead the serendipitous mapping of HC7N simultaneously with NH3 has allowed us to detect HC7N at low abundances in regions where it otherwise may not have been looked for. This result extends the known star-forming regions containing significant HC7N emission from typically quiescent regions, like the Taurus molecular cloud, to more complex, active environments. △ Less

Submitted 6 September, 2013; originally announced September 2013.

Comments: 19 pages, 13 figures, accepted to MNRAS. Version with full resolution figures available at http://www.dunlap.utoronto.ca/~friesen/Friesen_HC7N.pdf

The Protocluster G18.67+0.03: A Test Case for Class I Methanol Masers as Evolutionary Indicators for Massive Star Formation

Authors: C. J. Cyganowski, C. L. Brogan, T. R. Hunter, Q. Zhang, R. K. Friesen, R. Indebetouw, C. J. Chandler

Abstract: We present high angular resolution Submillimeter Array (SMA) and Karl G. Jansky Very Large Array (VLA) observations of the massive protocluster G18.67+0.03. Previously targeted in maser surveys of GLIMPSE Extended Green Objects (EGOs), this cluster contains three Class I methanol maser sources, providing a unique opportunity to test the proposed role of Class I masers as evolutionary indicators fo… ▽ More We present high angular resolution Submillimeter Array (SMA) and Karl G. Jansky Very Large Array (VLA) observations of the massive protocluster G18.67+0.03. Previously targeted in maser surveys of GLIMPSE Extended Green Objects (EGOs), this cluster contains three Class I methanol maser sources, providing a unique opportunity to test the proposed role of Class I masers as evolutionary indicators for massive star formation. The millimeter observations reveal bipolar molecular outflows, traced by 13CO(2-1) emission, associated with all three Class I maser sources. Two of these sources (including the EGO) are also associated with 6.7 GHz Class II methanol masers; the Class II masers are coincident with millimeter continuum cores that exhibit hot core line emission and drive active outflows, as indicated by the detection of SiO(5-4). In these cases, the Class I masers are coincident with outflow lobes, and appear as clear cases of excitation by active outflows. In contrast, the third Class I source is associated with an ultracompact HII region, and not with Class II masers. The lack of SiO emission suggests the 13CO outflow is a relic, consistent with its longer dynamical timescale. Our data show that massive young stellar objects associated only with Class I masers are not necessarily young, and provide the first unambiguous evidence that Class I masers may be excited by both young (hot core) and older (UC HII) MYSOs within the same protocluster. △ Less

Submitted 19 October, 2012; v1 submitted 11 October, 2012; originally announced October 2012.

Comments: Astrophysical Journal Letters, accepted. emulateapj, 7 pages including 4 figures and 1 table. Figures compressed. v2: coauthor affiliation updated, emulateapj version

Testing Triggered Star Formation in Six H II Regions

Authors: William J. Dirienzo, Rémy Indebetouw, Crystal Brogan, Claudia J. Cyganowski, Edward B. Churchwell, Rachel K. Friesen

Abstract: We investigated six H II regions with infrared, bright rimmed bubble or cometary morphology, in search of quantitative evidence for triggered star formation, both collect and collapse and radiatively driven implosion. We identified and classified 458 Young Stellar Objects (YSOs) in and around the H II regions. YSOs were determined by fitting a collection of radiative transfer model spectral energy… ▽ More We investigated six H II regions with infrared, bright rimmed bubble or cometary morphology, in search of quantitative evidence for triggered star formation, both collect and collapse and radiatively driven implosion. We identified and classified 458 Young Stellar Objects (YSOs) in and around the H II regions. YSOs were determined by fitting a collection of radiative transfer model spectral energy distributions (SEDs) to infrared photometry for a large sample of point sources. We determined areas where there exist enhanced populations of relatively unevolved YSOs on the bright rims of these regions, suggesting that star formation has been triggered there. We further investigated the physical properties of the regions by using radio continuum emission as a proxy for ionizing flux powering the H II regions, and 13CO (1-0) observations to measure masses and gravitational stability of molecular clumps. We used an analytical model of collect and collapse triggered star formation, as well as a simulation of radiatively driven implosion, and thus we compare the observed properties of the molecular gas with those predicted in the triggering scenarios. Notably, those regions in our sample that show evidence of cometary, or "blister," morphology are more likely to show evidence of triggering. Keywords: HII regions -- ISM: bubbles -- Stars: formation -- Stars: protostars △ Less

Submitted 17 December, 2012; v1 submitted 9 October, 2012; originally announced October 2012.

Comments: [v1] 58 pages, 14 figures, published in the Astronomical Journal; [v2] Title and minor typos corrected to match published version, and references updated

Journal ref: 2012AJ....144..173D

First Results from a 1.3 cm EVLA Survey of Massive Protostellar Objects: G35.03+0.35

Authors: C. L. Brogan, T. R. Hunter, C. J. Cyganowski, R. K. Friesen, C. J. Chandler, R. Indebetouw

Abstract: We have performed a 1.3 centimeter survey of 24 massive young stellar objects (MYSOs) using the Expanded Very Large Array (EVLA). The sources in the sample exhibit a broad range of massive star formation signposts including Infrared Dark Clouds (IRDCs), UCHII regions, and extended 4.5 micron emission in the form of Extended Green Objects (EGOs). In this work, we present results for G35.03+0.35 whi… ▽ More We have performed a 1.3 centimeter survey of 24 massive young stellar objects (MYSOs) using the Expanded Very Large Array (EVLA). The sources in the sample exhibit a broad range of massive star formation signposts including Infrared Dark Clouds (IRDCs), UCHII regions, and extended 4.5 micron emission in the form of Extended Green Objects (EGOs). In this work, we present results for G35.03+0.35 which exhibits all of these phenomena. We simultaneously image the 1.3 cm ammonia (1,1) through (6,6) inversion lines, four methanol transitions, two H recombination lines, plus continuum at 0.05 pc resolution. We find three areas of thermal ammonia emission, two within the EGO (designated the NE and SW cores) and one toward an adjacent IRDC. The NE core contains an UCHII region (CM1) and a candidate HCHII region (CM2). A region of non-thermal, likely masing ammonia (3,3) and (6,6) emission is coincident with an arc of 44 GHz methanol masers. We also detect two new 25 GHz Class I methanol masers. A complementary Submillimeter Array 1.3 mm continuum image shows that the distribution of dust emission is similar to the lower-lying ammonia lines, all peaking to the NW of CM2, indicating the likely presence of an additional MYSO in this protocluster. By modeling the ammonia and 1.3 mm continuum data, we obtain gas temperatures of 20-220 K and masses of 20-130 solar. The diversity of continuum emission properties and gas temperatures suggest that objects in a range of evolutionary states exist concurrently in this protocluster. △ Less

Submitted 5 June, 2011; originally announced June 2011.

Comments: To appear in Astrophysical Journal Letters Special Issue on the EVLA. 16 pages, 3 figures. Includes the complete version of Figure 3, which was unable to fit into the journal article due to the number of panels

Journal ref: The Astrophysical Journal Letters, Volume 739, Issue 1, article id. L16, 6 pp., September 2011

The Initial Conditions of Clustered Star Formation III. The Deuterium Fractionation of the Ophiuchus B2 Core

Authors: R. K. Friesen, J. Di Francesco, P. C. Myers, A. Belloche, Y. L. Shirley, T. L. Bourke, P. André

Abstract: We present N2D+ 3-2 (IRAM) and H2D+ 1_11 - 1_10 and N2H+ 4-3 (JCMT) maps of the small cluster-forming Ophiuchus B2 core in the nearby Ophiuchus molecular cloud. In conjunction with previously published N2H+ 1-0 observations, the N2D+ data reveal the deuterium fractionation in the high density gas across Oph B2. The average deuterium fractionation R_D = N(N2D+)/N(N2H+) ~ 0.03 over Oph B2, with seve… ▽ More We present N2D+ 3-2 (IRAM) and H2D+ 1_11 - 1_10 and N2H+ 4-3 (JCMT) maps of the small cluster-forming Ophiuchus B2 core in the nearby Ophiuchus molecular cloud. In conjunction with previously published N2H+ 1-0 observations, the N2D+ data reveal the deuterium fractionation in the high density gas across Oph B2. The average deuterium fractionation R_D = N(N2D+)/N(N2H+) ~ 0.03 over Oph B2, with several small scale R_D peaks and a maximum R_D = 0.1. The mean R_D is consistent with previous results in isolated starless and protostellar cores. The column density distributions of both H2D+ and N2D+ show no correlation with total H2 column density. We find, however, an anticorrelation in deuterium fractionation with proximity to the embedded protostars in Oph B2 to distances >= 0.04 pc. Destruction mechanisms for deuterated molecules require gas temperatures greater than those previously determined through NH3 observations of Oph B2 to proceed. We present temperatures calculated for the dense core gas through the equating of non-thermal line widths for molecules (i.e., N2D+ and H2D+) expected to trace the same core regions, but the observed complex line structures in B2 preclude finding a reasonable result in many locations. This method may, however, work well in isolated cores with less complicated velocity structures. Finally, we use R_D and the H2D+ column density across Oph B2 to set a lower limit on the ionization fraction across the core, finding a mean x_e, lim >= few x 10^{-8}. Our results show that care must be taken when using deuterated species as a probe of the physical conditions of dense gas in star-forming regions. △ Less

Submitted 23 June, 2010; originally announced June 2010.

Comments: ApJ accepted

The JCMT Legacy Survey of the Gould Belt: a first look at Serpens with HARP

Authors: S. F. Graves, J. S. Richer, J. V. Buckle, A. Duarte-Cabral, G. A. Fuller, M. R. Hogerheijde, J. E. Owen, C. Brunt, H. M. Butner, B. Cavanagh, A. Chrysostomou, E. I. Curtis, C. J. Davis, M. Etxaluze, J. Di Francesco, P. Friberg, R. K. Friesen, J. S. Greaves, J. Hatchell, D. Johnstone, B. Matthews, H. Matthews, C. D. Matzner, D. Nutter, J. M. C. Rawlings , et al. (10 additional authors not shown)

Abstract: The Gould Belt Legacy Survey (GBS) on the JCMT has observed a region of 260 square arcminutes in 12CO J=3--2 emission, and a 190 square arcminute subset of this in 13CO and C18O towards the Serpens molecular cloud. We examine the global velocity structure of the non-outflowing gas, and calculate excitation temperatures and opacities. The large scale mass and energetics of the region are evaluated,… ▽ More The Gould Belt Legacy Survey (GBS) on the JCMT has observed a region of 260 square arcminutes in 12CO J=3--2 emission, and a 190 square arcminute subset of this in 13CO and C18O towards the Serpens molecular cloud. We examine the global velocity structure of the non-outflowing gas, and calculate excitation temperatures and opacities. The large scale mass and energetics of the region are evaluated, with special consideration for high velocity gas. We find the cloud to have a mass of 203 solar masses, and to be gravitationally bound, and that the kinetic energy of the outflowing gas is approximately seventy percent of the turbulent kinetic energy of the cloud. We identify compact outflows towards some of the submillimetre Class 0/I sources in the region △ Less

Submitted 4 June, 2010; originally announced June 2010.

Comments: 17 Pages, accepted by MNRAS

The Initial Conditions of Clustered Star Formation. II. N2H+ Observations of the Ophiuchus B Core

Authors: R. K. Friesen, J. Di Francesco, Y. Shimajiri, S. Takakuwa

Abstract: We present a Nobeyama 45 m Radio Telescope map and Australia Telescope Compact Array pointed observations of N2H+ 1-0 emission towards the clustered, low mass star forming Oph B Core within the Ophiuchus molecular cloud. We compare these data with previously published results of high resolution NH3 (1,1) and (2,2) observations in Oph B. We use 3D Clumpfind to identify emission features in the si… ▽ More We present a Nobeyama 45 m Radio Telescope map and Australia Telescope Compact Array pointed observations of N2H+ 1-0 emission towards the clustered, low mass star forming Oph B Core within the Ophiuchus molecular cloud. We compare these data with previously published results of high resolution NH3 (1,1) and (2,2) observations in Oph B. We use 3D Clumpfind to identify emission features in the single-dish N2H+ map, and find that the N2H+ `clumps' match well similar features previously identified in NH3 (1,1) emission, but are frequently offset to clumps identified at similar resolution in 850 micron continuum emission. Wide line widths in the Oph B2 sub-Core indicate non-thermal motions dominate the Core kinematics, and remain transonic at densities n ~ 3 x 10^5 cm^-3 with large scatter and no trend with N(H2). Non-thermal motions in Oph B1 and B3 are subsonic with little variation, but also show no trend with H2 column density. Over all Oph B, non-thermal N2H+ line widths are substantially narrower than those traced by NH3, making it unlikely NH3 and N2H+ trace the same material, but the v_LSR of both species agree well. We find evidence for accretion in Oph B1 from the surrounding ambient gas. The NH3/N2H+ abundance ratio is larger towards starless Oph B1 than towards protostellar Oph B2, similar to recent observational results in other star-forming regions. Small-scale structure is found in the ATCA N2H+ 1-0 emission, where emission peaks are again offset from continuum emission. In particular, the ~1 M_Sun B2-MM8 clump is associated with a N2H+ emission minimum and surrounded by a broken ring-like N2H+ emission structure, suggestive of N2H+ depletion. We find a strong general trend of decreasing N2H+ abundance with increasing N(H2) in Oph B which matches that found for NH3. △ Less

Submitted 19 November, 2009; originally announced November 2009.

Comments: 55 pages (manuscript), 15 figures, ApJ accepted

Journal ref: Astrophys.J.708:1002-1024,2010

Submillimeter Observations of the Quiescent Core - Ophiuchus A-N6

Authors: A. Pon, R. Plume, R. K. Friesen, J. Di Francesco, B. Matthews, E. A. Bergin

Abstract: We have observed the Oph A-N6 prestellar core in the following transitions: N2D+ J=3 to 2, DCO+ J=3 to 2 and J=5 to 4, HCO+ J=3 to 2, CS J=5 to 4 and J=7 to 6, and H13CO+ J=3 to 2 and J=4 to 3, using the James Clerk Maxwell Telescope. We also observed the NH3 (1,1) and (2,2) inversion transitions towards the Oph A-N6 peak with the Green Bank Telescope. We have found that the N6 core is composed… ▽ More We have observed the Oph A-N6 prestellar core in the following transitions: N2D+ J=3 to 2, DCO+ J=3 to 2 and J=5 to 4, HCO+ J=3 to 2, CS J=5 to 4 and J=7 to 6, and H13CO+ J=3 to 2 and J=4 to 3, using the James Clerk Maxwell Telescope. We also observed the NH3 (1,1) and (2,2) inversion transitions towards the Oph A-N6 peak with the Green Bank Telescope. We have found that the N6 core is composed of shells of different chemical composition due to the freezing out of chemical species at different densities. The undepleted species N2D+ appears to trace the high-density interior of the core, DCO+ and H13CO+ trace an intermediate region, and CS traces the outermost edges of the core. A distinct blue-red spectral asymmetry, indicative of infall motion, is clearly detected in the HCO+ spectra, suggesting that N6 is undergoing gravitational collapse. This collapse was possibly initiated by a decrease in turbulent support suggested by the fact that the non-thermal contribution to the line widths is smaller for the molecular species closer to the center of the core. We also present a temperature profile for the core. These observations support the claim that the Oph A-N6 core is an extremely young prestellar core, which may have been recently cut off from MHD support and begun to collapse. △ Less

Submitted 31 March, 2009; originally announced April 2009.

Comments: 18 pages, 8 figures, accepted by ApJ

Journal ref: Astrophys.J.698:1914-1923,2009

The Initial Conditions of Clustered Star Formation I: NH3 Observations of Dense Cores in Ophiuchus

Authors: R. K. Friesen, J. Di Francesco, Y. L. Shirley, P. C. Myers

Abstract: We present combined interferometer and single dish telescope data of NH3 (J,K) = (1,1) and (2,2) emission towards the clustered star forming Ophiuchus B, C and F Cores at high spatial resolution (~1200 AU) using the Australia Telescope Compact Array, the Very Large Array, and the Green Bank Telescope. While the large scale features of the NH3 (1,1) integrated intensity appear similar to 850 micr… ▽ More We present combined interferometer and single dish telescope data of NH3 (J,K) = (1,1) and (2,2) emission towards the clustered star forming Ophiuchus B, C and F Cores at high spatial resolution (~1200 AU) using the Australia Telescope Compact Array, the Very Large Array, and the Green Bank Telescope. While the large scale features of the NH3 (1,1) integrated intensity appear similar to 850 micron continuum emission maps of the Cores, on 15" (1800 AU) scales we find significant discrepancies between the dense gas tracers in Oph B, but good correspondence in Oph C and F. Using the Clumpfind structure identifying algorithm, we identify 15 NH3 clumps in Oph B, and 3 each in Oph C and F. Only five of the Oph B NH3 clumps are coincident within 30" (3600 AU) of a submillimeter clump. We find v_LSR varies little across any of the Cores, and additionally varies by only ~1.5 km/s between them. The observed NH3 line widths within the Oph B and F Cores are generally large and often mildly supersonic, while Oph C is characterized by narrow line widths which decrease to nearly thermal values. We find several regions of localized narrow line emission (Δv < 0.4 km/s), some of which are associated with NH3 clumps. We derive the kinetic temperatures of the gas, and find they are remarkably constant across Oph B and F, with a warmer mean value (T_K = 15 K) than typically found in isolated regions and consistent with previous results in clustered regions. Oph C, however, has a mean T_K = 12 K, decreasing to a minimum T_K = 9.4 K towards the submillimeter continuum peak, similar to previous studies of isolated starless cores. There is no significant difference in temperature towards protostars embedded in the Cores. [Abridged] △ Less

Submitted 3 March, 2009; originally announced March 2009.

Comments: ApJ, accepted

Journal ref: Astrophys.J.697:1457-1480,2009

Fourier Transform Spectroscopy of the submillimetre continuum emission from hot molecular cores

Authors: R. K. Friesen, D. Johnstone, D. A. Naylor, G. R. Davis

Abstract: We have used a Fourier Transform Spectrometer on the James Clerk Maxwell Telescope to study the submillimetre continuum emission from dust in three hot molecular cores (HMC). The spectral index beta of the dust emission for these sources has been determined solely within the 30 GHz wide 350 GHz (850 micron) passband to an accuracy comparable to those determined through multi-wavelength observati… ▽ More We have used a Fourier Transform Spectrometer on the James Clerk Maxwell Telescope to study the submillimetre continuum emission from dust in three hot molecular cores (HMC). The spectral index beta of the dust emission for these sources has been determined solely within the 30 GHz wide 350 GHz (850 micron) passband to an accuracy comparable to those determined through multi-wavelength observations. We find an average beta = 1.6, in agreement with spectral indices determined from previous submillimetre observations of these sources and with those determined for HMC in general. The largest single source of uncertainty in these results is the variability of the atmosphere at 350 GHz, and with better sky subtraction techniques we show that the dust spectral index can clearly be determined within one passband to high accuracy with a submillimetre FTS. Using an imaging FTS on SCUBA-2, the next generation wide-field submillimetre camera currently under development to replace SCUBA at the JCMT in 2006, we calculate that at 350 GHz it will be possible to determine beta to +/- 0.1 for sources as faint as 400 mJy/beam and to +/- 0.3 for sources as faint as 140 mJy/beam. △ Less

Submitted 16 May, 2005; originally announced May 2005.

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

Journal ref: Mon.Not.Roy.Astron.Soc. 361 (2005) 460-468