Nearby Star-forming Regions Research Articles

Multi-line characterization of whole molecular clouds using stratified random sampling

We have developed a new technique to characterize the multi-line emission from molecular clouds using statistical sampling. Our method uses available extinction maps to select a relatively small sample of cloud positions that cover the full range of column densities in the cloud, and that can be observed with only a modest investment of telescope time. Here we present the first results of applying this technique to the Perseus molecular cloud, a nearby star-forming region that contains a population of both isolated and clustered protostars. We have used the IRAM 30m telescope to observe a sample of 100 random positions that cover two orders of magnitude in H2 column density. These positions have been observed over the full 3mm wavelength band together with selected portions of the 2 and 1mm bands. We find that the emission properties of most species are strongly correlated with the column density over the whole cloud, and that they can be reproduced using a relative simple radiative transfer model.

A Mid-infrared Study of Directly Imaged Planetary-mass Companions Using Archival Spitzer/IRAC Images

Abstract The atmospheres and accretion disks of planetary-mass and substellar companions provide an unprecedented look into planet and moon formation processes, most notably the frequency and lifetime of circumplanetary disks. In our ongoing effort to leverage the extraordinary sensitivity of the Spitzer/Infrared Array Camera (IRAC) at 3.6, 4.5, 5.8, and 8.0 μm to study wide planetary-mass and substellar companions near the diffraction limit, we present point-spread function fitting photometry of archival Spitzer/IRAC images for nine stars (G0 to M4+M7) in nearby star-forming regions or stellar associations that host companions at separations of ρ = 1.″17–12.″33. We detect all system primaries in all four IRAC channels and recover eight low-mass companions in at least one IRAC channel for our sample, five of which have not been resolved previously in IRAC images. We measure nonphotospheric [3.6]–[8.0] colors for four of the system companions (DH Tau B, 2M0441 B, SR 12 c, and ROXs 42B b), confirming or discovering the presence of circumstellar or circum(sub)stellar disks. We detect fluxes consistent with photospheric emission for four other companions (AB Pic b, CHXR 73 b, 1RXS J1609 b, and HD 203030 b) that are unlikely to host disks. Combined with past detections of accretion or disk indicators, we determine the global disk frequency of young (<15 Myr) wide companions with masses near the deuterium-burning limit to be 56% ± 12%.

Chemically Fresh Gas Inflows Detected in a Nearby High-mass Star-forming Region

Abstract We report the detection of a chemically fresh inflow that is feeding high-mass young-stellar-object (HMYSO) growth in the nearby high-mass star-forming region G352.63 made with both the Atacama Large Millimeter/submillimeter Array (ALMA) and the Submillimeter Array (SMA). High-quality images of the dust and molecular lines from both ALMA and SMA have consistently revealed a gravitationally controlled cold (∼10 K) gas inflow of chemically fresh molecules (e.g., CCH and HC3N) toward the central HMYSO and its surrounding dense gas structure, which has a possible torus- or disk-like morphology. The HMYSO is also observed to have an outflow, which is nearly perpendicular to the torus and its parental filament, and thus can be clearly separated from the inflows. These kinematic features provide observational evidence to support the conjecture that the infalling streamers in high-mass star-forming regions could proceed in a similar process to that observed in low-mass counterparts. The chemically fresh infalling streamers could also be involved in the disk or torus configuration, fragmentation, and accretion bursts that occur in both simulations and observations.

Small Protoplanetary Disks in the Orion Nebula Cluster and OMC1 with ALMA

