Bringing high spatial resolution to the far-infrared: A giant leap for astrophysics.

The far-infrared (FIR) regime is one of the wavelength ranges where no astronomical data with sub-arcsecond spatial resolution exist. None of the medium-term satellite projects like SPICA, Millimetron, or the Origins Space Telescope will resolve this malady. For many research areas, however, information at high spatial and spectral resolution in the FIR, taken from atomic fine-structure lines, from highly excited carbon monoxide (CO), light hydrides, and especially from water lines would open the door for transformative science. A main theme will be to trace the role of water in proto-planetary discs, to observationally advance our understanding of the planet formation process and, intimately related to that, the pathways to habitable planets and the emergence of life. Furthermore, key observations will zoom into the physics and chemistry of the star-formation process in our own Galaxy, as well as in external galaxies. The FIR provides unique tools to investigate in particular the energetics of heating, cooling, and shocks. The velocity-resolved data in these tracers will reveal the detailed dynamics engrained in these processes in a spatially resolved fashion, and will deliver the perfect synergy with ground-based molecular line data for the colder dense gas.

Multiplicity of Northern bright O-type stars with optical long baseline interferometry. Results of the pilot survey.

The study of the multiplicity of massive stars gives hints of their formation processes and their evolution path. Optical interferometry is mandatory to fulfill our knowledge of their multiplicity by probing the separation gap between 1 and 50 mas. We demonstrated the capability of the new interferometric instrument MIRC-X, located at the CHARA array, to study a large sample of more than 120 (H < 7.5) O-type stars. We observed 29 O-type star systems, including a couple of systems in average atmospheric conditions around a magnitude of H = 7.5. Out of these 29 systems, we detected 18 companions in 16 different systems, resulting in a multiplicity fraction f m = 16 / 29 = 0.55, and a companion fraction of f c = 18 / 29 = 0.62. We observed for the first time 11 of these detected companions. This study concludes that a large survey on more than 120 Northern O-type stars is possible with MIRC-X.

The cosmic ray ionization and γ-ray budgets of star-forming galaxies.

Cosmic rays in star-forming galaxies are a dominant source of both diffuse γ-ray emission and ionization in gas too deeply shielded for photons to penetrate. Though the cosmic rays responsible for γ-rays and ionization are of different energies, they are produced by the same star formation-driven sources, and thus galaxies' star formation rates, γ-ray luminosities, and ionization rates should all be linked. In this paper, we use up-to-date cross-section data to determine this relationship, finding that cosmic rays in a galaxy of star formation rate [Formula: see text] and gas depletion time t dep produce a maximum primary ionization rate ζ ≈ 1 × 10-16(t dep/Gyr)-1 s-1 and a maximum γ-ray luminosity [Formula: see text] erg s-1 in the 0.1-100 GeV band. These budgets imply either that the ionization rates measured in Milky Way molecular clouds include a significant contribution from local sources that elevate them above the Galactic mean, or that CR-driven ionization in the Milky Way is enhanced by sources not linked directly to star formation. Our results also imply that ionization rates in starburst systems are only moderately enhanced compared to those in the Milky Way. Finally, we point out that measurements of γ-ray luminosities can be used to place constraints on galactic ionization budgets in starburst galaxies that are nearly free of systematic uncertainties on the details of cosmic ray acceleration.

Environmental variation of the low-mass IMF.

We use a series of magnetohydrodynamic simulations including both radiative and protostellar outflow feedback to study environmental variation of the initial mass function (IMF). The simulations represent a carefully-controlled experiment whereby we keep all dimensionless parameters of the flow constant except for those related to feedback. We show that radiation feedback suppresses the formation of lower mass objects more effectively as the surface density increases, but this only partially compensates for the decreasing Jeans mass in denser environments. Similarly, we find that protostellar outflows are more effective at suppressing the formation of massive stars in higher surface density environments. The combined effect of these two trends is towards an IMF with a lower characteristic mass and a narrower overall mass range in high surface density environments. We discuss the implications for these findings for the interpretation of observational evidence of IMF variation in early type galaxies.

