Circumstellar discs are the precursors of planetary systems and develop shortly after their host star has formed. In their early stages these discs are immersed in an environment rich in gas and neighbouring stars, which can be hostile for their survival. There are several environmental processes that can affect the evolution of circumstellar discs, and external photoevaporation is arguably one of the most important ones. Theoretical and observational evidence point to circumstellar discs losing mass quickly when in the vicinity of massive, bright stars. External photoevaporation can then constrain the time and material available to form planets. The stellar density of the region seems to influence the extent of the effects of external photoevaporation. In this work we perform simulations of star clusters with a range of stellar densities. Stars with masses $\mathrm{M}_* \leq 1.9$ \MSun are surrounded by a disc, and stars with masses $\mathrm{M}_* > 1.9$ \MSun are considered as emitting radiation. Our results show that external photoevaporation is efficient in depleting disc masses and that the degree of its effect is related to stellar density. Dense clusters have $\sim10\%$ surviving discs by 2 Myr, whereas sparse ones have $\sim50\%$. Surviving discs in sparser regions also span a larger range of masses. We compare our results to observations of the Lupus clouds, the Orion Nebula Cluster, the Orion Molecular Cloud-2, Taurus, and NGC 2024, and find that the trends observed between region density and disc masses are similar to those in our simulations.

Strong winds and ultraviolet (UV) radiation from O-type stars disrupt and ionize their molecular core birthplaces, sweeping up material into parsec-size shells. Owing to dissociation by starlight, the thinnest shells are expected to host low molecular abundances and therefore little star formation. Here, we expand previous maps taken with the IRAM 30m telescope and present square-degree 12CO and 13CO (J=2-1) maps of the wind-driven "Veil bubble'' that surrounds the Trapezium cluster and its natal Orion molecular core (OMC). Although widespread and extended CO emission is largely absent from the Veil, we show that several CO "globules'' exist and are embedded in the [CII]158um-bright shell that confines the bubble. This includes the first detection of quiescent CO at negative LSR velocities in Orion. Given the harsh UV irradiation conditions in this translucent material, the detection of CO globules is surprising. These globules are small (R=7,100 AU), not massive (M=0.3M_Sun), and are moderately dense: n_ H=4x10^4 cm^-3 (median values). They are confined by the external pressure of the shell, P_ext/k~10^7 cm^-3 K, and are likely magnetically supported. They are either transient objects formed by instabilities or have detached from pre-existing molecular structures, sculpted by the passing shock associated with the expanding shell and by UV radiation from the Trapezium. Some represent the first stages in the formation of small pillars, others of isolated small globules. Although their masses do not suggest they will form stars, one globule matches the position of a known YSO. The lack of extended CO in the "Veil shell'' demonstrates that feedback from massive stars expels, agitates, and reprocesses most of the disrupted molecular cloud gas, thereby limiting the star-formation rate in the region. The presence of globules is a result of this feedback.

Sulphur is one of the most abundant elements in the Universe. Surprisingly, sulphuretted molecules are not as abundant as expected in the interstellar medium, and the identity of the main sulphur reservoir is still an open question. Our goal is to investigate the H$_{2}$S chemistry in dark clouds, as this stable molecule is a potential sulphur reservoir. Using millimeter observations of CS, SO, H$_{2}$S, and their isotopologues, we determine the physical conditions and H$_{2}$S abundances along the cores TMC 1-C, TMC 1-CP, and Barnard 1b. The gas-grain model Nautilus is then used to model the sulphur chemistry and explore the impact of photo-desorption and chemical desorption on the H$_2$S abundance. Our model shows that chemical desorption is the main source of gas-phase H$_2$S in dark cores. The measured H$_{2}$S abundance can only be fitted if we assume that the chemical desorption rate decreases by more than a factor of 10 when $n_{\rm H}>2\times10^{4}$. This change in the desorption rate is consistent with the formation of thick H$_2$O and CO ice mantles on grain surfaces. The observed SO and H$_2$S abundances are in good agreement with our predictions adopting an undepleted value of the sulphur abundance. However, the CS abundance is overestimated by a factor of $5-10$. Along the three cores, atomic S is predicted to be the main sulphur reservoir. We conclude that the gaseous H$_2$S abundance is well reproduced, assuming undepleted sulphur abundance and chemical desorption as the main source of H$_2$S. The behavior of the observed H$_{2}$S abundance suggests a changing desorption efficiency, which would probe the snowline in these cores. Our model, however, overestimates the observed gas-phase CS abundance. Given the uncertainty in the sulphur chemistry, our data are consistent with a cosmic elemental S abundance with an uncertainty of a factor of 10.