Abstract The Orion Nebula Cluster (ONC) is the nearest dense star-forming region at ∼400 pc away, making it an ideal target to study the impact of high stellar density and proximity to massive stars (the Trapezium) on protoplanetary disk evolution. The OMC1 molecular cloud is a region of high extinction situated behind the Trapezium in which actively forming stars are shielded from the Trapezium’s strong radiation. In this work, we survey disks at high resolution with Atacama Large Millimeter/submillimeter Array at three wavelengths with resolutions of 0.″095 (3 mm; Band 3), 0.″048 (1.3 mm; Band 6), and 0.″030 (0.85 mm; Band 7) centered on radio Source I. We detect 127 sources, including 15 new sources that have not previously been detected at any wavelength. 72 sources are spatially resolved at 3 mm, with sizes from ∼8–100 au. We classify 76 infrared-detected sources as foreground ONC disks and the remainder as embedded OMC1 disks. The two samples have similar disk sizes, but the OMC1 sources have a dense and centrally concentrated spatial distribution, indicating they may constitute a spatially distinct subcluster. We find smaller disk sizes and a lack of large (>75 au) disks in both our samples compared to other nearby star-forming regions, indicating that environmental disk truncation processes are significant. While photoevaporation from nearby massive Trapezium stars may account for the smaller disks in the ONC, the embedded sources in OMC1 are hidden from this radiation and thus must truncated by some other mechanism, possibly dynamical truncation or accretion-driven contraction.

Investigating the structure of star-forming regions using INDICATE

ABSTRACT The ability to make meaningful comparisons between theoretical and observational data of star-forming regions is key to understanding the star formation process. In this paper, we test the performance of INDICATE, a new method to quantify the clustering tendencies of individual stars in a region, on synthetic star-forming regions with substructured, and smooth, centrally concentrated distributions. INDICATE quantifies the amount of stellar affiliation of each individual star, and also determines whether this affiliation is above random expectation for the star-forming region in question. We show that INDICATE cannot be used to quantify the overall structure of a region due to a degeneracy when applied to regions with different geometries. We test the ability of INDICATE to detect differences in the local stellar surface density and its ability to detect and quantify mass segregation. We then compare it to other methods such as the mass segregation ratio ΛMSR, the local stellar surface density ratio ΣLDR, and the cumulative distribution of stellar positions. INDICATE detects significant differences in the clustering tendencies of the most massive stars when they are at the centre of a smooth, centrally concentrated distribution, corresponding to areas of greater stellar surface density. When applied to a subset of the 50 most massive stars, we show INDICATE can detect signals of mass segregation. We apply INDICATE to the following nearby star-forming regions: Taurus, ONC, NGC 1333, IC 348, and ρ Ophiuchi and find a diverse range of clustering tendencies in these regions.

Bridging the Gap between Protoplanetary and Debris Disks: Separate Evolution of Millimeter and Micrometer-sized Dust

Abstract The connection between the nature of a protoplanetary disk and that of a debris disk is not well understood. Dust evolution, planet formation, and disk dissipation likely play a role in the processes involved. We aim to reconcile both manifestations of dusty circumstellar disks through a study of optically thin Class III disks and how they correlate to younger and older disks. In this work, we collect literature and Atacama Large Millimeter/submillimeter Array archival millimeter fluxes for 85 disks (8%) of all Class III disks across nearby star-forming regions. We derive millimeter-dust masses M dust and compare these with Class II and debris disk samples in the context of excess infrared luminosity, accretion rate, and age. The mean M dust of Class III disks is 0.29 ± 0.19 M ⊕. We propose a new evolutionary scenario wherein radial drift is very efficient for nonstructured disks during the Class II phase resulting in a rapid M dust decrease. In addition, we find possible evidence for long infrared protoplanetary disk timescales, ∼8 Myr, consistent with overall slow disk evolution. In structured disks, the presence of dust traps allows for the formation of planetesimal belts at large radii, such as those observed in debris disks. We propose therefore that the planetesimal belts in debris disks are the result of dust traps in structured disks, whereas protoplanetary disks without dust traps decrease in dust mass through radial drift and are therefore undetectable as debris disks after the gas dissipation. These results provide a hypothesis for a novel view of disk evolution.

The Young Stellar Population, Distance, and Cloud–Cloud Collision Induced Star Formation Scenario of the Trifid Nebula