Identification of interstellar amino acetonitrile in the hot molecular core G10.47＋0.03: Possible glycine survey candidate for the future.

Amino acids are the essential keys that contribute to the study of the formation of life. The simplest amino acid, glycine (NH2CH2COOH), has been searched for a long time in the interstellar medium, but all surveys of glycine have failed. Since the detection of glycine in the interstellar medium was extremely difficult, we aimed to search for the precursor of glycine. After detailed searches of the individual prebiotic molecular species, we successfully identified the emission lines of possible glycine precursor molecule amino acetonitrile (NH2CH2CN) towards the hot molecular core G10.47 + 0.03 using the Atacama Large Millimeter/Submillimeter Array. The estimated statistical column density of amino acetonitrile was (9.10 ± 0.7) × 1015 cm-2 with rotational temperature (Trot) 122 ± 8.8 K. The estimated fractional abundance of amino acetonitrile was 7.01 × 10-8. We found that the estimated fractional abundance of NH2CH2CN fairly agrees with the theoretical value predicted by the three-phase warm-up model from Garrod (2013).

Hierarchical star formation in the Magellanic Clouds with the VMC survey.

The VISTA Magellanic Clouds Survey (VMC) is a near-infrared survey of the Magellanic system. The VMC data has been exploited to detect and study statistically correlated young groups of stars - also known as "young stellar structures" - in the Large and Small Magellanic Clouds (LMC and SMC). We showcase the ~ 3000 recently detected young stellar structures in the LMC and their similarity to the fractal interstellar medium. We discuss how their properties indicate their formation mechanisms and that there are no preferred scales of star formation in the LMC.

A scaling relation for the molecular cloud lifetime in Milky Way-like galaxies.

We study the time evolution of molecular clouds across three Milky Way-like isolated disc galaxy simulations at a temporal resolution of 1 Myr and at a range of spatial resolutions spanning two orders of magnitude in spatial scale from ∼10 pc up to ∼1 kpc. The cloud evolution networks generated at the highest spatial resolution contain a cumulative total of ∼80 000 separate molecular clouds in different galactic-dynamical environments. We find that clouds undergo mergers at a rate proportional to the crossing time between their centroids, but that their physical properties are largely insensitive to these interactions. Below the gas-disc scale height, the cloud lifetime τlife obeys a scaling relation of the form τlife∝ℓ-0.3 with the cloud size ℓ, consistent with over-densities that collapse, form stars, and are dispersed by stellar feedback. Above the disc scale height, these self-gravitating regions are no longer resolved, so the scaling relation flattens to a constant value of ∼13 Myr, consistent with the turbulent crossing time of the gas disc, as observed in nearby disc galaxies.

Magnetic field amplification in accretion discs around the first stars: implications for the primordial IMF.

Magnetic fields play an important role in the dynamics of present-day molecular clouds. Recent work has shown that magnetic fields are equally important for primordial clouds, which form the first stars in the Universe. While the primordial magnetic field strength on cosmic scales is largely unconstrained, theoretical models strongly suggest that a weak seed field existed in the early Universe. We study how the amplification of such a weak field can influence the evolution of accretion discs around first stars, and thus affect the primordial initial mass function (IMF). We perform a suite of 3D ideal magneto-hydrodynamic simulations with different initial field strengths and numerical resolutions. We find that, in simulations with sufficient spatial resolution to resolve the Jeans scale during the collapse, even initially weak magnetic fields grow exponentially to become dynamically important due to both the so-called small-scale turbulent dynamo and the large-scale mean-field dynamo. Capturing the small-scale dynamo action depends primarily on how well we resolve the Jeans length, while capturing the large-scale dynamo depends on the Jeans resolution as well as the maximum absolute resolution. Provided enough resolution, we find that fragmentation does not depend strongly on the initial field strength, because even weak fields grow to become strong. However, fragmentation in runs with magnetic fields differs significantly from those without magnetic fields. We conclude that the development of dynamically strong magnetic fields during the formation of the first stars is likely inevitable, and that these fields had a significant impact on the primordial IMF.

