Showing 1–50 of 267 results for author: Evans, N J

CORINOS II. JWST-MIRI detection of warm molecular gas from an embedded, disk-bearing protostar

Authors: Colette Salyk, Yao-Lun Yang, Klaus M. Pontoppidan, Jennifer B. Bergner, Yuki Okoda, Jaeyeong Kim, Neal J. Evans II, Ilsedore Cleeves, Ewine F. van Dishoeck, Robin T. Garrod, Joel D. Green

Abstract: We present James Webb Space Telescope (JWST) Mid-InfraRed Instrument (MIRI) observations of warm CO and H$_2$O gas in emission toward the low-mass protostar IRAS 15398-3359, observed as part of the CORINOS program. The CO is detected via the rovibrational fundamental band and hot band near 5 $μ$m, whereas the H$_2$O is detected in the rovibrational bending mode at 6-8 $μ$m. Rotational analysis ind… ▽ More We present James Webb Space Telescope (JWST) Mid-InfraRed Instrument (MIRI) observations of warm CO and H$_2$O gas in emission toward the low-mass protostar IRAS 15398-3359, observed as part of the CORINOS program. The CO is detected via the rovibrational fundamental band and hot band near 5 $μ$m, whereas the H$_2$O is detected in the rovibrational bending mode at 6-8 $μ$m. Rotational analysis indicates that the CO originates in a hot reservoir of $1551\pm135$ K, while the water is much cooler at $212\pm 2$ K. Neither the CO nor the H$_2$O line images are significantly spatially extended, constraining the emission to within $\sim$40 au of the protostar. The compactness and high temperature of the CO are consistent with an origin in the embedded protostellar disk, or a compact disk wind. In contrast, the water must arise from a cooler region and requires a larger emitting area (compared to CO) to produce the observed fluxes. The water may arise from a more extended part of the disk, or from the inner portion of the outflow cavity. Thus, the origin of the molecular emission observed with JWST remains ambiguous. Better constraints on the overall extinction, comparison with realistic disk models, and future kinematically-resolved observations may all help to pinpoint the true emitting reservoirs. △ Less

Submitted 21 July, 2024; originally announced July 2024.

Comments: 15 pages, 13 figures. Accepted for publication in ApJ

Does Metallicity affect Protoplanetary Disk Fraction? Answers from the Outer Milky Way

Authors: Sudeshna Patra, Jessy Jose, Neal J. Evans II