Abstract The Trifid Nebula is a young, nearby star-forming region where star formation is proposed to have been triggered by cloud–cloud collision (CCC), based on observations of molecular clouds. It offers a unique opportunity to test whether the CCC hypothesis is supported by the spatial distribution and star formation chronology of young stars. We present the first study of the optically visible pre-main sequence (PMS) population of the region using riHα imaging and Gaia astrometry. Combined with an analysis of young stellar objects (YSOs) using infrared imaging, we capture the spatial distribution and star formation chronology of the young stellar population. From the analysis, 15 Flat/Class I YSOs, 46 Class II YSOs, and 41 accreting PMS stars are identified (diskless/non-accreting sources are not included in the analysis). The distance based on Gaia parallaxes is ∼1250 pc, significantly closer than previously reported. The Class II YSOs and PMS stars (∼1.5 Myr old) are spread toward the edge of the molecular clouds. They are slightly younger than the estimated crossing time of ∼2.7 Myr and closer to the estimated dynamical age ∼0.85 Myr. Younger Class I YSOs are more concentrated spatially. There exists a cavity devoid of young stars where the two clouds overlap. This evidence suggests that the current generation of stars formed after the collision of two clouds ∼1 Myr ago, and this result can be corroborated using future spectroscopic studies.

The JCMT Transient Survey: Four-year Summary of Monitoring the Submillimeter Variability of Protostars

Abstract We present the four-year survey results of monthly submillimeter monitoring of eight nearby (<500 pc) star-forming regions by the JCMT Transient Survey. We apply the Lomb–Scargle Periodogram technique to search for and characterize variability on 295 submillimeter peaks brighter than 0.14 Jy beam−1, including 22 disk sources (Class II), 83 protostars (Class 0/I), and 190 starless sources. We uncover 18 secular variables, all of them protostars. No single-epoch burst or drop events and no inherently stochastic sources are observed. We classify the secular variables by their timescales into three groups: Periodic, Curved, and Linear. For the Curved and Periodic cases, the detectable fractional amplitude, with respect to mean peak brightness, is ∼4% for sources brighter than ∼0.5 Jy beam−1. Limiting our sample to only these bright sources, the observed variable fraction is 37% (16 out of 43). Considering source evolution, we find a similar fraction of bright variables for both Class 0 and Class I. Using an empirically motivated conversion from submillimeter variability to variation in mass accretion rate, six sources (7% of our full sample) are predicted to have years-long accretion events during which the excess mass accreted reaches more than 40% above the total quiescently accreted mass: two previously known eruptive Class I sources, V1647 Ori and EC 53 (V371 Ser), and four Class 0 sources, HOPS 356, HOPS 373, HOPS 383, and West 40. Considering the full protostellar ensemble, the importance of episodic accretion on few years timescale is negligible—only a few percent of the assembled mass. However, given that this accretion is dominated by events on the order of the observing time window, it remains uncertain as to whether the importance of episodic events will continue to rise with decades-long monitoring.

Quantifying Variability of Young Stellar Objects in the Mid-infrared Over 6 Years with the Near-Earth Object Wide-field Infrared Survey Explorer

Abstract Variability in young stellar objects (YSOs) can be caused by various time-dependent phenomena associated with star formation, including accretion rates, geometric changes in the circumstellar disks, stochastic hydromagnetic interactions between stellar surfaces and inner-disk edges, reconnections within the stellar magnetosphere, and hot/cold spots on stellar surfaces. We uncover and characterize ∼1700 variables from a sample of ∼5400 YSOs in nearby low-mass star-forming regions using mid-IR light curves obtained from the 6.5 yr Near-Earth Object Wide-field Infrared Survey Explorer All Sky Survey. The mid-IR variability traces a wide range of dynamical, physical, and geometrical phenomenon. We classify six types of YSO mid-IR variability based on their light curves: secular variability (linear, curved, and periodic) and stochastic variability (burst, drop, and irregular). YSOs in earlier evolutionary stages have higher fractions of variables and higher amplitudes for the variability, with the recurrence timescale of FUor-type outbursts (defined here as ΔW1 or ΔW2 > 1 mag followed by inspection of candidates) of ∼1000 yr in the early embedded protostellar phase. Known eruptive young stars and subluminous objects show fractions of variables similar to the fraction (∼55%) found in typical protostars, suggesting that these two distinct types are not distinct in variability over the 6.5 yr timescale. Along with brightness variability, we also find a diverse range of secular color variations, which can be attributed to a competitive interplay between the variable accretion luminosity of the central source and the variable extinction by material associated with the accretion process.

