Showing 1–50 of 69 results for author: Busquet, G

Multiple chemical tracers finally unveil the intricate NGC\,1333 IRAS\,4A outflow system. FAUST XVI

Authors: Layal Chahine, Cecilia Ceccarelli, Marta De Simone, Claire J. Chandler, Claudio Codella, Linda Podio, Ana López-Sepulcre, Nami Sakai, Laurent Loinard, Mathilde Bouvier, Paola Caselli, Charlotte Vastel, Eleonora Bianchi, Nicolás Cuello, Francesco Fontani, Doug Johnstone, Giovanni Sabatini, Tomoyuki Hanawa, Ziwei E. Zhang, Yuri Aikawa, Gemma Busquet, Emmanuel Caux, Aurore Durán, Eric Herbst, François Ménard , et al. (32 additional authors not shown)

Abstract: The exploration of outflows in protobinary systems presents a challenging yet crucial endeavour, offering valuable insights into the dynamic interplay between protostars and their evolution. In this study, we examine the morphology and dynamics of jets and outflows within the IRAS\,4A protobinary system. This analysis is based on ALMA observations of SiO(5--4), H$_2$CO(3$_{0,3}$--2$_{0,3}$), and H… ▽ More The exploration of outflows in protobinary systems presents a challenging yet crucial endeavour, offering valuable insights into the dynamic interplay between protostars and their evolution. In this study, we examine the morphology and dynamics of jets and outflows within the IRAS\,4A protobinary system. This analysis is based on ALMA observations of SiO(5--4), H$_2$CO(3$_{0,3}$--2$_{0,3}$), and HDCO(4$_{1,4}$--3$_{1,3}$) with a spatial resolution of $\sim$150\,au. Leveraging an astrochemical approach involving the use of diverse tracers beyond traditional ones has enabled the identification of novel features and a comprehensive understanding of the broader outflow dynamics. Our analysis reveals the presence of two jets in the redshifted emission, emanating from IRAS\,4A1 and IRAS\,4A2, respectively. Furthermore, we identify four distinct outflows in the region for the first time, with each protostar, 4A1 and 4A2, contributing to two of them. We characterise the morphology and orientation of each outflow, challenging previous suggestions of bends in their trajectories. The outflow cavities of IRAS\,4A1 exhibit extensions of 10$''$ and 13$''$ with position angles (PA) of 0$^{\circ}$ and -12$^{\circ}$, respectively, while those of IRAS\,4A2 are more extended, spanning 18$''$ and 25$''$ with PAs of 29$^{\circ}$ and 26$^{\circ}$. We propose that the misalignment of the cavities is due to a jet precession in each protostar, a notion supported by the observation that the more extended cavities of the same source exhibit lower velocities, indicating they may stem from older ejection events. △ Less

Submitted 21 May, 2024; originally announced May 2024.

ALMA-IMF XIII: N$_2$H$^+$ kinematic analysis on the intermediate protocluster G353.41

Authors: R. H. Álvarez-Gutiérrez, A. M. Stutz, N. Sandoval-Garrido, F. Louvet, F. Motte, R. Galván-Madrid, N. Cunningham, P. Sanhueza, M. Bonfand, S. Bontemps, A. Gusdorf, A. Ginsburg, T. Csengeri, S. D. Reyes, J. Salinas, T. Baug, L. Bronfman, G. Busquet, D. J. Díaz-González, M. Fernandez-Lopez, A. Guzmán, A. Koley, H. -L. Liu, F. A. Olguin, M. Valeille-Manet , et al. (1 additional authors not shown)

Abstract: The ALMA-IMF Large Program provides multi-tracer observations of 15 Galactic massive protoclusters at matched sensitivity and spatial resolution. We focus on the dense gas kinematics of the G353.41 protocluster traced by N$_2$H$^+$ (1$-$0), with an spatial resolution $\sim$0.02 pc. G353.41, at a distance of $\sim$2 kpc, has a mass of $\sim$2500 M$_{\odot}$ within $1.3\times1.3$ pc$^2$. We extract… ▽ More The ALMA-IMF Large Program provides multi-tracer observations of 15 Galactic massive protoclusters at matched sensitivity and spatial resolution. We focus on the dense gas kinematics of the G353.41 protocluster traced by N$_2$H$^+$ (1$-$0), with an spatial resolution $\sim$0.02 pc. G353.41, at a distance of $\sim$2 kpc, has a mass of $\sim$2500 M$_{\odot}$ within $1.3\times1.3$ pc$^2$. We extract the N$_2$H$^+$ isolated line component and we decompose it by fitting up to 3 Gaussian velocity components. This allows us to identify velocity structures that are impossible to identify in the traditional position-velocity diagram. We identify multiple velocity gradients on large ($\sim$1 pc) and small scales ($\sim$0.2 pc). We find good agreement between the N$_2$H$^+$ velocities and the previously reported DCN core velocities, suggesting that cores are kinematically coupled to the dense gas in which they form. We measure 9 converging ``V-shaped'' velocity gradients ($\sim20$ km/s/pc), located in filaments, which are sometimes associated with cores near their point of convergence. The average inflow timescale is $\sim67$ kyr, or about twice the free-fall time of cores in the same area ($\sim33$ kyr) but substantially shorter than protostar lifetime estimates ($\sim$0.5 Myr). We derive mass accretion rates in the range of $(0.35-8.77)\,\times\,10^{-4}$ M$_{\odot}$/yr. This feeding might lead to further filament collapse and formation of new cores. We suggest that the protocluster is collapsing on large scales, but the velocity signature of collapse is slow compared to pure free-fall. These data are consistent with a comparatively slow global protocluster contraction under gravity, and faster core formation within, suggesting the formation of multiple generations of stars over the protocluster lifetime. △ Less

Submitted 13 June, 2024; v1 submitted 10 April, 2024; originally announced April 2024.

Comments: Accepted in A&A, 22 pages, 21 figures, 1 interactive figure, 4 tables

Can radial temperature profiles be inferred using NH3 (1, 1) and (2, 2) observations?

Authors: Robert Estalella, Aina Palau, Gemma Busquet

Abstract: A number of works infer radial temperature profiles of envelopes surrounding young stellar objects using several rotational transitions in a pixel-by-pixel or azimuthally-averaged basis. However, in many cases the assumption that the rotational temperature is constant along the line of sight is made, while this is not the case when a partially resolved envelope, assumed to be spherically symmetric… ▽ More A number of works infer radial temperature profiles of envelopes surrounding young stellar objects using several rotational transitions in a pixel-by-pixel or azimuthally-averaged basis. However, in many cases the assumption that the rotational temperature is constant along the line of sight is made, while this is not the case when a partially resolved envelope, assumed to be spherically symmetric, is used to obtain values of temperature for different projected radii. This kind of analysis (homogeneous analysis) is intrinsically inconsistent. By using a spherical envelope model to interpret NH3 (1, 1) and (2, 2) observations, we tested how robust it is to infer radial temperature profiles of an envelope. The temperature and density of the model envelope are power laws of radius, but the density can be flat for an inner central part. The homogeneous analysis was applied to obtain radial temperature profiles, and resulted that for small projected radii, where the optical depth of the lines is high, the homogeneous temperature can be much higher than the actual envelope temperature. In general, for larger projected radii, both the temperature and the temperature power-law index can be underestimated by as much as 40%, and 0.15, respectively. We applied this study to the infrared dark cloud G14.225-0.506 for which the radial temperature profile was previously derived from the dust emission at submillimeter wavelengths and the spectral energy distribution. As expected, the homogeneous analysis underestimated both the temperature and the temperature power-law index. △ Less

Submitted 25 January, 2024; originally announced January 2024.

Comments: Accepted for publication in MNRAS

Radio survey of the stellar population in the infrared dark cloud G14.225-0.506

Authors: Elena Díaz-Márquez, Roger Grau, Gemma Busquet, Josep Miquel Girart, Álvaro Sánchez-Monge, Aina Palau, Matthew S. Povich, Nacho Añez-López, Hauyu Baobab Liu, Qizhou Zhang, Robert Estalella