Abstract: The role of metallicity in shaping protoplanetary disk evolution remains poorly comprehended. This study analyzes the disk fraction of 10 young (0.9-2.1 Myr) and low-metallicity (0.34-0.83 Z$_{\odot}$) clusters located in the outer Milky Way with Galactocentric distances between 10 and 13 kpc. Using $JHK$ data obtained from UKIDSS, the calculated disk fraction values for low-mass stars (0.2-2 M… ▽ More The role of metallicity in shaping protoplanetary disk evolution remains poorly comprehended. This study analyzes the disk fraction of 10 young (0.9-2.1 Myr) and low-metallicity (0.34-0.83 Z$_{\odot}$) clusters located in the outer Milky Way with Galactocentric distances between 10 and 13 kpc. Using $JHK$ data obtained from UKIDSS, the calculated disk fraction values for low-mass stars (0.2-2 M$_{\odot}$) ranged from 42% to 7%. To enhance the statistical reliability of our analysis, eight additional low-metallicity clusters are sourced from previous studies with metallicity range 0.25-0.85 Z$_{\odot}$ along with our sample, resulting in a total of 18 regions with low-metallicity. We find that low-metallicity clusters exhibit on average $2.6\pm0.2$ times lower disk fraction compared to solar-metallicity clusters in all the age bins we have. Within the age range we can probe, our study does not find evidence of faster disk decay in sub-solar metallicity regions compared to solar-metallicity regions. Furthermore, we observe a positive correlation between cluster disk fraction and metallicity for two different age groups of 0.3-1.4 and 1.4-2.5 Myr. We emphasize that both cluster age and metallicity significantly affect the fraction of stars with evidence of inner disks. △ Less

Submitted 8 May, 2024; originally announced May 2024.

Comments: 24 pages, 17 figures, 2 tables, accepted for publication in The Astrophysical Journal

JWST/MIRI detection of suprathermal OH rotational emissions: probing the dissociation of the water by Lyman alpha photons near the protostar HOPS 370

Authors: David A. Neufeld, P. Manoj, Himanshu Tyagi, Mayank Narang, Dan M. Watson, S. Thomas Megeath, Ewine F. Van Dishoeck, Robert A. Gutermuth, Thomas Stanke, Yao-Lun Yang, Adam E. Rubinstein, Guillem Anglada, Henrik Beuther, Alessio Caratti o Garatti, Neal J. Evans II, Samuel Federman, William J. Fischer, Joel Green, Pamela Klaassen, Leslie W. Looney, Mayra Osorio, Pooneh Nazari, John J. Tobin, Lukasz Tychoniec, Scott Wolk

Abstract: Using the MIRI/MRS spectrometer on JWST, we have detected pure rotational, suprathermal OH emissions from the vicinity of the intermediate-mass protostar HOPS 370 (OMC2/FIR3). These emissions are observed from shocked knots in a jet/outflow, and originate in states of rotational quantum number as high as 46 that possess excitation energies as large as $E_U/k = 4.65 \times 10^4$ K. The relative str… ▽ More Using the MIRI/MRS spectrometer on JWST, we have detected pure rotational, suprathermal OH emissions from the vicinity of the intermediate-mass protostar HOPS 370 (OMC2/FIR3). These emissions are observed from shocked knots in a jet/outflow, and originate in states of rotational quantum number as high as 46 that possess excitation energies as large as $E_U/k = 4.65 \times 10^4$ K. The relative strengths of the observed OH lines provide a powerful diagnostic of the ultraviolet radiation field in a heavily-extinguished region ($A_V \sim 10 - 20$) where direct UV observations are impossible. To high precision, the OH line strengths are consistent with a picture in which the suprathermal OH states are populated following the photodissociation of water in its $\tilde B - X$ band by ultraviolet radiation produced by fast ($\sim 80\,\rm km\,s^{-1}$) shocks along the jet. The observed dominance of emission from symmetric ($A^\prime$) OH states over that from antisymmetric ($A^{\prime\prime}$) states provides a distinctive signature of this particular population mechanism. Moreover, the variation of intensity with rotational quantum number suggests specifically that Ly$α$ radiation is responsible for the photodissociation of water, an alternative model with photodissociation by a 10$^4$ K blackbody being disfavored at a high level of significance. Using measurements of the Br$α$ flux to estimate the Ly$α$ production rate, we find that $\sim 4\%$ of the Ly$α$ photons are absorbed by water. Combined with direct measurements of water emissions in the $ν_2 = 1 -0$ band, the OH observations promise to provide key constraints on future models for the diffusion of Ly$α$ photons in the vicinity of a shock front. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 30 pages, 7 figures. Accepted for publication in ApJ Letters

The ALMA-QUARKS Survey: II. the ACA 1.3 mm continuum source catalog and the assembly of dense gas in massive star-forming clumps

Authors: Fengwei Xu, Ke Wang, Tie Liu, Lei Zhu, Guido Garay, Xunchuan Liu, Paul Goldsmith, Qizhou Zhang, Patricio Sanhueza, Shengli Qin, Jinhua He, Mika Juvela, Anandmayee Tej, Hongli Liu, Shanghuo Li, Kaho Morii, Siju Zhang, Jianwen Zhou, Amelia Stutz, Neal J. Evans, Kim Kee-Tae, Shengyuan Liu, Diego Mardones, Guangxing Li, Leonardo Bronfman , et al. (8 additional authors not shown)

Abstract: Leveraging the high resolution, high sensitivity, and wide frequency coverage of the Atacama Large Millimeter/submillimeter Array (ALMA), the QUARKS survey, standing for "Querying Underlying mechanisms of massive star formation with ALMA-Resolved gas Kinematics and Structures", is observing 139 massive star-forming clumps at ALMA Band 6 ($λ\sim$ 1.3 mm). This paper introduces the Atacama Compact A… ▽ More Leveraging the high resolution, high sensitivity, and wide frequency coverage of the Atacama Large Millimeter/submillimeter Array (ALMA), the QUARKS survey, standing for "Querying Underlying mechanisms of massive star formation with ALMA-Resolved gas Kinematics and Structures", is observing 139 massive star-forming clumps at ALMA Band 6 ($λ\sim$ 1.3 mm). This paper introduces the Atacama Compact Array (ACA) 7-m data. Combining multi-wavelength data, we provide the first edition of QUARKS atlas, offering insights into the multiscale and multiphase interstellar medium in high-mass star formation. The ACA 1.3 mm catalog includes 207 continuum sources that are called ACA sources. Their gas kinetic temperatures are estimated using three formaldehyde (H$_2$CO) transitions with a non-LTE radiation transfer model, and the mass and density are derived from a dust emission model. The ACA sources are massive (16-84 percentile values of 6-160 $M_{\odot}$), gravity-dominated ($M\propto R^{1.1}$) fragments within massive clumps, with supersonic turbulence ($\mathcal{M}>1$) and embedded star-forming protoclusters. We find a linear correlation between the masses of the fragments and the massive clumps, with a ratio of 6% between the two. When considering the fragments as representative of dense gas, the ratio indicates a dense gas fraction (DGF) of 6%, although with a wide scatter ranging from 1% to 10%. If we consider the QUARKS massive clumps to be what is observed at various scales, then the size-independent DGF indicates a self-similar fragmentation or collapsing mode in protocluster formation. With the ACA data over four orders of magnitude of luminosity-to-mass ratio ($L/M$), we find that the DGF increases significantly with $L/M$, which indicates clump evolutionary stage. We observed a limited fragmentation at the subclump scale, which can be explained by dynamic global collapse process. △ Less

Submitted 4 April, 2024; v1 submitted 2 April, 2024; originally announced April 2024.

Comments: 24 pages, 7 figures. Accepted for publication in Research in Astronomy and Astrophysics. QUARKS atlas link: https://drive.google.com/file/d/1KTqXxCDduYepvLd9kIvZVSSytK48OmfL/view?usp=sharing

JWST observations of $^{13}$CO$_{2}$ ice: Tracing the chemical environment and thermal history of ices in protostellar envelopes

Authors: Nashanty G. C. Brunken, Will R. M. Rocha, Ewine F. van Dishoeck, Robert Gutermuth, Himanshu Tyagi, Katerina Slavicinska, Pooneh Nazari, S. Thomas Megeath, Neal J. Evans II, Mayank Narang, P. Manoj, Adam E. Rubinstein, Dan M. Watson, Leslie W. Looney, Harold Linnartz, Alessio Caratti o Garatti, Henrik Beuther, Hendrik Linz, Pamela Klaassen, Charles A. Poteet, Samuel Federman, Guillem Anglada, Prabhani Atnagulov, Tyler L. Bourke, William J. Fischer , et al. (16 additional authors not shown)

Abstract: The structure and composition of simple ices can be modified during stellar evolution by protostellar heating. Key to understanding the involved processes are thermal and chemical tracers that can diagnose the history and environment of the ice. The 15.2 $μ$m bending mode of $^{12}$CO$_2$ has proven to be a valuable tracer of ice heating events but suffers from grain shape and size effects. A viab… ▽ More The structure and composition of simple ices can be modified during stellar evolution by protostellar heating. Key to understanding the involved processes are thermal and chemical tracers that can diagnose the history and environment of the ice. The 15.2 $μ$m bending mode of $^{12}$CO$_2$ has proven to be a valuable tracer of ice heating events but suffers from grain shape and size effects. A viable alternative tracer is the weaker $^{13}$CO$_2$ isotopologue band at 4.39 $μ$m which has now become accessible at high S/N with the $\textit{James Webb}$ Space Telescope (JWST). We present JWST NIRSpec observations of $^{13}$CO$_2$ ice in five deeply embedded Class 0 sources spanning a wide range in luminosities (0.2 - 10$^4$ L$_{\odot}$ ) taken as part of the Investigating Protostellar Accretion Across the Mass Spectrum (IPA) program. The band profiles vary significantly, with the most luminous sources showing a distinct narrow peak at 4.38 $μ$m. We first apply a phenomenological approach and show that a minimum of 3-4 Gaussian profiles are needed to fit the $^{13}$CO$_2$ absorption feature. We then combine these findings with laboratory data and show that a 15.2 $μ$m $^{12}$CO$_2$ band inspired five-component decomposition can be applied for the isotopologue band where each component is representative of CO$_2$ ice in a specific molecular environment. The final solution consists of cold mixtures of CO$_2$ with CH$_3$OH, H$_2$O and CO as well as segregated heated pure CO$_2$ ice. Our results are in agreement with previous studies of the $^{12}$CO$_2$ ice band, further confirming that $^{13}$CO$_{2}$ is a useful alternative tracer of protostellar heating events. We also propose an alternative solution consisting only of heated CO$_2$:CH$_3$OH and CO$_2$:H$_2$O ices and warm pure CO$_2$ ice for decomposing the ice profiles of the two most luminous sources in our sample. △ Less

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

Journal ref: A&A 685, A27 (2024)

Hunt for complex cyanides in protostellar ices with JWST: Tentative detection of CH$_3$CN and C$_2$H$_5$CN

Authors: P. Nazari, W. R. M. Rocha, A. E. Rubinstein, K. Slavicinska, M. G. Rachid, E. F. van Dishoeck, S. T. Megeath, R. Gutermuth, H. Tyagi, N. Brunken, M. Narang, P. Manoj, D. M. Watson, N. J. Evans II, S. Federman, J. Muzerolle Page, G. Anglada, H. Beuther, P. Klaassen, L. W. Looney, M. Osorio, T. Stanke, Y. -L. Yang

Abstract: Nitrogen-bearing complex organic molecules have been commonly detected in the gas phase but not yet in interstellar ices. This has led to the long-standing question of whether these molecules form in the gas phase or in ices. $\textit{James Webb}$ Space Telescope ($\textit{JWST}$) offers the sensitivity, spectral resolution, and wavelength coverage needed to detect them in ices and investigate whe… ▽ More Nitrogen-bearing complex organic molecules have been commonly detected in the gas phase but not yet in interstellar ices. This has led to the long-standing question of whether these molecules form in the gas phase or in ices. $\textit{James Webb}$ Space Telescope ($\textit{JWST}$) offers the sensitivity, spectral resolution, and wavelength coverage needed to detect them in ices and investigate whether their abundance ratios are similar in gas and ice. We report the first tentative detection of CH$_3$CN, C$_2$H$_5$CN, and the simple molecule, N$_2$O, based on the CN-stretch band in interstellar ices toward three (HOPS 153, HOPS 370, and IRAS 20126+4104) out of the five protostellar systems observed as part of the Investigating Protostellar Accretion (IPA) GO program with $\textit{JWST}$-NIRSpec. We also provide upper limits for the two other sources with smaller luminosities in the sample. We detect OCN$^-$ in the ices of all sources with typical CH$_3$CN/OCN$^-$ ratios of around 1. Ice and gas column density ratios of the nitrogen-bearing species with respect to each other are better matched than those with respect to methanol, which are a factor of ${\sim}5$ larger in the ices than the gas. We attribute the elevated ice column densities with respect to methanol to the difference in snowline locations of nitrogen-bearing molecules and of methanol, biasing the gas-phase observations toward fewer nitrogen-bearing molecules. Moreover, we find tentative evidence for enhancement of OCN$^-$, CH$_3$CN, and C$_2$H$_5$CN in warmer ices, although formation of these molecules likely starts along with methanol in the cold prestellar phase. Future surveys combining NIRSpec and MIRI, and additional laboratory spectroscopic measurements of C$_2$H$_5$CN ice, are necessary for robust detection and conclusions on the formation history of complex cyanides. △ Less

Submitted 15 January, 2024; originally announced January 2024.

Comments: Accepted for publication in A&A

Journal ref: A&A 686, A71 (2024)

IPA. Class 0 Protostars Viewed in CO Emission Using JWST/NIRSpec

Authors: Adam E. Rubinstein, Himanshu Tyagi, Pooneh Nazari, Robert Gutermuth, Samuel Federman, Mayank Narang, Will R. M. Rocha, Nashanty Brunken, Katie Slavicinska, Neal J. Evans II, Joel D. Green, Dan M. Watson, Henrik Beuther, Tyler Bourke, Alessio Caratti o Garatti, Lee Hartmann, Pamela Klaassen, Hendrik Linz, Leslie W. Looney, Puravankara Manoj, S. Thomas Megeath, James Muzerolle Page, Thomas Stanke, John J. Tobin, Ewine F. van Dishoeck , et al. (2 additional authors not shown)

Abstract: We investigate the bright CO fundamental emission in the central regions of five Class 0 protostars using the JWST's Near-Infrared Spectrograph (NIRSpec) and provide clues to what processes excite the gas. CO line emission images are extracted for a forest of $\sim$150 ro-vibrational transitions from two vibrational bands, $v=1-0$ and $v=2-1$. However, ${}^{13}$CO is not detected, and thus we can… ▽ More We investigate the bright CO fundamental emission in the central regions of five Class 0 protostars using the JWST's Near-Infrared Spectrograph (NIRSpec) and provide clues to what processes excite the gas. CO line emission images are extracted for a forest of $\sim$150 ro-vibrational transitions from two vibrational bands, $v=1-0$ and $v=2-1$. However, ${}^{13}$CO is not detected, and thus we can only statistically constrain the ${}^{12}$CO optical depth. Using noise measurements to determine upper limits to the ${}^{13}$CO emission, the flux ratio of ${}^{12}$CO/${}^{13}$CO indicates that the ${}^{12}$CO emission itself is not optically thick for ro-vibrational transitions with upper state rotational quantum number $J_u \geq 15$. We construct population diagrams to estimate the rotational temperature and number of molecules from extinction-corrected CO line fluxes assuming CO emission is optically thin. Two different temperature components are required for $v=1$ ($\sim600-1000$ K and $\sim1500-3500$ K), while one hotter component is required for $v=2$ ($\sim2000-6000$ K). The vibrational temperature is $\sim 900$ K among our sources and shows no trend with luminosity. Using vibrational temperatures and the inferred total amount of CO molecules for our sources, the total warm gas mass correlates strongly with luminosity ranging from $\sim$0.1 $\rm M_{Earth}$ for the low-mass protostars to $\sim$1 M$_{\rm sun}$ for the high-mass protostars. Interpreting the distribution of gas column densities and temperatures depends on radiative and chemical processes affecting CO. The presence of a $v=2$ population may indicate CO gas radiatively excited. Selective UV photodissociation of CO isotopologues around our high-mass sources may explain their depletion of ${}^{13}$CO. △ Less

Submitted 12 December, 2023; originally announced December 2023.

Comments: 30 pages, 6 figures, 5 tables, submitted to ApJ

The CO outflow components ejected by a recent accretion event in B335

Authors: Chul-Hwan Kim, Jeong-Eun Lee, Carlos Contreras Peña, Doug Johnstone, Gregory J. Herczeg, John J. Tobin, Neal J. Evans II

Abstract: Protostellar outflows often present a knotty appearance, providing evidence of sporadic accretion in stellar mass growth. To understand the direct relation between mass accretion and ejection, we analyze the contemporaneous accretion activity and associated ejection components in B335. B335 has brightened in the mid-IR by 2.5 mag since 2010, indicating increased luminosity, presumably due to incre… ▽ More Protostellar outflows often present a knotty appearance, providing evidence of sporadic accretion in stellar mass growth. To understand the direct relation between mass accretion and ejection, we analyze the contemporaneous accretion activity and associated ejection components in B335. B335 has brightened in the mid-IR by 2.5 mag since 2010, indicating increased luminosity, presumably due to increased mass accretion rate onto the protostar. ALMA observations of 12CO emission in the outflow reveal high-velocity emission, estimated to have been ejected 4.6 - 2 years before the ALMA observation and consistent with the jump in mid-IR brightness. The consistency in timing suggests that the detected high-velocity ejection components are directly linked to the most recent accretion activity. We calculated the kinetic energy, momentum, and force for the ejection component associated with the most recent accretion activity and found that at least, about 1.0% of accreted mass has been ejected. More accurate information on the jet inclination and the temperature of the ejected gas components will better constrain the ejected mass induced by the recently enhanced accretion event. △ Less

Submitted 10 December, 2023; originally announced December 2023.

Comments: 11 pages, 5 figures, and 1 table

Discovery of a collimated jet from the low luminosity protostar IRAS 16253$-$2429 in a quiescent accretion phase with the JWST

Authors: Mayank Narang, Manoj P., Himanshu Tyagi, Dan M. Watson, S. Thomas Megeath, Samuel Federman, Adam E. Rubinstein, Robert Gutermuth, Alessio Caratti o Garatti, Henrik Beuther, Tyler L. Bourke, Ewine F. Van Dishoeck, Neal J. Evans II, Guillem Anglada, Mayra Osorio, Thomas Stanke, James Muzerolle, Leslie W. Looney, Yao-Lun Yang, John J. Tobin, Pamela Klaassen, Nicole Karnath, Prabhani Atnagulov, Nashanty Brunken, William J. Fischer , et al. (14 additional authors not shown)

Abstract: Investigating Protostellar Accretion (IPA) is a JWST Cycle~1 GO program that uses NIRSpec IFU and MIRI MRS to obtain 2.9--28~$μ$m spectral cubes of young, deeply embedded protostars with luminosities of 0.2 to 10,000~L$_{\odot}$ and central masses of 0.15 to 12~M$_{\odot}$. In this Letter, we report the discovery of a highly collimated atomic jet from the Class~0 protostar IRAS~16253$-$2429, the l… ▽ More Investigating Protostellar Accretion (IPA) is a JWST Cycle~1 GO program that uses NIRSpec IFU and MIRI MRS to obtain 2.9--28~$μ$m spectral cubes of young, deeply embedded protostars with luminosities of 0.2 to 10,000~L$_{\odot}$ and central masses of 0.15 to 12~M$_{\odot}$. In this Letter, we report the discovery of a highly collimated atomic jet from the Class~0 protostar IRAS~16253$-$2429, the lowest luminosity source ($L_\mathrm{bol}$ = 0.2 $L_\odot$) in the IPA program. The collimated jet is detected in multiple [Fe~II] lines, [Ne~II], [Ni~II], and H~I lines, but not in molecular emission. The atomic jet has a velocity of about 169~$\pm$~15~km\,s$^{-1}$, after correcting for inclination. The width of the jet increases with distance from the central protostar from 23 to~60 au, corresponding to an opening angle of 2.6~$\pm$~0.5\arcdeg. By comparing the measured flux ratios of various fine structure lines to those predicted by simple shock models, we derive a shock {speed} of 54~km\,s$^{-1}$ and a preshock density of 2.0$\times10^{3}$~cm$^{-3}$ at the base of the jet. {From these quantities and using a suite of jet models and extinction laws we compute a mass loss rate between $0.4 -1.1\times10^{-10}~M_{\odot}$~yr~$^{-1}$.} The low mass loss rate is consistent with simultaneous measurements of low mass accretion rate ($2.4~\pm~0.8~\times~10^{-9}~M_{\odot}$~yr$^{-1}$) for IRAS~16253$-$2429 from JWST observations (Watson et al. in prep), indicating that the protostar is in a quiescent accretion phase. Our results demonstrate that very low-mass protostars can drive highly collimated, atomic jets, even during the quiescent phase. △ Less

Submitted 11 January, 2024; v1 submitted 21 October, 2023; originally announced October 2023.

Comments: Accepted to ApJL. Comments and feedback welcome

Investigating Protostellar Accretion-Driven Outflows Across the Mass Spectrum: JWST NIRSpec IFU 3-5~$μ$m Spectral Mapping of Five Young Protostars

Authors: Samuel Federman, S. Thomas Megeath, Adam E. Rubinstein, Robert Gutermuth, Mayank Narang, Himanshu Tyagi, P. Manoj, Guillem Anglada, Prabhani Atnagulov, Henrik Beuther, Tyler L. Bourke, Nashanty Brunken, Alessio Caratti o Garatti, Neal J. Evans II, William J. Fischer, Elise Furlan, Joel Green, Nolan Habel, Lee Hartmann, Nicole Karnath, Pamela Klaassen, Hendrik Linz, Leslie W. Looney, Mayra Osorio, James Muzerolle Page , et al. (13 additional authors not shown)

Abstract: Investigating Protostellar Accretion is a Cycle 1 JWST program using the NIRSpec+MIRI integral field units to obtain 2.9--28 $μ$m spectral cubes of five young protostars with luminosities of 0.2-10,000 L$_{\odot}$ in their primary accretion phase. This paper introduces the NIRSpec 2.9--5.3 $μ$m data of the inner 840-9000 au with spatial resolutions from 28-300 au. The spectra show rising continuum… ▽ More Investigating Protostellar Accretion is a Cycle 1 JWST program using the NIRSpec+MIRI integral field units to obtain 2.9--28 $μ$m spectral cubes of five young protostars with luminosities of 0.2-10,000 L$_{\odot}$ in their primary accretion phase. This paper introduces the NIRSpec 2.9--5.3 $μ$m data of the inner 840-9000 au with spatial resolutions from 28-300 au. The spectra show rising continuum emission; deep ice absorption; emission from H$_{2}$, H~I, and [Fe~II]; and the CO fundamental series in emission and absorption. Maps of the continuum emission show scattered light cavities for all five protostars. In the cavities, collimated jets are detected in [Fe~II] for the four $< 320$~L$_{\odot}$ protostars, two of which are additionally traced in Br-$α$. Knots of [Fe~II] emission are detected toward the most luminous protostar, and knots of [FeII] emission with dynamical times of $< 30$~yrs are found in the jets of the others. While only one jet is traced in H$_2$, knots of H$_2$ and CO are detected in the jets of four protostars. H$_2$ is seen extending through the cavities, showing that they are filled by warm molecular gas. Bright H$_2$ emission is seen along the walls of a single cavity, while in three cavities narrow shells of H$_2$ emission are found, one of which has an [Fe~II] knot at its apex. These data show cavities containing collimated jets traced in atomic/ionic gas surrounded by warm molecular gas in a wide-angle wind and/or gas accelerated by bow shocks in the jets. △ Less

Submitted 24 April, 2024; v1 submitted 5 October, 2023; originally announced October 2023.

Comments: 26 pages, 12 figures

Journal ref: ApJ 966 41 (2024)

The ALMA Survey of Star Formation and Evolution in Massive Protoclusters with Blue Profiles (ASSEMBLE): Core Growth, Cluster Contraction, and Primordial Mass Segregation

Authors: Fengwei Xu, Ke Wang, Tie Liu, Mengyao Tang, Neal J. Evans II, Aina Palau, Kaho Morii, Jinhua He, Patricio Sanhueza, Hong-Li Liu, Amelia Stutz, Qizhou Zhang, Xi Chen, Pak Shing Li, Gilberto C. Gómez, Enrique Vázquez-Semadeni, Shanghuo Li, Xiaofeng Mai, Xing Lu, Meizhu Liu, Li Chen, Chuanshou Li, Hongqiong Shi, Zhiyuan Ren, Di Li , et al. (18 additional authors not shown)

Abstract: The ALMA Survey of Star Formation and Evolution in Massive Protoclusters with Blue Profiles (ASSEMBLE) aims to investigate the process of mass assembly and its connection to high-mass star formation theories in protoclusters in a dynamic view. We observed 11 massive (Mclump>1000 Msun), luminous (Lbol>10,000 Lsun), and blue-profile (infall signature) clumps by ALMA with resolution of 2200-5500 au a… ▽ More The ALMA Survey of Star Formation and Evolution in Massive Protoclusters with Blue Profiles (ASSEMBLE) aims to investigate the process of mass assembly and its connection to high-mass star formation theories in protoclusters in a dynamic view. We observed 11 massive (Mclump>1000 Msun), luminous (Lbol>10,000 Lsun), and blue-profile (infall signature) clumps by ALMA with resolution of 2200-5500 au at 350 GHz (870 um) in continuum and line emission. 248 dense cores were identified, including 106 cores showing protostellar signatures and 142 prestellar core candidates. Compared to early-stage infrared dark clouds (IRDCs) by ASHES, the core mass and surface density within the ASSEMBLE clumps exhibited significant increment, suggesting concurrent core accretion during the evolution of the clumps. The maximum mass of prestellar cores was found to be 2 times larger than that in IRDCs, indicating evolved protoclusters have the potential to harbor massive prestellar cores. The mass relation between clumps and their most massive core (MMCs) is observed in ASSEMBLE but not in IRDCs, which is suggested to be regulated by multiscale mass accretion. The mass correlation between the core clusters and their MMCs has a steeper slope compared to that observed in stellar clusters, which can be due to fragmentation of the MMC and stellar multiplicity. We observe a decrease in core separation and an increase in central concentration as protoclusters evolve. We confirm primordial mass segregation in the ASSEMBLE protoclusters, possibly resulting from gravitational concentration and/or gas accretion. △ Less

Submitted 26 September, 2023; originally announced September 2023.

Comments: 37 pages, 13 figures, 5 tables; accepted for publication in ApJS

Spectral survey of a Hot core with an Eruptive Accretion in S255IR NIRS3 (SHEA): The discovery of class I and class II millimeter methanol maser transitions

Authors: Giseon Baek, Jeong-Eun Lee, Neal J. Evans II, Tomoya Hirota, Yuri Aikawa, Ji-hyun Kang, Jungha Kim, Jes K. Jørgensen

Abstract: We report the detection of the millimeter CH$_3$OH masers including a new detection of class I (11$_{0,11}$-10$_{1,10}$A) and class II (6$_{1,5}$-5$_{2,4}$E) maser transitions toward the high-mass protostar S255IR NIRS3 in post-burst phase. The CH$_3$OH emissions were detected as a mixture of maser and thermal characteristics. We examine the detected transitions using an excitation diagram and LTE… ▽ More We report the detection of the millimeter CH$_3$OH masers including a new detection of class I (11$_{0,11}$-10$_{1,10}$A) and class II (6$_{1,5}$-5$_{2,4}$E) maser transitions toward the high-mass protostar S255IR NIRS3 in post-burst phase. The CH$_3$OH emissions were detected as a mixture of maser and thermal characteristics. We examine the detected transitions using an excitation diagram and LTE model spectra and compare the observed properties with those of thermal lines. Class II CH$_3$OH maser transitions showed distinctive intensity and velocity distributions from those of thermal transitions. Bright distinct emission components in addition to the fragmented and arc-shaped emissions are only detected in class I CH$_3$OH maser transitions toward southern and western directions from the protostellar position, implying the presence of the slow outflow shocks. △ Less

Submitted 11 August, 2023; originally announced August 2023.

Comments: 9 pages, 4 Figures, 1 Table, Accepted for publication in ApJL

Report number: ApJL 954 L25

Journal ref: The Astrophysical Journal Letters, 2023, Volume 954, Number 1

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

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

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

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

Comments: 28 pages, 11 figures

Triple spiral arms of a triple protostar system imaged in molecular lines

Authors: Jeong-Eun Lee, Tomoaki Matsumoto, Hyun-Jeong Kim, Seokho Lee, Daniel Harsono, Jaehan Bae, Neal J. Evans II, Shu-ichiro Inutsuka, Minho Choi, Ken'ichi Tatematsu, Jae-Joon Lee, Dan Jaffe

Abstract: Most stars form in multiple star systems. For a better understanding of their formation processes, it is important to resolve the individual protostellar components and the surrounding envelope and disk material at the earliest possible formation epoch because the formation history can be lost in a few orbital timescales. Here we present the ALMA observational results of a young multiple protostel… ▽ More Most stars form in multiple star systems. For a better understanding of their formation processes, it is important to resolve the individual protostellar components and the surrounding envelope and disk material at the earliest possible formation epoch because the formation history can be lost in a few orbital timescales. Here we present the ALMA observational results of a young multiple protostellar system, IRAS 04239+2436, where three well-developed large spiral arms were detected in the shocked SO emission. Along the most conspicuous arm, the accretion streamer was also detected in the SO$_2$ emission. The observational results are complemented by numerical magneto-hydrodynamic simulations, where those large arms only appear in magnetically weakened clouds. The numerical simulations also suggest that the large triple spiral arms are the result of gravitational interactions between compact triple protostars and the turbulent infalling envelope. △ Less

Submitted 10 June, 2023; originally announced June 2023.

Comments: Accepted for publication in ApJ

High-resolution APEX/LAsMA $^{12}$CO and $^{13}$CO (3-2) observation of the G333 giant molecular cloud complex : I. Evidence for gravitational acceleration in hub-filament systems

Authors: J. W. Zhou, F. Wyrowski, S. Neupane, J. S. Urquhart, N. J. Evans II, E. Vázquez-Semadeni, K. M. Menten, Y. Gong, T. Liu

Abstract: Hub-filament systems are suggested to be the birth cradles of high-mass stars and clusters. We apply the FILFINDER algorithm to the integrated intensity maps of the 13CO (3-2) line to identify filaments in the G333 complex, and extract the velocity and intensity along the filament skeleton from moment maps. Clear velocity and density fluctuations are seen along the filaments, allowing us to fit ve… ▽ More Hub-filament systems are suggested to be the birth cradles of high-mass stars and clusters. We apply the FILFINDER algorithm to the integrated intensity maps of the 13CO (3-2) line to identify filaments in the G333 complex, and extract the velocity and intensity along the filament skeleton from moment maps. Clear velocity and density fluctuations are seen along the filaments, allowing us to fit velocity gradients around the intensity peaks. The velocity gradients fitted to the LAsMA data and ALMA data agree with each other over the scales covered by ALMA observations in the ATOMS survey. Changes of velocity gradient with scale indicate a ''funnel'' structure of the velocity field in PPV space, indicative of a smooth, continuously increasing velocity gradient from large to small scales, and thus consistent with gravitational acceleration. The typical velocity gradient corresponding to a 1 pc scale is ~1.6km/s/pc. Assuming free-fall, we estimate a kinematic mass within 1 pc of ~1190 M$_\odot$, which is consistent with typical masses of clumps in the ATLASGAL survey. We find direct evidence for gravitational acceleration from comparison of the observed accelerations to those predicted by free-fall onto dense hubs. On large scales, we find that the inflow may be driven by the larger scale structure, consistent with hierarchical structure in the molecular cloud and gas inflow from large to small scales. The hub-filament structures at different scales may be organized into a hierarchical system extending up to the largest scales probed, through the coupling of gravitational centers at different scales. We argue that the ''funnel'' structure in PPV space can be an effective probe for the gravitational collapse motions in molecular clouds. The large scale gas inflow is driven by gravity, implying that the molecular clouds in G333 complex may be in the state of global gravitational collapse. △ Less

Submitted 21 May, 2023; originally announced May 2023.

Journal ref: A&A 676, A69 (2023)

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

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

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

Submitted 14 February, 2023; originally announced February 2023.

Comments: Accepted for publication in ApJ

First Detection of Radio Recombination Lines of Ions Heavier than Helium

Authors: Xunchuan Liu, Tie Liu, Zhiqiang Shen, Paul F. Goldsmith, Neal J. Evans II, Sheng-Li Qin, Qiuyi Luo, Yu Cheng, Sheng-Yuan Liu, Fengyao Zhu, Ken'ichi Tatematsu, Meizhu Liu, Dongting Yang, Chuanshou Li, Li Cen, Juan Li, Xing Lu, Qilao Gu, Rongbing Zhao, Bing Li, Yajun Wu, Weiye Zhong, Zhang Zhao, Jinqing Wang, Qinghui Liu , et al. (10 additional authors not shown)

Abstract: We report the first detection of radio recombination lines (RRLs) of ions heavier than helium. In a highly sensitive multi-band (12--50 GHz) line survey toward Orion KL with the TianMa 65-m Radio Telescope (TMRT), we successfully detected more than fifteen unblended $α$ lines of RRLs of singly ionized species (XII) recombined from XIII. The Ka-band (35--50 GHz) spectrum also shows tentative signal… ▽ More We report the first detection of radio recombination lines (RRLs) of ions heavier than helium. In a highly sensitive multi-band (12--50 GHz) line survey toward Orion KL with the TianMa 65-m Radio Telescope (TMRT), we successfully detected more than fifteen unblended $α$ lines of RRLs of singly ionized species (XII) recombined from XIII. The Ka-band (35--50 GHz) spectrum also shows tentative signals of $β$ lines of ions. The detected lines can be successfully crossmatched with the the rest frequencies of RRLs of CII and/or OII. This finding greatly expands the connotation of ion RRLs, since before this work only two blended lines (105$α$ and 121$α$) of HeII had been reported. Our detected lines can be fitted simultaneously under assumption of local thermodynamic equilibrium (LTE). An abundance of CIII and OIII of 8.8$\times$10$^{-4}$ is obtained, avoiding the complexities of optical/infrared observations and the blending of RRLs of atoms. It is consistent with but approaches the upper bound of the value (10$^{-4}$--$10^{-3}$) estimated from optical/infrared observations. The effects of dielectronic recombination may contribute to enhancing the level populations even at large $n$. We expect future observations using radio interferometers could break the degeneracy between C and O, and help to reveal the ionization structure and dynamical evolution of various ionized regions. △ Less

Submitted 7 February, 2023; originally announced February 2023.

Comments: 3 figures, 2 tables, accepted by A&A Letter

Journal ref: A&A 671, L1 (2023)

Models of Rotating Infall for the B335 Protostar

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

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

Submitted 7 December, 2022; originally announced December 2022.

Comments: 22 pages, 20 figures, in press ApJ

CORINOS I: JWST/MIRI Spectroscopy and Imaging of a Class 0 protostar IRAS 15398-3359

Authors: Yao-Lun Yang, Joel D. Green, Klaus M. Pontoppidan, Jennifer B. Bergner, L. Ilsedore Cleeves, Neal J. Evans II, Robin T. Garrod, Mihwa Jin, Chul Hwan Kim, Jaeyeong Kim, Jeong-Eun Lee, Nami Sakai, Christopher N. Shingledecker, Brielle Shope, John J. Tobin, Ewine van Dishoeck

Abstract: The origin of complex organic molecules (COMs) in young Class 0 protostars has been one of the major questions in astrochemistry and star formation. While COMs are thought to form on icy dust grains via gas-grain chemistry, observational constraints on their formation pathways have been limited to gas-phase detection. Sensitive mid-infrared spectroscopy with JWST enables unprecedented investigatio… ▽ More The origin of complex organic molecules (COMs) in young Class 0 protostars has been one of the major questions in astrochemistry and star formation. While COMs are thought to form on icy dust grains via gas-grain chemistry, observational constraints on their formation pathways have been limited to gas-phase detection. Sensitive mid-infrared spectroscopy with JWST enables unprecedented investigation of COM formation by measuring their ice absorption features. We present an overview of JWST/MIRI MRS spectroscopy and imaging of a young Class 0 protostar, IRAS 15398-3359, and identify several major solid-state absorption features in the 4.9-28 $μ$m wavelength range. These can be attributed to common ice species, such as H$_2$O, CH$_3$OH, NH$_3$, and CH$_4$, and may have contributions from more complex organic species, such as C$_2$H$_5$OH and CH$_3$CHO. The MRS spectra show many weaker emission lines at 6-8 $μ$m, which are due to warm CO gas and water vapor, possibly from a young embedded disk previously unseen. Finally, we detect emission lines from [Fe II], [Ne II], [S I], and H$_2$, tracing a bipolar jet and outflow cavities. MIRI imaging serendipitously covers the south-western (blue-shifted) outflow lobe of IRAS 15398-3359, showing four shell-like structures similar to the outflows traced by molecular emission at sub-mm wavelengths. This overview analysis highlights the vast potential of JWST/MIRI observations and previews scientific discoveries in the coming years. △ Less

Submitted 14 November, 2022; v1 submitted 22 August, 2022; originally announced August 2022.

Comments: Accepted to ApJL. The methods of data reduction and spectral extraction are updated

Tracers of Dense Gas in the Outer Galaxy

Authors: Sudeshna Patra, Neal J. Evans II, Kee-Tae Kim, Mark Heyer, Jens Kauffmann, Jessy Jose, Manash R. Samal, Swagat R. Das

Abstract: We have mapped HCN and HCO$^{+}$ (J = 1 $\to$ 0) line emission toward a sample of seven star-forming regions (with 12 + log[O/H] range from 8.34 to 8.69) in the outer Milky Way (Galactocentric distance > 9.5 kpc), using the 14-meter radio telescope of the Taeduk Radio Astronomy Observatory (TRAO). We compare these two molecular lines with other conventional tracers of dense gas, millimeter-wave co… ▽ More We have mapped HCN and HCO$^{+}$ (J = 1 $\to$ 0) line emission toward a sample of seven star-forming regions (with 12 + log[O/H] range from 8.34 to 8.69) in the outer Milky Way (Galactocentric distance > 9.5 kpc), using the 14-meter radio telescope of the Taeduk Radio Astronomy Observatory (TRAO). We compare these two molecular lines with other conventional tracers of dense gas, millimeter-wave continuum emission from dust and extinction thresholds ($A_{V} \geq 8$ mag), inferred from the $^{13}$CO line data. HCN and HCO$^{+}$ correlate better with the millimeter emission than with the extinction criterion. A significant amount of luminosity comes from regions below the extinction criterion and outside the millimeter clump for all the clouds. The average fraction of HCN luminosity from within the regions with $A_{V} \geq 8$ mag is $0.343\pm0.225$; for the regions of millimeter emission, it is $0.478\pm0.149$. Based on a comparison with column density maps from Herschel, HCN and HCO$^{+}$ trace dense gas in high column density regions better than does $^{13}$CO. HCO$^{+}$ is less concentrated than HCN for outer Galaxy targets, in contrast with the inner Galaxy sample, suggesting that metallicity may affect the interpretation of tracers of dense gas. The conversion factor between the dense gas mass ($M_{dense}$) and line luminosities of HCN and HCO$^{+}$, when integrated over the whole cloud, is comparable with factors used in extragalactic studies. △ Less

Submitted 23 July, 2022; originally announced July 2022.

Comments: 20 pages, 8 figures, 5 tables, accepted for publication in The Astronomical Journal

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- XI. From inflow to infall in hub-filament systems

Authors: Jian-Wen Zhou, Tie Liu, Neal J. Evans II, Guido Garay, Paul F. Goldsmith, Gilberto C. Gomez, Enrique Vazquez-Semadeni, Hong-Li Liu, Amelia M. Stutz, Ke Wang, Mika Juvela, Jinhua He, Di Li, Leonardo Bronfman, Xunchuan Liu, Feng-Wei Xu, Anandmayee Tej, L. K. Dewangan, Shanghuo Li, Siju Zhang, Chao Zhang, Zhiyuan Ren, Kenichi Tatematsu, Pak Shing Li, Chang Won Lee , et al. (15 additional authors not shown)

Abstract: We investigate the presence of hub-filament systems in a large sample of 146 active proto-clusters, using H$^{13}$CO$^{+}$ J=1-0 molecular line data obtained from the ATOMS survey. We find that filaments are ubiquitous in proto-clusters, and hub-filament systems are very common from dense core scales ($\sim$0.1 pc) to clump/cloud scales ($\sim$1-10 pc). The proportion of proto-clusters containing… ▽ More We investigate the presence of hub-filament systems in a large sample of 146 active proto-clusters, using H$^{13}$CO$^{+}$ J=1-0 molecular line data obtained from the ATOMS survey. We find that filaments are ubiquitous in proto-clusters, and hub-filament systems are very common from dense core scales ($\sim$0.1 pc) to clump/cloud scales ($\sim$1-10 pc). The proportion of proto-clusters containing hub-filament systems decreases with increasing dust temperature ($T_d$) and luminosity-to-mass ratios ($L/M$) of clumps, indicating that stellar feedback from H{\sc ii} regions gradually destroys the hub-filament systems as proto-clusters evolve. Clear velocity gradients are seen along the longest filaments with a mean velocity gradient of 8.71 km s$^{-1}$pc$^{-1}$ and a median velocity gradient of 5.54 km s$^{-1}$pc$^{-1}$. We find that velocity gradients are small for filament lengths larger than $\sim$1~pc, probably hinting at the existence of inertial inflows, although we cannot determine whether the latter are driven by large-scale turbulence or large-scale gravitational contraction. In contrast, velocity gradients below $\sim$1~pc dramatically increase as filament lengths decrease, indicating that the gravity of the hubs or cores starts to dominate gas infall at small scales. We suggest that self-similar hub-filament systems and filamentary accretion at all scales may play a key role in high-mass star formation. △ Less

Submitted 16 June, 2022; originally announced June 2022.

Comments: 16 pages

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

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

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

Submitted 13 April, 2022; originally announced April 2022.

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

Slow Star Formation in the Milky Way: Theory Meets Observations

Authors: Neal J. Evans II, Jeong-Gyu Kim, Eve C. Ostriker

Abstract: The observed star formation rate of the Milky Way can be explained by applying a metallicity-dependent factor to convert CO luminosity to molecular gas mass and a star formation efficiency per free-fall time that depends on the virial parameter of a molecular cloud. These procedures also predict the trend of star formation rate surface density with Galactocentric radius. The efficiency per free-fa… ▽ More The observed star formation rate of the Milky Way can be explained by applying a metallicity-dependent factor to convert CO luminosity to molecular gas mass and a star formation efficiency per free-fall time that depends on the virial parameter of a molecular cloud. These procedures also predict the trend of star formation rate surface density with Galactocentric radius. The efficiency per free-fall time variation with virial parameter plays the major role in bringing theory into agreement with observations for the total star formation rate, while the metallicity dependence of the CO luminosity to mass conversion is most notable in the variation with Galactocentric radius. Application of these changes resolves a factor of over 100 discrepancy between observed and theoretical star formation rates that has been known for nearly 50 years. △ Less

Submitted 5 April, 2022; originally announced April 2022.

Comments: In press, Astrophysical Journal Letters 11 pages, 6 figures

Nobeyama Survey of Inward Motions toward Cores in Orion Identified by SCUBA-2

Authors: Ken'ichi Tatematsu, You-Ting Yeh, Naomi Hirano, Sheng-Yuan Liu, Tie Liu, Somnath Dutta, Dipen Sahu, Neal J. Evans II, Mika Juvela, Hee-Weon Yi, Jeong-Eun Lee, Patricio Sanhueza, Shanghuo Li, David Eden, Gwanjeong Kim, Chin-Fei Lee, Yuefang Wu, Kee-Tae Kim, L. Viktor T'oth, Minho Choi, Miju Kang, Mark A. Thompson, Gary A. Fuller, 20 Di Li, Ke Wang , et al. (6 additional authors not shown)

Abstract: In this study, 36 cores (30 starless and 6 protostellar) identified in Orion were surveyed to search for inward motions. We used the Nobeyama 45 m radio telescope, and mapped the cores in the $J = 1\rightarrow0$ transitions of HCO$^+$, H$^{13}$CO$^+$, N$_2$H$^+$, HNC, and HN$^{13}$C. The asymmetry parameter $δV$, which was the ratio of the difference between the HCO$^+$ and H$^{13}$CO$^+$ peak vel… ▽ More In this study, 36 cores (30 starless and 6 protostellar) identified in Orion were surveyed to search for inward motions. We used the Nobeyama 45 m radio telescope, and mapped the cores in the $J = 1\rightarrow0$ transitions of HCO$^+$, H$^{13}$CO$^+$, N$_2$H$^+$, HNC, and HN$^{13}$C. The asymmetry parameter $δV$, which was the ratio of the difference between the HCO$^+$ and H$^{13}$CO$^+$ peak velocities to the H$^{13}$CO$^+$ line width, was biased toward negative values, suggesting that inward motions were more dominant than outward motions. Three starless cores (10% of all starless cores surveyed) were identified as cores with blue-skewed line profiles (asymmetric profiles with more intense blue-shifted emission), and another two starless cores (7%) were identified as candidate blue-skewed line profiles. The peak velocity difference between HCO$^+$ and H$^{13}$CO$^+$ of them was up to 0.9 km s$^{-1}$, suggesting that some inward motions exceeded the speed of sound for the quiescent gas ($\sim10-17$ K). The mean of $δV$ of the five aforementioned starless cores was derived to be $-$0.5$\pm$0.3. One core, G211.16$-$19.33North3, observed using the ALMA ACA in DCO$^+$ $J = 3\rightarrow2$ exhibited blue-skewed features. Velocity offset in the blue-skewed line profile with a dip in the DCO$^+$ $J = 3\rightarrow2$ line was larger ($\sim 0.5$ km s$^{-1}$) than that in HCO$^+$ $J = 1\rightarrow0$ ($\sim 0.2$ km s$^{-1}$), which may represent gravitational acceleration of inward motions. It seems that this core is at the last stage in the starless phase, judging from the chemical evolution factor version 2.0 (CEF2.0). △ Less

Submitted 24 March, 2022; originally announced March 2022.

Comments: 45 pages, 28 figures, ApJ in press

Journal ref: Astrophysical J. 931 (2022) 33

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions-IX. A pilot study towards IRDC G034.43+00.24 on multi-scale structures and gas kinematics

Authors: Hong-Li Liu, Anandmayee Tej, Tie Liu, Paul F. Goldsmith, Amelia Stutz, Mika Juvela, Sheng-Li Qin, Feng-Wei Xu, Leonardo Bronfman, Neal J. Evans, Anindya Saha, Namitha Issac, Ken'ichi Tatematsu, Ke Wang, Shanghuo Li, Siju Zhang, Tapas Baug, Lokesh Dewangan, Yue-Fang Wu, Yong Zhang, Chang Won Lee, Xun-Chuan Liu, Jianwen Zhou, Archana Soam

Abstract: We present a comprehensive study of the gas kinematics associated with density structures at different spatial scales in the filamentary infrared dark cloud, G034.43+00.24 (G34). This study makes use of the H13CO+ (1-0) molecular line data from the ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS) survey, which has spatial and velocity resolution of 0.04 pc and 0.2 km/s, r… ▽ More We present a comprehensive study of the gas kinematics associated with density structures at different spatial scales in the filamentary infrared dark cloud, G034.43+00.24 (G34). This study makes use of the H13CO+ (1-0) molecular line data from the ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS) survey, which has spatial and velocity resolution of 0.04 pc and 0.2 km/s, respectively. Several tens of dendrogram structures have been extracted in the position-position-velocity space of H13CO+, which include 21 small-scale leaves and 20 larger-scale branches. Overall, their gas motions are supersonic but they exhibit the interesting behavior where leaves tend to be less dynamically supersonic than the branches. For the larger-scale, branch structures, the observed velocity-size relation (i.e., velocity variation/dispersion versus size) are seen to follow the Larson scaling exponent while the smaller-scale, leaf structures show a systematic deviation and display a steeper slope. We argue that the origin of the observed kinematics of the branch structures is likely to be a combination of turbulence and gravity-driven ordered gas flows. In comparison, gravity-driven chaotic gas motion is likely at the level of small-scale leaf structures. The results presented in our previous paper and this current follow-up study suggest that the main driving mechanism for mass accretion/inflow observed in G34 varies at different spatial scales. We therefore conclude that a scale-dependent combined effect of turbulence and gravity is essential to explain the star-formation processes in G34. △ Less