Multicomponent, multiwavelength benchmarks for source- and filament-extraction methods

Modern multiwavelength observations of star-forming regions that reveal complex, highly structured molecular clouds require adequate extraction methods that provide both complete detection of the structures and their accurate measurements. The omnipresence of filamentary structures and their physical connection to prestellar cores make it necessary to use methods that are able to disentangle and extract both sources and filaments. It is fundamentally important to test all extraction methods on standard benchmarks to compare their detection and measurement qualities and fully understand their capabilities before their scientific applications. A recent publication described getsf, the new method for source and filament extraction that employs the separation of the structural components, a successor to getsources, getfilaments, and getimages (collectively referred to as getold). This new paper describes a detailed benchmarking of both getsf and getold using two multicomponent, multiwavelength benchmarks resembling the Herschel observations of the nearby star-forming regions. Each benchmark consists of simulated images at six Herschel wavelengths and one additional derived surface density image with a 13″ resolution. The structural components of the benchmark images include a background cloud, a dense filament, hundreds of starless and protostellar cores, and instrumental noise. Five variants of benchmark images of different complexity are derived from the two benchmarks and are used to perform the source and filament extractions with getsf and getold. A formalism for evaluating source detection and measurement qualities is presented, allowing quantitative comparisons of different extraction methods in terms of their completeness, reliability, and goodness, as well as the detection and measurement accuracies and the overall quality. A detailed analysis of the benchmarking results shows that getsf has better source and filament extraction qualities than getold and that the best choice of the images for source detection with getsf is the high-resolution surface density, alone or with the other available Herschel images. The benchmarks explored in this paper provide the standard tests for calibrating existing and future source- and filament-extraction methods to choose the best tool for astrophysical studies.

The Per-Tau Shell: A Giant Star-forming Spherical Shell Revealed by 3D Dust Observations

Abstract A major question in the field of star formation is how molecular clouds form out of the diffuse interstellar medium (ISM). Recent advances in 3D dust mapping are revolutionizing our view of the structure of the ISM. Using the highest-resolution 3D dust map to date, we explore the structure of a nearby star-forming region, which includes the well-known Perseus and Taurus molecular clouds. We reveal an extended near-spherical shell, 156 pc in diameter (hereafter called the “Per-Tau Shell”), in which the Perseus and Taurus clouds are embedded. We also find a large ring structure at the location of Taurus (hereafter called the “Tau Ring”). We discuss a formation scenario for the Per-Tau Shell, in which previous stellar and supernova feedback events formed a large expanding shell, where the swept-up ISM has condensed to form both the shell and the Perseus and Taurus molecular clouds within it. We present auxiliary observations of H i, Hα, 26Al, and X-rays that further support this scenario, and estimate the Per-Tau Shell’s age to be ≈6–22 Myrs. The Per-Tau shell offers the first 3D observational view of a phenomenon long-hypothesized theoretically, molecular cloud formation and star formation triggered by previous stellar and supernova feedback.

A Solar System formation analogue in the Ophiuchus star-forming complex

Anomalies among the daughter nuclei of the extinct short-lived radionuclides (SLRs) in the calcium-aluminum-rich inclusions (CAIs) indicate that the Solar System must have been born near a source of the SLRs so that they could be incorporated before they decayed away. $\gamma$-rays from one such living SLR, $^{26}$Al, are detected in only a few nearby star-forming regions. Here we employ multi-wavelength observations to demonstrate that one such region, Ophiuchus, containing many pre-stellar cores that may serve as analogs for the emerging Solar System, is inundated with $^{26}$Al from the neighboring Upper-Scorpius association, and so may provide concrete guidance for how SLR enrichment proceeded in the Solar System complementary to the meteoritics. We demonstrate via Bayesian forward modeling drawing on a wide range of observational and theoretical results that this $^{26}$Al likely 1) arises from supernova explosions, 2) arises from multiple stars, 3) has enriched the gas prior to the formation of the cores, and 4) gives rise to a broad distribution of core enrichment spanning about two orders of magnitude. This means that if the spread in CAI ages is small, as it is in the Solar System, protoplanetary disks must suffer a global heating event.