The life cycle of the Central Molecular Zone - II. Distribution of atomic and molecular gas tracers.

We use the hydrodynamical simulation of our inner Galaxy presented in Armillotta et al. to study the gas distribution and kinematics within the Central Molecular Zone (CMZ). We use a resolution high enough to capture the gas emitting in dense molecular tracers such as NH3 and HCN, and simulate a time window of 50 Myr, long enough to capture phases during which the CMZ experiences both quiescent and intense star formation. We then post-process the simulated CMZ to calculate its spatially dependent chemical and thermal state, producing synthetic emission data cubes and maps of both H i and the molecular gas tracers CO, NH3, and HCN. We show that, as viewed from Earth, gas in the CMZ is distributed mainly in two parallel and elongated features extending from positive longitudes and velocities to negative longitudes and velocities. The molecular gas emission within these two streams is not uniform, and it is mostly associated with the region where gas flowing towards the Galactic Centre through the dust lanes collides with gas orbiting within the ring. Our simulated data cubes reproduce a number of features found in the observed CMZ. However, some discrepancies emerge when we use our results to interpret the position of individual molecular clouds. Finally, we show that, when the CMZ is near a period of intense star formation, the ring is mostly fragmented as a consequence of supernova feedback, and the bulk of the emission comes from star-forming molecular clouds. This correlation between morphology and star formation rate should be detectable in observations of extragalactic CMZs.

REFINING THE CENSUS OF THE UPPER SCORPIUS ASSOCIATION WITH GAIA.

We have refined the census of stars and brown dwarfs in the Upper Sco association (~ 10 Myr, ~145 pc) by 1) updating the selection of candidate members from our previous survey to include the high-precision astrometry from the second data release of Gaia, 2) obtaining spectra of a few hundred candidate members to measure their spectral types and verify their youth, and 3) assessing the membership (largely with Gaia astrometry) of 2020 stars toward Upper Sco that show evidence of youth in this work and previous studies. We arrive at a catalog of 1761 objects that are adopted as members of Upper Sco. The distribution of spectral types among the adopted members is similar to those in other nearby star-forming regions, indicating a similar initial mass function. In previous studies, we have compiled mid-infrared photometry from WISE and the Spitzer Space Telescope for members of Upper Sco and used those data to identify the stars that show evidence of circumstellar disks; we present the same analysis for our new catalog of members. As in earlier work, we find that the fraction of members with disks increases with lower stellar masses, ranging from ≲ 10% for > 1 M ⊙ to ~ 22% for 0.01-0.3 M ⊙. Finally, we have estimated the relative ages of Upper Sco and other young associations using their sequences of low-mass stars in M G RP versus G BP - G RP. This comparison indicates that Upper Sco is a factor of two younger than the ß Pic association (21-24 Myr) according to both non-magnetic and magnetic evolutionary models.

A SURVEY FOR NEW STARS AND BROWN DWARFS IN THE OPHIUCHUS STAR-FORMING COMPLEX.

We have performed a survey for new members of the Ophiuchus cloud complex using high-precision astrometry from the second data release of Gaia, proper motions measured with multi-epoch images from the Spitzer Space Telescope, and color-magnitude diagrams constructed with photometry from various sources. Through spectroscopy of candidates selected with those data, we have identified 155 new young stars. Based on available measurements of kinematics, we classify 102, 47, and six of those stars as members of Ophiuchus, Upper Sco, and other populations in Sco-Cen, respectively. We have also assessed the membership of all other stars in the vicinity of Ophiuchus that have spectroscopic evidence of youth from previous studies, arriving at a catalog of 373 adopted members of the cloud complex. For those adopted members, we have compiled mid-IR photometry from Spitzer and the Wide-field Infrared Survey Explorer and have used mid-IR colors to identify and classify circumstellar disks. We find that 210 of the members show evidence of disks, including 48 disks that are in advanced stages of evolution. Finally, we have estimated the relative median ages of the populations near the Ophiuchus clouds and the surrounding Upper Sco association using absolute K-band magnitudes (M K ) based on Gaia parallaxes. If we adopt an age 10 Myr for Upper Sco, then the relative values of M K imply median ages of ~ 2 Myr for L1689 and embedded stars in L1688, 3-4 Myr for low-extinction stars near L1688, and ~ 6 Myr for the group containing ρ Oph.