Submitted 22 February, 2022; originally announced February 2022.

Comments: 11 pages, 6 figures, and 1 table. To appear in MNRAS

Atomic Shocks in the Outflow of L1551 IRS 5 Identified with SOFIA-upGREAT Observations of [OI]

Authors: Yao-Lun Yang, Neal J. Evans II, Agata Karska, Lars E. Kristensen, Rebeca Aladro, Jon P. Ramsey, Joel D. Green, Jeong-Eun Lee

Abstract: We present velocity resolved SOFIA/upGREAT observations of [OI] and [CII] lines toward a Class I protostar, L1551 IRS 5, and its outflows. The SOFIA observations detect [OI] emission toward only the protostar and [CII] emission toward the protostar and the red-shifted outflow. The [OI] emission has a width of $\sim$100 km s$^{-1}$ only in the blue-shifted velocity, suggesting an origin in shocked… ▽ More We present velocity resolved SOFIA/upGREAT observations of [OI] and [CII] lines toward a Class I protostar, L1551 IRS 5, and its outflows. The SOFIA observations detect [OI] emission toward only the protostar and [CII] emission toward the protostar and the red-shifted outflow. The [OI] emission has a width of $\sim$100 km s$^{-1}$ only in the blue-shifted velocity, suggesting an origin in shocked gas. The [CII] lines are narrow, consistent with an origin in a photodissociation region. Differential dust extinction from the envelope due to the inclination of the outflows is the most likely cause of the missing red-shifted [OI] emission. Fitting the [OI] line profile with two Gaussian components, we find one component at the source velocity with a width of $\sim$20 km s$^{-1}$ and another extremely broad component at -30 km s$^{-1}$ with a width of 87.5 km s$^{-1}$, the latter of which has not been seen in L1551 IRS 5. The kinematics of these two components resemble cavity shocks in molecular outflows and spot shocks in jets. Radiative transfer calculations of the [OI], high-J CO, and H$_2$O lines in the cavity shocks indicate that [OI] dominates the oxygen budget, making up more than 70% of the total gaseous oxygen abundance and suggesting [O]/[H] of $\sim$1.5$\times$10$^{-4}$. Attributing the extremely broad [OI] component to atomic winds, we estimate the intrinsic mass loss rate of (1.3$\pm$0.8)$\times$ 10$^{-6}$ M$_{\odot}$ yr$^{-1}$. The intrinsic mass loss rates derived from low-J CO, [OI], and HI are similar, supporting the model of momentum-conserving outflows, where the atomic wind carries most momentum and drives the molecular outflows. △ Less