Aims. To demonstrate that `INDICATE' is a powerful spatial analysis tool which when combined with kinematic data from Gaia DR2 can be used to robustly probe star formation history. Methods. We compared the dynamic & spatial distributions of young stellar objects (YSOs) at various evolutionary stages in NGC 2264 using Gaia DR2 proper motion data and INDICATE. Results. The dynamic & spatial behaviours of YSOs at different evolutionary stages are distinct. Dynamically, Class II YSOs predominately have non-random trajectories that are consistent with known substructures, whereas Class III YSOs have random trajectories with no clear expansion or contraction patterns. Spatially, there is a correlation between the evolutionary stage and source concentration: 69.4% of Class 0/I, 27.9% of Class II, and 7.7% of Class III objects are found to be clustered. The proportion of YSOs clustered with objects of the same class also follows this trend. Class 0/I objects are both found to be more tightly clustered with the general populous/objects of the same class than Class IIs and IIIs by a factor of 1.2/4.1 and 1.9/6.6, respectively. An exception to these findings is within 0.05deg of S Mon where Class III objects mimic the behaviours of Class II sources across the wider cluster region. Our results suggest (i) current YSOs distributions are a result of dynamical evolution, (ii) prolonged star formation has been occurring sequentially, and (iii) stellar feedback from S Mon is causing YSOs to appear as more evolved sources. Conclusions. Designed to provide a quantitative measure of clustering behaviours, INDICATE is a powerful tool with which to perform rigorous spatial analyses. Our findings are consistent with what is known about NGC 2264, effectively demonstrating that when combined with kinematic data from Gaia DR2 INDICATE can be used to robustly study the star formation history of a cluster.

Context. The gas kinetic temperature (TK) determines the physical and chemical evolution of the Interestellar Medium (ISM). However, obtaining reliable TK estimates usually requires expensive observations including the combination of multi-line analysis and dedicated radiative transfer calculations. Aims. This work explores the use of HCN and HNC observations, and particularly its I(HCN)/I(HNC) intensity ratio of their J=1-0 lines, as direct probe of the gas kinetic temperature in the molecular ISM. Methods. We obtained a new set of large-scale observations of both HCN and HNC (1-0) lines along the Integral Shape Filament (ISF) in Orion. In combination with ancillary gas and dust temperature measurements, we find a systematic temperature dependence of the observed I(HCN)/I(HNC) intensity ratio across our maps. Additional comparisons with chemical models demonstrate that these observed I(HCN)/I(HNC) variations are driven by the effective destruction and isomerization mechanisms of HNC under low energy barriers. Results. The observed variations of I(HCN)/I(HNC) with TK can be described with a two-part linear function. This empirical calibration is then used to create a temperature map of the entire ISF. Comparisons with similar dust temperature measurements in this cloud, as well as in other regions and galactic surveys, validate this simple technique to obtain direct estimates of the gas kinetic temperature in a wide range of physical conditions and scales with an optimal working range between 15 K and 40 K. Conclusions. Both observations and models demonstrate the strong sensitivity of the I(HCN)/I(HNC) ratio to the gas kinetic temperature. Since these lines are easily obtained in observations of local and extragalactic sources, our results highlight the potential use of this observable as new chemical thermometer for the ISM.

Protoplanetary discs (PPDs) in the Orion Nebula Cluster (ONC) are irradiated by UV fields from the massive star $\theta^1$C. This drives thermal winds, inducing mass loss rates of up to $\dot{M}_\mathrm{wind}\sim 10^{-7}\,M_\odot$/yr in the `proplyds' (ionised PPDs) close to the centre. For the mean age of the ONC and reasonable initial PPD masses, such mass loss rates imply that discs should have been dispersed. However, ~80% of stars still exhibit a NIR excess, suggesting that significant circumstellar mass remains. This `proplyd lifetime problem' has persisted since the discovery of photoevaporating discs in the core of the ONC by O'Dell & Wen (1994). In this work, we demonstrate how an extended period of star formation can solve this problem. Coupling N-body calculations and a viscous disc evolution model, we obtain high disc fractions at the present day. This is partly due to the migration of older stars outwards, and younger stars inwards such that the most strongly irradiated PPDs are also the youngest. We show how the disc mass distribution can be used to test the recent claims in the literature for multiple stellar populations in the ONC. Our model also explains the recent finding that host mass and PPD mass are only weakly correlated, in contrast with other regions of similar age. We conclude that the status of the ONC as the archetype for understanding the influence of environment on planet formation is undeserved; the complex star formation history (involving star formation episodes within ~0.8 Myr of the present day) results in confusing signatures in the PPD population.

The mass evolution of protoplanetary disks is driven by both internal processes and external factors, such as photoevaporation. Disentangling these two effects, however, has remained difficult. We measure the dust masses of a sample of 132 disks in the Orion Molecular Cloud (OMC)-2 region, and compare them to (i) externally photoevaporated disks in the Trapezium cluster, and (ii) disks in nearby low-mass star forming regions (SFRs). This allows us to test if initial disk properties are the same in high- and low-mass SFRs, and enables a direct measurement of the effect of external photoevaporation on disks. A ~ $20' \times 4'$ mosaic of 3 mm continuum observations from the Atacama Large Millimeter/submillimeter Array (ALMA) was used to measure the fluxes of 132 disks and 35 protostars >0.5 pc away from the Trapezium. We identify and characterize a sample of 34 point sources not included in the Spitzer catalog on which the sample is based. Of the disks, 37 (28%) are detected, with masses ranging from 7-270 M_e. The detection rate for protostars is higher at 69%. Disks near the Trapezium are found to be less massive by a factor $0.18^{+0.18}_{-0.11}$, implying a mass loss rate of $8\times10^{-8}$ M_sun/yr. Our observations allow us to distinguish the impact of time and environment on disk evolution in a single SFR. The disk mass distribution in OMC-2 is statistically indistinguishable from that in nearby low-mass SFRs, like Lupus and Taurus. We conclude that age is the main factor determining the evolution of these disks. This result is robust with respect to assumptions of dust temperature, sample incompleteness and biases. The difference between the OMC-2 and Trapezium cluster samples is consistent with mass loss driven by far-ultraviolet radiation near the Trapezium. Together, this implies that in isolation, disk formation and evolution proceed similarly, regardless of cloud mass.