Effects of grain alignment efficiency on synthetic dust polarization observations of molecular clouds.

It is well known that the polarized continuum emission from magnetically aligned dust grains is determined to a large extent by local magnetic field structure. However, the observed significant anticorrelation between polarization fraction and column density may be strongly affected, perhaps even dominated by variations in grain alignment efficiency with local conditions, in contrast to standard assumptions of a spatially homogeneous grain alignment efficiency. Here we introduce a generic way to incorporate heterogeneous grain alignment into synthetic polarization observations of molecular clouds (MCs), through a simple model where the grain alignment efficiency depends on the local gas density as a power law. We justify the model using results derived from radiative torque alignment theory. The effects of power-law heterogeneous alignment models on synthetic observations of simulated MCs are presented. We find that the polarization fraction-column density correlation can be brought into agreement with observationally determined values through heterogeneous alignment, though there remains degeneracy with the relative strength of cloud-scale magnetized turbulence and the mean magnetic field orientation relative to the observer. We also find that the dispersion in polarization angles-polarization fraction correlation remains robustly correlated despite the simultaneous changes to both observables in the presence of heterogeneous alignment.

Abundances of sulphur molecules in the Horsehead nebula First NS+ detection in a photodissociation region.

CONTEXT: Sulphur is one of the most abundant elements in the Universe (S/H∼1.3×10 -5 ) and plays a crucial role in biological systems on Earth. The understanding of its chemistry is therefore of major importance. AIMS: Our goal is to complete the inventory of S-bearing molecules and their abundances in the prototypical photodissociation region (PDR) the Horsehead nebula to gain insight into sulphur chemistry in UV irradiated regions. Based on the WHISPER (Wide-band High-resolution Iram-30m Surveys at two positions with Emir Receivers) millimeter (mm) line survey, our goal is to provide an improved and more accurate description of sulphur species and their abundances towards the core and PDR positions in the Horsehead. METHODS: The Monte Carlo Markov Chain (MCMC) methodology and the molecular excitation and radiative transfer code RADEX were used to explore the parameter space and determine physical conditions and beam-averaged molecular abundances. RESULTS: A total of 13 S-bearing species (CS, SO, SO2, OCS, H2CS - both ortho and para - HDCS, C2S, HCS+, SO+, H2S, S2H, NS and NS+) have been detected in the two targeted positions. This is the first detection of SO+ in the Horsehead and the first detection of NS+ in any PDR. We find a differentiated chemical behaviour between C-S and O-S bearing species within the nebula. The C-S bearing species C2S and o-H2CS present fractional abundances a factor of > two higher in the core than in the PDR. In contrast, the O-S bearing molecules SO, SO2, and OCS present similar abundances towards both positions. A few molecules, SO+, NS, and NS+, are more abundant towards the PDR than towards the core, and could be considered as PDR tracers. CONCLUSIONS: This is the first complete study of S-bearing species towards a PDR. Our study shows that CS, SO, and H2S are the most abundant S-bearing molecules in the PDR with abundances of ∼ a few 10-9. We recall that SH, SH+, S, and S+ are not observable at the wavelengths covered by the WHISPER survey. At the spatial scale of our observations, the total abundance of S atoms locked in the detected species is < 10-8, only ∼0.1% of the cosmic sulphur abundance.