Submitted 22 November, 2021; originally announced November 2021.

Comments: 19 pages, 14 figures, accepted for publication in The Astrophysical Journal

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- IV. Radio Recombination Lines and evolution of star formation efficiencies

Authors: C. Zhang, Neal J. Evans II, T. Liu, J. -W. Wu, Ke Wang, H. -L. Liu, F. -Y. Zhu, Z. -Y. Ren, L. K. Dewangan, Chang Won Lee, Shanghuo Li, L. Bronfman, A. Tej, D. Li

Abstract: We report detection of radio recombination line (RRL) H$_{40α}$ toward 75 sources, with data obtained from ACA observations in the ATOMS survey of 146 active Galactic star forming regions. We calculated ionized gas mass and star formation rate with H40U line emission. The mass of ionized gas is significantly smaller than molecular gas mass, indicating that ionized gas is negligible in the star for… ▽ More We report detection of radio recombination line (RRL) H$_{40α}$ toward 75 sources, with data obtained from ACA observations in the ATOMS survey of 146 active Galactic star forming regions. We calculated ionized gas mass and star formation rate with H40U line emission. The mass of ionized gas is significantly smaller than molecular gas mass, indicating that ionized gas is negligible in the star forming clumps of the ATOMS sample. The star formation rate (SFR$_{{\rm H}_{40α}}$) estimated with RRL H$_{40α}$ agrees well with that (SFR$_{\rm L_{bol}}$) calculated with the total bolometric luminosity (L$_{\rm bol}$) when SFR $\gtrsim 5 {\rm M_\odot My}r^{-1}$, suggesting that millimeter RRLs could well sample the upper part of the initial mass function (IMF) and thus be good tracers for SFR. We also study the relationships between L$_{\rm bol}$ and the molecular line luminosities (L0mol) of CS J=2-1 and HC$_3$N J=11-10 for all the 146 ATOMS sources. The Lbol-L0mol correlations of both the CS J=2-1 and HC3N J=11-10 lines appear approximately linear and these transitions have success in predicting L$_{\rm bol}$ similar to that of more commonly used transitions. The L$_{\rm bol}$-to-L$_{\rm mol}$ ratios or SFR-to-mass ratios (star formation efficiency; SFE) do not change with galactocentric distances (R$_{\rm GC}$). Sources with H$_{40α}$ emission (or H$_{\rm II}$ regions) show higher L$_{\rm bol}$-to-L$_{\rm mol}$ than those without H$_{40α}$ emission, which may be an evolutionary effect. △ Less