The past two decades have seen extensive surveys of the far-infrared to submillimeter continuum emission in the plane of our Galaxy. We line out prospects for the coming decade for corresponding molecular and atomic line surveys which are needed to fully understand the formation of the dense structures that give birth to clusters and stars out of the diffuse interstellar medium. We propose to work towards Galaxy wide surveys in mid-J CO lines to trace shocks from colliding clouds, Galaxy-wide surveys for atomic Carbon lines in order to get a detailed understanding of the relation of atomic and molecular gas in clouds, and to perform extensive surveys of the structure of the dense parts of molecular clouds to understand the importance of filaments/fibers over the full range of Galactic environments and to study how dense cloud cores are formed from the filaments. This work will require a large (50m) Single Dish submillimeter telescope equipped with massively multipixel spectrometer arrays, such as envisaged by the AtLAST project.

GEMS is an IRAM 30m Large Program whose aim is determining the elemental depletions and the ionization fraction in a set of prototypical star-forming regions. This paper presents the first results from the prototypical dark cloud TMC 1. Extensive millimeter observations have been carried out with the IRAM 30m telescope (3mm and 2mm) and the 40m Yebes telescope (1.3cm and 7mm) to determine the fractional abundances of CO, HCO+, HCN, CS, SO, HCS+, and N2H+ in three cuts which intersect the dense filament at the well-known positions TMC 1-CP, TMC 1-NH3, and TMC 1-C, covering a visual extinction range from Av~3 to ~20mag. Two phases with differentiated chemistry can be distinguished: i) the translucent envelope with molecular hydrogen densities of (1-5)x10$^3$ cm$^{-3}$; and ii) the dense phase, located at Av>10mag, with molecular hydrogen densities >10$^4$ cm$^{-3}$. Observations and modeling show that the gas phase abundances of C and O progressively decrease along the C+/C/CO transition zone where C/H~8x10$^{-5}$ and C/O~0.8-1, until the beginning of the dense phase at Av~10 mag. This is consistent with the grain temperatures being below the CO evaporation temperature in this region. In the case of sulfur, a strong depletion should occur before the translucent phase where we estimate a S/H~(0.4 - 2.2) x10$^{-6}$, an abundance ~7-40 times lower than the solar value. A second strong depletion must be present during the formation of the thick icy mantles to achieve the values of S/H measured in the dense cold cores (S/H~8x10$^{-8}$). Based on our chemical modeling, we constrain the value of $\zeta_{\rm H_2}$ to ~(0.5 - 1.8) x10$^{-16}$ s$^{-1}$ in the translucent cloud.

Recent observations of lensed sources have shown that the faintest ($M_{\mathrm{UV}} \approx -15\,\mathrm{mag}$) galaxies observed at z=6-8 appear to be extremely compact. Some of them have inferred sizes of less than 40 pc for stellar masses between $10^6$ and $10^7\,\mathrm{M}_{\odot}$, comparable to individual super star clusters or star cluster complexes at low redshift. High-redshift, low-mass galaxies are expected to show a clumpy, irregular morphology and if star clusters form in each of these well-separated clumps, the observed galaxy size would be much larger than the size of an individual star forming region. As supernova explosions impact the galaxy with a minimum delay time that exceeds the time required to form a massive star cluster, other processes are required to explain the absence of additional massive star forming regions. In this work we investigate whether the radiation of a young massive star cluster can suppress the formation of other detectable clusters within the same galaxy already before supernova feedback can affect the galaxy. We find that in low-mass ($M_{200} \lesssim 10^{10}\,\mathrm{M}_{\odot}$) haloes, the radiation from a compact star forming region with an initial mass of $10^{7}\,\mathrm{M}_{\odot}$ can keep gas clumps with Jeans masses larger than $\approx 10^{7}\,\mathrm{M}_{\odot}$ warm and ionized throughout the galaxy. In this picture, the small intrinsic sizes measured in the faintest $z=6-8$ galaxies are a natural consequence of the strong radiation field that stabilises massive gas clumps. A prediction of this mechanism is that the escape fraction for ionizing radiation is high for the extremely compact, high-z sources.

The binary T Tauri system JW 566 in the Orion Molecular Cloud underwent an energetic, short-lived flare observed at submillimetre wavelengths by the SCUBA-2 instrument on 26 November 2016 (UT). The emission faded by nearly 50% during the 31 minute integration. The simultaneous source fluxes averaged over the observation are 500 +/- 107 mJy/beam at 450 microns and 466 +/- 47 mJy/beam at 850 microns. The 850 micron flux corresponds to a radio luminosity of $L_{\nu}=8\times10^{19}$ erg/s/Hz, approximately one order of magnitude brighter (in terms of $\nu L_{\nu}$) than that of a flare of the young star GMR-A, detected in Orion in 2003 at 3mm. The event may be the most luminous known flare associated with a young stellar object and is also the first coronal flare discovered at sub-mm wavelengths. The spectral index between 450 microns and 850 microns of $\alpha = 0.11$ is broadly consistent with non-thermal emission. The brightness temperature was in excess of $6\times10^{4}$ K. We interpret this event to be a magnetic reconnection that energised charged particles to emit gyrosynchrotron/synchrotron radiation.

We use the $\mathit{Gaia}$ DR2 distances of about 700 mid-infrared selected young stellar objects in the benchmark giant molecular cloud Orion A to infer its 3D shape and orientation. We find that Orion A is not the fairly straight filamentary cloud that we see in (2D) projection, but instead a cometary-like cloud oriented toward the Galactic plane, with two distinct components: a denser and enhanced star-forming (bent) Head, and a lower density and star-formation quieter $\sim$75 pc long Tail. The true extent of Orion A is not the projected $\sim$40 pc but $\sim$90 pc, making it by far the largest molecular cloud in the local neighborhood. Its aspect ratio ($\sim$30:1) and high column-density fraction ($\sim45\%$) make it similar to large-scale Milky Way filaments ("bones"), despite its distance to the galactic mid-plane being an order of magnitude larger than typically found for these structures.