Gas phase Elemental abundances in Molecular cloudS (GEMS): I. The prototypical dark cloud TMC 1.

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 (3 mm and 2 mm) and the 40m Yebes telescope (1.3 cm and 7 mm) 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 A V ~ 3 to ~20 mag. Two phases with differentiated chemistry can be distinguished: i) the translucent envelope with molecular hydrogen densities of 1-5×103 cm-3; and ii) the dense phase, located at A V > 10 mag, with molecular hydrogen densities >104 cm-3. Observations and modeling show that the gas phase abundances of C and O progressively decrease along the C+/C/CO transition zone (A V ~ 3 mag) where C/H ~ 8×10-5 and C/O~0.8-1, until the beginning of the dense phase at A V ~ 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) ×10-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 ~8×10-8). Based on our chemical modeling, we constrain the value of ζ H2 to ~ (0.5 - 1.8) ×10-16 s-1 in the translucent cloud.

A Census of Star Formation in the Outer Galaxy: The SMOG Field.

In this paper we undertake a study of the 21 deg2 SMOG field, a Spitzer cryogenic mission Legacy program to map a region of the outer Milky Way toward the Perseus and outer spiral arms with the IRAC and MIPS instruments. We identify 4648 YSOs across the field. Using the DBSCAN method, we identify 68 clusters or aggregations of YSOs in the region, having eight or more members. We identify 1197 Class I objects, 2632 Class II objects, and 819 Class III objects, of which 45 are candidate transition disk objects, utilizing the MIPS 24 photometry. The ratio of YSOs identified as members of clusters was 2872/4648, or 62%. The ratios of Class I to Class II YSOs in the clusters are broadly consistent with those found in the inner Galactic and nearby Gould Belt young star formation regions. The clustering properties indicate that the protostars may be more tightly bound to their natal sites than the Class II YSOs, and the Class III YSOs are generally widely distributed. We further perform an analysis of the WISE data of the SMOG field to determine how the lower resolution and sensitivity of WISE affects the identification of YSOs as compared to Spitzer: we identify 931 YSOs using combined WISE and 2MASS photometry, or 20% (931/4648) of the total number identified with Spitzer. Performing the same clustering analysis finds 31 clusters that reliably trace the larger associations identified with the Spitzer data. Twelve of the clusters identified have previously measured distances from the WISE H II survey. SEDFitter modeling of these YSOs is reported, leading to an estimation of the initial mass function in the aggregate of these clusters that approximates that found in the inner Galaxy, implying that the processes behind stellar mass distribution during star formation are not widely affected by the lower density and metallicity of the outer Galaxy.

Thermal Jeans Fragmentation within ~ 1000 AU in OMC-1S.

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.

Revisiting the case of R Monocerotis: Is CO removed at R < 20 au?⋆.

CONTEXT: To our knowledge, R Mon is the only B0 star in which a gaseous Keplerian disk has been detected. However, there is some controversy about the spectral type of R Mon. Some authors propose that it could be a later B8e star, where disks are more common. AIMS: Our goal is to re-evaluate the R Mon spectral type and characterize its protoplanetary disk. METHODS: The spectral type of R Mon has been re-evaluated using the available continuum data and UVES emission lines. We used a power-law disk model to fit previous 12CO 1→0 and 2→1 interferometric observations and the PACS CO data to investigate the disk structure. Interferometric detections of 13CO J=1→0, HCO+ 1→0, and CN 1→0 lines using the IRAM Plateau de Bure Interferometer (PdBI) are presented. The HCN 1→0 line was not detected. RESULTS: Our analysis confirms that R Mon is a B0 star. The disk model compatible with the 12CO 1→0 and 2→1 interferometric observations falls short of predicting the observed fluxes of the 1431 CO lines suggest the existence of a region empty of CO at R≲20 au in the proto-planetary disk. The intense emission of the HCO+ and CN lines shows the strong influence of UV photons on gas chemistry. CONCLUSIONS: The observations gathered in this paper are consistent with the presence of a transition disk with a cavity of R in ≳20 au around R Mon. This size is similar to the photoevaporation radius that supports the interpretation that UV photoevaporation is main disk dispersal mechanism in massive stars.