Abstract: The IRDC G14.225-0.506 is associated with a network of filaments, which result in two different dense hubs, as well as with several signposts of star formation activity. The aim of this work is to study the cm continuum emission to characterize the stellar population in G14.2. We performed deep (~1.5-3 microJy) radio continuum observations at 6 and 3.6 cm using the VLA in the A configuration (~0.3… ▽ More The IRDC G14.225-0.506 is associated with a network of filaments, which result in two different dense hubs, as well as with several signposts of star formation activity. The aim of this work is to study the cm continuum emission to characterize the stellar population in G14.2. We performed deep (~1.5-3 microJy) radio continuum observations at 6 and 3.6 cm using the VLA in the A configuration (~0.3''). We have also made use of observations taken during different days to study the presence of variability at short timescales. We detected a total of 66 sources, 32 in the northern region G14.2-N and 34 in the southern region G14.2-S. Ten of the sources are found to be variable. Based on their spectral index, the emission in G14.2-N is mainly dominated by non-thermal sources while G14.2-S contains more thermal emitters. Approximately 75% of the sources present a counterpart at other wavelengths. In the inner 0.4~pc region around the center of each hub, the number of IR sources in G14.2-N is larger than in G14.2-S by a factor of 4. We also studied the relation between the radio luminosity and the bolometric luminosity, finding that the thermal emission of the studied sources is compatible with thermal radio jets. For our sources with X-ray counterparts, the non-thermal emitters follow a Güdel-Benz relation with k = 0.03. We found similar levels of fragmentation between G14.2-N and G14.2-S, suggesting that both regions are most likely twin hubs. The non-thermal emission found in the less evolved objects suggests that G14.2-N may be composed of more massive YSOs as well as being in a more advanced evolutionary stage, consistent with the filament-halo gradient in age and mass from previous works. Our results confirm a wider evolutionary sequence starting in G14.2-S as the youngest part, followed by G14.2-N, and ending with the most evolved region M17. △ Less

Submitted 21 November, 2023; originally announced November 2023.

Comments: 28 pages, 22 figures. Accepted for publication in A&A

The role of turbulence in high-mass star formation: Subsonic and transonic turbulence are ubiquitously found at early stages

Authors: Chao Wang, Ke Wang, Feng-Wei Xu, Patricio Sanhueza, Hauyu Baobab Liu, Qizhou Zhang, Xing Lu, F. Fontani, Paola Caselli, Gemma Busquet, Jonathan C. Tan, Di Li, J. M. Jackson, Thushara Pillai, Paul T. P. Ho, Andrés E. Guzmán, Nannan Yue

Abstract: Context. Traditionally, supersonic turbulence is considered to be one of the most likely mechanisms to slow down the gravitational collapse in dense clumps, thereby enabling the formation of massive stars. However, several recent studies have raised differing points of view based on observations carried out with sufficiently high spatial and spectral resolution. These studies call for a re-evaluat… ▽ More Context. Traditionally, supersonic turbulence is considered to be one of the most likely mechanisms to slow down the gravitational collapse in dense clumps, thereby enabling the formation of massive stars. However, several recent studies have raised differing points of view based on observations carried out with sufficiently high spatial and spectral resolution. These studies call for a re-evaluation of the role turbulence plays in massive star-forming regions. Aims. Our aim is to study the gas properties, especially the turbulence, in a sample of massive star-forming regions with sufficient spatial and spectral resolution, which can both resolve the core fragmentation and the thermal line width. Methods. We observed NH3 metastable lines with the Very Large Array (VLA) to assess the intrinsic turbulence. Results. Analysis of the turbulence distribution histogram for 32 identified NH3 cores reveals the presence of three distinct components. Furthermore, our results suggest that (1) sub- and transonic turbulence is a prevalent (21 of 32) feature of massive star-forming regions and those cold regions are at early evolutionary stage. This investigation indicates that turbulence alone is insufficient to provide the necessary internal pressure required for massive star formation, necessitating further exploration of alternative candidates; and (2) studies of seven multi-core systems indicate that the cores within each system mainly share similar gas properties and masses. However, two of the systems are characterized by the presence of exceptionally cold and dense cores that are situated at the spatial center of each system. Our findings support the hub-filament model as an explanation for this observed distribution △ Less

Submitted 7 February, 2024; v1 submitted 27 October, 2023; originally announced October 2023.

Comments: 34 pages, 15 figures, 4 tables. Accepted for publication on A&A

Disk and Envelope Streamers of the GGD27-MM1 Massive Protostar

Authors: M. Fernández-López, J. M. Girart, J. A. López-Vázquez, R. Estalella, G. Busquet, S. Curiel, N. Añez-López

Abstract: We present new Atacama Large (sub)Millimeter Array 0.98 mm observations of the continuum emission and several molecular lines toward the high-mass protostellar system GGD27-MM1, driving the HH 80-81 radio-jet. The detailed analysis of the continuum and the CH$_3$CN molecular emission allows us to separate the contributions from the dust content of the disk (extending up to 190 au), the molecular c… ▽ More We present new Atacama Large (sub)Millimeter Array 0.98 mm observations of the continuum emission and several molecular lines toward the high-mass protostellar system GGD27-MM1, driving the HH 80-81 radio-jet. The detailed analysis of the continuum and the CH$_3$CN molecular emission allows us to separate the contributions from the dust content of the disk (extending up to 190 au), the molecular content of the disk (extending from 140 to 360 au), and the content of the envelope, revealing the presence of several possible accretion streamers (also seen in other molecular tracers, such as CH$_3$OH). We analyze the physical properties of the system, producing temperature and column density maps, and radial profiles for the disk and the envelope. We qualitatively reproduce the trajectories and line-of-sight velocities of the possible streamers using a theoretical model approach. An ad-hoc model of a flared disk comprising a hot dust disk embedded in cold gas fits the H$_2$S emission, which revealed the molecular disk as crescent-shape with a prominent central absorption. Another fit to the central absorption spectrum suggests that the absorption is probably caused by different external cold layers from the envelope or the accretion streamers. Finally, the analysis of the rotation pattern of the different molecular transitions in the molecular disk, suggests that there is an inner zone devoid of molecular content. △ Less

Submitted 12 July, 2023; originally announced July 2023.

Comments: 36 pages, 28 figures, 6 tables. Accepted by the ApJ

ALMA-IMF. VII. First release of the full spectral line cubes: Core kinematics traced by DCN J=(3-2)

Authors: N. Cunningham, A. Ginsburg, R. Galván-Madrid, F. Motte, T. Csengeri, A. M. Stutz, M. Fernández-López, R. H. Álvarez-Gutiérrez, M. Armante, T. Baug, M. Bonfand, S. Bontemps, J. Braine, N. Brouillet, G. Busquet, D. J. Díaz-González, J. Di Francesco, A. Gusdorf, F. Herpin, H. Liu, A. López-Sepulcre, F. Louvet, X. Lu, L. Maud, T. Nony , et al. (8 additional authors not shown)

Abstract: ALMA-IMF is an Atacama Large Millimeter/submillimeter Array (ALMA) Large Program designed to measure the core mass function (CMF) of 15 protoclusters chosen to span their early evolutionary stages. It further aims to understand their kinematics, chemistry, and the impact of gas inflow, accretion, and dynamics on the CMF. We present here the first release of the ALMA-IMF line data cubes (DR1), prod… ▽ More ALMA-IMF is an Atacama Large Millimeter/submillimeter Array (ALMA) Large Program designed to measure the core mass function (CMF) of 15 protoclusters chosen to span their early evolutionary stages. It further aims to understand their kinematics, chemistry, and the impact of gas inflow, accretion, and dynamics on the CMF. We present here the first release of the ALMA-IMF line data cubes (DR1), produced from the combination of two ALMA 12m-array configurations. The data include 12 spectral windows, with eight at 1.3mm and four at 3mm. The broad spectral coverage of ALMA-IMF (~6.7 GHz bandwidth coverage per field) hosts a wealth of simple atomic, molecular, ionised, and complex organic molecular lines. We describe the line cube calibration done by ALMA and the subsequent calibration and imaging we performed. We discuss our choice of calibration parameters and optimisation of the cleaning parameters, and we demonstrate the utility and necessity of additional processing compared to the ALMA archive pipeline. As a demonstration of the scientific potential of these data, we present a first analysis of the DCN (3-2) line. We find that DCN traces a diversity of morphologies and complex velocity structures, which tend to be more filamentary and widespread in evolved regions and are more compact in the young and intermediate-stage protoclusters. Furthermore, we used the DCN (3-2) emission as a tracer of the gas associated with 595 continuum cores across the 15 protoclusters, providing the first estimates of the core systemic velocities and linewidths within the sample. We find that DCN (3-2) is detected towards a higher percentage of cores in evolved regions than the young and intermediate-stage protoclusters and is likely a more complete tracer of the core population in more evolved protoclusters. The full ALMA 12m-array cubes for the ALMA-IMF Large Program are provided with this DR1 release. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: 75 pages (21 main body; 54 appendix), 37 figures. The ALMA-IMF DR1 line release is hosted at https://dataverse.harvard.edu/dataverse/alma-imf

Journal ref: A&A 678, A194 (2023)

Matter ejections behind the highs and lows of the transitional millisecond pulsar PSR J1023+0038

Authors: M. C. Baglio, F. Coti Zelati, S. Campana, G. Busquet, P. D'Avanzo, S. Giarratana, M. Giroletti, F. Ambrosino, S. Crespi, A. Miraval Zanon, X. Hou, D. Li, J. Li, P. Wang, D. M. Russell, D. F. Torres, K. Alabarta, P. Casella, S. Covino, D. M. Bramich, D. de Martino, M. Méndez, S. E. Motta, A. Papitto, P. Saikia , et al. (1 additional authors not shown)

Abstract: Transitional millisecond pulsars are an emerging class of sources that link low-mass X-ray binaries to millisecond radio pulsars in binary systems. These pulsars alternate between a radio pulsar state and an active low-luminosity X-ray disc state. During the active state, these sources exhibit two distinct emission modes (high and low) that alternate unpredictably, abruptly, and incessantly. X-ray… ▽ More Transitional millisecond pulsars are an emerging class of sources that link low-mass X-ray binaries to millisecond radio pulsars in binary systems. These pulsars alternate between a radio pulsar state and an active low-luminosity X-ray disc state. During the active state, these sources exhibit two distinct emission modes (high and low) that alternate unpredictably, abruptly, and incessantly. X-ray to optical pulsations are observed only during the high mode. The root cause of this puzzling behaviour remains elusive. This paper presents the results of the most extensive multi-wavelength campaign ever conducted on the transitional pulsar prototype, PSR J1023+0038, covering from the radio to X-rays. The campaign was carried out over two nights in June 2021 and involved 12 different telescopes and instruments, including XMM-Newton, HST, VLT/FORS2 (in polarimetric mode), ALMA, VLA, and FAST. By modelling the broadband spectral energy distributions in both emission modes, we show that the mode switches are caused by changes in the innermost region of the accretion disc. These changes trigger the emission of discrete mass ejections, which occur on top of a compact jet, as testified by the detection of at least one short-duration millimetre flare with ALMA at the high-to-low mode switch. The pulsar is subsequently re-enshrouded, completing our picture of the mode switches. △ Less

Submitted 28 August, 2023; v1 submitted 23 May, 2023; originally announced May 2023.

Comments: 25 pages, 12 figures, 9 tables. Accepted for publication on Astronomy and Astrophysics

Journal ref: A&A 677, A30 (2023)

FAUST VIII. The protostellar disk of VLA 1623-2417 W and its streamers imaged by ALMA

Authors: S. Mercimek, L. Podio, C. Codella, L. Chahine, A. López-Sepulcre, S. Ohashi, L. Loinard, D. Johnstone, F. Menard, N. Cuello, P. Caselli, J. Zamponi, Y. Aikawa, E. Bianchi, G. Busquet, J. E. Pineda, M. Bouvier, M. De Simone, Y. Zhang, N. Sakai, C. J. Chandler, C. Ceccarelli, F. Alves, A. Durán, D. Fedele , et al. (3 additional authors not shown)

Abstract: More than 50% of solar-mass stars form in multiple systems. It is therefore crucial to investigate how multiplicity affects the star and planet formation processes at the protostellar stage. We report continuum and C$^{18}$O (2-1) observations of the VLA 1623-2417 protostellar system at 50 au angular resolution as part of the ALMA Large Program FAUST. The 1.3 mm continuum probes the disks of VLA 1… ▽ More More than 50% of solar-mass stars form in multiple systems. It is therefore crucial to investigate how multiplicity affects the star and planet formation processes at the protostellar stage. We report continuum and C$^{18}$O (2-1) observations of the VLA 1623-2417 protostellar system at 50 au angular resolution as part of the ALMA Large Program FAUST. The 1.3 mm continuum probes the disks of VLA 1623A, B, and W, and the circumbinary disk of the A1+A2 binary. The C$^{18}$O emission reveals, for the first time, the gas in the disk-envelope of VLA 1623W. We estimate the dynamical mass of VLA 1623W, $M_{\rm dyn}=0.45\pm0.08$ M$_{\odot}$, and the mass of its disk, $M_{\rm disk}\sim6\times10^{-3}$ M$_{\odot}$. C$^{18}$O also reveals streamers that extend up to 1000 au, spatially and kinematically connecting the envelope and outflow cavities of the A1+A2+B system with the disk of VLA 1623W. The presence of the streamers, as well as the spatial ($\sim$1300 au) and velocity ($\sim$2.2 km/s) offset of VLA 1623W suggest that either sources W and A+B formed in different cores, interacting between them, or that source W has been ejected from the VLA 1623 multiple system during its formation. In the latter case, the streamers may funnel material from the envelope and cavities of VLA 1623AB onto VLA 1623W, thus concurring to set its final mass and chemical content. △ Less

Submitted 28 March, 2023; originally announced March 2023.

Comments: 10 pages, 5 figures, accepted by MNRAS

ALMA-IMF. V. Prestellar and protostellar core populations in the W43 cloud complex

Authors: T. Nony, R. Galvan-Madrid, F. Motte, Y. Pouteau, N. Cunningham, F. Louvet, A. M. Stutz, B. Lefloch, S. Bontemps, N. Brouillet, A. Ginsburg, I. Joncour, F. Herpin, P. Sanhueza, T. Csengeri, A. P. M. Towner, M. Bonfand, M. Fernández-López, T. Baug, L. Bronfman, G. Busquet, J. Di Francesco, A. Gusdorf, X. Lu, F. Olguin , et al. (2 additional authors not shown)

Abstract: The origin of the stellar initial mass function (IMF) and its relation with the core mass function (CMF) are actively debated issues with important implications in astrophysics. Recent observations in the W43 molecular complex of top-heavy CMFs, with an excess of high-mass cores compared to the canonical mass distribution, raise questions about our understanding of the star formation processes and… ▽ More The origin of the stellar initial mass function (IMF) and its relation with the core mass function (CMF) are actively debated issues with important implications in astrophysics. Recent observations in the W43 molecular complex of top-heavy CMFs, with an excess of high-mass cores compared to the canonical mass distribution, raise questions about our understanding of the star formation processes and their evolution in space and time. We aim to compare populations of protostellar and prestellar cores in three regions imaged in the ALMA-IMF Large Program. We created an homogeneous core catalogue in W43, combining a new core extraction in W43-MM1 with the catalogue of W43-MM2&MM3 presented in a previous work. Our detailed search for protostellar outflows enabled us to identify between 23 and 30 protostellar cores out of 127 cores in W43-MM1 and between 42 and 51 protostellar cores out of 205 cores in W43-MM2&MM3. Cores with neither outflows nor hot core emission are classified as prestellar candidates. We found a similar fraction of cores which are protostellar in the two regions, about 35%. This fraction strongly varies in mass, from 15-20% at low mass, between 0.8 and 3$M_{\odot} $ up to about 80% above 16$M_{\odot}$. Protostellar cores are found to be, on average, more massive and smaller in size than prestellar cores. Our analysis also revealed that the high-mass slope of the prestellar CMF in W43, $α=-1.46_{-0.19}^{+0.12}$, is consistent with the Salpeter slope, and thus the top-heavy form measured for the global CMF, $α=-0.96$, is due to the protostellar core population. Our results could be explained by clump-fed models in which cores grow in mass, especially during the protostellar phase, through inflow from their environment. The difference between the slopes of the prestellar and protostellar CMFs moreover implies that high-mass cores grow more in mass than low-mass cores. △ Less

Submitted 15 March, 2023; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: 28 pages, 18 figures. Accepted for publication in A&A. Title and appendices updated

Journal ref: A&A 674, A75 (2023)

ALMA-IMF VI -- Investigating the origin of stellar masses: Core mass function evolution in the W43-MM2&MM3 mini-starburst

Authors: Y. Pouteau, F. Motte, T. Nony, M. Gonzalez, I. Joncour, J. -F. Robitaille, G. Busquet, R. Galvan-Madrid, A. Gusdorf, P. Hennebelle, A. Ginsburg, T. Csengeri, P. Sanhueza, P. Dell'Ova, A. M. Stutz, A. P. M. Towner, N. Cunningham, F. Louvet, A. Men'shchikov, M. Fernandez-Lopez, N. Schneider, M. Armante, J. Bally, T. Baug, M. Bonfand , et al. (13 additional authors not shown)

Abstract: Among the most central open questions regarding the initial mass function (IMF) of stars is the impact of environment on the shape of the core mass function (CMF) and thus potentially on the IMF. The ALMA-IMF Large Program aims to investigate the variations in the core distributions with cloud characteristics, as diagnostic observables of the formation process and evolution of clouds. The present… ▽ More Among the most central open questions regarding the initial mass function (IMF) of stars is the impact of environment on the shape of the core mass function (CMF) and thus potentially on the IMF. The ALMA-IMF Large Program aims to investigate the variations in the core distributions with cloud characteristics, as diagnostic observables of the formation process and evolution of clouds. The present study focuses on the W43-MM2&MM3 mini-starburst, whose CMF has recently been found to be top-heavy with respect to the Salpeter slope. W43-MM2&MM3 harbors a rich cluster that contains a statistically significant number of cores, which was previously characterized in Paper III. We applied a multi-scale decomposition technique to the ALMA 1.3 mm and 3 mm continuum images to define six subregions. For each subregion we characterized the high column density probability distribution function, n-PDF, and the shape of the cloud gas using the 1.3 mm image. Using the core catalog, we investigate correlations between the CMF and cloud and core properties. We classify the subregions into different stages of evolution, from quiescent to burst to post-burst, based on the surface number density of cores, number of outflows, and UCHii presence. The high-mass end of the subregion CMFs varies from being close to the Salpeter slope (quiescent) to top-heavy (burst and post-burst). Moreover, the second tail of the n-PDF varies from steep, to flat like observed for the high mass star-forming clouds. We found that subregions with flat second n-PDF tails display top-heavy CMFs. The CMF may evolve from Salpeter to top-heavy throughout the star formation process from the quiescent to the burst phase. This scenario raises the question of if the CMF might revert again to Salpeter as the cloud approaches the end of its star formation stage, a hypothesis that remains to be tested. △ Less

Submitted 20 February, 2023; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: Accepted for publication in A&A (december, 5th 2022) after language editing; 25 pages, 14 figures

Journal ref: A&A 674, A76 (2023)

ALMA-IMF IV -- A comparative study of the main hot cores in W43-MM1: detection, temperature and molecular composition

Authors: N. Brouillet, D. Despois, J. Molet, T. Nony, F. Motte, A. Gusdorf, F. Louvet, S. Bontemps, F. Herpin, M. Bonfand, T. Csengeri, A. Ginsburg, N. Cunningham, R. Galvan-Madrid, L. Maud, G. Busquet, L. Bronfman, M. Fernandez-Lopez, D. L. Jeff, B. Lefloch, Y. Pouteau, P. Sanhueza, A. M. Stutz, M. Valeille-Manet

Abstract: W43-MM1 is a young region, very rich in terms of high-mass star formation. We aim to systematically identify the massive cores which contain a hot core and compare their molecular composition. We used ALMA high-spatial resolution (2500 au) data of W43-MM1 to identify line-rich protostellar cores and make a comparative study of their temperature and molecular composition. The identification of hot… ▽ More W43-MM1 is a young region, very rich in terms of high-mass star formation. We aim to systematically identify the massive cores which contain a hot core and compare their molecular composition. We used ALMA high-spatial resolution (2500 au) data of W43-MM1 to identify line-rich protostellar cores and make a comparative study of their temperature and molecular composition. The identification of hot cores is based on both the spatial distribution of the complex organic molecules and the contribution of molecular lines relative to the continuum intensity. We rely on the analysis of CH3CN and CH3CCH to estimate the temperatures of the selected cores. Finally, we rescale the spectra of the different hot cores based on their CH3OCHO line intensities to directly compare the detections and line intensities of the other species. W43-MM1 turns out to be a region rich in massive hot cores with at least 1 less massive and 7 massive hot cores. The excitation temperature of CH3CN is of the same order for all of them (120-160 K). There is a factor of up to 30 difference in the intensity of the complex organic molecules (COMs) lines. However the molecular emission of the hot cores appears to be the same within a factor 2-3. This points towards both a similar chemical composition and excitation of most of the COMs over these massive cores, which span about an order of magnitude in core mass. In contrast, CH3CCH emission is found to preferentially trace more the envelope, with a temperature ranging from 50 K to 90 K. Core 1, the most massive hot core of W43-MM1, shows a richer line spectrum than the other cores. In core 2, the emission of O-bearing molecules does not peak at the dust continuum core center; the blue and red shifted emission correspond to the outflow lobes, suggesting a formation via the sublimation of the ice mantles through shocks or UV irradiation on the walls of the cavity. △ Less

Submitted 7 July, 2022; originally announced July 2022.

Comments: 19 pages, 17 figures. Accepted for publication in A&A

Journal ref: A&A 665, A140 (2022)

FAUST VI. VLA 1623--2417 B: a new laboratory for astrochemistry around protostars on 50 au scale

Authors: C. Codella, A. López-Sepulcre, S. Ohashi, C. J. Chandler, M. De Simone, L. Podio, C. Ceccarelli, N. Sakai, F. Alves, A. Durán, D. Fedele, L. Loinard, S. Mercimek, N. Murillo, E. Bianchi, M. Bouvier, G. Busquet, P. Caselli, F. Dulieu, S. Feng, T. Hanawa, D. Johnstone, B. Lefloch, L. T. Maud, G. Moellenbrock , et al. (3 additional authors not shown)

Abstract: The ALMA interferometer, with its unprecedented combination of high-sensitivity and high-angular resolution, allows for (sub-)mm wavelength mapping of protostellar systems at Solar System scales. Astrochemistry has benefited from imaging interstellar complex organic molecules in these jet-disk systems. Here we report the first detection of methanol (CH3OH) and methyl formate (HCOOCH3) emission tow… ▽ More The ALMA interferometer, with its unprecedented combination of high-sensitivity and high-angular resolution, allows for (sub-)mm wavelength mapping of protostellar systems at Solar System scales. Astrochemistry has benefited from imaging interstellar complex organic molecules in these jet-disk systems. Here we report the first detection of methanol (CH3OH) and methyl formate (HCOOCH3) emission towards the triple protostellar system VLA1623-2417 A1+A2+B, obtained in the context of the ALMA Large Program FAUST. Compact methanol emission is detected in lines from Eu = 45 K up to 61 K and 537 K towards components A1 and B, respectively. LVG analysis of the CH3OH lines towards VLA1623-2417 B indicates a size of 0.11-0.34 arcsec (14-45 au), a column density N(CH3OH) = 10^16-10^17 cm-2, kinetic temperature > 170 K, and volume density > 10^8 cm-3. An LTE approach is used for VLA1623-2417 A1, given the limited Eu range, and yields Trot < 135 K. The methanol emission around both VLA1623-2417 A1 and B shows velocity gradients along the main axis of each disk. Although the axial geometry of the two disks is similar, the observed velocity gradients are reversed. The CH3OH spectra from B shows two broad (4-5 km s-1) peaks, which are red- and blue-shifted by about 6-7 km s-1 from the systemic velocity. Assuming a chemically enriched ring within the accretion disk, close to the centrifugal barrier, its radius is calculated to be 33 au. The methanol spectra towards A1 are somewhat narrower (about 4 km s-1), implying a radius of 12-24 au. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: MNRAS, in press

Hot methanol in the [BHB2007] 11 protobinary system: hot corino versus shock origin? : FAUST V

Authors: C. Vastel, F. Alves, C. Ceccarelli, M. Bouvier, I. Jimenez-Serra, T. Sakai, P. Caselli, L. Evans, F. Fontani, R. Le Gal, C. J. Chandler, B. Svoboda, L. Maud, C. Codella, N. Sakai, A. Lopez-Sepulcre, G. Moellenbrock, Y. Aikawa, N. Balucani, E. Bianchi, G. Busquet, E. Caux, S. Charnley, N. Cuello, M. De Simone , et al. (41 additional authors not shown)

Abstract: Methanol is a ubiquitous species commonly found in the molecular interstellar medium. It is also a crucial seed species for the building-up of the chemical complexity in star forming regions. Thus, understanding how its abundance evolves during the star formation process and whether it enriches the emerging planetary system is of paramount importance. We used new data from the ALMA Large Program F… ▽ More Methanol is a ubiquitous species commonly found in the molecular interstellar medium. It is also a crucial seed species for the building-up of the chemical complexity in star forming regions. Thus, understanding how its abundance evolves during the star formation process and whether it enriches the emerging planetary system is of paramount importance. We used new data from the ALMA Large Program FAUST (Fifty AU STudy of the chemistry in the disk/envelope system of Solar-like protostars) to study the methanol line emission towards the [BHB2007] 11 protobinary system (sources A and B), where a complex structure of filaments connecting the two sources with a larger circumbinary disk has been previously detected. Twelve methanol lines have been detected with upper energies in the range [45-537] K along with one 13CH3OH transition. The methanol emission is compact and encompasses both protostars, separated by only 28 au and presents three velocity components, not spatially resolved by our observations, associated with three different spatial regions, with two of them close to 11B and the third one associated with 11A. A non-LTE radiative transfer analysis of the methanol lines concludes that the gas is hot and dense and highly enriched in methanol with an abundance as high as 1e-5. Using previous continuum data, we show that dust opacity can potentially completely absorb the methanol line emission from the two binary objects. Although we cannot firmly exclude other possibilities, we suggest that the detected hot methanol is resulting from the shocked gas from the incoming filaments streaming towards [BHB2007] 11 A and B, respectively. Higher spatial resolution observations are necessary to confirm this hypothesis. △ Less

Submitted 21 June, 2022; originally announced June 2022.

Comments: 17 pages, accepted in A&A

Journal ref: A&A 664, A171 (2022)

ALMA-IMF III -- Investigating the origin of stellar masses: Top-heavy core mass function in the W43-MM2&MM3 mini-starburst

Authors: Y. Pouteau, F. Motte, T. Nony, R. Galván-Madrid, A. Men'shchikov, S. Bontemps, J. -F. Robitaille, F. Louvet, A. Ginsburg, F. Herpin, A. López-Sepulcre, P. Dell'Ova, A. Gusdorf, P. Sanhueza, A. M. Stutz, N. Brouillet, B. Thomasson, M. Armante, T. Baug, G. Busquet, T. Csengeri, N. Cunningham, M. Fernández-López, H. -L. Liu, F. Olguin , et al. (13 additional authors not shown)

Abstract: The ALMA-IMF Large Program observed the W43-MM2-MM3 ridge, whose 1.3mm and 3mm ALMA 12m array continuum images reach a 2500au spatial resolution. We used both the best-sensitivity and the line-free ALMA-IMF images, reduced the noise with the multi-resolution segmentation technique MnGSeg, and derived the most complete and most robust core catalog possible. Using two different extraction software p… ▽ More The ALMA-IMF Large Program observed the W43-MM2-MM3 ridge, whose 1.3mm and 3mm ALMA 12m array continuum images reach a 2500au spatial resolution. We used both the best-sensitivity and the line-free ALMA-IMF images, reduced the noise with the multi-resolution segmentation technique MnGSeg, and derived the most complete and most robust core catalog possible. Using two different extraction software packages, getsf and GExt2D, we identified 200 compact sources, whose 100 common sources have on average fluxes consistent to within 30%. We filtered sources with non-negligible free-free contamination and corrected fluxes from line contamination, resulting in a W43-MM2-MM3 catalog of 205 getsf cores. With a median deconvolved FWHM size of 3400au, core masses range from 0.1Msun to 70Msun and the getsf catalog is 90% complete down to 0.8Msun. The high-mass end of the core mass function (CMF) of W43-MM2-MM3 is top-heavy compared to the canonical IMF. Fitting the cumulative CMF with a single power law of the form N(>logM)\propto M^a, we measured a=-0.95\pm0.04, compared to the canonical a=-1.35 Salpeter IMF slope. The slope of the CMF is robust with respect to map processing, extraction software package, and reasonable variations in the assumptions taken to estimate core masses. We explore several assumptions on how cores transfer their mass to stars and sub-fragment to predict the IMF resulting from the W43-MM2-MM3 CMF. In stark contrast to the commonly accepted paradigm, our result argues against the universality of the CMF shape. More robust functions of the star-formation efficiency and core sub-fragmentation are required to better predict the resulting IMF, here suggested to remain top-heavy at the end of the star-formation phase. If confirmed, the IMFs emerging from starburst events could inherit their top-heavy shape from their parental CMFs, challenging the IMF universality. △ Less

Submitted 28 April, 2022; v1 submitted 7 March, 2022; originally announced March 2022.

Comments: 28 pages, 16 figures, 6 tables, A&A accepted on 02/08/22, revised following the language referee report

Journal ref: A&A 664, A26 (2022)

FAUST III. Misaligned rotations of the envelope, outflow, and disks in the multiple protostellar system of VLA 1623$-$2417

Authors: Satoshi Ohashi, Claudio Codella, Nami Sakai, Claire J. Chandler, Cecilia Ceccarelli, Felipe Alves, Davide Fedele, Tomoyuki Hanawa, Aurora Durán, Cécile Favre, Ana López-Sepulcre, Laurent Loinard, Seyma Mercimek, Nadia M. Murillo, Linda Podio, Yichen Zhang, Yuri Aikawa, Nadia Balucani, Eleonora Bianchi, Mathilde Bouvier, Gemma Busquet, Paola Caselli, Emmanuel Caux, Steven Charnley, Spandan Choudhury , et al. (47 additional authors not shown)

Abstract: We report a study of the low-mass Class-0 multiple system VLA 1623AB in the Ophiuchus star-forming region, using H$^{13}$CO$^+$ ($J=3-2$), CS ($J=5-4$), and CCH ($N=3-2$) lines as part of the ALMA Large Program FAUST. The analysis of the velocity fields revealed the rotation motion in the envelope and the velocity gradients in the outflows (about 2000 au down to 50 au). We further investigated the… ▽ More We report a study of the low-mass Class-0 multiple system VLA 1623AB in the Ophiuchus star-forming region, using H$^{13}$CO$^+$ ($J=3-2$), CS ($J=5-4$), and CCH ($N=3-2$) lines as part of the ALMA Large Program FAUST. The analysis of the velocity fields revealed the rotation motion in the envelope and the velocity gradients in the outflows (about 2000 au down to 50 au). We further investigated the rotation of the circum-binary VLA 1623A disk as well as the VLA 1623B disk. We found that the minor axis of the circum-binary disk of VLA 1623A is misaligned by about 12 degrees with respect to the large-scale outflow and the rotation axis of the envelope. In contrast, the minor axis of the circum-binary disk is parallel to the large-scale magnetic field according to previous dust polarization observations, suggesting that the misalignment may be caused by the different directions of the envelope rotation and the magnetic field. If the velocity gradient of the outflow is caused by rotation, the outflow has a constant angular momentum and the launching radius is estimated to be $5-16$ au, although it cannot be ruled out that the velocity gradient is driven by entrainments of the two high-velocity outflows. Furthermore, we detected for the first time a velocity gradient associated with rotation toward the VLA 16293B disk. The velocity gradient is opposite to the one from the large-scale envelope, outflow, and circum-binary disk. The origin of its opposite gradient is also discussed. △ Less

Submitted 18 January, 2022; originally announced January 2022.

Comments: 27 pages, 21 figures, 2 Tables, Accepted for publication in ApJ

A train of shocks at 3000 au scale? Exploring the clash of an expanding bubble into the NGC 1333 IRAS 4 region. SOLIS XIV

Authors: Marta De Simone, Claudio Codella, Cecilia Ceccarelli, Ana López-Sepulcre, Roberto Neri, Pedro Ruben Rivera-Ortiz, Gemma Busquet, Paola Caselli, Eleonora Bianchi, Francesco Fontani, Bertrand Lefloch, Yoko Oya, Jaime E. Pineda

Abstract: There is evidence that the star formation process is linked to the intricate net of filaments in molecular clouds, which may be also due to gas compression from external triggers. We studied the southern region of the Perseus NGC 1333 molecular cloud, known to be heavily shaped by similar external triggers, to shed light on the process that perturbed the filament where the Class 0 IRAS4 protostars… ▽ More There is evidence that the star formation process is linked to the intricate net of filaments in molecular clouds, which may be also due to gas compression from external triggers. We studied the southern region of the Perseus NGC 1333 molecular cloud, known to be heavily shaped by similar external triggers, to shed light on the process that perturbed the filament where the Class 0 IRAS4 protostars lie. We use new IRAM-NOEMA observations of SiO and CH3OH, both known to trace violent events as shocks, toward IRAS 4A as part of the Large Program Seeds Of Life in Space (SOLIS). We detected three parallel elongated ($>$6000 au) structures, called fingers, with narrow line profiles (~1.5 $km s^{-1}$) peaked at the cloud systemic velocity, tracing gas with high density (5-20 $10^5 cm^{-3}$) and high temperature (80-160 K). They are chemically different, with the northern finger traced by both SiO and CH3OH ([CH3OH]/[SiO]~160-300), while the other two only by SiO ([CH3OH]/[SiO]$<$ 40). Among various possibilities, a train of three shocks, distanced by $>$5000 yr, would be consistent with the observations if a substantial fraction of silicon, frozen onto the grain mantles, is released by the shocks.We suggest that the shock train is due to an expanding gas bubble, coming behind NGC 1333 from the southwest and clashing against the filament, where IRAS 4A lies. Finally, we propose a solution to the two-decades long debate on the nature and origin of the widespread narrow SiO emission observed in the south part of NGC 1333, namely that it is due to unresolved trains of shocks. △ Less

Submitted 18 February, 2022; v1 submitted 10 January, 2022; originally announced January 2022.

Comments: accepted to MNRAS, stac083

ALMA-IMF II -- investigating the origin of stellar masses: Continuum Images and Data Processing

Authors: A. Ginsburg, T. Csengeri, R. Galván-Madrid, N. Cunningham, R. H. Álvarez-Gutiérrez, T. Baug, M. Bonfand, S. Bontemps, G. Busquet, D. J. Díaz-González, M. Fernández-López, A. Guzmán, F. Herpin, H. Liu, A. López-Sepulcre, F. Louvet, L. Maud, F. Motte, F. Nakamura, T. Nony, F. A. Olguin, Y. Pouteau, P. Sanhueza, A. M. Stutz, A. P. M. Towner , et al. (27 additional authors not shown)

Abstract: We present the first data release of the ALMA-IMF Large Program, which covers the 12m-array continuum calibration and imaging. The ALMA-IMF Large Program is a survey of fifteen dense molecular cloud regions spanning a range of evolutionary stages that aims to measure the core mass function (CMF). We describe the data acquisition and calibration done by the Atacama Large Millimeter/submillimeter Ar… ▽ More We present the first data release of the ALMA-IMF Large Program, which covers the 12m-array continuum calibration and imaging. The ALMA-IMF Large Program is a survey of fifteen dense molecular cloud regions spanning a range of evolutionary stages that aims to measure the core mass function (CMF). We describe the data acquisition and calibration done by the Atacama Large Millimeter/submillimeter Array (ALMA) observatory and the subsequent calibration and imaging we performed. The image products are combinations of multiple 12m array configurations created from a selection of the observed bandwidth using multi-term, multi-frequency synthesis imaging and deconvolution. The data products are self-calibrated and exhibit substantial noise improvements over the images produced from the delivered data. We compare different choices of continuum selection, calibration parameters, and image weighting parameters, demonstrating the utility and necessity of our additional processing work. Two variants of continuum selection are used and will be distributed: the "best-sensitivity" data, which include the full bandwidth, including bright emission lines that contaminate the continuum, and "cleanest", which select portions of the spectrum that are unaffected by line emission. We present a preliminary analysis of the spectral indices of the continuum data, showing that the ALMA products are able to clearly distinguish free-free emission from dust emission, and that in some cases we are able to identify optically thick emission sources. The data products are made public with this release. △ Less

Submitted 13 May, 2023; v1 submitted 15 December, 2021; originally announced December 2021.

Comments: Data released on Zenodo at https://zenodo.org/record/5702966. Accepted for publication in A&A

Journal ref: A&A 662, A9 (2022)

ALMA-IMF I -- Investigating the origin of stellar masses: Introduction to the Large Program and first results

Authors: F. Motte, S. Bontemps, T. Csengeri, Y. Pouteau, F. Louvet, A. M. Stutz, N. Cunningham, A. López-Sepulcre, N. Brouillet, R. Galván-Madrid, A. Ginsburg, L. Maud, A. Men'shchikov, F. Nakamura, T. Nony, P. Sanhueza, R. H. Álvarez-Gutiérrez, M. Armante, T. Baug, M. Bonfand, G. Busquet, E. Chapillon, D. Díaz-González, M. Fernández-López, A. E. Guzmán , et al. (39 additional authors not shown)

Abstract: The ALMA-IMF Large Program imaged a total noncontiguous area of 53pc2, covering 15 extreme, nearby protoclusters of the Milky Way. They were selected to span relevant early protocluster evolutionary stages. Our 1.3mm and 3mm observations provide continuum images that are homogeneously sensitive to point-like cores with masses of 0.2 and 0.6Msun, respectively, with a matched spatial resolution of 2… ▽ More The ALMA-IMF Large Program imaged a total noncontiguous area of 53pc2, covering 15 extreme, nearby protoclusters of the Milky Way. They were selected to span relevant early protocluster evolutionary stages. Our 1.3mm and 3mm observations provide continuum images that are homogeneously sensitive to point-like cores with masses of 0.2 and 0.6Msun, respectively, with a matched spatial resolution of 2000au. We also detect lines that probe the protocluster structure, kinematics, chemistry, and feedback over scales from clouds to filaments to cores. We classify ALMA-IMF protoclusters as Young, Intermediate, or Evolved based on the amount of dense gas in the cloud that has potentially been impacted by HII regions. The ALMA-IMF catalog contains 700 cores that span a mass range of 0.15-250Msun at a typical size of 2100au. We show that this core sample has no significant distance bias and can be used to build core mass functions at similar physical scales. Significant gas motions, which we highlight here in the G353.41 region, are traced down to core scales and can be used to look for inflowing gas streamers and to quantify the impact of the possible associated core mass growth on the shape of the CMF with time. Our first analysis does not reveal any significant evolution of the matter concentration from clouds to cores or from the youngest to more evolved protoclusters, indicating that cloud dynamical evolution and stellar feedback have for the moment only had a slight effect on the structure of high-density gas in our sample. Furthermore, the first-look analysis of the line richness toward bright cores indicates that the survey encompasses several tens of hot cores, of which we highlight the most massive in the G351.77 cloud. Their homogeneous characterization can be used to constrain the emerging molecular complexity in protostars of high to intermediate masses. △ Less

Submitted 15 December, 2021; originally announced December 2021.

Comments: 27 pages, 9 figures, accepted for publication by A&A

Journal ref: A&A 662, A8 (2022)

HCN/HNC chemistry in shocks: a study of L1157-B1 with ASAI

Authors: B. Lefloch, G. Busquet, S. Viti, C. Vastel, E. Mendoza, M. Benedettini, C. Codella, L. Podio, A. Schutzer, P. R. Rivera-Ortiz, J. R. D. Lépine, R. Bachiller

Abstract: HCN and its isomer HNC play an important role in molecular cloud chemistry and the formation of more complex molecules. We investigate here the impact of protostellar shocks on the HCN and HNC abundances from high-sensitivity IRAM 30m observations of the prototypical shock region L1157-B1 and the envelope of the associated Class 0 protostar, as a proxy for the pre-shock gas. The isotopologues H… ▽ More HCN and its isomer HNC play an important role in molecular cloud chemistry and the formation of more complex molecules. We investigate here the impact of protostellar shocks on the HCN and HNC abundances from high-sensitivity IRAM 30m observations of the prototypical shock region L1157-B1 and the envelope of the associated Class 0 protostar, as a proxy for the pre-shock gas. The isotopologues H$^{12}$CN, HN$^{12}$C, H$^{13}$CN, HN$^{13}$C, HC$^{15}$N, H$^{15}$NC, DCN and DNC were all detected towards both regions. Abundances and excitation conditions were obtained from radiative transfer analysis of molecular line emission under the assumption of Local Thermodynamical Equilibrium. In the pre-shock gas, the abundances of the HCN and HNC isotopologues are similar to those encountered in dark clouds, with a HCN/HNC abundance ratio $\approx 1$ for all isotopologues. A strong D-enrichment (D/H$\approx 0.06$) is measured in the pre-shock gas. There is no evidence of $^{15}$N fractionation neither in the quiescent nor in the shocked gas. At the passage of the shock, the HCN and HNC abundances increase in the gas phase in different manners so that the HCN/HNC relative abundance ratio increases by a factor 20. The gas-grain chemical and shock model UCLCHEM allows us to reproduce the observed trends for a C-type shock with pre-shock density $n$(H)= $10^5$cm$^{-3}$ and shock velocity $V_s= 40$km/s. We conclude that the HCN/HNC variations across the shock are mainly caused by the sputtering of the grain mantle material in relation with the history of the grain ices. △ Less

Submitted 22 July, 2021; originally announced July 2021.

Comments: 14 pages, 7 figures. Accepted for publication in MNRAS

Seed of Life in Space (SOLIS) XI. First measurement of nitrogen fractionation in shocked clumps of the L1157 protostellar outflow

Authors: M. Benedettini, S. Viti, C. Codella, C. Ceccarelli, R. Neri, A. Lopez-Sepulcre, E. Bianchi, G. Busquet, P. Caselli, F. Fontani, B. Lefloch, L. Podio, S. Spezzano, C. Vastel

Abstract: The isotopic ratio of nitrogen presents a wide range of values in the Solar System and in star forming system whose origin is still unclear. Chemical reactions in the gas phase are one of the possible processes that could modify the $^{14}$N/$^{15}$N ratio. We aim at investigating if and how the passage of a shock wave in the interstellar medium, can affect the relative fraction of nitrogen isotop… ▽ More The isotopic ratio of nitrogen presents a wide range of values in the Solar System and in star forming system whose origin is still unclear. Chemical reactions in the gas phase are one of the possible processes that could modify the $^{14}$N/$^{15}$N ratio. We aim at investigating if and how the passage of a shock wave in the interstellar medium, can affect the relative fraction of nitrogen isotopes. The ideal place for such a study is the L1157 outflow, where several shocked clumps are present. We present the first measurement of the $^{14}$N/$^{15}$N ratio in the two shocked clumps, B1 and B0, of the protostellar outflow L1157, derived from the interferomteric maps of the H$^{13}$CN(1-0) and the HC$^{15}$N(1-0) lines. In B1, we find that the H$^{13}$CN(1-0) and HC$^{15}$N(1-0) emission traces the front of the clump, with averaged column density of $N$(H$^{13}$CN) $\sim$ 7$\times$10$^{12}$ cm$^{-2}$ and $N$(HC$^{15}$N) $\sim$ 2$\times$10$^{12}$ cm$^{-2}$. In this region the ratio H$^{13}$CN(1-0)/HC$^{15}$N(1-0) is quite uniform with an average value of $\sim$ 5$\pm$1. The same average value is also measured in the smaller clump B0e. Assuming the standard $^{12}$C/$^{13}$C = 68, we obtain $^{14}$N/$^{15}$N = 340$\pm$70, similar to those usually found in prestellar cores and protostars. We analysed the prediction of a chemical shock model for several shock conditions and we found that the nitrogen and carbon fractionations do not vary much for the first period after the shock. The observed H$^{13}$CN/HC$^{15}$N can be reproduced by a non-dissociative, C-type shock with parameters in agreement with previous modelling of L1157-B1. Both observations and chemical models indicate that the shock propagation does not affect the nitrogen isotopic ratio that remains similar to that measured in lower temperature gas in prestellar cores and in protostellar envelopes. △ Less

Submitted 2 April, 2021; v1 submitted 30 November, 2020; originally announced November 2020.

Comments: 7 pages, 6 figures

Journal ref: A&A 645, A91 (2021)

VLA ammonia observations of L1287: analysis of the "Guitar" core and two filaments

Authors: Inma Sepúlveda, Robert Estalella, Guillem Anglada, Rosario López, Angels Riera, Gemma Busquet, Aina Palau, José M. Torrelles, Luis F. Rodríguez

Abstract: The present work aims at studying the dense gas of the molecular cloud LDN 1287 (L1287), which harbors a double FU Ori system, an energetic molecular outflow and a still-forming cluster of deeply embedded low-mass, young stellar objects, showing a high level of fragmentation. We present optical Halpha and [SII], and VLA NH$_3$ (1,1) and (2,2) observations with an angular resolution of ~3.5''. The… ▽ More The present work aims at studying the dense gas of the molecular cloud LDN 1287 (L1287), which harbors a double FU Ori system, an energetic molecular outflow and a still-forming cluster of deeply embedded low-mass, young stellar objects, showing a high level of fragmentation. We present optical Halpha and [SII], and VLA NH$_3$ (1,1) and (2,2) observations with an angular resolution of ~3.5''. The observed NH$_3$ spectra have been analyzed with the HfS tool, fitting simultaneously three different velocity components. The NH$_3$ emission from L1287 comes from four different structures: a core associated with RNO 1, a guitar-shaped core (the "Guitar") and two interlaced filaments (the Blue and Red Filaments) roughly centered towards the binary FU Ori system RNO 1B/C and its associated cluster. Regarding the Guitar Core, there are clear signatures of gas infall onto a central mass that has been estimated to be ~2.1 $M_\odot$. Regarding the two filaments, they have radii ~0.03 pc, masses per unit length ~50 $M_\odot$ pc$^{-1}$, and are near isothermal equilibrium. A central cavity, probably related with the outflow, and also traced by the Halpha and [SII] emission, is identified, with several young stellar objects near its inner walls. Both filaments show clear signs of perturbation by the high-velocity gas of the outflows driven by one or several young stellar objects of the cluster. The Blue and Red filaments are coherent in velocity and have nearly subsonic gas motions, except at the position of the embedded sources. Velocity gradients across the Blue Filament can be interpreted either as infalling material onto the filament or rotation. Velocity gradients along the filaments are interpreted as infall motions towards a gravitational well at the intersection of both filaments. △ Less

Submitted 3 November, 2020; originally announced November 2020.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 644, A128 (2020)

The role of the magnetic field in the fragmentation process: the case of G14.225-0.506

Authors: N. Añez-López, G. Busquet, P. M. Koch, J. M. Girart, H. B. Liu, F. Santos, N. L. Chapman, G. Novak, A. Palau, P. T. P. Ho, Q. Zhang

Abstract: B-fields are predicted to play a role in the formation of filamentary structures and their fragmentation process. We aim at investigating the role of the B-field in the process of core fragmentation toward the hub-filament systems in the IRDC G14.2, which present different fragmentation level. We performed observations of the thermal dust polarization at 350 μm using the CSO toward the hubs. We ap… ▽ More B-fields are predicted to play a role in the formation of filamentary structures and their fragmentation process. We aim at investigating the role of the B-field in the process of core fragmentation toward the hub-filament systems in the IRDC G14.2, which present different fragmentation level. We performed observations of the thermal dust polarization at 350 μm using the CSO toward the hubs. We applied the polarization--intensity-gradient method to estimate the significance of the B-field over the G-force. The B-field in Hub-N shows a uniform structure along the E-W orientation, perpendicular to the major axis of the hub-filament system. The I-gradient in Hub-N displays a local minimum coinciding with the dust core MM1a detected with interferometric observations. The B-field orientation is perturbed when approaching the dust core. Hub-S shows 2 local minima, reflecting the bimodal distribution of the B-field. In Hub-N, both E and W of the hub-filament system, the I-gradient and the B-field are parallel whereas they tend to be perpendicular when penetrating the filaments and hub. The analysis of the δ- and Σ B-maps indicate that, the B-field cannot prevent the collapse, suggesting that the B-field is initially dragged by the infalling motion and aligned with it, or is channeling material toward the central ridge from both sides. Values of Σ B > 1 are found toward a N-S ridge encompassing the dust emission peak, indicating that in this region B-field dominates over G-force, or that with the current angular resolution we cannot resolve an hypothetical more complex structure. We estimated the B-field strength, the MtF ratio and the A-M number, and found differences between the 2 hubs. The different levels of fragmentation observed in these 2 hubs could arise from the differences in the B-field properties rather than from different intensity of the G-field. △ Less

Submitted 6 November, 2020; v1 submitted 26 October, 2020; originally announced October 2020.

Comments: 14 pages, 9 figures

Journal ref: A&A 644, A52 (2020)

Simultaneous evolution of the virial parameter and star formation rate in molecular clumps undergoing global hierarchical collapse

Authors: Vianey Camacho, Enrique Vázquez-Semadeni, Aina Palau, Gemma Busquet, Manuel Zamora-Avilés

Abstract: We compare dense clumps and cores in a numerical simulation of molecular clouds (MCs) undergoing global hierarchical collapse (GHC) to observations in two MCs at different evolutionary stages, the Pipe and the G14.225 clouds, to test the ability of the GHC scenario to follow the early evolution of the energy budget and star formation activity of these structures. In the simulation, we select a reg… ▽ More We compare dense clumps and cores in a numerical simulation of molecular clouds (MCs) undergoing global hierarchical collapse (GHC) to observations in two MCs at different evolutionary stages, the Pipe and the G14.225 clouds, to test the ability of the GHC scenario to follow the early evolution of the energy budget and star formation activity of these structures. In the simulation, we select a region that contains cores of sizes and densities similar to the Pipe cores, and find that it evolves through accretion, developing substructure similar to that of G14.225 cloud after $\sim 1.6$ Myr. Within this region, we follow the evolution of the Larson ratio $\mathcal{L} \equiv σ_{\rm v}/R^{1/2}$, where $σ_{\rm v}$ is the velocity dispersion and $R$ is the size, the virial parameter $α$, and the star formation activity of the cores/clumps. In the simulation, we find that as the region evolves: $i)$ its clumps have $\mathcal{L}$ and $α$ values first consistent with those of the Pipe substructures and later with those of G14.225; $ii)$ the individual cores first exhibit a decrease in $α$ followed by an increase when star formation begins; $iii)$ collectively, the ensemble of cores/clumps reproduces the observed trend of lower $α$ for higher-mass objects, and $iv)$ the star formation rate and star formation efficiency increase monotonically. We suggest that this evolution is due to the simultaneous loss of externally-driven compressive kinetic energy and increase of the self-gravity-driven motions. We conclude that the GHC scenario provides a realistic description of the evolution of the energy budget of the clouds' substructure at early times, which occurs simultaneously with an evolution of the star formation activity. △ Less