Multi-wavelength observations in the sub-mm regime provide information on the distribution of both the dust column density and the effective dust temperature in molecular clouds. In this study, we created high-resolution and high-dynamic-range maps of the Pipe nebula region and explored the value of dust-temperature measurements in particular towards the dense cores embedded in the cloud. The maps are based on data from the Herschel and Planck satellites, and calibrated with a near-infrared extinction map based on 2MASS observations. We have considered a sample of previously defined cores and found that the majority of core regions contain at least one local temperature minimum. Moreover, we observed an anti-correlation between column density and temperature. The slope of this anti-correlation is dependent on the region boundaries and can be used as a metric to distinguish dense from diffuse areas in the cloud if systematic effects are addressed appropriately. Employing dust-temperature data thus allows us to draw conclusions on the thermodynamically dominant processes in this sample of cores: external heating by the interstellar radiation field and shielding by the surrounding medium. In addition, we have taken a first step towards a physically motivated core definition by recognising that the column-density-temperature anti-correlation is sensitive to the core boundaries. Dust-temperature maps therefore clearly contain valuable information about the physical state of the observed medium.

We study the kinematics of the dense gas of starless and protostellar cores traced by the N2D+(2-1), N2H+(1-0), DCO+(2-1), and H13CO+(1-0) transitions along the L1495 filament and the kinematic links between the cores and the surrounding molecular cloud. We measure velocity dispersions, local and total velocity gradients and estimate the specific angular momenta of 13 dense cores in the four transitions using the on-the-fly observations with the IRAM 30 m antenna. To study a possible connection to the filament gas, we use the fit results of the C18O(1-0) survey performed by Hacar et al. (2013). All cores show similar properties along the 10 pc-long filament. N2D+(2-1) shows the most centrally concentrated structure, followed by N2H+(1-0) and DCO+(2-1), which show similar spatial extent, and H13CO+(1-0). The non-thermal contribution to the velocity dispersion increases from higher to lower density tracers. The change of magnitude and direction of the total velocity gradients depending on the tracer used indicates that internal motions change at different depths within the cloud. N2D+ and N2H+ show smaller gradients than the lower density tracers DCO+ and H13CO+, implying a loss of specific angular momentum at small scales. At the level of cloud-core transition, the core's external envelope traced by DCO+ and H13CO+ is spinning up, consistent with conservation of angular momentum during core contraction. C18O traces the more extended cloud material whose kinematics is not affected by the presence of dense cores. The decrease in specific angular momentum towards the centres of the cores shows the importance of local magnetic fields to the small scale dynamics of the cores. The random distributions of angles between the total velocity gradient and large scale magnetic field suggests that the magnetic fields may become important only in the high density gas within dense cores.

Abridged. Are all filaments bundles of fibers? To address this question, we have investigated the gas organization within the paradigmatic Integral Shape Filament (ISF). We combined two new ALMA Cycle 3 mosaics with previous IRAM 30m observations to produce a high-dynamic range N$_2$H$^+$(1-0) emission map of the ISF tracing its high-density material and velocity structure down to scales of 0.009 pc. From the analysis of the gas kinematics, we identify a total of 55 dense fibers in the central region of the ISF. Independently of their location, these fibers are characterized by transonic internal motions, lengths of ~0.15 pc, and masses per-unit-length close to those expected in hydrostatic equilibrium. The ISF fibers are spatially organized forming a dense bundle with multiple hub-like associations likely shaped by the local gravitational potential. Within this complex network, the ISF fibers show a compact radial emission profile with a median FWHM of 0.035 pc systematically narrower than the previously proposed universal 0.1 pc filament width. Our ALMA observations reveal complex bundles of fibers in the ISF, suggesting strong similarities between the internal substructure of this massive filament and previously studied lower-mass objects. The fibers show identical dynamic properties in both low- and high-mass regions, and their widespread detection suggests a preferred organizational mechanism of gas in which the physical fiber dimensions (width and length) are self-regulated depending on their intrinsic gas density. Combined with previous works, we identify a systematic increase of the surface density of fibers as a function of the total mass per-unit-length in filamentary clouds. Based on this empirical correlation, we propose a unified star-formation scenario where the observed differences between low- and high-mass clouds emerge naturally from the initial concentration of fibers.

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

We present a solution for the observed core fragmentation of filaments in the Taurus L1517 dark cloud which previously could not be explained (Hacar et. al 2011). Core fragmentation is a vital step for the formation of stars. Observations suggest a connection to the filamentary structure of the cloud gas but it remains unclear which process is responsible. We show that the gravitational instability process of an isothermal cylinder can account for the exhibited fragmentation under the assumption that the perturbation grows on the dominant wavelength. We use numerical simulations with the code RAMSES, estimate observed column densities and line-of-sight velocities and compare them to the observations. A critical factor for the observed fragmentation is that cores grow by redistributing mass within the filament and thus the density between the cores decreases over the fragmentation process. This often leads to wrong dominant wavelength estimates as it is strongly dependent on the initial central density. We argue that non-linear effects also play an important role on the evolution of the fragmentation. Once the density perturbation grows above the critical line-mass, non-linearity leads to an enhancement of the central core density in comparison to the analytical prediction. Choosing the correct initial conditions with perturbation strengths of around 20\%, leads to inclination corrected line-of-sight velocities and central core densities within the observational measurement error in a realistic evolution time.