The dense warm ionized medium in the inner Galaxy

Context. Ionized interstellar gas is an important component of the interstellar medium and its lifecycle. The recent evidence for a widely distributed highly ionized warm interstellar gas with a density intermediate between the warm ionized medium (WIM) and compact H II regions suggests that there is a major gap in our understanding of the interstellar gas. Aims. Our goal is to investigate the properties of the dense WIM in the Milky Way using spectrally resolved SOFIA GREAT [N II] 205 μm fine-structure lines and Green Bank Telescope hydrogen radio recombination lines (RRL) data, supplemented by spectrally unresolved Herschel PACS [N II] 122μm data, and spectrally resolved 12CO. Methods. We observed eight lines of sight (LOS) in the 20° < l < 30° region in the Galactic plane. We analyzed spectrally resolved lines of [N II] at 205 μm and RRL observations, along with the spectrally unresolved Herschel PACS 122 μm emission, using excitation and radiative transfer models to determine the physical parameters of the dense WIM. We derived the kinetic temperature, as well as the thermal and turbulent velocity dispersions from the [N II] and RRL linewidths. Results. The regions with [N II] 205 μm emission are characterized by electron densities, n(e) ~ 10−35 cm−3, temperatures range from 3400 to 8500 K, and nitrogen column densities N(N+) ~ 7 × 1016 to 3 × 1017 cm−2. The ionized hydrogen column densities range from 6 × 1020 to 1.7 × 1021 cm−2 and the fractional nitrogen ion abundance x(N+) ~ 1.1 × 10−4 to 3.0 × 10−4, implying an enhanced nitrogen abundance at a distance ~4.3 kpc from the Galactic Center. The [N II] 205 μm emission lines coincide with CO emission, although often with an offset in velocity, which suggests that the dense warm ionized gas is located in, or near, star-forming regions, which themselves are associated with molecular gas. Conclusions. These dense ionized regions are found to contribute ≳50% of the observed [C II] intensity along these LOS. The kinetic temperatures we derive are too low to explain the presence of N+ resulting from electron collisional ionization and/or proton charge transfer of atomic nitrogen. Rather, these regions most likely are ionized by extreme ultraviolet (EUV) radiation from nearby star-forming regions or as a result of EUV leakage through a clumpy and porous interstellar medium.

X-Ray Superflares from Pre-main-sequence Stars: Flare Energetics and Frequency

Abstract Solar-type stars exhibit their highest levels of magnetic activity during their early convective pre-main-sequence (PMS) phase of evolution. The most powerful PMS flares, superflares and megaflares, have peak X-ray luminosities of log ( L X ) = 30.5 – 34.0 erg s−1 and total energies log ( E X ) = 34 – 38 erg. Among >24,000 X-ray-selected young (t ≲ 5 Myr) members of 40 nearby star-forming regions from our earlier Chandra MYStIX and SFiNCs surveys, we identify and analyze a well-defined sample of 1086 X-ray superflares and megaflares, the largest sample ever studied. Most are considerably more powerful than optical/X-ray superflares detected on main-sequence stars. This study presents energy estimates of these X-ray flares and the properties of their host stars. These events are produced by young stars of all masses over evolutionary stages ranging from protostars to diskless stars, with the occurrence rate positively correlated with stellar mass. Flare properties are indistinguishable for disk-bearing and diskless stars indicating star–disk magnetic fields are not involved. A slope of α ≃ 2 in the flare energy distributions dN / dE X ∝ E X − α is consistent with that of optical/X-ray flaring from older stars and the Sun. Megaflares ( log ( E X ) > 36.2 erg) from solar-mass stars have an occurrence rate of 1.7 − 0.6 + 1.0 flares/star/year and contribute at least 10%–20% to the total PMS X-ray energetics. These explosive events may have important astrophysical effects on protoplanetary disk photoevaporation, ionization of disk gas, production of spallogenic radionuclides in disk solids, and hydrodynamic escape of young planetary atmospheres. Our following paper details plasma and magnetic loop modeling of the >50 brightest X-ray megaflares.