Submitted 23 September, 2020; v1 submitted 4 February, 2020; originally announced February 2020.

Comments: Accepted for publication in ApJ

Modeling the accretion disk around the high-mass protostar GGD 27-MM1

Authors: N. Añez-Lopez, M. Osorio, G. Busquet, J. M. Girart, E. Macıas, C. Carrasco-Gonzalez, S. Curiel, R. Estalella, M. Fernandez-Lopez, R. Galvan-Madrid, J. Kwon, J. M. Torrelles

Abstract: Recent high-angular resolution (40 mas) ALMA observations at 1.14 mm resolve a compact (R~200 au) flattened dust structure perpendicular to the HH 80-81 jet emanating from the GGD 27-MM1 high-mass protostar, making it a robust candidate for a true accretion disk. The jet/disk system (HH 80-81 / GGD 27-MM1) resemble those found in association with low- and intermediate-mass protostars. We present r… ▽ More Recent high-angular resolution (40 mas) ALMA observations at 1.14 mm resolve a compact (R~200 au) flattened dust structure perpendicular to the HH 80-81 jet emanating from the GGD 27-MM1 high-mass protostar, making it a robust candidate for a true accretion disk. The jet/disk system (HH 80-81 / GGD 27-MM1) resemble those found in association with low- and intermediate-mass protostars. We present radiative transfer models that fit the 1.14 mm ALMA dust image of this disk which allow us to obtain its physical parameters and predict its density and temperature structure. Our results indicate that this accretion disk is compact (Rdisk~170 au) and massive (5Msun), about 20% of the stellar mass of 20 Msun. We estimate the total dynamical mass of the star-disk system from the molecular line emission finding a range between 21 and 30 Msun, which is consistent with our model. We fit the density and temperature structures found by our model with power law functions. These results suggest that accretion disks around massive stars are more massive and hotter than their low-mass siblings, but they still are quite stable. We also compare the temperature distribution in the GGD 27-MM1 disk with that found in low- and intermediate-mass stars and discuss possible implications on the water snow line. We have also carried about a study of the distance based on Gaia DR2 data and the population of young stellar objects (YSOs) in this region, and from the extinction maps. We conclude that the source distance is in within 1.2 and 1.4 kpc, closer than what was derived in previous studies (1.7 kpc). △ Less