We have investigated the global dynamical state of the Integral Shaped Filament in the Orion A cloud using new N$_2$H$^+$ (1-0) large-scale, IRAM30m observations. Our analysis of its internal gas dynamics reveals the presence of accelerated motions towards the Orion Nebula Cluster, showing a characteristic blue-shifted profile centred at the position of the OMC-1 South region. The properties of these observed gas motions (profile, extension, and magnitude) are consistent with the expected accelerations for the gravitational collapse of the OMC-1 region and explain both the physical and kinematic structure of this cloud.

Are the initial conditions for clustered star formation the same as for non-clustered star formation? To investigate the initial gas properties in young proto-clusters we carried out a comprehensive and high-sensitivity study of the internal structure, density, temperature, and kinematics of the dense gas content of the NGC1333 region in Perseus, one of the nearest and best studied embedded clusters. The analysis of the gas velocities in the Position-Position-Velocity space reveals an intricate underlying gas organization both in space and velocity. We identified a total of 14 velocity-coherent, (tran-)sonic structures within NGC1333, with similar physical and kinematic properties than those quiescent, star-forming (aka fertile) fibers previously identified in low-mass star-forming clouds. These fibers are arranged in a complex spatial network, build-up the observed total column density, and contain the dense cores and protostars in this cloud. Our results demonstrate that the presence of fibers is not restricted to low-mass clouds but can be extended to regions of increasing mass and complexity. We propose that the observational dichotomy between clustered and non-clustered star-forming regions might be naturally explained by the distinct spatial density of fertile fibers in these environments.

This paper continues our study of the foreground population to the Orion molecular clouds. The goal is to characterize the foreground population north of NGC 1981 and to investigate the star formation history in the large Orion star-forming region. We focus on a region covering about 25 square degrees, centered on the $\epsilon$ Orionis supergiant (HD 37128, B0\,Ia) and covering the Orion Belt asterism. We used a combination of optical (SDSS) and near-infrared (2MASS) data, informed by X-ray (\textit{XMM-Newton}) and mid-infrared (WISE) data, to construct a suite of color-color and color-magnitude diagrams for all available sources. We then applied a new statistical multiband technique to isolate a previously unknown stellar population in this region. We identify a rich and well-defined stellar population in the surveyed region that has about 2\,000 objects that are mostly M stars. We infer the age for this new population to be at least 5\, Myr and likely $\sim10$\,Myr and estimate a total of about 2\,500 members, assuming a normal IMF. This new population, which we call the Orion Belt population, is essentially extinction-free, disk-free, and its spatial distribution is roughly centered near $\epsilon$ Ori, although substructure is clearly present. The Orion Belt population is likely the low-mass counterpart to the Ori OB Ib subgroup. Although our results do not rule out Blaauw's sequential star formation scenario for Orion, we argue that the recently proposed blue streams scenario provides a better framework on which one can explain the Orion star formation region as a whole. We speculate that the Orion Belt population could represent the evolved counterpart of an Orion nebula-like cluster.

(Abridged) Many of the observed CO line profiles exhibit broad linewidths that greatly exceed the thermal broadening expected within molecular clouds. These suprathermal CO linewidths are assumed to be originated from the presence of unresolved supersonic motions inside clouds. Typically overlooked in the literature, in this paper we aim to quantify the impact of the opacity broadening effects on the current interpretation of the CO suprathermal line profiles. Without any additional contributions to the gas velocity field, a large fraction of the apparently supersonic (${\cal M}\sim$2-3) linewidths measured in both $^{12}$CO and $^{13}$CO (J=1-0) lines can be explained by the saturation of their corresponding sonic-like, optically-thin C$^{18}$O counterparts assuming standard isotopic fractionation. Combined with the presence of multiple components detected in our C$^{18}$O spectra, these opacity effects seem to be also responsible of the highly supersonic linewidths (${\cal M}>$8-10) detected in the broadest $^{12}$CO and $^{13}$CO spectra in Taurus. Our results demonstrate that most of the suprathermal $^{12}$CO and $^{13}$CO linewidths could be primarily created by a combination of opacity broadening effects and multiple gas velocity components blended in these saturated emission lines. Once corrected by their corresponding optical depth, each of these gas components present transonic intrinsic linewidths consistently traced by the three CO isotopologues within a factor of 2. Highly correlated and velocity-coherent at large scales, the largest and highly supersonic velocity differences inside clouds are generated by the relative motions between individual gas components. This highly discretized structure of the molecular gas traced in CO suggest that the gas dynamics inside molecular clouds could be better described by the properties of a fully-resolved macroscopic turbulence.

We compare APOGEE radial velocities (RVs) of young stars in the Orion A cloud with CO line gas emission and find a correlation between the two at large-scales, in agreement with previous studies. However, at smaller scales we find evidence for the presence of substructure in the stellar velocity field. Using a Friends-of-Friends approach we identify 37 stellar groups with almost identical RVs. These groups are not randomly distributed but form elongated chains or strings of stars with five or more members with low velocity dispersion, across lengths of 1-1.5~pc. The similarity between the kinematic properties of the APOGEE strings and the internal velocity field of the chains of dense cores and fibers recently identified in the dense ISM is striking and suggests that for most of the Orion A cloud, young stars keep memory of the parental gas substructure where they originated.