A Stellar Mass Dependence of Structured Disks: A Possible Link with Exoplanet Demographics

Abstract Gaps in protoplanetary disks have long been hailed as signposts of planet formation. However, a direct link between exoplanets and disks remains hard to identify. We present a large sample study of ALMA disk surveys of nearby star-forming regions to disentangle this connection. All disks are classified as either structured (transition, ring, extended) or nonstructured (compact) disks. Although low-resolution observations may not identify large-scale substructure, we assume that an extended disk must contain substructure from a dust evolution argument. A comparison across ages reveals that structured disks retain high dust masses up to at least 10 Myr, whereas the dust mass of compact, nonstructured disks decreases over time. This can be understood if the dust mass evolves primarily by radial drift, unless drift is prevented by pressure bumps. We identify a stellar mass dependence of the fraction of structured disks. We propose a scenario linking this dependence with that of giant exoplanet occurrence rates. We show that there are enough exoplanets to account for the observed disk structures if transitional disks are created by exoplanets more massive than Jupiter and ring disks by exoplanets more massive than Neptune, under the assumption that most of those planets eventually migrate inwards. On the other hand, the known anticorrelation between transiting super-Earths and stellar mass implies those planets must form in the disks without observed structure, consistent with formation through pebble accretion in drift-dominated disks. These findings support an evolutionary scenario where the early formation of giant planets determines the disk’s dust evolution and its observational appearance.

Compact disks. An explanation to faint CO emission in Lupus disks

Context. ALMA disk surveys have shown that a large fraction of observed protoplanetary disks in nearby star-forming regions (SFRs) are fainter than expected in CO isotopolog emission. Disks not detected in 13CO line emission are also faint and often unresolved in the continuum emission at an angular resolution of around 0.2 arcsec. Aims. Focusing on the Lupus SFR, the aim of this work is to investigate whether this population comprises radially extended and low-mass disks – as commonly assumed so far – or intrinsically radially compact disks, an interpretation that we propose in this paper. The latter scenario was already proposed for individual sources or small samples of disks, while this work targets a large population of disks in a single young SFR for which statistical arguments can be made. Methods. We ran a new grid of physical–chemical models of compact disks with the physical–chemical code DALI in order to cover a region of the parameter space that has not been explored before with this code. We compared these models with 12CO and 13CO ALMA observations of faint disks in the Lupus SFR, and report the simulated integrated continuum and CO isotopolog fluxes of the new grid of compact models. Results. Lupus disks that are not detected in 13CO emission and have faint or undetected 12CO emission are consistent with compact disk models. For disks with a limited radial extent, the emission of CO isotopologs is mostly optically thick and scales with the surface area, that is, it is fainter for smaller objects. The fraction of compact disks is potentially between roughly 50% and 60% of the entire Lupus sample. Deeper observations of 12CO and 13CO at a moderate angular resolution will allow us to distinguish whether faint disks are intrinsically compact or extended but faint, without the need to resolve them. If the fainter end of the disk population observed by ALMA disk surveys is consistent with such objects being very compact, this will either create a tension with viscous spreading or require MHD winds or external processes to truncate the disks.

A Low-mass Cold and Quiescent Core Population in a Massive Star Protocluster