Submitted 27 November, 2019; originally announced November 2019.

Comments: 20 pages, 15 figures

Filamentary Accretion Flows in the Infrared Dark Cloud G14.225-0.506 Revealed by ALMA

Authors: Huei-Ru Vivien Chen, Qizhou Zhang, M. C. H. Wright, Gemma Busquet, Yuxin Lin, Hauyu Baobab Liu, F. A. Olguin, Patricio Sanhueza, Fumitaka Nakamura, Aina Palau, Satoshi Ohashi, Ken'ichi Tatematsu, Li-Wen Liao

Abstract: Filaments are ubiquitous structures in molecular clouds and play an important role in the mass assembly of stars. We present results of dynamical stability analyses for filaments in the infrared dark cloud G14.225$-$0.506, where a delayed onset of massive star formation was reported in the two hubs at the convergence of multiple filaments of parsec length. Full-synthesis imaging is performed with… ▽ More Filaments are ubiquitous structures in molecular clouds and play an important role in the mass assembly of stars. We present results of dynamical stability analyses for filaments in the infrared dark cloud G14.225$-$0.506, where a delayed onset of massive star formation was reported in the two hubs at the convergence of multiple filaments of parsec length. Full-synthesis imaging is performed with the Atacama Large Millimeter/submillimeter Array (ALMA) to map the $\mathrm{N_2H^+} \; (1-0)$ emission in two hub-filament systems with a spatial resolution of $\sim 0.034 \; \mathrm{pc}$. Kinematics are derived from sophisticated spectral fitting algorithm that accounts for line blending, large optical depth, and multiple velocity components. We identify five velocity coherent filaments and derive their velocity gradients with principal component analysis. The mass accretion rates along the filaments are up to $10^{-4} \; \mathrm{M_\odot \, \mathrm{yr^{-1}}}$ and are significant enough to affect the hub dynamics within one free-fall time ($\sim 10^5 \; \mathrm{yr}$). The $\mathrm{N_2H^+}$ filaments are in equilibrium with virial parameter $α_\mathrm{vir} \sim 1.2$. We compare $α_\mathrm{vir}$ measured in the $\mathrm{N_2H^+}$ filaments, $\mathrm{NH_3}$ filaments, $870 \; μ\mathrm{m}$ dense clumps, and $3 \; \mathrm{mm}$ dense cores. The decreasing trend in $α_\mathrm{vir}$ with decreasing spatial scales persists, suggesting an increasingly important role of gravity at small scales. Meanwhile, $α_\mathrm{vir}$ also decreases with decreasing non-thermal motions. In combination with the absence of high-mass protostars and massive cores, our results are consistent with the global hierarchical collapse scenario. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 31 pages, 21 figures, 2 movies (plus 2 more movies in the published version), accepted for publication in The Astrophysical Journal