Orion A hosts the nearest massive star factory, thus offering a unique opportunity to resolve the processes connected with the formation of both low- and high-mass stars. Here we present the most detailed and sensitive near-infrared (NIR) observations of the entire molecular cloud to date. With the unique combination of high image quality, survey coverage, and sensitivity, our NIR survey of Orion A aims at establishing a solid empirical foundation for further studies of this important cloud. In this first paper we present the observations, data reduction, and source catalog generation. To demonstrate the data quality, we present a first application of our catalog to estimate the number of stars currently forming inside Orion A and to verify the existence of a more evolved young foreground population. We used the European Southern Observatory's (ESO) Visible and Infrared Survey Telescope for Astronomy (VISTA) to survey the entire Orion A molecular cloud in the NIR $J, H$, and $K_S$ bands, covering a total of $\sim$18.3 deg$^2$. We implemented all data reduction recipes independently of the ESO pipeline. Estimates of the young populations toward Orion A are derived via the $K_S$-band luminosity function. Our catalog (799995 sources) increases the source counts compared to the Two Micron All Sky Survey by about an order of magnitude. The 90% completeness limits are 20.4, 19.9, and 19.0 mag in $J, H$, and $K_S$, respectively. The reduced images have 20% better resolution on average compared to pipeline products. We find between 2300 and 3000 embedded objects in Orion A and confirm that there is an extended foreground population above the Galactic field, in agreement with previous work. The Orion A VISTA catalog represents the most detailed NIR view of the nearest massive star-forming region and provides a fundamental basis for future studies of star formation processes toward Orion.

Filaments play a central role in the molecular clouds' evolution, but their internal dynamical properties remain poorly characterized. To further explore the physical state of these structures, we have investigated the kinematic properties of the Musca cloud. We have sampled the main axis of this filamentary cloud in $^{13}$CO and C$^{18}$O (2--1) lines using APEX observations. The different line profiles in Musca shows that this cloud presents a continuous and quiescent velocity field along its $\sim$6.5 pc of length. With an internal gas kinematics dominated by thermal motions (i.e., $\sigma_{NT}/c_s\lesssim1$) and large-scale velocity gradients, these results reveal Musca as the longest velocity-coherent, sonic-like object identified so far in the ISM. The transonic properties of Musca present a clear departure from the predicted supersonic velocity dispersions expected in the Larson's velocity dispersion-size relationship, and constitute the first observational evidence of a filament fully decoupled from the turbulent regime over multi-parsec scales.

Aims: We resolve the length-scales for filament formation and fragmentation (res. <=0.1pc), in particular the Jeans length and cylinder fragmentation scale. Methods: We observed the prototypical high-mass star-forming filament IRDC18223 with the Plateau de Bure Interferometer (PdBI) in the 3.2mm continuum and N2H+(1-0) line emission in a ten field mosaic at a spatial resolution of ~4'' (~14000AU). Results: The dust continuum emission resolves the filament into a chain of at least 12 relatively regularly spaced cores. The mean separation between cores is ~0.40(+-0.18)pc. While this is approximately consistent with the fragmentation of an infinite, isothermal, gravitationally bound gas cylinder, a high mass-to-length ratio of M/l~1000M_sun/pc requires additional turbulent and/or magnetic support against radial collapse of the filament. The N2H+(1-0) data reveal a velocity gradient perpendicular to the main filament. Although rotation of the filament cannot be excluded, the data are also consistent with the main filament being comprised of several velocity-coherent sub-filaments. Furthermore, this velocity gradient perpendicular to the filament resembles recent results toward Serpens south that are interpreted as signatures of filament formation within magnetized and turbulent sheet-like structures. Lower-density gas tracers ([CI] and C18O) reveal a similar red/blueshifted velocity structure on scales around 60'' east and west of the IRDC18223 filament. This may tentatively be interpreted as a signature of the large-scale cloud and the smaller-scale filament being kinematically coupled. We do not identify a velocity gradient along the axis of the filament. This may either be due to no significant gas flows along the filamentary axis, but it may partly also be caused by a low inclination angle of the filament with respect to the plane of the sky that could minimize such signature.

We present the first 7.5'x11.5' velocity-resolved map of the [CII]158um line toward the Orion molecular cloud-1 (OMC-1) taken with the Herschel/HIFI instrument. In combination with far-infrared (FIR) photometric images and velocity-resolved maps of the H41alpha hydrogen recombination and CO J=2-1 lines, this data set provides an unprecedented view of the intricate small-scale kinematics of the ionized/PDR/molecular gas interfaces and of the radiative feedback from massive stars. The main contribution to the [CII] luminosity (~85%) is from the extended, FUV-illuminated face of the cloud G_0>500, n_H>5x10^3 cm^-3) and from dense PDRs (G_0~10^4, n_H~10^5 cm^-3) at the interface between OMC-1 and the HII region surrounding the Trapezium cluster. Around 15% of the [CII] emission arises from a different gas component without CO counterpart. The [CII] excitation, PDR gas turbulence, line opacity (from [13CII]) and role of the geometry of the illuminating stars with respect to the cloud are investigated. We construct maps of the [CII]/FIR and FIR/M_Gas ratios and show that [CII]/FIR decreases from the extended cloud component (10^-2-10^-3) to the more opaque star-forming cores (10^-3-10^-4). The lowest values are reminiscent of the "[CII] deficit" seen in local ultra-luminous IR galaxies hosting vigorous star formation. Spatial correlation analysis shows that the decreasing [CII]/FIR ratio correlates better with the column density of dust through the molecular cloud than with FIR/M_Gas. We conclude that the [CII] emitting column relative to the total dust column along each line of sight is responsible for the observed [CII]/FIR variations through the cloud.