Abstract Pre-stellar cores represent the initial conditions of star formation. Although these initial conditions in nearby low-mass star-forming regions have been investigated in detail, such initial conditions remain vastly unexplored for massive star-forming regions. We report the detection of a cluster of low-mass starless and pre-stellar core candidates in a massive star protocluster-forming cloud, NGC 6334S. With the Atacama Large Millimeter/submillimeter Array (ALMA) observations at a ∼0.02 pc spatial resolution, we identified 17 low-mass starless core candidates that do not show any evidence of protostellar activity. These candidates present small velocity dispersions, high fractional abundances of NH2D, high NH3 deuterium fractionations, and are completely dark in the infrared wavelengths from 3.6 up to 70 μm. Turbulence is significantly dissipated and the gas kinematics are dominated by thermal motions toward these candidates. Nine out of the 17 cores are gravitationally bound, and therefore are identified as pre-stellar core candidates. The embedded cores of NGC 6334S show a wide diversity in masses and evolutionary stages.

Ultraviolet spectra of extreme nearby star-forming regions: Evidence for an overabundance of very massive stars

ABSTRACT As deep spectroscopic campaigns extend to higher redshifts and lower stellar masses, the interpretation of galaxy spectra depends increasingly upon models for very young stellar populations. Here we present new HST/COS ultraviolet spectroscopy of seven nearby (<120 Mpc) star-forming regions hosting very young stellar populations (∼4–20 Myr) with optical Wolf–Rayet stellar wind signatures, ideal laboratories in which to benchmark these stellar models. We detect nebular C iii] in all seven, but at equivalent widths uniformly <10 Å. This suggests that even for very young stellar populations, the highest equivalent width C iii] emission at ≥15 Å is reserved for inefficiently cooled gas at metallicities at or below that of the SMC. The spectra also reveal strong C iv P-Cygni profiles and broad He ii emission formed in the winds of massive stars, including some of the most prominent He ii stellar wind lines ever detected in integrated spectra. We find that the latest stellar population synthesis prescriptions with improved treatment of massive stars nearly reproduce the entire range of stellar He ii wind strengths observed here. However, we find that these models cannot simultaneously match the strongest wind features alongside the optical nebular line constraints. This discrepancy can be naturally explained by an overabundance of very massive stars produced by a high incidence of binary mass transfer and mergers occurring on short ≲10 Myr time-scales, suggesting these processes may be crucial for understanding systems dominated by young stars both nearby and in the early Universe.

A neural network-based methodology to select young stellar object candidates from IR surveys

Context. Observed young stellar objects (YSOs) are used to study star formation and characterize star-forming regions. For this purpose, YSO candidate catalogs are compiled from various surveys, especially in the infrared (IR), and simple selection schemes in color-magnitude diagrams (CMDs) are often used to identify and classify YSOs. Aims. We propose a methodology for YSO classification through machine learning (ML) using Spitzer IR data. We detail our approach in order to ensure reproducibility and provide an in-depth example on how to efficiently apply ML to an astrophysical classification. Methods. We used feedforward artificial neural networks (ANNs) that use the four IRAC bands (3.6, 4.5, 5.8, and 8 μm) and the 24 μm MIPS band from Spitzer to classify point source objects into CI and CII YSO candidates or as contaminants. We focused on nearby (≲1 kpc) star-forming regions including Orion and NGC 2264, and assessed the generalization capacity of our network from one region to another. Results. We found that ANNs can be efficiently applied to YSO classification with a contained number of neurons (∼25). Knowledge gathered on one star-forming region has shown to be partly efficient for prediction in new regions. The best generalization capacity was achieved using a combination of several star-forming regions to train the network. Carefully rebalancing the training proportions was necessary to achieve good results. We observed that the predicted YSOs are mainly contaminated by under-constrained rare subclasses like Shocks and polycyclic aromatic hydrocarbons (PAHs), or by the vastly dominant other kinds of stars (mostly on the main sequence). We achieved above 90% and 97% recovery rate for CI and CII YSOs, respectively, with a precision above 80% and 90% for our most general results. We took advantage of the great flexibility of ANNs to define, for each object, an effective membership probability to each output class. Using a threshold in this probability was found to efficiently improve the classification results at a reasonable cost of object exclusion. With this additional selection, we reached 90% and 97% precision on CI and CII YSOs, respectively, for more than half of them. Our catalog of YSO candidates in Orion (365 CI, 2381 CII) and NGC 2264 (101 CI, 469 CII) predicted by our final ANN, along with the class membership probability for each object, is publicly available at the CDS. Conclusions. Compared to usual CMD selection schemes, ANNs provide a possibility to quantitatively study the properties and quality of the classification. Although some further improvement may be achieved by using more powerful ML methods, we established that the result quality depends mostly on the training set construction. Improvements in YSO identification with IR surveys using ML would require larger and more reliable training catalogs, either by taking advantage of current and future surveys from various facilities like VLA, ALMA, or Chandra, or by synthesizing such catalogs from simulations.