Discovery of the elusive radical NCO and confirmation of H2NCO+ in space.

The isocyanate radical (NCO) is the simplest molecule containing the backbone of the peptide bond, C(=O)-N. This bond has a prebiotic interest since is the one linking two amino acids to form large chains of proteins. It is also present in some organic molecules observed in space such as HNCO, NH2CHO and CH3NCO. In this letter we report the first detection in space of NCO towards the dense core L483. We also report the identification of the ion H2NCO+, definitively confirming its presence in space, and observations of HNCO, HOCN, and HCNO in the same source. For NCO, we derive a column density of 2.2×1012 cm-2, which means that it is only ~5 times less abundant than HNCO. We find that H2NCO+, HOCN and HCNO have abundances relative to HNCO of 1/400, 1/80, and 1/160, respectively. Both NCO and H2NCO+ are involved in the production of HNCO and several of its isomers. We have updated our previous chemical models involving NCO and the production of the CHNO isomers. Taking into account the uncertainties in the model, the observed abundances are reproduced relatively well. Indeed, the detection of NCO and H2NCO+ in L483 supports the chemical pathways to the formation of the detected CHNO isomers. Sensitive observations of NCO in sources where other molecules containing the C(=O)-N subunit have been detected could help in elucidating its role in prebiotic chemistry in space.

ALMA observations of the young protostellar system Barnard 1b: signatures of an incipient hot corino in B1b-S.

The Barnard 1b core shows signatures of being at the earliest stages of low-mass star formation, with two extremely young and deeply embedded protostellar objects. Hence, this core is an ideal target to study the structure and chemistry of the first objects formed in the collapse of prestellar cores. We present ALMA Band 6 spectral line observations at ~0.6″ of angular resolution towards Barnard 1b. We have extracted the spectra towards both protostars, and used a Local Thermodynamic Equilibrium (LTE) model to reproduce the observed line profiles. B1b-S shows rich and complex spectra, with emission from high energy transitions of complex molecules, such as CH3OCOH and CH3CHO, including vibrational level transitions. We have tentatively detected for the first time in this source emission from NH2CN, NH2CHO, CH3CH2OH, CH2OHCHO, CH3CH2OCOH and both aGg' and gGg' conformers of (CH2OH)2. This is the first detection of ethyl formate (CH3CH2OCOH) towards a low-mass star forming region. On the other hand, the spectra of the FHSC candidate B1b-N are free of COMs emission. In order to fit the observed line profiles in B1b-S, we used a source model with two components: an inner hot and compact component (200 K, 0.35″) and an outer and colder one (60 K, 0.6″). The resulting COM abundances in B1b-S range from 10-13 for NH2CN and NH2CHO, up to 10-9 for CH3OCOH. Our ALMA Band 6 observations reveal the presence of a compact and hot component in B1b-S, with moderate abundances of complex organics. These results indicate that a hot corino is being formed in this very young Class 0 source.

NEW YOUNG STARS AND BROWN DWARFS IN THE UPPER SCORPIUS ASSOCIATION.

To improve the census of the Upper Sco association (~11 Myr, ~145 pc), we have identified candidate members using parallaxes, proper motions, and color-magnitude diagrams from several wide-field imaging surveys and have obtained optical and infrared spectra of several hundred candidates to measure their spectral types and assess their membership. We also have performed spectroscopy on a smaller sample of previously known or suspected members to refine their spectral types and evidence of membership. We have classified 530 targets as members of Upper Sco, 377 of which lack previous spectroscopy. Our new compilation of all known members of the association contains 1631 objects. Although the census of Upper Sco has expanded significantly over the last decade, there remain hundreds of candidates that lack spectroscopy. The precise parallaxes and proper motions from the second data release of Gaia should extend down to substellar masses in Upper Sco, which will greatly facilitate the identification of the undiscovered members.