Submitted 29 September, 2021; originally announced October 2021.

TIMES II: Investigating the Relation Between Turbulence and Star-forming Environments in Molecular Clouds

Authors: Hyeong-Sik Yun, Jeong-Eun Lee, Neal J. Evans II, Stella S. R. Offner, Mark H. Heyer, Jungyeon Cho, Brandt A. L. Gaches, Yao-Lun Yang, How-Huan Chen, Yunhee Choi, Yong-Hee Lee, Giseon Baek, Minho Choi, Jongsoo Kim, Hyunwoo Kang, Seokho Lee, Ken'ichi Tatematsu

Abstract: We investigate the effect of star formation on turbulence in the Orion A and Ophiuchus clouds using principal component analysis (PCA). We measure the properties of turbulence by applying PCA on the spectral maps in $^{13}$CO, C$^{18}$O, HCO$^+$ $J=$1$-$0, and CS $J=$2$-$1. First, the scaling relations derived from PCA of the $^{13}$CO maps show that the velocity difference ($δv$) for a given spat… ▽ More We investigate the effect of star formation on turbulence in the Orion A and Ophiuchus clouds using principal component analysis (PCA). We measure the properties of turbulence by applying PCA on the spectral maps in $^{13}$CO, C$^{18}$O, HCO$^+$ $J=$1$-$0, and CS $J=$2$-$1. First, the scaling relations derived from PCA of the $^{13}$CO maps show that the velocity difference ($δv$) for a given spatial scale ($L$) is the highest in the integral shaped filament (ISF) and L1688, where the most active star formation occurs in the two clouds. The $δv$ increases with the number density and total bolometric luminosity of the protostars in the sub-regions. Second, in the ISF and L1688 regions, the $δv$ of C$^{18}$O, HCO$^+$, and CS are generally higher than that of $^{13}$CO, which implies that the dense gas is more turbulent than the diffuse gas in the star-forming regions; stars form in dense gas, and dynamical activities associated with star formation, such as jets and outflows, can provide energy into the surrounding gas to enhance turbulent motions. △ Less

Submitted 28 July, 2021; originally announced July 2021.

Comments: 29 pages, 17 figures, 7 tables, accepted for publication in ApJ

Planck Galactic Cold Clumps at High Galactic Latitude-A Study with CO Lines

Authors: Fengwei Xu, Yuefang Wu, Tie Liu, Xunchuan Liu, Chao Zhang, Jarken Esimbek, Sheng-Li Qin, Di Li, Ke Wang, Jinghua Yuan, Fanyi Meng, Tianwei Zhang, David Eden, K. Tatematsu, Neal J. Evans, Paul. F. Goldsmith, Qizhou Zhang, C. Henkel, Hee-Weon Yi, Jeong-Eun Lee, Mika Saajasto, Gwangeong Kim, Mika Juvela, Dipen Sahu, Shin-Ying Hsu , et al. (6 additional authors not shown)

Abstract: Gas at high Galactic latitude is a relatively little-noticed component of the interstellar medium. In an effort to address this, forty-one Planck Galactic Cold Clumps at high Galactic latitude (HGal; $|b|>25^{\circ}$) were observed in $^{12}$CO, $^{13}$CO and C$^{18}$O J=1-0 lines, using the Purple Mountain Observatory 13.7-m telescope. $^{12}$CO (1-0) and $^{13}$CO (1-0) emission was detected in… ▽ More Gas at high Galactic latitude is a relatively little-noticed component of the interstellar medium. In an effort to address this, forty-one Planck Galactic Cold Clumps at high Galactic latitude (HGal; $|b|>25^{\circ}$) were observed in $^{12}$CO, $^{13}$CO and C$^{18}$O J=1-0 lines, using the Purple Mountain Observatory 13.7-m telescope. $^{12}$CO (1-0) and $^{13}$CO (1-0) emission was detected in all clumps while C$^{18}$O (1-0) emission was only seen in sixteen clumps. The highest and average latitudes are $71.4^{\circ}$ and $37.8^{\circ}$, respectively. Fifty-one velocity components were obtained and then each was identified as a single clump. Thirty-three clumps were further mapped at 1$^\prime$ resolution and 54 dense cores were extracted. Among dense cores, the average excitation temperature $T_{\mathrm{ex}}$ of $^{12}$CO is 10.3 K. The average line widths of thermal and non-thermal velocity dispersions are $0.19$ km s$^{-1}$ and $0.46$ km s$^{-1}$ respectively, suggesting that these cores are dominated by turbulence. Distances of the HGal clumps given by Gaia dust reddening are about $120-360$ pc. The ratio of $X_{13}$/$X_{18}$ is significantly higher than that in the solar neighbourhood, implying that HGal gas has a different star formation history compared to the gas in the Galactic disk. HGal cores with sizes from $0.01-0.1$ pc show no notable Larson's relation and the turbulence remains supersonic down to a scale of slightly below $0.1$ pc. None of the HGal cores which bear masses from 0.01-1 $M_{\odot}$ are gravitationally bound and all appear to be confined by outer pressure. △ Less

Submitted 20 July, 2021; v1 submitted 17 July, 2021; originally announced July 2021.

Comments: 35 pages, 13 figures

Which Molecular Cloud Structures Are Bound?

Authors: Neal J. Evans II, Mark Heyer. Marc-Antoine Miville-Deschênes, Quang Nguyen-Luong. Manuel Merello

Abstract: We analyze surveys of molecular cloud structures defined by tracers ranging from CO $J = 1-0$ through $^{13}$CO $J = 1-0$ to dust emission together with NH$_3$ data. The mean value of the virial parameter and the fraction of mass in bound structures depends on the method used to identify structures. Generally, the virial parameter decreases and the fraction of mass in bound structures increases wi… ▽ More We analyze surveys of molecular cloud structures defined by tracers ranging from CO $J = 1-0$ through $^{13}$CO $J = 1-0$ to dust emission together with NH$_3$ data. The mean value of the virial parameter and the fraction of mass in bound structures depends on the method used to identify structures. Generally, the virial parameter decreases and the fraction of mass in bound structures increases with the effective density of the tracer, the surface density and mass of the structures, and the distance from the center of a galaxy. For the most complete surveys of structures in the Galaxy defined by CO $J = 1-0$, the fraction of mass that is in bound structures is 0.19. For catalogs of other galaxies based on CO $J = 2-1$, the fraction is 0.35. These results offer substantial alleviation of the fundamental problem of slow star formation. If only clouds found to be bound are counted and they are assumed to collapse in a free-fall time at their mean cloud density, the sum over all clouds in a complete survey of the Galaxy yields a predicted star formation rate of 46 solar masses per year, a factor of 6.5 less than if all clouds are bound. △ Less

Submitted 12 July, 2021; originally announced July 2021.

Comments: 26 pages, 14 figures, accepted to ApJ

TIMES I: a Systematic Observation in Multiple Molecular Lines Toward the Orion A and Ophiuchus Clouds

Authors: Hyeong-Sik Yun, Jeong-Eun Lee, Yunhee Choi, Neal J. Evans II, Stella S. R. Offner, Mark H. Heyer, Brandt A. L. Gaches, Yong-Hee Lee, Giseon Baek, Minho Choi, Hyunwoo Kang, Seokho Lee, Ken'ichi Tatematsu, Yao-Lun Yang, How-Huan Chen, Youngung Lee, Jae Hoon Jung, Changhoon Lee, Jungyeon Cho

Abstract: We have used the Taeduk Radio Astronomy Observatory to observe the Orion A and Ophiuchus clouds in the $J=$1$-$0 lines of $^{13}$CO, C$^{18}$O, HCN, HCO$^+$, and N$_2$H$^+$ and the $J=$2$-$1 line of CS. The fully sampled maps with uniform noise levels are used to create moment maps. The variations of the line intensity and velocity dispersion with total column density, derived from dust emission m… ▽ More We have used the Taeduk Radio Astronomy Observatory to observe the Orion A and Ophiuchus clouds in the $J=$1$-$0 lines of $^{13}$CO, C$^{18}$O, HCN, HCO$^+$, and N$_2$H$^+$ and the $J=$2$-$1 line of CS. The fully sampled maps with uniform noise levels are used to create moment maps. The variations of the line intensity and velocity dispersion with total column density, derived from dust emission maps, are presented and compared to previous work. The CS line traces dust column density over more than one order of magnitude, and the N$_2$H$^+$ line best traces the highest column density regime ($\log(N_\mathrm{H_2}$) $>$ 22.8). Line luminosities, integrated over the cloud, are compared to those seen in other galaxies. The HCO$^+$-to-HCN luminosity ratio in the Orion A cloud is similar to that of starburst galaxies, while that in the Ophiuchus cloud is in between those of active galactic nuclei and starburst galaxies. △ Less