Multi-scale view of star formation in IRAS 21078+5211: from clump fragmentation to disk wind

Context. Star formation (SF) is a multi-scale process in which the mode of fragmentation of the collapsing clump on scales of 0.1–1 pc determines the mass reservoir and affects the accretion process of the individual protostars on scales of 10–100 au. Aims. We want to investigate the nearby (located at 1.63 ± 0.05 kpc) high-mass star-forming region IRAS 21078+5211 at linear scales from ~1 pc down to ~10 au. Methods. We combine the data of two recent programs: the NOrthern Extended Millimeter Array large project CORE and the Protostellar Outflows at the EarliesT Stages (POETS) survey. The former provides images of the 1 mm dust continuum and molecular line emissions with a linear resolution of ≈600 au covering a field of view up to ≈0.5 pc. The latter targets the ionized gas and 22 GHz water masers, mapping linear scales from a few 103 au down to a few astronomical units. Results. In IRAS 21078+5211, a highly fragmented cluster (size ~0.1 pc) of molecular cores is observed, located at the density peak of an elongated (size ~1 pc) molecular cloud. A small (≈1 km s−1 per 0.1 pc) LSR velocity (VLSR) gradient is detected across the major axis of the molecular cloud. Assuming we are observing a mass flow from the harboring cloud to the cluster, we derive a mass infall rate of ≈10−4 M⊙ yr−1. The most massive cores (labeled 1, 2, and 3) are found at the center of the cluster, and these are the only ones that present a signature of protostellar activity in terms of emission from high-excitation molecular lines or a molecular outflow. The masses of the young stellar objects (YSOs) inside these three cores are estimated in the range 1–6 M⊙. We reveal an extended (size ~0.1 pc), bipolar collimated molecular outflow emerging from core 1. We believe this is powered by the compact (size ≲1000 au) radio jet discovered in the POETS survey, ejected by a YSO embedded in core 1 (named YSO-1), since the molecular outflow and the radio jet are almost parallel and have a comparable momentum rate. By means of high-excitation lines, we find a large (≈14 km s−1 over 500 au) VLSR gradient at the position of YSO-1, oriented approximately perpendicular to the radio jet. Assuming this is an edge-on, rotating disk and fitting a Keplerian rotation pattern, we determine the YSO-1 mass to be 5.6 ± 2.0 M⊙. The water masers observed in the POETS survey emerge within 100–300 au from YSO-1 and are unique tracers of the jet kinematics. Their three-dimensional (3D) velocity pattern reveals that the gas flows along, and rotates about, the jet axis. We show that the 3D maser velocities are fully consistent with the magneto-centrifugal disk-wind models predicting a cylindrical rotating jet. Under this hypothesis, we determine the jet radius to be ≈ 16 au and the corresponding launching radius and terminal velocity to be ≈ 2.2 au and ≈ 200 km s−1, respectively. Conclusions. Complementing high-angular resolution, centimeter and millimeter interferometric observations in thermal tracers with Very Long Baseline Interferometry of molecular masers, is invaluable in studying high-mass SF. The combination of these twodatasets allows us to connect the events that we see at large scales, as clump fragmentation and mass flows, with the physical processes identified at small scales, specifically, accretion and ejection in disk-jet systems.