Unveiling a cluster of protostellar disks around the massive protostar GGD27 MM1

Authors: G. Busquet, J. M. Girart, R. Estalella, M. Fernández-López, R. Galván-Madrid, G. Anglada, C. Carrasco-González, N. Añez-López, S. Curiel, M. Osorio, L. F. Rodríguez, J. M. Torrelles

Abstract: We used ALMA to observe the star-forming region GGD27 at 1.14 mm with an unprecedented angular resolution, 40 mas (56 au) and sensitivity (0.002 Msun). We detected a cluster of 25 continuum sources, most of which are likely tracing disks around Class 0/I protostars. Excluding the two most massive objects, disks masses are in the range 0.003-0.05 Msun. The analysis of the cluster properties indicat… ▽ More We used ALMA to observe the star-forming region GGD27 at 1.14 mm with an unprecedented angular resolution, 40 mas (56 au) and sensitivity (0.002 Msun). We detected a cluster of 25 continuum sources, most of which are likely tracing disks around Class 0/I protostars. Excluding the two most massive objects, disks masses are in the range 0.003-0.05 Msun. The analysis of the cluster properties indicates that GGD27 displays moderate subclustering. This result combined with the dynamical timescale of the radio jet (10000 years) suggests the youthfulness of the cluster. The lack of disk mass segregation signatures may support this too. We found a clear paucity of disks with Rdisk >100 au. The median value of the radius is 34 au, smaller than the median of 92 au for Taurus but comparable to the value found in Ophiuchus and in the Orion Nebula Cluster. In GGD27 there is no evidence of a distance-dependent disk mass distribution (i. e., disk mass depletion due to external photoevaporation), most likely due to the cluster youth. There is a clear deficit of disks for distances <0.02 pc. Only for distances >0.04 pc stars can form larger and more massive disks, suggesting that dynamical interactions far from the cluster center are weaker, although the small disks found could be the result of disk truncation. This work demonstrates the potential to characterize disks from low-mass YSOs in distant and massive (still deeply embedded) clustered environments. △ Less

Submitted 20 February, 2019; originally announced February 2019.

Comments: Accepted for publication in Astronomy & Astrophysics Letters

Journal ref: A&A 623, L8 (2019)

Extreme fragmentation and complex kinematics at the center of the L1287 cloud

Authors: Carmen Juárez, Hauyu-Baobab Liu, Josep M. Girart, Aina Palau, Gemma Busquet, Roberto Galván-Madrid, Naomi Hirano, Yuxin Lin

Abstract: The filamentary infrared dark cloud L1287 is actively forming a dense cluster of low-mass YSOs at its inner $\sim$0.1 pc region. To help understand the origin of this low-mass YSO cluster, the present work aims at resolving the gas structures and kinematics. We have performed $\sim$1$"$ angular resolution ($\sim$930 AU) SMA observations at $\sim$1.3 mm. From a $\sim$2$"$ resolution 1.3 mm continuu… ▽ More The filamentary infrared dark cloud L1287 is actively forming a dense cluster of low-mass YSOs at its inner $\sim$0.1 pc region. To help understand the origin of this low-mass YSO cluster, the present work aims at resolving the gas structures and kinematics. We have performed $\sim$1$"$ angular resolution ($\sim$930 AU) SMA observations at $\sim$1.3 mm. From a $\sim$2$"$ resolution 1.3 mm continuum image we identified six dense cores, namely SMA1-6 with masses in the range of $\sim0.4-4$ M$_\odot$. From a $\sim$1$"$ resolution 1.3 mm continuum image, we find a high fragmentation level, with 14 compact millimeter sources within 0.1 pc (two of them associated with the known accretion outburst YSOs RNO 1C and RNO 1B). The dense gas tracer DCN (3--2) traces well the dust continuum emission and shows a clear velocity gradient along the NW-SE direction centered at SMA3. There is another velocity gradient with opposite direction around the most luminous YSO IRAS 00338+6312. The fragmentation within 0.1 pc in L1287 is very high compared to other regions at the same spatial scales. The incoherent motions of dense gas flows are sometimes interpreted by being influenced by (proto)stellar feedback (e.g., outflows), which is not yet ruled out in this particular target source. The directions of the velocity gradients traced by DCN are approximately perpendicular to those of the dominant CO outflow(s). Therefore, we alternatively hypothesize that the velocity gradients revealed by DCN trace the convergence from the $\gtrsim$0.1 pc scales infalling motion towards the rotational motions around the more compact ($\sim0.02$ pc) sources. This global molecular gas converging flow may feed the formation of the dense low-mass YSO cluster. IRAS 00338+6312 is the most likely powering source of the dominant CO outflow. A compact blue-shifted outflow from RNO 1C is also identified. △ Less

Submitted 20 November, 2018; originally announced November 2018.

Comments: 14 pages, 11 figures. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 621, A140 (2019)

Distribution of Serpens South protostars revealed with ALMA

Authors: Adele L. Plunkett, Manuel Fernández-López, Héctor G. Arce, Gemma Busquet, Diego Mardones, Michael M. Dunham

Abstract: Aims: We investigated the masses and spatial distributions of pre-stellar and protostellar candidates in the young, low-mass star forming region Serpens South, where active star formation is known to occur along a predominant filamentary structure. Previous observations used to study these distributions have been limited by two important observational factors: (1) sensitivity limits that leave the… ▽ More Aims: We investigated the masses and spatial distributions of pre-stellar and protostellar candidates in the young, low-mass star forming region Serpens South, where active star formation is known to occur along a predominant filamentary structure. Previous observations used to study these distributions have been limited by two important observational factors: (1) sensitivity limits that leave the lowest-mass sources undetected, or (2) resolution limits that cannot distinguish binaries and/or cluster members in close proximity. Methods: Recent millimeter-wavelength interferometry observations can now uncover faint and/or compact sources in order to study a more complete population of protostars, especially in nearby ($D<500$ pc) clusters. Here we present ALMA observations of 1 mm (Band 6) continuum in a $3 \times 2$ arcminutes region at the center of Serpens South. Our angular resolution of $\sim1$ arcsec is equivalent to $\sim400$ au, corresponding to scales of envelopes and/or disks of protostellar sources. Results: We detect 52 sources with 1 mm continuum, and we measure masses of $0.002 - 0.9$ solar masses corresponding to gas and dust in the disk and/or envelope of the protostellar system. For the deeply embedded (youngest) sources with no IR counterparts, we find evidence of mass segregation and clustering according to: the Minimum Spanning Tree method, distribution of projected separations between unique sources, and concentration of higher-mass sources near to the dense gas at the cluster center. Conclusions: The mass segregation of the mm sources is likely primordial rather than dynamical given the young age of this cluster, compared with segregation time. This is the first case to show this for mm sources in a low-mass protostellar cluster environment. △ Less

Submitted 6 April, 2018; originally announced April 2018.

Comments: Accepted to A&A, FITS data to be available at CDS

Journal ref: A&A 615, A9 (2018)

Resolving the polarized dust emission of the disk around the massive star powering the HH~80-81 radio jet

Authors: J. M. Girart, M. Fernandez-Lopez, Z. -Y. Li, H. Yang, R. Estalella, G. Anglada, N. Añez-Lopez, G. Busquet, C. Carrasco-Gonzalez, S. Curiel, R. Galvan-Madrid, J. F. Gomez, I. de Gregorio-Monsalvo, I. Jimenez-Serra, R. Krasnopolsky, J. Marti, M. Osorio, M. Padovani, R. Rao, L. F. Rodriguez, J. M. Torrelles