Filamentary structures are common in molecular clouds. Explaining how they fragment to dense cores is a missing step in understanding their role in star formation. We perform a case study of whether low-mass filaments are close-to hydrostatic prior to their fragmentation, and whether their fragmentation agrees with gravitational fragmentation models. For this, we study the 6.5 pc long Musca molecular cloud that is an ideal candidate for a filament at an early stage of fragmentation. We employ dust extinction mapping in conjunction with near-infrared data from the NEWFIRM instrument, and 870 um dust continuum emission data from the LABOCA instrument, to estimate column densities. We use the data to identify fragments from the cloud and to determine the radial density distribution of its filamentary part. We compare the cloud's morphology with 13CO and C18O line emission observed with the APEX/SHeFI instrument. The Musca cloud is pronouncedly fragmented at its ends, but harbours a remarkably well-defined, 1.6 pc long filament in its Center region. The line mass of the filament is 21-31 Ms pc^-1 and FWHM 0.07 pc. Its radial profile can be fitted with a Plummer profile that has the power-index of 2.6 \pm 11%, flatter than that of an infinite hydrostatic filament. The profile can also be fitted with a hydrostatic cylinder truncated by external pressure. These models imply a central density of 5-10 x 10^4 cm^-3. The fragments in the cloud have a mean separation of 0.4 pc, in agreement with gravitational fragmentation. These properties, together with the subsonic and velocity-coherent nature of the cloud, suggest a scenario in which an initially hydrostatic cloud is currently gravitationally fragmenting. The fragmentation has started a few tenths of a Myr ago from the cloud ends, leaving its center yet relatively non-fragmented, possibly because of gravitational focusing in a finite geometry.

(Abridged) We study the kinematics of the dense gas in the Taurus L1495/B213 filamentary region to investigate the mechanism of core formation. We use observations of N2H+(1-0) and C18O(2-1) carried out with the IRAM 30m telescope. We find that the dense cores in L1495/B213 are significantly clustered in linear chain-like groups about 0.5pc long. The internal motions in these chains are mostly subsonic and the velocity is continuous, indicating that turbulence dissipation in the cloud has occurred at the scale of the chains and not at the smaller scale of the individual cores. The chains also present an approximately constant abundance of N2H+ and radial intensity profiles that can be modeled with a density law that follows a softened power law. A simple analysis of the spacing between the cores using an isothermal cylinder model indicates that the cores have likely formed by gravitational fragmentation of velocity-coherent filaments. Combining our analysis of the cores with our previous study of the large-scale C18O emission from the cloud, we propose a two-step scenario of core formation in L1495/B213. In this scenario, named "fray and fragment," L1495/B213 originated from the supersonic collision of two flows. The collision produced a network of intertwined subsonic filaments or fibers ("fray" step). Some of these fibers accumulated enough mass to become gravitationally unstable and fragment into chains of closely-spaced cores. This scenario may also apply to other regions of star formation.

The physical properties of the so-called Ostriker isothermal, non-rotating filament have been classically used as benchmark to interpret the stability of the filaments observed in nearby clouds. However, such static picture seems to contrast with the more dynamical state observed in different filaments. In order to explore the physical conditions of filaments under realistic conditions, in this work we theoretically investigate how the equilibrium structure of a filament changes in a rotating configuration. To do so, we solve the hydrostatic equilibrium equation assuming both uniform and differential rotations independently. We obtain a new set of equilibrium solutions for rotating and pressure truncated filaments. These new equilibrium solutions are found to present both radial and projected column density profiles shallower than their Ostriker-like counterparts. Moreover, and for rotational periods similar to those found in the observations, the centrifugal forces present in these filaments are also able to sustain large amounts of mass (larger than the mass attained by the Ostriker filament) without being necessary unstable. Our results indicate that further analysis on the physical state of star-forming filaments should take into account rotational effects as stabilizing agents against gravity

(abridged) [...] Methods: In a continued study of the molecular core population of the Pipe Nebula, we present a molecular-line survey of 52 cores. Previous research has shown a variety of different chemical evolutionary stages among the cores. Using the Mopra radio telescope, we observed the ground rotational transitions of HCO+, H13CO+, HCN, H13CN, HNC, and N2H+. These data are complemented with near-infrared extinction maps to constrain the column densities, effective dust temperatures derived from Herschel data, and NH3-based gas kinetic temperatures. Results: The target cores are located across the nebula, span visual extinctions between 5 and 67 mag, and effective dust temperatures (averaged along the lines of sight) between 13 and 19 K. The extinction-normalized integrated line intensities, a proxy for the abundance in constant excitation conditions of optically thin lines, vary within an order of magnitude for a given molecule. The effective dust temperatures and gas kinetic temperatures are correlated, but the effective dust temperatures are consistently higher than the gas kinetic temperatures. Combining the molecular line and temperature data, we find that N2H+ is only detected toward the coldest and densest cores while other lines show no correlation with these core properties. Conclusions: Within this large sample, N2H+ is the only species to exclusively trace the coldest and densest cores, in agreement with chemical considerations. In contrast, the common high-density tracers HCN and HNC are present in a majority of cores, demonstrating the utility of these molecules to characterize cores over a large range of extinctions. The correlation between the effective dust temperatures and the gas kinetic temperatures suggests that the former are dominated by dust that is both dense and thermodynamically coupled to the dense gas traced by NH3. [...]