Submitted 8 June, 2021; originally announced June 2021.

Comments: 62 pages, 53 figures, 6 tables, accepted for publication in ApJS

Molecular Cloud Cores with High Deuterium Fractions: Nobeyama Mapping Survey

Authors: Ken'ichi Tatematsu, Gwanjeong Kim, Tie Liu, Neal J. Evans II, Hee-Weon Yi, Jeong-Eun Lee, Yuefang Wu, Naomi Hirano, Sheng-Yuan Liu, Somnath Dutta, Dipen Sahu, Patricio Sanhueza, Kee-Tae Kim, Mika Juvela, L. Viktor T'oth, Orsolya Feh'er, Jinhua He, J. X. Ge, Siyi Feng, Minho Choi, Miju Kang, Mark A. Thompson, Gary A. Fuller, Di Li, Isabelle Ristorcelli , et al. (13 additional authors not shown)

Abstract: We present the results of on-the-fly mapping observations of 44 fields containing 107 SCUBA-2 cores in the emission lines of molecules, N$_2$H$^+$, HC$_3$N, and CCS at 82$-$94 GHz using the Nobeyama 45-m telescope. This study aimed at investigating the physical properties of cores that show high deuterium fractions and might be close to the onset of star formation. We found that the distributions… ▽ More We present the results of on-the-fly mapping observations of 44 fields containing 107 SCUBA-2 cores in the emission lines of molecules, N$_2$H$^+$, HC$_3$N, and CCS at 82$-$94 GHz using the Nobeyama 45-m telescope. This study aimed at investigating the physical properties of cores that show high deuterium fractions and might be close to the onset of star formation. We found that the distributions of the N$_2$H$^+$ and HC$_3$N line emissions are approximately similar to that of 850-$μ$m dust continuum emission, whereas the CCS line emission is often undetected or is distributed in a clumpy structure surrounding the peak position of the 850-$μ$m dust continuum emission. Occasionally (12%), we observe the CCS emission which is an early-type gas tracer toward the young stellar object, probably due to local high excitation. Evolution toward star formation does not immediately affect nonthermal velocity dispersion. △ Less

Submitted 7 June, 2021; originally announced June 2021.

Comments: 45 pages, 15 figures, ApJS, in press

Journal ref: Astrophys.J.Suppl. 256 (2021) 25

ATOMS:ALMA Three-millimeter Observations of Massive Star-forming regions -- III :Catalogues of candidate hot molecular cores and Hyper/Ultra compact HII regions

Authors: Hong-Li Liu, Tie Liu, Neal J. Evans, Ke Wang, Guido Garay, Sheng-Li Qin, Shanghuo Li, Amelia Stutz, Paul F. Goldsmith, Sheng-Yuan Liu, Anandmayee Tej, Qizhou Zhang, Mika Juvela, Di Li, Jun-Zhi Wang, Leonardo Bronfman, Zhiyuan Ren, Yue-Fang Wu, Kee-Tae Kim, Chang-Won Lee, Kenichi Tatematsu, Maria. R. Cunningham, Xun-Chuan Liu, Jing-Wen Wu, Tomoya Hirota , et al. (29 additional authors not shown)

Abstract: We have identified 453 compact dense cores in 3 mm continuum emission maps in the ATOMS (ALMA Three-millimeter Observations of Massive Star-forming regions) survey, and compiled three catalogues of high-mass star forming cores. One catalogue, referred to as H/UC-HII catalogue, includes 89 cores that enshroud hyper/ultra compact (H/UC) HII regions as characterized by associated compact H40alpha emi… ▽ More We have identified 453 compact dense cores in 3 mm continuum emission maps in the ATOMS (ALMA Three-millimeter Observations of Massive Star-forming regions) survey, and compiled three catalogues of high-mass star forming cores. One catalogue, referred to as H/UC-HII catalogue, includes 89 cores that enshroud hyper/ultra compact (H/UC) HII regions as characterized by associated compact H40alpha emission. A second catalogue, referred to as pure s-cHMC, includes 32 candidate Hot Molecular Cores (HMCs) showing rich spectra (N>20lines) of complex organic molecules (COMs) but not associated with H/UC-HII regions. The third catalogue, referred to as pure w-cHMC, includes 58 candidate HMCs with relatively low levels of COM richness and not associated with H/UC-HII regions. These three catalogues of dense cores provide an important foundation for future studies of the early stages of high-mass star formation across the Milky Way. We also find that nearly half of H/UC-HII cores are candidate HMCs. From the number counts of COM-containing and H/UC-HII cores, we suggest that the duration of high-mass protostellar cores showing chemically rich features is at least comparable to the lifetime of H/UC-HII regions. For cores in the H/UC-HII catalogue, the width of the H40alpha line increases as the core size decreases, suggesting that the non-thermal dynamical and/or pressure line-broadening mechanisms dominate on the smaller scales of the H/UC-HII cores. △ Less

Submitted 7 May, 2021; originally announced May 2021.

Comments: 17 pages, five tables, and 11 figures. Accepted for publication at MNRAS

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

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

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

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

Comments: Published in The Astrophysical Journal Letters (ApJL)

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

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

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

Submitted 27 October, 2020; originally announced October 2020.

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

Molecular Cloud Cores with High Deuterium Fraction: Nobeyama Single-Pointing Survey

Authors: Gwanjeong Kim, Kenichi Tatematsu, Tie Liu, Miss Hee-Weon Yi, Jinhua He, Naomi Hirano, Sheng-Yuan Liu, Minho Choi, Patricio Sanhueza, L. Viktor Toth, Neal J. Evans, Siyi Feng, Mika Juvela, Kee-Tae Kim, Charlotte Vastel, Jeong-Eun Lee, Quang Nguyn-Lu'o'ng, Miju Kang, Isabelle Ristorcelli, O. Fehér, Yuefang Wu, Satoshi Ohashi, Ke Wang, Ryo Kandori, Tomoya Hirota , et al. (7 additional authors not shown)

Abstract: We present the results of a single-pointing survey of 207 dense cores embedded in Planck Galactic Cold Clumps distributed in five different environments ($λ$ Orionis, Orion A, B, Galactic plane, and high latitudes) to identify dense cores on the verge of star formation for the study of the initial conditions of star formation. We observed these cores in eight molecular lines at 76-94 GHz using the… ▽ More We present the results of a single-pointing survey of 207 dense cores embedded in Planck Galactic Cold Clumps distributed in five different environments ($λ$ Orionis, Orion A, B, Galactic plane, and high latitudes) to identify dense cores on the verge of star formation for the study of the initial conditions of star formation. We observed these cores in eight molecular lines at 76-94 GHz using the Nobeyama 45-m telescope. We find that early-type molecules (e.g., CCS) have low detection rates and that late-type molecules (e.g., N$_2$H$^+$, c-C$_3$H$_2$) and deuterated molecules (e.g., N$_2$D$^+$, DNC) have high detection rates, suggesting that most of the cores are chemically evolved. The deuterium fraction (D/H) is found to decrease with increasing distance, indicating that it suffers from differential beam dilution between the D/H pair of lines for distant cores ($>$1 kpc). For $λ$ Orionis, Orion A, and B located at similar distances, D/H is not significantly different, suggesting that there is no systematic difference in the observed chemical properties among these three regions. We identify at least eight high D/H cores in the Orion region and two at high latitudes, which are most likely to be close to the onset of star formation. There is no clear evidence of the evolutionary change in turbulence during the starless phase, suggesting that the dissipation of turbulence is not a major mechanism for the beginning of star formation as judged from observations with a beam size of 0.04 pc. △ Less

Submitted 23 July, 2020; originally announced July 2020.

Comments: 55 pages, 12 Figures, 8 Tables, accepted for publication in ApJS

LEGO II: A 3 mm molecular line study covering 100 pc of one of the most actively star-forming portions within the Milky Way Disc

Authors: A. T. Barnes, J. Kauffmann, F. Bigiel, N. Brinkmann, D. Colombo, A. E Guzmán, W. J. Kim, L. Szűcs, V. Wakelam, S. Aalto, T. Albertsson, N. J. Evans II, S. C. O. Glover, P. F. Goldsmith, C. Kramer, K. Menten, Y. Nishimura, S. Viti, Y. Watanabe, A. Weiss, M. Wienen, H. Wiesemeyer, F. Wyrowski

Abstract: The current generation of (sub)mm-telescopes has allowed molecular line emission to become a major tool for studying the physical, kinematic, and chemical properties of extragalactic systems, yet exploiting these observations requires a detailed understanding of where emission lines originate within the Milky Way. In this paper, we present 60$^{\prime\prime}$ ($\sim$3pc) resolution observations of… ▽ More The current generation of (sub)mm-telescopes has allowed molecular line emission to become a major tool for studying the physical, kinematic, and chemical properties of extragalactic systems, yet exploiting these observations requires a detailed understanding of where emission lines originate within the Milky Way. In this paper, we present 60$^{\prime\prime}$ ($\sim$3pc) resolution observations of many 3mm-band molecular lines across a large map of the W49 massive star-forming region ($\sim$100$\times$100pc at 11kpc), which were taken as part of the "LEGO" IRAM-30m large project. We find that the spatial extent or brightness of the molecular line transitions are not well correlated with their critical densities, highlighting abundance and optical depth must be considered when estimating line emission characteristics. We explore how the total emission and emission efficiency (i.e. line brightness per H$_{2}$ column density) of the line emission vary as a function of molecular hydrogen column density and dust temperature. We find that there is not a single region of this parameter space responsible for the brightest and most efficiently emitting gas for all species. For example, we find that the HCN transition shows high emission efficiency at high column density ($10^{22}$cm$^{-2}$) and moderate temperatures (35K), whilst e.g. N$_2$H$^+$ emits most efficiently towards lower temperatures ($10^{22}$cm$^{-2}$; <20K). We determine $X_{\mathrm{CO} (1-0)} \sim 0.3 \times 10^{20} \mathrm{cm^{-2}(Kkms^{-1})^{-1}}$, and $α_{\mathrm{HCN} (1-0)} \sim 30\mathrm{M_\odot(Kkms^{-1}pc^2)^{-1}}$, which both differ significantly from the commonly adopted values. In all, these results suggest caution should be taken when interpreting molecular line emission. △ Less

Submitted 3 September, 2020; v1 submitted 21 July, 2020; originally announced July 2020.

Comments: 33+2 pages, 17+3 figures, 4+4 tables. Tables available online in a machine-readable format. Accepted for publication in MNRAS

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- I. Survey description and a first look at G9.62+0.19

Authors: Tie Liu, Neal J. Evans, Kee-Tae Kim, Paul F. Goldsmith, Sheng-Yuan Liu, Qizhou Zhang, Kenichi Tatematsu, Ke Wang, Mika Juvela, Leonardo Bronfman, Maria. R. Cunningham, Guido Garay, Tomoya Hirota, Jeong-Eun Lee, Sung-Ju Kang, Di Li, Pak-Shing Li, Diego Mardones, Sheng-Li Qin, Isabelle Ristorcelli, Anandmayee Tej, L. Viktor Toth, Jing-Wen Wu, Yue-Fang Wu, Hee-weon Yi , et al. (20 additional authors not shown)

Abstract: The "ATOMS," standing for {\it ALMA Three-millimeter Observations of Massive Star-forming regions}, survey has observed 146 active star forming regions with ALMA Band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary str… ▽ More The "ATOMS," standing for {\it ALMA Three-millimeter Observations of Massive Star-forming regions}, survey has observed 146 active star forming regions with ALMA Band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the "ATOMS" survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS $J = 2-1$, HCO$^+$ $J = 1-0$ and HCN $J = 1-0$, are found to show extended gas emission in low density regions within the clump; less than 25\% of their emission is from dense cores. SO, CH$_3$OH, H$^{13}$CN and HC$_3$N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over $\sim$1 pc), which may be caused by slow shocks from large--scale colliding flows or H{\sc ii} regions. Stellar feedback from an expanding H{\sc ii} region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analyzed with other survey data, e.g., "MALT90", "Orion B", "EMPIRE", "ALMA\_IMF", and "ALMAGAL", to deepen our understandings of "dense gas" star formation scaling relations and massive proto-cluster formation. △ Less

Submitted 8 July, 2020; v1 submitted 2 June, 2020; originally announced June 2020.

Comments: published on MNRAS

Journal ref: MNRAS, 496, 2790 (2020)

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- II. Compact objects in ACA observations and star formation scaling relations

Authors: Tie Liu, Neal J. Evans, Kee-Tae Kim, Pail F. Goldsmith, Sheng-Yuan Liu, Qizhou Zhang, Kenichi Tatematsu, Ke Wang, Mika Juvela, Leonardo Bronfman, Maria. R. Cunningham, Guido Garay, Tomoya Hirota, Jeong-Eun Lee, Sung-Ju Kang, Di Li, Pak-Shing Li, Diego Mardones, Sheng-Li Qin, Isabelle Ristorcelli, Anandmayee Tej, L. Viktor Toth, Jing-Wen Wu, Yue-Fang Wu, Hee-weon Yi , et al. (19 additional authors not shown)

