STATCONT is a python-based tool designed to determine the continuum emission level in spectral data, in particular for sources with a line-rich spectrum. The tool inspects the intensity distribution ...of a given spectrum and automatically determines the continuum level by using different statistical approaches. The different methods included in STATCONT are tested against synthetic data. We conclude that the sigma-clipping algorithm provides the most accurate continuum level determination, together with information on the uncertainty in its determination. This uncertainty can be used to correct the final continuum emission level, resulting in the here called ‘corrected sigma-clipping method’ or c-SCM. The c-SCM has been tested against more than 750 different synthetic spectra reproducing typical conditions found towards astronomical sources. The continuum level is determined with a discrepancy of less than 1% in 50% of the cases, and less than 5% in 90% of the cases, provided at least 10% of the channels are line free. The main products of STATCONT are the continuum emission level, together with a conservative value of its uncertainty, and datacubes containing only spectral line emission, i.e., continuum-subtracted datacubes. STATCONT also includes the option to estimate the spectral index, when different files covering different frequency ranges are provided.
We present synthetic ALMA observations of Keplerian, protostellar discs in the Class 0 stage studying the emission of molecular tracers like 13CO, C18O, HCO+, H13CO+, N2H+, and H2CO. We model the ...emission of discs around low- and intermediate-mass protostars. We show that under optimal observing conditions ALMA is able to detect the discs already in the earliest stage of protostellar evolution, although the emission is often concentrated to the innermost 50 au. Therefore, a resolution of a few 0.1 arcsec might be too low to detect Keplerian discs around Class 0 objects. We also demonstrate that under optimal conditions for edge-on discs Keplerian rotation signatures are recognisable, from which protostellar masses can be inferred. For this we here introduce a new approach, which allows us to determine protostellar masses with higher fidelity than before. Furthermore, we show that it is possible to reveal Keplerian rotation even for strongly inclined discs and that ALMA should be able to detect possible signs of fragmentation in face-on discs. In order to give some guidance for future ALMA observations, we investigate the influence of varying observing conditions and source distances. We show that it is possible to probe Keplerian rotation in inclined discs with an observing time of 2 h and a resolution of 0.1 arcsec, even in the case of moderate weather conditions. Furthermore, we demonstrate that under optimal conditions, Keplerian discs around intermediate-mass protostars should be detectable up to kpc distances.
We have studied four complex organic molecules (COMs), the oxygen-bearing methyl formate (CH
3
OCHO) and dimethyl ether (CH
3
OCH
3
) as well as the nitrogen-bearing formamide (NH
2
CHO) and ethyl ...cyanide (C
2
H
5
CN), towards a large sample of 39 high-mass star-forming regions representing different evolutionary stages, from early to evolved phases. We aim to identify potential correlations and chemical links between the molecules and to trace their evolutionary sequence through the star formation process. We analysed spectra obtained at 3, 2, and 0.9 mm with the IRAM-30m telescope. We derived the main physical parameters for each species by fitting the molecular lines. We compared them and evaluated their evolution while also taking several other interstellar environments into account. We report detections in 20 sources, revealing a clear dust absorption effect on column densities. Derived abundances range between ~ 10
−10
−10
−7
for CH
3
OCHO and CH
3
OCH
3
, ~ 10
−12
−10
−10
for NH
2
CHO, and ~ 10
−11
−10
−9
for C
2
H
5
CN. The abundances of CH
3
OCHO, CH
3
OCH
3
, and C
2
H
5
CN are very strongly correlated (
r
≥ 0.92) across ~ 4 orders of magnitude. We note that CH
3
OCHO and CH
3
OCH
3
show the strongest correlations in most parameters, and a nearly constant ratio (~ 1) over a remarkable ~ 9 orders of magnitude in luminosity for the following wide variety of sources: pre-stellar to evolved cores, low- to high-mass objects, shocks, Galactic clouds, and comets. This indicates that COMs chemistry is likely early developed and then preserved through evolved phases. Moreover, the molecular abundances clearly increase with evolution, covering ~ 6 orders of magnitude in the luminosity/mass ratio. We consider CH
3
OCHO and CH
3
OCH
3
to be most likely chemically linked. They could, for example, share a common precursor, or be formed one from the other. Based on correlations, ratios, and the evolutionary trend, we propose a general scenario for all COMs, involving a formation in the cold, earliest phases of star formation and a following increasing desorption with the progressive thermal and shock-induced heating of the evolving core.
Context. Submillimeter Array (SMA) 870 μm polarization observations of the hot molecular core G31.41+0.31 revealed one of the clearest examples up to date of an hourglass-shaped magnetic field ...morphology in a high-mass star-forming region. Aims. To better establish the role that the magnetic field plays in the collapse of G31.41+0.31, we carried out Atacama Large Millimeter/ submillimeter Array (ALMA) observations of the polarized dust continuum emission at 1.3 mm with an angular resolution four times higher than that of the previous (sub)millimeter observations to achieve an unprecedented image of the magnetic field morphology. Methods. We used ALMA to perform full polarization observations at 233 GHz (Band 6). The resulting synthesized beam is 0′′.28×0′′.20 which, at the distance of the source, corresponds to a spatial resolution of ~875 au. Results. The observations resolve the structure of the magnetic field in G31.41+0.31 and allow us to study the field in detail. The polarized emission in the Main core of G31.41+0.41is successfully fit with a semi-analytical magnetostatic model of a toroid supported by magnetic fields. The best fit model suggests that the magnetic field is well represented by a poloidal field with a possible contribution of a toroidal component of ~10% of the poloidal component, oriented southeast to northwest at approximately −44° and with an inclination of approximately −45°. The magnetic field is oriented perpendicular to the northeast to southwest velocity gradient detected in this core on scales from 103 to 104 au. This supports the hypothesis that the velocity gradient is due to rotation of the core and suggests that such a rotation has little effect on the magnetic field. The strength of the magnetic field estimated in the central region of the core with the Davis–Chandrasekhar-Fermi method is ~8–13 mG and implies that the mass-to-flux ratio in this region is slightly supercritical. Conclusions. The magnetic field in G31.41+0.31 maintains an hourglass-shaped morphology down to scales of <1000 au. Despite the magnetic field being important in G31.41+0.31, it is not enough to prevent fragmentation and collapse of the core, as demonstrated by the presence of at least four sources embedded in the center of the core.
Aims. We describe the assignment of a previously unidentified interstellar absorption line to ArH+ and discuss its relevance in the context of hydride absorption in diffuse gas with a low H2 ...fraction. The confidence of the assignment to ArH+ is discussed, and the column densities are determined toward several lines of sight. The results are then discussed in the framework of chemical models, with the aim of explaining the observed column densities. Methods. We fitted the spectral lines with multiple velocity components, and determined column densities from the line-to-continuum ratio. The column densities of ArH+ were compared to those of other species, tracing interstellar medium (ISM) components with different H2 abundances. We constructed chemical models that take UV radiation and cosmic ray ionization into account. Results. Thanks to the detection of two isotopologues, 36ArH+ and 38ArH+, we are confident about the carrier assignment to ArH+. NeH+ is not detected with a limit of NeH+/ArH+ ≤ 0.1. The derived column densities agree well with the predictions of chemical models. ArH+ is a unique tracer of gas with a fractional H2 abundance of 10-4 − 10-3 and shows little correlation to H2O+, which traces gas with a fractional H2 abundance of ≈0.1. Conclusions. A careful analysis of variations in the ArH+, OH+, H2O+, and HF column densities promises to be a faithful tracer of the distribution of the H2 fractional abundance by providing unique information on a poorly known phase in the cycle of interstellar matter and on its transition from atomic diffuse gas to dense molecular gas traced by CO emission. Abundances of these species put strong observational constraints upon magnetohydrodynamical (MHD)simulations of the interstellar medium, and potentially could evolve into a tool characterizing the ISM. Paradoxically, the ArH+ molecule is a better tracer of almost purely atomic hydrogen gas than Hi itself, since Hi can also be present in gas with a significant molecular content, but ArH+ singles out gas that is >99.9% atomic.
High-mass stars and star clusters commonly form within hub-filament systems. Monoceros R2 (hereafter Mon R2), at a distance of 830 pc, harbors one of the closest such systems, making it an excellent ...target for case studies.
We investigate the morphology, stability and dynamical properties of the Mon R2 hub-filament system.
We employ observations of the
CO and C
O 1→0 and 2→1 lines obtained with the IRAM-30m telescope. We also use H
column density maps derived from
dust emission observations.
We identified the filamentary network in Mon R2 with the DisPerSE algorithm and characterized the individual filaments as either main (converging into the hub) or secondary (converging to a main filament) filaments. The main filaments have line masses of 30-100
pc
and show signs of fragmentation, while the secondary filaments have line masses of 12-60
pc
and show fragmentation only sporadically. In the context of Ostriker's hydrostatic filament model, the main filaments are thermally supercritical. If non-thermal motions are included, most of them are trans-critical. Most of the secondary filaments are roughly transcritical regardless of whether non-thermal motions are included or not. From the morphology and kinematics of the main filaments, we estimate a mass accretion rate of 10
-10
yr
into the central hub. The secondary filaments accrete into the main filaments with a rate of 0.1-0.4×10
yr
. The main filaments extend into the central hub. Their velocity gradients increase towards the hub, suggesting acceleration of the gas.We estimate that with the observed infall velocity, the mass-doubling time of the hub is ~ 2:5 Myr, ten times larger than the free-fall time, suggesting a dynamically old region. These timescales are comparable with the chemical age of the Hii region. Inside the hub, the main filaments show a ring- or a spiral-like morphology that exhibits rotation and infall motions. One possible explanation for the morphology is that gas is falling into the central cluster following a spiral-like pattern.
Context. Circumstellar discs around massive stars could mediate the accretion onto the star from the infalling envelope, and could minimize the effects of radiation pressure. Despite such a crucial ...role, only a few convincing candidates have been provided for discs around deeply embedded O-type (proto)stars. Aims. In order to establish whether disc-mediated accretion is the formation mechanism for the most massive stars, we have searched for circumstellar, rotating discs around a limited sample of six luminous (>105L⊙) young stellar objects. These objects were selected on the basis of their IR and radio properties in order to maximize the likelihood of association with disc+jet systems. Methods. We used ALMA with ~0.̋2 resolution to observe a large number of molecular lines typical of hot molecular cores. In this paper we limit our analysis to two disc tracers (methyl cyanide, CH3CN, and its isotopologue, 13CH3CN), and an outflow tracer (silicon monoxide, SiO). Results. We reveal many cores, although their number depends dramatically on the target. We focus on the cores that present prominent molecular line emission. In six of these a velocity gradient is seen across the core,three of which show evidence of Keplerian-like rotation. The SiO data reveal clear but poorly collimated bipolar outflow signatures towards two objects only. This can be explained if real jets are rare (perhaps short-lived) in very massive objects and/or if stellar multiplicity significantly affects the outflow structure.For all cores with velocity gradients, the velocity field is analysed through position–velocity plots to establish whether the gas is undergoing rotation with νrot ∝ R− α, as expected for Keplerian-like discs. Conclusions. Our results suggest that in three objects we are observing rotation in circumstellar discs, with three more tentative cases, and one core where no evidence for rotation is found. In all cases but one, we find that the gas mass is less than the mass of any embedded O-type star, consistent with the (putative) discs undergoing Keplerian-like rotation. With the caveat of low number statistics, we conclude that the disc detection rate could be sensitive to the evolutionary stage of the young stellar object. In young, deeply embedded sources, the evidence for discs could be weak because of confusion with the surrounding envelope, while in the most evolved sources the molecular component of the disc could have already been dispersed. Only in those objects that are at an intermediate stage of the evolution would the molecular disc be sufficiently prominent and relatively less embedded to be detectable by mm/submm observations.
As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41 +0.31 with ALMA at 1.4 mm with an angular resolution of ~0.′′22 (~1700 ...au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is associated with two free-free continuum sources. The Main core looks featureless and homogeneous in dust continuum emission and does not present any hint of fragmentation. Each transition of CH3CN and CH3OCHO, both ground and vibrationally excited, as well as those of CH3CN isotopologues, shows a clear velocity gradient along the NE–SW direction, with velocity linearly increasing with distance from the center, consistent with solid-body rotation. However, when comparing the velocity field of transitions with different upper level energies, the rotation velocity increases with increasing energy of the transition, which suggests that the rotation speeds up toward the center. Spectral lines towardtoward the dust continuum peak show an inverse P-Cygni profile that supports the existence of infall in the core. The infall velocity increases with the energy of the transition suggesting that the infall is accelerating toward the center of the core, consistent with gravitational collapse. Despite the monolithic appearance of the Main core, the presence of red-shifted absorption, the existence of two embedded free-free sources at the center, and the rotational spin-up are consistent with an unstable core undergoing fragmentation with infall and differential rotation due to conservation of angular momentum. Therefore, the most likely explanation for the monolithic morphology is that the large opacity of the dust emission prevents the detection of any inhomogeneity in the core.
Context. Theory predicts, and observations confirm, that the column density ratio of a molecule containing D to its counterpart containing H can be used as an evolutionary tracer in the low-mass star ...formation process. Aims. Since it remains unclear if the high-mass star formation process is a scaled-up version of the low-mass one, we investigated whether the relation between deuteration and evolution can be applied to the high-mass regime. Methods. With the IRAM-30 m telescope, we observed rotational transitions of N2D+ and N2H+ and derived the deuterated fraction in 27 cores within massive star-forming regions understood to represent different evolutionary stages of the massive-star formation process. Results. The abundance of N2D+ is higher at the pre-stellar/cluster stage, then drops during the formation of the protostellar object(s) as in the low-mass regime, remaining relatively constant during the ultra-compact HII region phase. The objects with the highest fractional abundance of N2D+ are starless cores with properties very similar to typical pre-stellar cores of lower mass. The abundance of N2D+ is lower in objects with higher gas temperatures as in the low-mass case but does not seem to depend on gas turbulence. Conclusions. Our results indicate that the N2D+-to-N2H+ column density ratio can be used as an evolutionary indicator in both low- and high-mass star formation, and that the physical conditions influencing the abundance of deuterated species likely evolve similarly during the processes that lead to the formation of both low- and high-mass stars.
Context.
To better understand the formation of high-mass stars, it is fundamental to investigate how matter accretes onto young massive stars, how it is ejected, and how all this differs from the ...low-mass case. The massive protocluster G31.41+0.31 is the ideal target to study all these processes because observations at millimeter and centimeter wavelengths have resolved the emission of the Main core into at least four massive dust continuum sources, named A, B, C, and D, within 1″ or 0.018 pc, and have identified signatures of infall and several outflows associated with the core.
Aims.
We study the interplay between infall and outflow in G31.41+0.31, at a spatial resolution of a few 100 au, by investigating their properties and their possible impact on the core.
Methods.
We carried out molecular line observations of typical high-density tracers, such as CH
3
CN or H
2
CO, and shock and outflow tracers, such as SiO, with ALMA at 1.4 mm that achieved an angular resolution of ~0.′′09 (~340 au).
Results.
The observations reveal inverse P Cygni profiles in CH
3
CN and H
2
CO toward the four sources embedded in the Main core, suggesting that all of them are undergoing collapse. The infall rates estimated from the redshifted absorption are on the order of ~10
−2
M
⊙
yr
−1
. The individual infall rates imply that the accretion timescale of the Main core is an order of magnitude smaller than its rotation timescale. This confirms that rotating toroids such as the G31 Main core are non-equilibrium, transient collapsing structures that need to be constantly replenished with fresh material from a large-scale reservoir. For sources B, C, and D the infall could be accelerating inside the sources, while for source A the presence of a second emission component complicates the interpretation. The SiO observations have revealed the presence of at least six outflows in the G31.41+0.31 star-forming region, and suggest that each of the four sources embedded in the Main core drives a molecular outflow. The outflow rates are on the order of ~10
−5
–10
−4
M
⊙
yr
−1
, depending on the SiO abundance. The feedback of the outflows appears sufficient to sustain the turbulence in the core and to eventually disrupt the core or prevent its further collapse. The mass accretion rates onto the individual sources, estimated from the highest values of the outflow mass loss rates, are on the order of 10
−4
M
⊙
yr
−1
. The difference of about two orders of magnitude between the accretion rates and the infall rates suggests that the central stars and the jets and outflows might not be efficient in removing disk material, which could lead to gravitational instabilities, fragmentation of the disk, and episodic accretion events.
Conclusions.
Infall and outflows are simultaneously present in all four sources embedded in the high-mass star-forming core G31.41+0.31. This indicates that these sources are still actively accreting and have not reached their final mass yet, which confirms the youth of this massive core.