The physical properties of the so-called Ostriker isothermal filament (Ostriker 1964) have been classically used as benchmark to interpret the stability of the filaments observed in nearby clouds. However, recent continuum studies have shown that the internal structure of the filaments depart from the isothermality, typically exhibiting radially increasing temperature gradients. The presence of internal temperature gradients within filaments suggests that the equilibrium configuration of these objects should be therefore revisited. The main goal of this work is to theoretically explore how the equilibrium structure of a filament changes in a non-isothermal configuration. We solve the hydrostatic equilibrium equation assuming temperature gradients similar to those derived from observations. We obtain a new set of equilibrium solutions for non-isothermal filaments with both linear and asymptotically constant temperature gradients. Our results show that, for sufficiently large internal temperature gradients, a non-isothermal filament could present significantly larger masses per unit length and shallower density profiles than the isothermal filament without collapsing by its own gravity. We conclude that filaments can reach an equilibrium configuration under non-isothermal conditions. Detailed studies of both the internal mass distribution and temperature gradients within filaments are then needed in order to judge the physical state of filaments.

Context. A small group of bipolar protostellar outflows display strong emission from shock-tracer molecules such as SiO and CH3OH, and are generally referred to as "chemically active." The best-studied outflow from this group is the one in L 1157. Aims. We study the molecular emission from the bipolar outflow powered by the very young stellar object HH 114 MMS and compare its chemical composition with that of the L1157 outflow. Methods. We have used the IRAM 30m radio telescope to observe a number of transitions from CO, SiO, CH3OH, SO, CS, HCN, and HCO+ toward the HH 114 MMS outflow. The observations consist of maps and a two-position molecular survey. Results. The HH 114 MMS outflow presents strong emission from a number of shock-tracer molecules that dominate the appearance of the maps around the central source. The abundance of these molecules is comparable to the abundance in L 1157. Conclusions. The outflow from HH 114 MMS is a spectacular new case of a chemically active outflow.

(Abridged) Context. Core condensation is a critical step in the star-formation process, but is still poorly characterized observationally. Aims. We have studied the 10 pc-long L1495/B213 complex in Taurus to investigate how dense cores have condensed out of the lower-density cloud material. Results. From the N$_2$H$^+$ emission, we identify 19 dense cores, some starless and some protostellar. They are not distributed uniformly, but tend to cluster with relative separations on the order of 0.25 pc. From the C$^{18}$O emission, we identify multiple velocity components in the gas. We have characterized them by fitting gaussians to the spectra, and by studying the distribution of the fits in position-position-velocity space. In this space, the C$^{18}$O components appear as velocity-coherent structures, and we have identified them automatically using a dedicated algorithm (FIVe: Friends In Velocity). Using this algorithm, we have identified 35 filamentary components with typical lengths of 0.5 pc, sonic internal velocity dispersions, and mass-per-unit-length close to the stability threshold of isothermal cylinders at 10 K. Core formation seems to have occurred inside the filamentary components via fragmentation, with a small number of fertile components with larger mass-per-unit-length being responsible for most cores in the cloud. At large scales, the filamentary components appear grouped into families, which we refer to as bundles. Conclusions. Core formation in L1495/B213 has proceeded by hierarchical fragmentation. The cloud fragmented first into several pc-scale regions. Each of these regions later fragmented into velocity-coherent filaments of about 0.5 pc in length. Finally, a small number of these filaments fragmented quasi-statically and produced the individual dense cores we see today.

Context. Low-mass star-forming cores differ from their surrounding molecular cloud in turbulence, shape, and density structure. Aims. We aim to understand how dense cores form out of the less dense cloud material by studying the connection between these two regimes. Methods. We observed the L1517 dark cloud in C18O(1-0), N2H+(1-0), and SO(JN=32-21) with the FCRAO 14m telescope, and in the 1.2mm dust continuum with the IRAM 30m telescope. Results. Most of the gas in the cloud lies in four filaments that have typical lengths of 0.5 pc. Five starless cores are embedded in these filaments and have chemical compositions indicative of different evolutionary stages. The filaments have radial profiles of C18O(1-0) emission with a central flattened region and a power-law tail, and can be fitted approximately as isothermal, pressure-supported cylinders. The filaments, in addition, are extremely quiescent. They have subsonic internal motions and are coherent in velocity over their whole length. The large-scale motions in the filaments can be used to predict the velocity inside the cores, indicating that core formation has not decoupled the dense gas kinematically from its parental material. In two filaments, these large-scale motions consist of oscillations in the velocity centroid, and a simple kinematic model suggests that they may be related to core-forming flows. Conclusions. Core formation in L1517 seems to have occurred in two steps. First, the subsonic, velocity-coherent filaments have condensed out of the more turbulent ambient cloud. Then, the cores fragmented quasi-statically and inherited the kinematics of the filaments. Turbulence dissipation has therefore occurred mostly on scales on the order of 0.5 pc or larger, and seems to have played a small role in the formation of the individual cores.

(Abridged) We present a molecular survey of the outflows powered by L1448-mm and IRAS 04166+2706, two sources with prominent wing and extremely high velocity (EHV) components in their CO spectra. The molecular composition of the two outflows presents systematic changes with velocity that we analyze by dividing the outflow in three chemical regimes, two of them associated with the wing component and the other the EHV gas. The analysis of the two wing regimes shows that species like H2CO and CH3OH favor the low-velocity gas, while SiO and HCN are more abundant in the fastest gas. We also find that the EHV regime is relatively rich in O-bearing species, as is not only detected in CO and SiO (already reported elsewhere), but also in SO, CH3OH, and H2CO (newly reported here), with a tentative detection in HCO+. At the same time, the EHV regime is relatively poor in C-bearing molecules like CS and HCN. We suggest that this difference in composition arises from a lower C/O ratio in the EHV gas. The different chemical compositions of the wing and EHV regimes suggest that these two outflow components have different physical origins. The wing component is better explained by shocked ambient gas, although none of the existing shock models explains all observed features. The peculiar composition of the EHV gas may reflect its origin as a dense wind from the protostar or its surrounding disk.