Abstract: We report studies of the relationships between the total bolometric luminosity ($L_{\rm bol}$ or $L_{\rm TIR}$) and the molecular line luminosities of $J=1-0$ transitions of H$^{13}$CN, H$^{13}$CO$^+$, HCN, and HCO$^+$ with data obtained from ACA observations in the "ATOMS" survey of 146 active Galactic star forming regions. The correlations between $L_{\rm bol}$ and molecular line luminosities… ▽ More We report studies of the relationships between the total bolometric luminosity ($L_{\rm bol}$ or $L_{\rm TIR}$) and the molecular line luminosities of $J=1-0$ transitions of H$^{13}$CN, H$^{13}$CO$^+$, HCN, and HCO$^+$ with data obtained from ACA observations in the "ATOMS" survey of 146 active Galactic star forming regions. The correlations between $L_{\rm bol}$ and molecular line luminosities $L'_{\rm mol}$ of the four transitions all appear to be approximately linear. Line emission of isotopologues shows as large scatters in $L_{\rm bol}$-$L'_{\rm mol}$ relations as their main line emission. The log($L_{\rm bol}$/$L'_{\rm mol}$) for different molecular line tracers have similar distributions. The $L_{\rm bol}$-to-$L'_{\rm mol}$ ratios do not change with galactocentric distances ($R_{\rm GC}$) and clump masses ($M_{\rm clump}$). The molecular line luminosity ratios (HCN-to-HCO$^+$, H$^{13}$CN-to-H$^{13}$CO$^+$, HCN-to-H$^{13}$CN and HCO$^+$-to-H$^{13}$CO$^+$) all appear constant against $L_{\rm bol}$, dust temperature ($T_{\rm d}$), $M_{\rm clump}$ and $R_{\rm GC}$. Our studies suggest that both the main lines and isotopologue lines are good tracers of the total masses of dense gas in Galactic molecular clumps. The large optical depths of main lines do not affect the interpretation of the slopes in star formation relations. We find that the mean star formation efficiency (SFE) of massive Galactic clumps in the "ATOMS" survey is reasonably consistent with other measures of the SFE for dense gas, even those using very different tracers or examining very different spatial scales. △ Less

Submitted 8 July, 2020; v1 submitted 25 May, 2020; originally announced May 2020.

Comments: Published on MNRAS. The full tables are included in Tables.pdf or Tables.tex files, which can be downloaded from source files

Journal ref: MNRAS, 496, 2821 (2020)

Cloud structures in M17 SWex : Possible cloud-cloud collision

Authors: Shinichi W. Kinoshita, Fumitaka Nakamura, Quang Nguyen-Luong, Benjamin Wu, Tomomi Shimoikura, Koji Sugitani, Kazuhito Dobashi, Hideaki Takemura, Patricio Sanhueza, Kee-Tae Kim, Hyunwoo Kang, Neal J. Evans, Glenn J. White, Cassandra Fallscheer

Abstract: Using wide-field $^{13}$CO ($J=1-0$) data taken with the Nobeyama 45-m telescope, we investigate cloud structures of the infrared dark cloud complex in M17 with SCIMES. In total, we identified 118 clouds that contain 11 large clouds with radii larger than 1 pc. The clouds are mainly distributed in the two representative velocity ranges of 10 $-$ 20 km s$^{-1}$ and 30 $-$ 40 km s$^{-1}$. By compari… ▽ More Using wide-field $^{13}$CO ($J=1-0$) data taken with the Nobeyama 45-m telescope, we investigate cloud structures of the infrared dark cloud complex in M17 with SCIMES. In total, we identified 118 clouds that contain 11 large clouds with radii larger than 1 pc. The clouds are mainly distributed in the two representative velocity ranges of 10 $-$ 20 km s$^{-1}$ and 30 $-$ 40 km s$^{-1}$. By comparing with the ATLASGAL catalog, we found that the majority of the $^{13}$CO clouds with 10 $-$ 20 km s$^{-1}$ and 30 $-$ 40 km s$^{-1}$ are likely located at distances of 2 kpc (Sagittarius arm) and 3 kpc (Scutum arm), respectively. Analyzing the spatial configuration of the identified clouds and their velocity structures, we attempt to reveal the origin of the cloud structure in this region. Here we discuss three possibilities: (1) overlapping with different velocities, (2) cloud oscillation, and (3) cloud-cloud collision. From the position-velocity diagrams, we found spatially-extended faint emission between $\sim$ 20 km s$^{-1}$ and $\sim$ 35 km s$^{-1}$, which is mainly distributed in the spatially-overlapped areas of the clouds. We also found that in some areas where clouds with different velocities overlapped, the magnetic field orientation changes abruptly. The distribution of the diffuse emission in the position-position-velocity space and the bending magnetic fields appear to favor the cloud-cloud collision scenario compared to other scenarios. In the cloud-cloud collision scenario, we propose that two $\sim$35 km s$^{-1}$ foreground clouds are colliding with clouds at $\sim$20 km s$^{-1}$ with a relative velocity of 15 km s$^{-1}$. These clouds may be substructures of two larger clouds having velocities of $\sim$ 35 km s$^{-1}$ ($\gtrsim 10^3 $ M$_{\odot}$) and $\sim$ 20 km s$^{-1}$ ($\gtrsim 10^4 $ M$_{\odot}$), respectively. △ Less

Submitted 23 May, 2020; v1 submitted 21 May, 2020; originally announced May 2020.

Comments: 34 pages, 13 figures, Accepted for publication in PASJ

ALMA ACA and Nobeyama observations of two Orion cores in deuterated molecular lines

Authors: Ken'ichi Tatematsu, Tie Liu, Gwanjeong Kim, Hee-Weon Yi, Jeong-Eun Lee, Naomi Hirano, Sheng-Yuan Liu, Satoshi Ohashi, Patricio Sanhueza, James Di Francesco, Neal J. Evans II, Gary A. Fuller, Ryo Kandori, Minho Choi, Miju Kang, Siyi Feng, Tomoya Hirota, Takeshi Sakai, Xing Lu, Quang Nguyen Lu'o'ng, Mark A. Thompson, Yuefang Wu, Di Li, Kee-Tae Kim, Ke Wang , et al. (3 additional authors not shown)

Abstract: We mapped two molecular cloud cores in the Orion A cloud with the ALMA ACA 7-m Array and with the Nobeyama 45-m radio telescope. These cores have bright N$_2$D$^+$ emission in single-pointing observations with the Nobeyama 45-m radio telescope, have relatively high deuterium fraction, and are thought to be close to the onset of star formation. One is a star-forming core, and the other is starless.… ▽ More We mapped two molecular cloud cores in the Orion A cloud with the ALMA ACA 7-m Array and with the Nobeyama 45-m radio telescope. These cores have bright N$_2$D$^+$ emission in single-pointing observations with the Nobeyama 45-m radio telescope, have relatively high deuterium fraction, and are thought to be close to the onset of star formation. One is a star-forming core, and the other is starless. These cores are located along filaments observed in N$_2$H$^+$, and show narrow linewidths of 0.41 km s$^{-1}$ and 0.45 km s$^{-1}$ in N$_2$D$^+$, respectively, with the Nobeyama 45-m telescope. Both cores were detected with the ALMA ACA 7m Array in the continuum and molecular lines at Band 6. The starless core G211 shows clumpy structure with several sub-cores, which in turn show chemical differences. Also, the sub-cores in G211 have internal motions that are almost purely thermal. The starless sub-core G211D, in particular, shows a hint of the inverse P Cygni profile, suggesting infall motion. The star-forming core G210 shows an interesting spatial feature of two N$_2$D$^+$ peaks of similar intensity and radial velocity located symmetrically with respect to the single dust continuum peak. One interpretation is that the two N$_2$D$^+$ peaks represent an edge-on pseudo-disk. The CO outflow lobes, however, are not directed perpendicular to the line connecting both N$_2$D$^+$ peaks. △ Less

Submitted 26 April, 2020; originally announced April 2020.

Comments: 31 pages, 17 figures, 9 tables, to appear in ApJ (accepted April 23, 2020)

Journal ref: Astrophys.J. 895 (2020) 119

Star Formation Occurs in Dense Gas, but What Does "Dense" Mean?

Authors: Neal J. Evans II, Kee-Tae Kim, Jingwen Wu, Zhang Chao, Mark Heyer, Tie Liu, Quang Nguyen-Luong, Jens Kauffmann

Abstract: We report results of a project to map HCN and HCO+ J = 1-0 emission toward a sample of molecular clouds in the inner Galaxy, all containing dense clumps that are actively engaged in star formation. We compare these two molecular line tracers with millimeter continuum emission and extinction, as inferred from 13CO, as tracers of dense gas in molecular clouds. The fraction of the line luminosity fro… ▽ More We report results of a project to map HCN and HCO+ J = 1-0 emission toward a sample of molecular clouds in the inner Galaxy, all containing dense clumps that are actively engaged in star formation. We compare these two molecular line tracers with millimeter continuum emission and extinction, as inferred from 13CO, as tracers of dense gas in molecular clouds. The fraction of the line luminosity from each tracer that comes from the dense gas, as measured by AV > 8 mag, varies substantially from cloud to cloud. In all cases, a substantial fraction (in most cases, the majority) of the total luminosity arises in gas below the AV > 8 mag threshold and outside the region of strong mm continuum emission. Measurements of the luminosity of HCN toward other galaxies will likely be dominated by such gas at lower surface density. Substantial, even dominant, contributions to the total line luminosity can arise in gas with densities typical of the cloud as a whole (densities about 100 per cubic cm). Defining the dense clump from the HCN or HCO+ emission itself, similarly to previous studies, leads to a wide range of clump properties, with some being considerably larger and less dense than in previous studies. HCN and HCO+ have similar ability to trace dense gas for the clouds in this sample. For the two clouds with low virial parameters, the 13CO is definitely a worse tracer of the dense gas, but for the other four, it is equally good (or bad) at tracing dense gas. △ Less

Submitted 14 April, 2020; originally announced April 2020.

Comments: Accepted by the Astrophysical Journal,

Constraining the Infalling Envelope Models of Embedded Protostars: BHR 71 and its Hot Corino

Authors: Y. -L. Yang, N. J. Evans II, A. Smith, J. -E. Lee, J. J. Tobin, S. Terebey, H. Calcutt, J. K. Jorgensen, J. D. Green, T. L. Bourke

Abstract: The collapse of the protostellar envelope results in the growth of the protostar and the development of a protoplanetary disk, playing a critical role during the early stages of star formation. Characterizing the gas infall in the envelope constrains the dynamical models of star formation. We present unambiguous signatures of infall, probed by optically thick molecular lines, toward an isolated em… ▽ More The collapse of the protostellar envelope results in the growth of the protostar and the development of a protoplanetary disk, playing a critical role during the early stages of star formation. Characterizing the gas infall in the envelope constrains the dynamical models of star formation. We present unambiguous signatures of infall, probed by optically thick molecular lines, toward an isolated embedded protostar, BHR 71 IRS1. The three dimensional radiative transfer calculations indicate that a slowly rotating infalling envelope model following the "inside-out" collapse reproduces the observations of both HCO$^{+}$ $J=4\rightarrow3$ and CS $J=7\rightarrow6$ lines, and the low velocity emission of the HCN $J=4\rightarrow3$ line. The envelope has a model-derived age of 12000$\pm$3000 years after the initial collapse. The envelope model underestimates the high velocity emission at the HCN $J=4\rightarrow3$ and H$^{13}$CN $J=4\rightarrow3$ lines, where outflows or a Keplerian disk may contribute. The ALMA observations serendipitously discover the emission of complex organic molecules (COMs) concentrated within a radius of 100 au, indicating that BHR 71 IRS1 harbors a hot corino. Eight species of COMs are identified, including CH$_{3}$OH and CH$_{3}$OCHO, along with H$_{2}$CS, SO$_{2}$ and HCN $v_{2}=1$. The emission of methyl formate and $^{13}$C-methanol shows a clear velocity gradient within a radius of 50 au, hinting at an unresolved Keplerian rotating disk. △ Less

Submitted 4 February, 2020; originally announced February 2020.

Comments: 42 pages, 38 figures, accepted for publication in ApJ

Large-scale molecular gas distribution in the M17 cloud complex: dense gas conditions of massive star formation?

Authors: Quang Nguyen-Luong, Fumitaka Nakamura, Koji Sugitani, Tomomi Shimoikura, Kazuhito Dobashi, Shinichi W. Kinoshita, Kee-Tae Kim, Hynwoo Kang, Patricio Sanhueza, Neal J. Evans II, Glenn J. White

Abstract: The non-uniform distribution of gas and protostars in molecular clouds is caused by combinations of various physical processes that are difficult to separate. We explore this non-uniform distribution in the M17 molecular cloud complex that hosts massive star formation activity using the $^{12}$CO ($J=1-0$) and $^{13}$CO ($J=1-0$) emission lines obtained with the Nobeyama 45m telescope. Differences… ▽ More The non-uniform distribution of gas and protostars in molecular clouds is caused by combinations of various physical processes that are difficult to separate. We explore this non-uniform distribution in the M17 molecular cloud complex that hosts massive star formation activity using the $^{12}$CO ($J=1-0$) and $^{13}$CO ($J=1-0$) emission lines obtained with the Nobeyama 45m telescope. Differences in clump properties such as mass, size, and gravitational boundedness reflect the different evolutionary stages of the M17-H{\scriptsize II} and M17-IRDC clouds. Clumps in the M17-H{\scriptsize II} cloud are denser, more compact, and more gravitationally bound than those in M17-IRDC. While M17-H{\scriptsize II} hosts a large fraction of very dense gas (27\%) that has column density larger than the threshold of $\sim$ 1 g cm$^{-2}$ theoretically predicted for massive star formation, this very dense gas is deficient in M17-IRDC (0.46\%). Our HCO$^+$ ($J=1-0$) and HCN ($J=1-0$) observations with the TRAO 14m telescope, {\nlqb trace all gas with column density higher than $3\times 10^{22}$ cm$^{-2}$}, confirm the deficiency of high density ($\gtrsim 10^5$ cm$^{-3}$) gas in M17-IRDC. Although M17-IRDC is massive enough to potentially form massive stars, its deficiency of very dense gas and gravitationally bound clumps can explain the current lack of massive star formation. △ Less

Submitted 26 January, 2020; originally announced January 2020.

Comments: We welcome all comments to make the paper better. Please send comments to my email. Thank you

Magnetic fields in the infrared dark cloud G34.43+0.24

Authors: Archana Soam, Tie Liu, B-G Andersson, Chang Won Lee, Junhao Liu, Mika Juvela, Pak Shing Li, Paul F. Goldsmith, Qizhou Zhang, Patrick M. Koch, Kee-Tae Kim, Keping Qiu, Neal J. Evans II, SCOPE collaborations

Abstract: We present the B-fields mapped in IRDC G34.43+0.24 using 850\,$μ$m polarized dust emission observed with the POL-2 instrument at JCMT. We examine the magnetic field geometries and strengths in the northern, central, and southern regions of the filament. The overall field geometry is ordered and aligned closely perpendicular to the filament's main axis, particularly in regions containing the centra… ▽ More We present the B-fields mapped in IRDC G34.43+0.24 using 850\,$μ$m polarized dust emission observed with the POL-2 instrument at JCMT. We examine the magnetic field geometries and strengths in the northern, central, and southern regions of the filament. The overall field geometry is ordered and aligned closely perpendicular to the filament's main axis, particularly in regions containing the central clumps MM1 and MM2, whereas MM3 in the north has field orientations aligned with its major axis. The overall field orientations are uniform at large (POL-2 at 14$\arcsec$ and SHARP at 10$\arcsec$) to small scales (TADPOL at 2.5$\arcsec$ and SMA at 1.5$\arcsec$) in the MM1 and MM2 regions. SHARP/CSO observations in MM3 at 350\,$μ$m from Tang et al. show a similar trend as seen in our POL-2 observations. TADPOL observations demonstrate a well-defined field geometry in MM1/MM2 consistent with MHD simulations of accreting filaments. We obtained a plane-of-sky magnetic field strength of 470$\pm$190\,$μ$G, 100$\pm$40\,$μ$G, and 60$\pm$34\,$μ$G in the central, northern and southern regions of G34, respectively, using the updated Davis-Chandrasekhar-Fermi relation. The estimated value of field strength, combined with column density and velocity dispersion values available in the literature, suggests G34 to be marginally critical with criticality parameter $\rm λ$ values 0.8$\pm$0.4, 1.1$\pm$0.8, and 0.9$\pm$0.5 in the central, northern, and southern regions, respectively. The turbulent motions in G34 are sub-Alfvénic with Alfvénic Mach numbers of 0.34$\pm$0.13, 0.53$\pm$0.30, and 0.49$\pm$0.26 in the three regions. The observed aligned B-fields in G34.43+0.24 are consistent with theoretical models suggesting that B-fields play an important role in guiding the contraction of the cloud driven by gravity. △ Less

Submitted 9 August, 2019; originally announced August 2019.

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

The Disk Gas Mass and the Far-IR Revolution

Authors: Edwin A. Bergin, Klaus M. Pontoppidan, Charles M. Bradford, L. Ilsedore Cleeves, Neal J. Evans, Maryvonne Gerin, Paul F. Goldsmith, Quentin Kral, Gary J. Melnick, Melissa McClure, Karin Oberg, Thomas L. Roellig, Edward Wright, Richard Teague, Jonathan P. Williams, Ke Zhang

Abstract: The gaseous mass of protoplanetary disks is a fundamental quantity in planet formation. The presence of gas is necessary to assemble planetesimals, it determines timescales of giant planet birth, and it is an unknown factor for a wide range of properties of planet formation, from chemical abundances (X/H) to the mass efficiency of planet formation. The gas mass obtained from traditional tracers, s… ▽ More The gaseous mass of protoplanetary disks is a fundamental quantity in planet formation. The presence of gas is necessary to assemble planetesimals, it determines timescales of giant planet birth, and it is an unknown factor for a wide range of properties of planet formation, from chemical abundances (X/H) to the mass efficiency of planet formation. The gas mass obtained from traditional tracers, such as dust thermal continuum and CO isotopologues, are now known to have significant (1 - 2 orders of magnitude) discrepancies. Emission from the isotopologue of H2, hydrogen deuteride (HD), offers an alternative measurement of the disk gas mass. Of all of the regions of the spectrum, the far-infrared stands out in that orders of magnitude gains in sensitivity can be gleaned by cooling a large aperture telescope to 8 K. Such a facility can open up a vast new area of the spectrum to exploration. One of the primary benefits of this far-infrared revolution would be the ability to survey hundreds of planet-forming disks in HD emission to derive their gaseous masses. For the first time, we will have statistics on the gas mass as a function of evolution, tracing birth to dispersal as a function of stellar spectral type. These measurements have broad implications for our understanding of the time scale during which gas is available to form giant planets, the dynamical evolution of the seeds of terrestrial worlds, and the resulting chemical composition of pre-planetary embryos carrying the elements needed for life. Measurements of the ground-state line of HD requires a space-based observatory operating in the far-infrared at 112 microns. △ Less

Submitted 20 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey

Variability in the Assembly of Protostellar Systems

Authors: Joel D. Green, Yao-Lun Yang, Tom Megeath, Doug Johnstone, John Tobin, Sarah Sadavoy, Klaus Pontoppidan, Stella Offner, Neal J. Evans, Dan M. Watson, Jennifer Hatchell, Ian Stephens, Zhi-Yun Li, Jacob White, Robert A. Gutermuth, Will Fischer, Agata Karska, Jens Kauffmann, Mike Dunham, Hector Arce

Abstract: Understanding the collapse of clouds and the formation of protoplanetary disks is essential to understanding the formation of stars and planets. Infall and accretion, the mass-aggregation processes that occur at envelope and disk scales, drive the dynamical evolution of protostars. While the observations of protostars at different stages constrain their evolutionary tracks, the impact of variabili… ▽ More Understanding the collapse of clouds and the formation of protoplanetary disks is essential to understanding the formation of stars and planets. Infall and accretion, the mass-aggregation processes that occur at envelope and disk scales, drive the dynamical evolution of protostars. While the observations of protostars at different stages constrain their evolutionary tracks, the impact of variability due to accretion bursts on dynamical and chemical evolution of the source is largely unknown. The lasting effects on protostellar envelopes and disks are tracked through multi-wavelength and time domain observational campaigns, requiring deep X-ray, infrared, and radio imaging and spectroscopy, at a sufficient level of spatial detail to distinguish contributions from the various substructures (i.e., envelope from disk from star from outflow). Protostellar models derived from these campaigns will illuminate the initial chemical state of protoplanetary disks during the epoch of giant planet formation. Insight from individual star formation in the Milky Way is also necessary to understand star formation rates in extragalactic sources. This cannot be achieved with ground-based observatories and is not covered by currently approved instrumentation. Requirements: High (v < 10 km/s for survey; v < 1 km/s for followup) spectral resolution capabilities with relatively rapid response times in the IR (3-500 um), X-ray (0.1-10 keV), and radio (cm) are critical to follow the course of accretion and outflow during an outburst. Complementary, AU-scale radio observations are needed to probe the disk accretion zone, and 10 AU-scale to probe chemical and kinematic structures of the disk-forming regions, and track changes in the dust, ice, and gas within protostellar envelopes. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: Submitted to the US Astronomy and Astrophysics Decadal Survey (Astro2020)

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

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

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

Submitted 26 February, 2019; originally announced February 2019.

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

ALMA Observations of the massive molecular outflow G331.512-0.103 II: physical properties, kinematics, and geometry modeling

Authors: Carlos Hervías-Caimapo, Manuel Merello, Leonardo Bronfman, Lars Åke-Nyman, Guido Garay, Nadia Lo, Neal J. Evans II, Cristian López-Calderón, Edgar Mendoza

Abstract: We present observations and analysis of the massive molecular outflow G331.512-0.103, obtained with ALMA band 7, continuing the work from Merello et al. (2013). Several lines were identified in the observed bandwidth, consisting of two groups: lines with narrow profiles, tracing the emission from the core ambient medium; and lines with broad velocity wings, tracing the outflow and shocked gas emis… ▽ More We present observations and analysis of the massive molecular outflow G331.512-0.103, obtained with ALMA band 7, continuing the work from Merello et al. (2013). Several lines were identified in the observed bandwidth, consisting of two groups: lines with narrow profiles, tracing the emission from the core ambient medium; and lines with broad velocity wings, tracing the outflow and shocked gas emission. The physical and chemical conditions, such as density, temperature, and fractional abundances are calculated. The ambient medium, or core, has a mean density of $\sim 5\times 10^6$ cm$^{-3}$ and a temperature of $\sim 70$ K. The SiO and SO$_2$ emission trace the very dense and hot part of the shocked outflow, with values of $n_{\rm H_2}\sim10^9$ cm$^{-3}$ and $T \sim 160-200$ K. The interpretation of the molecular emission suggests an expanding cavity geometry powered by stellar winds from a new-born UCHII region, alongside a massive and high-velocity molecular outflow. This scenario, along with the estimated physical conditions, is modeled using the 3D geometry radiative transfer code MOLLIE for the SiO(J$=8-7$) molecular line. The main features of the outflow and the expanding shell are reproduced by the model. △ Less

Submitted 23 December, 2018; originally announced December 2018.

Comments: 19 pages, 17 figures, 2 appendices. Accepted for publication by The Astrophysical Journal

Thermal balance and comparison of gas and dust properties of dense clumps in the Hi-GAL survey

Authors: Manuel Merello, Sergio Molinari, Kazi L. J. Rygl, Neal J. Evans II, Davide Elia, Eugenio Schisano, Alessio Traficante, Yancy Shirley, Brian Svoboda, Paul F. Goldsmith

Abstract: We present a comparative study of physical properties derived from gas and dust emission in a sample of 1068 dense Galactic clumps. The sources are selected from the crossmatch of the Hi-GAL survey with 16 catalogues of NH$_3$ line emission in its lowest inversion (1,1) and (2,2) transitions. The sample covers a large range in masses and bolometric luminosities, with surface densities above… ▽ More We present a comparative study of physical properties derived from gas and dust emission in a sample of 1068 dense Galactic clumps. The sources are selected from the crossmatch of the Hi-GAL survey with 16 catalogues of NH$_3$ line emission in its lowest inversion (1,1) and (2,2) transitions. The sample covers a large range in masses and bolometric luminosities, with surface densities above $Σ=0.1$ g cm$^{-2}$ and with low virial parameters $α<1$. The comparison between dust and gas properties shows an overall agreement between $T_{\textit{kin}}$ and $T_{\textit{dust}}$ at volumetric densities $n\gtrsim1.2\times10^{4}$ cm$^{-3}$, and a median fractional abundance $χ$(NH$_3$)$=1.46\times10^{-8}$. While the protostellar clumps in the sample have small differences between $T_{\textit{kin}}$ and $T_{\textit{dust}}$, prestellar clumps have a median ratio $T_{\textit{kin}}/T_{\textit{dust}}=1.24$, suggesting that these sources are thermally decoupled. A correlation is found between the evolutionary tracer $L/M$ and the parameters $T_{\textit{kin}}/T_{\textit{dust}}$ and $χ$(NH$_3$) in prestellar sources and protostellar clumps with $L/M<1$ L$_\odot$ M$_\odot^{-1}$. In addition, a weak correlation is found between non-thermal velocity dispersion and the $L/M$ parameter, possibly indicating an increase of turbulence with protostellar evolution in the interior of clumps. Finally, different processes are discussed to explain the differences between gas and dust temperatures in prestellar candidates, and the origin of non-thermal motions observed in the clumps. △ Less

Submitted 14 December, 2018; originally announced December 2018.

Comments: 28 pages, 22 figures, accepted for publication in MNRAS