Abstract: Here we present deep (16 mumJy), very high (40 mas) angular resolution 1.14 mm, polarimetric, Atacama Large Millimeter/submillimeter Array (ALMA) observations towards the massive protostar driving the HH 80-81 radio jet. The observations clearly resolve the disk oriented perpendicular to the radio jet, with a radius of ~0.171 arcsec (~291 au at 1.7 kpc distance). The continuum brightness temperatu… ▽ More Here we present deep (16 mumJy), very high (40 mas) angular resolution 1.14 mm, polarimetric, Atacama Large Millimeter/submillimeter Array (ALMA) observations towards the massive protostar driving the HH 80-81 radio jet. The observations clearly resolve the disk oriented perpendicular to the radio jet, with a radius of ~0.171 arcsec (~291 au at 1.7 kpc distance). The continuum brightness temperature, the intensity profile, and the polarization properties clearly indicate that the disk is optically thick for a radius of R<170 au. The linear polarization of the dust emission is detected almost all along the disk and its properties suggest that dust polarization is produced mainly by self-scattering. However, the polarization pattern presents a clear differentiation between the inner (optically thick) part of the disk and the outer (optically thin) region of the disk, with a sharp transition that occurs at a radius of 0.1 arcsec (~170 au). The polarization characteristics of the inner disk suggest that dust settling has not occurred yet with a maximum dust grain size between 50 and 500 mum. The outer part of the disk has a clear azimuthal pattern but with a significantly higher polarization fraction compared to the inner disk. This pattern is broadly consistent with self-scattering of a radiation field that is beamed radially outward, as expected in the optically thin outer region, although contribution from non-spherical grains aligned with respect to the radiative flux cannot be excluded. △ Less

Submitted 16 March, 2018; originally announced March 2018.

Comments: Accepted for publication in the Astrophysical Journal Letters

The circumestellar disk of the B0 protostar powering the HH 80-81 radio jet

Authors: J. M. Girart, R. Estalella, M. Fernández-López, S. Curiel, P Frau, R. Galvan-Madrid, R. Rao, G. Busquet, C. Juárez

Abstract: We present subarcsecond angular resolution observations carried out with the Submillimeter Array (SMA) at 880 $μ$m centered at the B0-type protostar GGD27~MM1, the driving source of the parsec scale HH 80-81 jet. We constrain its polarized continuum emission to be $\lesssim0.8\%$ at this wavelength. Its submm spectrum is dominated by sulfur-bearing species tracing a rotating disk--like structure (… ▽ More We present subarcsecond angular resolution observations carried out with the Submillimeter Array (SMA) at 880 $μ$m centered at the B0-type protostar GGD27~MM1, the driving source of the parsec scale HH 80-81 jet. We constrain its polarized continuum emission to be $\lesssim0.8\%$ at this wavelength. Its submm spectrum is dominated by sulfur-bearing species tracing a rotating disk--like structure (SO and SO$_2$ isotopologues mainly), but also shows HCN-bearing and CH$_3$OH lines, which trace the disk and the outflow cavity walls excavated by the HH 80-81 jet. The presence of many sulfurated lines could indicate the presence of shocked gas at the disk's centrifugal barrier or that MM1 is a hot core at an evolved stage. The resolved SO$_2$ emission traces very well the disk kinematics and we fit the SMA observations using a thin-disk Keplerian model, which gives the inclination (47$^{\circ}$), the inner ($\lesssim170$ AU) and outer ($\sim950-1300$~AU) radii and the disk's rotation velocity (3.4 km s$^{-1}$ at a putative radius of 1700 AU). We roughly estimate a protostellar dynamical mass of 4-18\msun. MM2 and WMC cores show, comparatively, an almost empty spectra suggesting that they are associated with extended emission detected in previous low-angular resolution observations, and therefore indicating youth (MM2) or the presence of a less massive object (WMC). △ Less

Submitted 20 July, 2017; originally announced July 2017.

Comments: 11 pages, 8 figures. Accepted for publication in the Astrophysical Journal

The L1157-B1 astrochemical laboratory: testing the origin of DCN

Authors: G. Busquet, F. Fontani, S. Viti, C. Codella, B. Lefloch, M. Benedettini, C. Ceccarelli

Abstract: L1157-B1 is the brightest shocked region of the large-scale molecular outflow, considered the prototype of chemically rich outflows, being the ideal laboratory to study how shocks affect the molecular gas. Several deuterated molecules have been previously detected with the IRAM 30m, most of them formed on grain mantles and then released into the gas phase due to the shock. We aim to observationall… ▽ More L1157-B1 is the brightest shocked region of the large-scale molecular outflow, considered the prototype of chemically rich outflows, being the ideal laboratory to study how shocks affect the molecular gas. Several deuterated molecules have been previously detected with the IRAM 30m, most of them formed on grain mantles and then released into the gas phase due to the shock. We aim to observationally investigate the role of the different chemical processes at work that lead to formation the of DCN and test the predictions of the chemical models for its formation. We performed high-angular resolution observations with NOEMA of the DCN(2-1) and H13CN(2-1) lines to compute the deuterated fraction, Dfrac(HCN). We detected emission of DCN(2-1) and H13CN(2-1) arising from L1157-B1 shock. Dfrac(HCN) is ~4x10$^{-3}$ and given the uncertainties, we did not find significant variations across the bow-shock. Contrary to HDCO, whose emission delineates the region of impact between the jet and the ambient material, DCN is more widespread and not limited to the impact region. This is consistent with the idea that gas-phase chemistry is playing a major role in the deuteration of HCN in the head of the bow-shock, where HDCO is undetected as it is a product of grain-surface chemistry. The spectra of DCN and H13CN match the spectral signature of the outflow cavity walls, suggesting that their emission result from shocked gas. The analysis of the time dependent gas-grain chemical model UCL-CHEM coupled with a C-type shock model shows that the observed Dfrac(HCN) is reached during the post-shock phase, matching the dynamical timescale of the shock. Our results indicate that the presence of DCN in L1157-B1 is a combination of gas-phase chemistry that produces the widespread DCN emission, dominating in the head of the bow-shock, and sputtering from grain mantles toward the jet impact region. △ Less

Submitted 14 July, 2017; v1 submitted 27 June, 2017; originally announced June 2017.

Comments: Accepted for publication in A&A. 7 pages, 5 Figures, 1 Table

Journal ref: A&A 604, A20 (2017)

Thermal Jeans fragmentation within 1000 AU in OMC-1S

Authors: Aina Palau, Luis A. Zapata, Carlos G. Roman-Zuniga, Alvaro Sanchez-Monge, Robert Estalella, Gemma Busquet, Josep M. Girart, Asuncion Fuente, Benoit Commercon

Abstract: We present subarcsecond 1.3 mm continuum ALMA observations towards the Orion Molecular Cloud 1 South (OMC-1S) region, down to a spatial resolution of 74 AU, which reveal a total of 31 continuum sources. We also present subarcsecond 7 mm continuum VLA observations of the same region, which allow to further study fragmentation down to a spatial resolution of 40 AU. By applying a Mean Surface Density… ▽ More We present subarcsecond 1.3 mm continuum ALMA observations towards the Orion Molecular Cloud 1 South (OMC-1S) region, down to a spatial resolution of 74 AU, which reveal a total of 31 continuum sources. We also present subarcsecond 7 mm continuum VLA observations of the same region, which allow to further study fragmentation down to a spatial resolution of 40 AU. By applying a Mean Surface Density of Companions method we find a characteristic spatial scale at ~560 AU, and we use this spatial scale to define the boundary of 19 `cores' in OMC-1S as groupings of millimeter sources. We find an additional characteristic spatial scale at ~2900 AU, which is the typical scale of the filaments in OMC-1S, suggesting a two-level fragmentation process. We measured the fragmentation level within each core and find a higher fragmentation towards the southern filament. In addition, the cores of the southern filament are also the densest (within 1100 AU) cores in OMC-1S. This is fully consistent with previous studies of fragmentation at spatial scales one order of magnitude larger, and suggests that fragmentation down to 40 AU seems to be governed by thermal Jeans processes in OMC-1S. △ Less

Submitted 20 February, 2018; v1 submitted 14 June, 2017; originally announced June 2017.

Comments: Accepted to ApJ

Cloud structure of three Galactic infrared dark star-forming regions from combining ground and space based bolometric observations

Authors: Yuxin Lin, Hauyu Baobab Liu, James E. Dale, Di Li, Gemma Busquet, Zhi-Yu Zhang, Adam Ginsburg, Roberto Galvan-Madrid, Attila Kovacs, Eric Koch, Lei Qian, Ke Wang, Steve Longmore, Huei-Ru Chen, Daniel Walker

Abstract: We have modified the iterative procedure introduced by Lin et al. (2016), to systematically combine the submm images taken from ground based (e.g., CSO, JCMT, APEX) and space (e.g., Herschel, Planck) telescopes. We applied the updated procedure to observations of three well studied Infrared Dark Clouds (IRDCs): G11.11-0.12, G14.225-0.506 and G28.34+0.06, and then performed single-component, modifi… ▽ More We have modified the iterative procedure introduced by Lin et al. (2016), to systematically combine the submm images taken from ground based (e.g., CSO, JCMT, APEX) and space (e.g., Herschel, Planck) telescopes. We applied the updated procedure to observations of three well studied Infrared Dark Clouds (IRDCs): G11.11-0.12, G14.225-0.506 and G28.34+0.06, and then performed single-component, modified black-body fits to derive $\sim$10$"$ resolution dust temperature and column density maps. The derived column density maps show that these three IRDCs exhibit complex filamentary structures embedding with rich clumps/cores. We compared the column density probability distribution functions (N-PDFs) and two-point correlation (2PT) functions of the column density field between these IRDCs with several OB cluster-forming regions. Based on the observed correlation and measurements, and complementary hydrodynamical simulations for a 10$^{4}$ $\rm M_{\odot}$ molecular cloud, we hypothesize that cloud evolution can be better characterized by the evolution of the (column) density distribution function and the relative power of dense structures as a function of spatial scales, rather than merely based on the presence of star-forming activity. Based on the small analyzed sample, we propose four evolutionary stages, namely: {\it cloud integration, stellar assembly, cloud pre-dispersal and dispersed-cloud.} The initial {\it cloud integration} stage and the final {\it dispersed cloud} stage may be distinguished from the two intermediate stages by a steeper than $-$4 power-law index of the N-PDF. The {\it cloud integration} stage and the subsequent {\it stellar assembly} stage are further distinguished from each other by the larger luminosity-to-mass ratio ($>$40 $\rm L_{\odot}/M_{\odot}$) of the latter. △ Less

Submitted 21 April, 2017; originally announced April 2017.

Comments: 14 pages, 13 figures; accepted for publication in ApJ

The shocked gas of the BHR71 outflow observed by Herschel: indirect evidence for an atomic jet

Authors: M. Benedettini, A. Gusdorf, B. Nisini, B. Lefloch, S. Anderl, G. Busquet, C. Ceccarelli, C. Codella, S. Leurini, L. Podio

Abstract: In the BHR71 region, two low-mass protostars drive two distinguishable outflows. They constitute an ideal laboratory to investigate the effects of shock chemistry and the mechanisms that led to their formation. We aim to define the morphology of the warm gas component of the BHR 71 outflow and at modelling its shocked component. We present the first far infrared Herschel images of the BHR71 outflo… ▽ More In the BHR71 region, two low-mass protostars drive two distinguishable outflows. They constitute an ideal laboratory to investigate the effects of shock chemistry and the mechanisms that led to their formation. We aim to define the morphology of the warm gas component of the BHR 71 outflow and at modelling its shocked component. We present the first far infrared Herschel images of the BHR71 outflow in the CO(14-13), H$_2$O (2$_{21}$-1$_{10}$), H$_2$O (2$_{12}$-1$_{01}$) and [OI] 145 $μ$m, lines, revealing the presence of several knots of warm, shocked gas associated with fast outflowing gas. In two of these knots we performed a detailed study of the physical conditions by comparing a large set of transitions from several molecules to a grid of shock models. Herschel lines ratios in the outflow knots are quite similar, showing that the excitation conditions of the fast moving gas do not change significantly within the first $\sim$ 0.068 pc of the outflow, apart at the extremity of the southern blue-shifted lobe that is expanding outside the molecular cloud. Rotational diagram, spectral line profile and LVG analysis of the CO lines in knot A show the presence of two gas components: one extended, cold ($T\sim$80 K) and dense ($n$(H$_2$) = 3$\times$10$^5$-4$\times$10$^6$ cm$^{-3}$) and another compact (18 arcsec), warm ($T$ = 1700-2200 K) with slightly lower density ($n$(H$_2$) = (2-6)$\times$10$^4$ cm$^{-3}$). In the two brightest knots (where we performed shock modelling) we found that H$_2$ and CO are well fitted with non-stationary (young) shocks. These models, however, significantly underestimate the observed fluxes of [OI] and OH lines, but are not too far off those of H$_2$O, calling for an additional, possibly dissociative, J-type shock component. Our modelling indirectly suggests that an additional shock component exists, possibly a remnant of the primary jet △ Less

Submitted 16 January, 2017; originally announced January 2017.

Comments: 16 pages, 15 figures

Dense core properties in the Infrared Dark cloud G14.225-0.506 revealed by ALMA

Authors: Satoshi Ohashi, Patricio Sanhueza, Huei-Ru Vivien Chen, Qizhou Zhang, Gemma Busquet, Fumitaka Nakamura, Aina Palau, Ken'ichi Tatematsu

Abstract: We have performed a dense core survey toward the Infrared Dark Cloud G14.225-0.506 at 3 mm continuum emission with the Atacama Large Millimeter/Submillimeter Array (ALMA). This survey covers the two hub-filament systems with an angular resolution of $\sim3$\arcsec ($\sim0.03$ pc). We identified 48 dense cores. Twenty out of the 48 cores are protostellar due to their association with young stellar… ▽ More We have performed a dense core survey toward the Infrared Dark Cloud G14.225-0.506 at 3 mm continuum emission with the Atacama Large Millimeter/Submillimeter Array (ALMA). This survey covers the two hub-filament systems with an angular resolution of $\sim3$\arcsec ($\sim0.03$ pc). We identified 48 dense cores. Twenty out of the 48 cores are protostellar due to their association with young stellar objects (YSOs) and/or X-ray point-sources, while the other 28 cores are likely prestellar and unrelated with known IR or X-ray emission. Using APEX 870 $μ$m continuum emission, we also identified the 18 clumps hosting these cores. Through virial analysis using the ALMA N$_2$H$^+$ and VLA/Effelsberg NH$_3$ molecular line data, we found a decreasing trend in the virial parameter with decreasing scales from filaments to clumps, and then to cores. The virial parameters of $0.1-1.3$ in cores, indicate that cores are likely undergoing dynamical collapse. The cumulative Core Mass Function (CMF) for the prestellar cores candidates has a power law index of $α=1.6$, with masses ranging from 1.5 to 22 $M_\odot$. We find no massive prestellar or protostellar cores. Previous studies suggest that massive O-tpye stars have not been produced yet in this region. Therefore, high-mass stars should be formed in the prestellar cores by accreting a significant amount of gas from the surrounding medium. Another possibility is that low-mass YSOs become massive by accreting from their parent cores that are fed by filaments. These two possibilities might be consistent with the scenario of global hierarchical collapse. △ Less

Submitted 26 October, 2016; originally announced October 2016.

Comments: 12 pages, 7 figures, 2 tables, accepted for publication in ApJ

Magnetically Dominated Parallel Interstellar Filaments at the Infrared Dark Cloud G14.225-0.506

Authors: Fabio P. Santos, Gemma Busquet, Gabriel A. P. Franco, Josep M. Girart, Qizhou Zhang

Abstract: The G14.225-0.506 infrared dark cloud (IRDC G14.2) displays a remarkable complex of parallel dense molecular filaments projected on the plane of the sky. Previous dust emission and molecular-line studies have speculated whether magnetic fields could have played an important role in the formation of such long-shaped structures, which are hosts to numerous young stellar sources. In this work we have… ▽ More The G14.225-0.506 infrared dark cloud (IRDC G14.2) displays a remarkable complex of parallel dense molecular filaments projected on the plane of the sky. Previous dust emission and molecular-line studies have speculated whether magnetic fields could have played an important role in the formation of such long-shaped structures, which are hosts to numerous young stellar sources. In this work we have conducted a vast polarimetric survey at optical and near-infrared wavelengths in order to study the morphology of magnetic field lines in IRDC G14.2 through the observation of background stars. The orientation of interstellar polarization, which traces magnetic field lines, is perpendicular to most of the filamentary features within the cloud. Additionally, the larger-scale molecular cloud as a whole exhibits an elongated shape also perpendicular to magnetic fields. Estimates of magnetic field strengths indicate values in the range $320 - 550\,μ$G, which allows sub-alfvénic conditions, but does not prevent the gravitational collapse of hub-filament structures, which in general are close to the critical state. These characteristics suggest that magnetic fields played the main role in regulating the collapse from large to small scales, leading to the formation of series of parallel elongated structures. The morphology is also consistent with numerical simulations that show how gravitational instabilities develop under strong magnetic fields. Finally, the results corroborate the hypothesis that a strong support from internal magnetic fields might explain why the cloud seems to be contracting on a time scale 2-3 times larger than what is expected from a free-fall collapse. △ Less

Submitted 27 September, 2016; v1 submitted 26 September, 2016; originally announced September 2016.

Comments: Accepted for publication at The Astrophysical Journal

What is controlling the fragmentation process in the Infrared Dark Cloud G14.225-0.506? Differet level of fragmentation in twin hubs

Authors: G. Busquet, R. Estalella, A. Palau, H. B. Liu, Q. Zhang, J. M. Girart, I. de Gregorio-Monsalvo, T. Pillai, G. Anglada, P. T. P. Ho

Abstract: We present observations of the 1.3 mm continuum emission toward hub-N and hub-S of the infrared dark cloud G14.225-0.506 carried out with the Submillimeter Array, together with observations of the dust emission at 870 and 350 microns obtained with APEX and CSO telescopes. The large scale dust emission of both hubs consists of a single peaked clump elongated in the direction of the associated filam… ▽ More We present observations of the 1.3 mm continuum emission toward hub-N and hub-S of the infrared dark cloud G14.225-0.506 carried out with the Submillimeter Array, together with observations of the dust emission at 870 and 350 microns obtained with APEX and CSO telescopes. The large scale dust emission of both hubs consists of a single peaked clump elongated in the direction of the associated filament. At small scales, the SMA images reveal that both hubs fragment into several dust condensations. The fragmentation level was assessed under the same conditions and we found that hub-N presents 4 fragments while hub-S is more fragmented, with 13 fragments identified. We studied the density structure by means of a simultaneous fit of the radial intensity profile at 870 and 350 microns and the spectral energy distribution adopting a Plummer-like function to describe the density structure. The parameters inferred from the model are remarkably similar in both hubs, suggesting that density structure could not be responsible in determining the fragmentation level. We estimated several physical parameters such as the level of turbulence and the magnetic field strength, and we found no significant differences between these hubs. The Jeans analysis indicates that the observed fragmentation is more consistent with thermal Jeans fragmentation compared with a scenario that turbulent support is included. The lower fragmentation level observed in hub-N could be explained in terms of stronger UV radiation effects from a nearby HII region, evolutionary effects, and/or stronger magnetic fields at small scales, a scenario that should be further investigated. △ Less

Submitted 8 February, 2016; originally announced February 2016.

Comments: 14 pages, 5 figures. Accepted for publication in ApJ

Ongoing star formation in the proto-cluster IRAS 22134+5834

Authors: Yuan Wang, Marc Audard, Francesco Fontani, Álvaro Sánchez-Monge, Gemma Busquet, Aina Palau, Henrik Beuther, Jonathan C. Tan, Robert Estalella, Andrea Isella, Frederic Gueth, Izaskun Jiménez-Serra

Abstract: IRAS 22134+5834 was observed in the centimeter with (E)VLA, 3~mm with CARMA, 2~mm with PdBI, and 1.3~mm with SMA, to study the continuum emission as well as the molecular lines, that trace different physical conditions of the gas to study the influence of massive YSOs on nearby starless cores, and the possible implications in the clustered star formation process. The multi-wavelength centimeter co… ▽ More IRAS 22134+5834 was observed in the centimeter with (E)VLA, 3~mm with CARMA, 2~mm with PdBI, and 1.3~mm with SMA, to study the continuum emission as well as the molecular lines, that trace different physical conditions of the gas to study the influence of massive YSOs on nearby starless cores, and the possible implications in the clustered star formation process. The multi-wavelength centimeter continuum observations revealed two radio sources within the cluster, VLA1 and VLA2. VLA1 is considered to be an optically thin UCHII region with a size of 0.01~pc and sits at the edge of the near-infrared (NIR) cluster. The flux of ionizing photons of the VLA1 corresponds to a B1 ZAMS star. VLA2 is associated with an infrared point source and has a negative spectral index. We resolved six millimeter continuum cores at 2~mm, MM2 is associated with the UCHII region VLA1, and other dense cores are distributed around the UCH{\sc ii} region. Two high-mass starless clumps (HMSC), HMSC-E (east) and HMSC-W (west), are detected around the NIR cluster with N$_2$H$^+$(1--0) and NH$_3$ emission, and show different physical and chemical properties. Two N$_2$D$^+$ cores are detected on an NH$_3$ filament close to the UCHII region, with a projected separation of $\sim$8000~AU at the assumed distance of 2.6~kpc. The kinematic properties of the molecular line emission confirm the expansion of the UCHII region and that the molecular cloud around the near infrared (NIR) cluster is also expanding. Our multi-wavelength study has revealed different generations of star formation in IRAS 22134+5834. The formed intermediate- to massive stars show strong impact on nearby starless clumps. We propose that while the stellar wind from the UCHII region and the NIR cluster drives the large scale bubble, the starless clumps and HMPOs formed at the edge of the cluster. △ Less

Submitted 19 November, 2015; v1 submitted 8 October, 2015; originally announced October 2015.

Comments: 22 pages, 17 figures, Accepted by A&A in Oct. 2015

Journal ref: A&A 587, A69 (2016)

The Spanish Square Kilometre Array White Book

Authors: J. A. Acosta-Pulido, I. Agudo, A. Alberdi, J. Alcolea, E. J. Alfaro, A. Alonso-Herrero, G. Anglada, P. Arnalte-Mur, Y. Ascasibar, B. Ascaso, R. Azulay, R. Bachiller, A. Baez-Rubio, E. Battaner, J. Blasco, C. B. Brook, V. Bujarrabal, G. Busquet, M. D. Caballero-Garcia, C. Carrasco-Gonzalez, J. Casares, A. J. Castro-Tirado, L. Colina, F. Colomer, I. de Gregorio-Monsalvo , et al. (94 additional authors not shown)

Abstract: The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on… ▽ More The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project. △ Less

Submitted 17 June, 2015; v1 submitted 10 June, 2015; originally announced June 2015.

Comments: Editors: Miguel A. Perez-Torres (Editor-in-chief), L. Verdes-Montenegro, J. C. Guirado, A. Alberdi, J. Martin-Pintado, R. Bachiller, D. Herranz, J. M. Girart, J. Gorgas, C. Hernandez-Monteagudo, S. Migliari and J. M. Rodriguez Espinosa; 289 pages, published in June 2015 by Sociedad Española de Astronomía

Gravity or turbulence? -III. Evidence of pure thermal Jeans fragmentation at ~0.1 pc scale

Authors: Aina Palau, Javier Ballesteros-Paredes, Enrique Vazquez-Semadeni, Alvaro Sanchez-Monge, Robert Estalella, S. Michael Fall, Luis A. Zapata, Vianey Camacho, Laura Gomez, Raul Naranjo-Romero, Gemma Busquet, Francesco Fontani

Abstract: We combine previously published interferometric and single-dish data of relatively nearby massive dense cores that are actively forming stars to test whether their `fragmentation level' is controlled by turbulent or thermal support. We find no clear correlation between the fragmentation level and velocity dispersion, nor between the observed number of fragments and the number of fragments expected… ▽ More We combine previously published interferometric and single-dish data of relatively nearby massive dense cores that are actively forming stars to test whether their `fragmentation level' is controlled by turbulent or thermal support. We find no clear correlation between the fragmentation level and velocity dispersion, nor between the observed number of fragments and the number of fragments expected when the gravitationally unstable mass is calculated including various prescriptions for `turbulent support'. On the other hand, the best correlation is found for the case of pure thermal Jeans fragmentation, for which we infer a core formation efficiency around 13 per cent, consistent with previous works. We conclude that the dominant factor determining the fragmentation level of star-forming massive dense cores at 0.1 pc scale seems to be thermal Jeans fragmentation. △ Less

Submitted 13 October, 2015; v1 submitted 28 April, 2015; originally announced April 2015.

Comments: accepted in MNRAS

Far-infrared line spectra of Seyfert galaxies from the Herschel-PACS Spectrometer

Authors: Luigi Spinoglio, Miguel Pereira-Santaella, Kalliopi M. Dasyra, Luca Calzoletti, Matthew A. Malkan, Silvia Tommasin, Gemma Busquet

Abstract: We present spectroscopic observations of FIR fine-structure lines of 26 Seyfert galaxies obtained with the Herschel-PACS spectrometer. These observations are complemented by spectroscopy with Spitzer-IRS and Herschel-SPIRE. The ratios of the OIII, NII, SIII and NeV lines have been used to determine electron densities in the ionised gas regions. The CI lines, observed with SPIRE, have been used to… ▽ More We present spectroscopic observations of FIR fine-structure lines of 26 Seyfert galaxies obtained with the Herschel-PACS spectrometer. These observations are complemented by spectroscopy with Spitzer-IRS and Herschel-SPIRE. The ratios of the OIII, NII, SIII and NeV lines have been used to determine electron densities in the ionised gas regions. The CI lines, observed with SPIRE, have been used to measure the densities in the neutral gas, while the OI lines provide a measure of the gas temperature, at densities below 10000 cm-3. Using the OI145/63um and SIII33/18um line ratios we find an anti-correlation of the temperature with the gas density. Using various fine-structure line ratios, we find that density stratification is common in these active galaxies. On average, the electron densities increase with the ionisation potential of the ions producing the NII, SIII and NeV emission. The infrared emission lines arise partly in the Narrow Line Region (NLR) photoionised by the AGN central engine, partly in HII regions photo ionised by hot stars and partly in neutral gas in photo-dissociated regions (PDRs). We attempt to separate the contributions to the line emission produced in these different regions by comparing our emission line ratios to empirical and theoretical values. In particular, we tried to separate the contribution of AGN and star formation by using a combination of Spitzer and Herschel lines, and we found that, besides the well known mid-IR line ratios, the mixed mid-IR/far-IR line ratio of OIII88um/OIV26um can reliably discriminate the two emission regimes, while the far-IR line ratio of CII157um/OI63um is only able to mildly separate the two regimes. By comparing the observed CII157um/NII205um ratio with photoionisation models, we also found that most of the CII emission in the galaxies we examined is due to PDRs. △ Less

Submitted 5 November, 2014; originally announced November 2014.

Comments: Accepted for publication in ApJ Main Journal on the 5th November 2014, 25 pages, 16 figures

The density structure of the L1157 molecular outflow

Authors: A. I. Gómez-Ruiz, C. Codella, B. Lefloch, M. Benedettini, G. Busquet, C. Ceccarelli, B. Nisini, L. Podio, S. Viti

Abstract: We present a multiline CS survey towards the brightest bow-shock B1 in the prototypical chemically active protostellar outflow L1157. We made use of (sub-)mm data obtained in the framework of the Chemical HErschel Surveys of Star forming regions (CHESS) and Astrochemical Surveys at IRAM (ASAI) key science programs. We detected $^{12}$C$^{32}$S, $^{12}$C$^{34}$S, $^{13}$C$^{32}$S, and $^{12}$C… ▽ More We present a multiline CS survey towards the brightest bow-shock B1 in the prototypical chemically active protostellar outflow L1157. We made use of (sub-)mm data obtained in the framework of the Chemical HErschel Surveys of Star forming regions (CHESS) and Astrochemical Surveys at IRAM (ASAI) key science programs. We detected $^{12}$C$^{32}$S, $^{12}$C$^{34}$S, $^{13}$C$^{32}$S, and $^{12}$C$^{33}$S emissions, for a total of 18 transitions, with $E_{\rm u}$ up to $\sim$ 180 K. The unprecedented sensitivity of the survey allows us to carefully analyse the line profiles, revealing high-velocity emission, up to 20 km s$^{-1}$ with respect to the systemic. The profiles can be well fitted by a combination of two exponential laws that are remarkably similar to what previously found using CO. These components have been related to the cavity walls produced by the $\sim$ 2000 yr B1 shock and the older ($\sim$ 4000 yr) B2 shock, respectively. The combination of low- and high-excitation CS emission was used to properly sample the different physical components expected in a shocked region. Our CS observations show that this molecule is highlighting the dense, $n_{\rm H_2}$ = 1--5 $\times$ 10$^{5}$ cm$^{-3}$, cavity walls produced by the episodic outflow in L1157. In addition, the highest excitation (E$_u$ $\geq$ 130 K) CS lines provide us with the signature of denser (1--5 $\times$ 10$^{6}$ cm$^{-3}$) gas, associated with a molecular reformation zone of a dissociative J-type shock, which is expected to arise where the precessing jet impacting the molecular cavities. The CS fractional abundance increases up to $\sim$ 10$^{-7}$ in all the kinematical components. This value is consistent with what previously found for prototypical protostars and it is in agreement with the prediction of the abundances obtained via the chemical code Astrochem. △ Less

Submitted 30 October, 2014; originally announced October 2014.

Comments: 11 pages, 10 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

Deuteration and evolution in the massive star formation process: the role of surface chemistry

Authors: F. Fontani, G. Busquet, Aina Palau, P. Caselli, A., Sanchez-Monge, J. C. Tan, M. Audard

Abstract: An ever growing number of observational and theoretical evidence suggests that the deuterated fraction (column density ratio between a species containing D and its hydrogenated counterpart, Dfrac) is an evolutionary indicator both in the low- and the high-mass star formation process. However, the role of surface chemistry in these studies has not been quantified from an observational point of view… ▽ More An ever growing number of observational and theoretical evidence suggests that the deuterated fraction (column density ratio between a species containing D and its hydrogenated counterpart, Dfrac) is an evolutionary indicator both in the low- and the high-mass star formation process. However, the role of surface chemistry in these studies has not been quantified from an observational point of view. In order to compare how the deuterated fractions of species formed only in the gas and partially or uniquely on grain surfaces evolve with time, we observed rotational transitions of CH3OH, 13CH3OH, CH2DOH, CH3OD at 3 and 1.3~mm, and of NH2D at 3~mm with the IRAM-30m telescope, and the inversion transitions (1,1) and (2,2) of NH3 with the GBT, towards most of the cores already observed by Fontani et al.~(2011, 2014) in N2H+, N2D+, HNC, DNC. NH2D is detected in all but two cores, regardless of the evolutionary stage. Dfrac(NH3) is on average above 0.1, and does not change significantly from the earliest to the most evolved phases, although the highest average value is found in the protostellar phase (~0.3). Few lines of CH2DOH and CH3OD are clearly detected, and only towards protostellar cores or externally heated starless cores. This work clearly confirms an expected different evolutionary trend of the species formed exclusively in the gas (N2D+ and N2H+) and those formed partially (NH2D and NH3) or totally (CH2DOH and CH3OH) on grain mantles. The study also reinforces the idea that Dfrac(N2H+) is the best tracer of massive starless cores, while high values of Dfrac(CH3OH) seem rather good tracers of the early protostellar phases, at which the evaporation/sputtering of the grain mantles is most efficient. △ Less

Submitted 27 October, 2014; originally announced October 2014.

Comments: 14 pages, 7 figures, 2 Appendices

The identification of filaments on far infrared and submillimiter images. Morphology, physical conditions and relation with star formation of filamentary structure

Authors: E. Schisano, K. L. Rygl, S. Molinari, G. Busquet, D. Elia, M. Pestalozzi, D. Polychroni, N. Billot, A. Noriega-Crespo, S. Carey, R. Paladini, T. J. T. Moore, R. Plume, S. C. O. Glover, E. Vazquez-Semadeni

Abstract: Observations of molecular clouds reveal a complex structure, with gas and dust often arranged in filamentary rather than spherical geometries. The associations of pre- and proto- stellar cores with the filaments suggest a direct link with the process of star formation. Any study of the properties of such filaments requires a representative samples from different enviroments and so an unbiased dete… ▽ More Observations of molecular clouds reveal a complex structure, with gas and dust often arranged in filamentary rather than spherical geometries. The associations of pre- and proto- stellar cores with the filaments suggest a direct link with the process of star formation. Any study of the properties of such filaments requires a representative samples from different enviroments and so an unbiased detection method. We developed such an approach using the Hessian matrix of a surface-brightness distribution to identify filaments and determine their physical and morphological properties. After testing the method on simulated, but realistic filaments, we apply the algorithms to column-density maps computed from Herschel observations of the Galactic Plane obtained by the Hi-GAL project. We identified ~500 filaments, in the longitude range of l=216.5 to l=225.5, with lengths from ~1 pc up to ~30 pc and widths between 0.1 pc and 2.5 pc. Average column densities are between 10^20 cm^-2 and 10^22 cm^-2. Filaments include the majority of dense material with N_H_2 > 6x10^21cm^-2. We find that the pre- and proto-stellar compact sources already identified in the same region are mostly associated with filaments. However, surface densities in excess of the expected critical values for high-mass star formation are only found on the filaments, indicating that these structures are necessary to channel material into the clumps. Furthermore, we analyze the gravitational stability of filaments and discuss their relationship with star formation. △ Less

Submitted 17 June, 2014; originally announced June 2014.

Comments: 32 pages, 23 Figures. Accepted for publication in ApJ

Fragmentation of massive dense cores down to ~1000 AU: Relation between fragmentation and density structure

Authors: Aina Palau, R. Estalella, J. M. Girart, A. Fuente, F. Fontani, B. Commercon, G. Busquet, S. Bontemps, A. Sanchez-Monge, L. A. Zapata, Q. Zhang, P. Hennebelle, J. Di Francesco

Abstract: In order to shed light on the main physical processes controlling fragmentation of massive dense cores, we present a uniform study of the density structure of 19 massive dense cores, selected to be at similar evolutionary stages, for which their relative fragmentation level was assessed in a previous work. We inferred the density structure of the 19 cores through a simultaneous fit of the radial i… ▽ More In order to shed light on the main physical processes controlling fragmentation of massive dense cores, we present a uniform study of the density structure of 19 massive dense cores, selected to be at similar evolutionary stages, for which their relative fragmentation level was assessed in a previous work. We inferred the density structure of the 19 cores through a simultaneous fit of the radial intensity profiles at 450 and 850 micron (or 1.2 mm in two cases) and the Spectral Energy Distribution, assuming spherical symmetry and that the density and temperature of the cores decrease with radius following power-laws. We find a weak (inverse) trend of fragmentation level and density power-law index, with steeper density profiles tending to show lower fragmentation, and vice versa. In addition, we find a trend of fragmentation increasing with density within a given radius, which arises from a combination of flat density profile and high central density and is consistent with Jeans fragmentation. We considered the effects of rotational-to-gravitational energy ratio, non-thermal velocity dispersion, and turbulence mode on the density structure of the cores, and found that compressive turbulence seems to yield higher central densities. Finally, a possible explanation for the origin of cores with concentrated density profiles, which are the cores showing no fragmentation, could be related with a strong magnetic field, consistent with the outcome of radiation magnetohydrodynamic simulations. △ Less

Submitted 31 January, 2014; originally announced January 2014.

Comments: accepted for publication in ApJ; a version with all (high resolution) figures is available at: http://www.ice.csic.es/files/palau/Palau14.pdf

The CHESS survey of the L1157-B1 bow-shock: high and low excitation water vapor

Authors: G. Busquet, B. Lefloch, M. Benedettini, C. Ceccarelli, C. Codella, S. Cabrit, B. Nisini, S. Viti, A. I. Gómez-Ruiz, A. Gusdorf, A. M. di Giorgio, L. Wiesenfeld

Abstract: Molecular outflows powered by young protostars strongly affect the kinematics and chemistry of the natal molecular cloud through strong shocks resulting in substantial modifications of the abundance of several species. As part of the "Chemical Herschel Surveys of Star forming regions" guaranteed time key program, we aim at investigating the physical and chemical conditions of H20 in the brightest… ▽ More Molecular outflows powered by young protostars strongly affect the kinematics and chemistry of the natal molecular cloud through strong shocks resulting in substantial modifications of the abundance of several species. As part of the "Chemical Herschel Surveys of Star forming regions" guaranteed time key program, we aim at investigating the physical and chemical conditions of H20 in the brightest shock region B1 of the L1157 molecular outflow. We observed several ortho- and para-H2O transitions using HIFI and PACS instruments on board Herschel, providing a detailed picture of the kinematics and spatial distribution of the gas. We performed a LVG analysis to derive the physical conditions of H2O shocked material, and ultimately obtain its abundance. We detected 13 H2O lines probing a wide range of excitation conditions. PACS maps reveal that H2O traces weak and extended emission associated with the outflow identified also with HIFI in the o-H2O line at 556.9 GHz, and a compact (~10") bright, higher-excitation region. The LVG analysis of H2O lines in the bow-shock show the presence of two gas components with different excitation conditions: a warm (Tkin~200-300 K) and dense (n(H2)~(1-3)x10^6 cm-3) component with an assumed extent of 10" and a compact (~2"-5") and hot, tenuous (Tkin~900-1400 K, n(H2)~10^3-10^4 cm-3) gas component, which is needed to account for the line fluxes of high Eu transitions. The fractional abundance of the warm and hot H2O gas components is estimated to be (0.7-2)x10^{-6} and (1-3)x10^{-4}, respectively. Finally, we identified an additional component in absorption in the HIFI spectra of H2O lines connecting with the ground state level, probably arising from the photodesorption of icy mantles of a water-enriched layer at the edges of the cloud. △ Less

Submitted 12 November, 2013; originally announced November 2013.

Comments: Accepted for publication in A&A. 12 pages, 9 figures, 4 tables

Two Mass Distributions in the L 1641 Molecular Clouds: The Herschel connection of Dense Cores and Filaments in Orion A

Authors: D. Polychroni, E. Schisano, D. Elia, A. Roy, S. Molinari, P. Martin, Ph. Andre, D. Turrini, K. L. J. Rygl, M. Benedettini, G. Busquet, A. M. di Giorgio, M. Pestalozzi, S. Pezzuto, D. Arzoumanian, S. Bontemps, J. Di Francesco, M. Hennemann, T. Hill, V. Konyves, A. Menshchikov, F. Motte, Q. Nguyen-Luong, N. Peretto, N. Schneider , et al. (1 additional authors not shown)

Abstract: We present the Herschel Gould Belt survey maps of the L1641 molecular clouds in Orion A. We extracted both the filaments and dense cores in the region. We identified which of dense sources are proto- or pre-stellar, and studied their association with the identified filaments. We find that although most (71%) of the pre-stellar sources are located on filaments there is still a significant fraction… ▽ More We present the Herschel Gould Belt survey maps of the L1641 molecular clouds in Orion A. We extracted both the filaments and dense cores in the region. We identified which of dense sources are proto- or pre-stellar, and studied their association with the identified filaments. We find that although most (71%) of the pre-stellar sources are located on filaments there is still a significant fraction of sources not associated with such structures. We find that these two populations (on and off the identified filaments) have distinctly different mass distributions. The mass distribution of the sources on the filaments is found to peak at 4 Solar masses and drives the shape of the CMF at higher masses, which we fit with a power law of the form dN/dlogM \propto M^{-1.4+/-0.4}. The mass distribution of the sources off the filaments, on the other hand, peaks at 0.8 Solar masses and leads to a flattening of the CMF at masses lower than ~4 Solar masses. We postulate that this difference between the mass distributions is due to the higher proportion of gas that is available in the filaments, rather than in the diffuse cloud. △ Less

Submitted 5 October, 2013; v1 submitted 9 September, 2013; originally announced September 2013.

Comments: Accepted for publication in the Astrophysical Journal Letters; 9 pages, 4 figures, 1 table

The B1 shock in the L1157 outflow as seen at high spatial resolution

Authors: M. Benedettini, S. Viti, C. Codella, F. Gueth, A. I. Gomez-Ruiz, R. Bachiller, M. T. Beltran, G. Busquet, C. Ceccarelli, B. Lefloch

Abstract: We present high spatial resolution (750 AU at 250 pc) maps of the B1 shock in the blue lobe of the L1157 outflow in four lines: CS (3-2), CH3OH (3_K-2_K), HC3N (16-15) and p-H2CO (2_02-3_01). The combined analysis of the morphology and spectral profiles has shown that the highest velocity gas is confined in a few compact (~ 5 arcsec) bullets while the lowest velocity gas traces the wall of the gas… ▽ More We present high spatial resolution (750 AU at 250 pc) maps of the B1 shock in the blue lobe of the L1157 outflow in four lines: CS (3-2), CH3OH (3_K-2_K), HC3N (16-15) and p-H2CO (2_02-3_01). The combined analysis of the morphology and spectral profiles has shown that the highest velocity gas is confined in a few compact (~ 5 arcsec) bullets while the lowest velocity gas traces the wall of the gas cavity excavated by the shock expansion. A large velocity gradient model applied to the CS (3-2) and (2-1) lines provides an upper limit of 10^6 cm^-3 to the averaged gas density in B1 and a range of 5x10^3< n(H2)< 5x10^5 cm^-3 for the density of the high velocity bullets. The origin of the bullets is still uncertain: they could be the result of local instabilities produced by the interaction of the jet with the ambient medium or could be clump already present in the ambient medium that are excited and accelerated by the expanding outflow. The column densities of the observed species can be reproduced qualitatively by the presence in B1 of a C-type shock and only models where the gas reaches temperatures of at least 4000 K can reproduce the observed HC3N column density. △ Less

Submitted 2 September, 2013; originally announced September 2013.

Comments: 13 pages, 12 figures

Broad N2H+ emission towards the protostellar shock L1157-B1

Authors: C. Codella, S. Viti, C. Ceccarelli, B. Lefloch, M. Benedettini, G. Busquet, P. Caselli, F. Fontani, A. Gómez-Ruiz, L. Podio, M. Vasta

Abstract: We present the first detection of N2H+ towards a low-mass protostellar outflow, namely the L1157-B1 shock, at about 0.1 pc from the protostellar cocoon. The detection was obtained with the IRAM 30-m antenna. We observed emission at 93 GHz due to the J = 1-0 hyperfine lines. The analysis of the emission coupled with the HIFI CHESS multiline CO observations leads to the conclusion that the observed… ▽ More We present the first detection of N2H+ towards a low-mass protostellar outflow, namely the L1157-B1 shock, at about 0.1 pc from the protostellar cocoon. The detection was obtained with the IRAM 30-m antenna. We observed emission at 93 GHz due to the J = 1-0 hyperfine lines. The analysis of the emission coupled with the HIFI CHESS multiline CO observations leads to the conclusion that the observed N2H+(1-0) line originates from the dense (> 10^5 cm-3) gas associated with the large (20-25 arcsec) cavities opened by the protostellar wind. We find a N2H+ column density of few 10^12 cm-2 corresponding to an abundance of (2-8) 10^-9. The N2H+ abundance can be matched by a model of quiescent gas evolved for more than 10^4 yr, i.e. for more than the shock kinematical age (about 2000 yr). Modelling of C-shocks confirms that the abundance of N2H+ is not increased by the passage of the shock. In summary, N2H+ is a fossil record of the pre-shock gas, formed when the density of the gas was around 10^4 cm-3, and then further compressed and accelerated by the shock. △ Less

Submitted 29 August, 2013; originally announced August 2013.

Comments: ApJ, in press