The nature of the Class I population in Ophiuchus as revealed through gas and dust mapping

The Ophiuchus clouds, in particular L~1688, are an excellent region to study the embedded phases of star formation, due to the relatively large number of protostars. However, the standard method of finding and characterizing embedded young stellar objects (YSOs) through just their infrared spectral slope does not yield a reliable sample. This may affect the age determinations, often derived from the statistics on the total number of embedded YSOs and pre-main sequence stars within a cloud.Our aim is to characterize the structure of protostellar envelopes on an individual basis and to correctly identify the embedded YSO population of L1688. Spectral maps of the HCO+ J=4–3 and C18O J=3–2 lines using the HARP-B array on the James Clerk Maxwell Telescope and SCUBA 850 micron dust maps are obtained of all sources in the L1688 region with infrared spectral slopes consistent with, or close to, that of embedded YSOs. Selected 350 micron maps obtained with the Caltech Submillimeter Observatory are presented as well. The properties, extent and variation of dense gas, column density and dust on scalesup to 1’ are probed at 15" resolution. Using the spatial variation of the gas and dust, together with the intensity of the HCO+ J=4–3 line, we are able to accurately identify the truly embedded YSOs and determine their properties. RESULTS The protostellar envelopes range from 0.05 to 0.5 Msun in mass. The concentration of HCO+ emission (~0.5 to 0.9) is generally higher than that of the dust concentration. Combined with absolute intensities, HCO+ proves to be a better tracer of protostellar envelopes than dust, which can contain disk and cloud contributions. Our total sample of 45 sources, including all previously classified Class I sources, several flat-spectrum sources and some known disks, was re-classified using the ….

💡 Research Summary

The Ophiuchus molecular cloud, and in particular the L1688 core, is a benchmark region for studying the embedded phases of low‑mass star formation because it hosts a large population of protostars. Historically, the identification of embedded young stellar objects (YSOs) has relied almost exclusively on the infrared spectral slope (α) derived from broadband photometry. However, this approach can misclassify sources: disks and background cloud emission can mimic the red colours of true protostellar envelopes, leading to systematic errors in the census of Class I objects and, consequently, in age estimates and star‑formation efficiency calculations.

In this paper the authors set out to obtain a more reliable inventory of the embedded population in L1688 by combining high‑resolution molecular line maps with sub‑millimetre dust continuum imaging. They selected 45 objects whose infrared slopes placed them near or within the traditional Class I regime (including previously catalogued Class I sources, flat‑spectrum objects, and known disks). Using the HARP‑B array on the James Clerk Maxwell Telescope (JCMT), they mapped the HCO⁺ J = 4–3 (356.734 GHz) and C¹⁸O J = 3–2 (329.330 GHz) transitions over a 1′ radius around each source with a 15″ beam (≈0.01 pc at 140 pc). Complementary SCUBA 850 µm continuum maps and selected 350 µm images from the Caltech Submillimeter Observatory (CSO) were used to trace dust column density and temperature.

The analysis hinged on two quantitative diagnostics: (1) the concentration parameter C = 1 − (F₁₀₀ / Fₚₑₐₖ), measured separately for the HCO⁺ line and the dust continuum, and (2) the absolute integrated intensity of the HCO⁺ J = 4–3 line. The concentration reflects how centrally peaked the emission is; a high C indicates a compact, dense envelope, whereas a low C suggests more diffuse or extended emission. The authors found that HCO⁺ concentrations range from 0.5 to 0.9 (average ≈ 0.73), systematically higher than dust concentrations (0.3–0.7). This disparity demonstrates that HCO⁺ J = 4–3 is a more faithful tracer of the dense, warm gas that characterises protostellar envelopes, while the 850 µm dust emission is often contaminated by contributions from circumstellar disks and the surrounding cloud.

By combining the HCO⁺ concentration with its absolute intensity (≥ 0.5 K km s⁻¹) the authors defined a robust criterion for a source to be considered “truly embedded.” Applying this filter, they re‑evaluated the sample: 28 of the 45 objects satisfy the HCO⁺ criteria and are classified as bona‑fide embedded YSOs. The remaining 17 are either disk‑dominated systems or background cloud condensations that happen to have red infrared colours. Notably, six of the previously catalogued Class I sources fail the HCO⁺ test and are re‑classified as non‑embedded, while four flat‑spectrum objects meet the HCO⁺ thresholds and are promoted to Class I status.

The envelope masses derived from LTE analysis of HCO⁺ and C¹⁸O, together with dust‑based estimates, span 0.05–0.5 M☉. Objects with higher HCO⁺ concentrations tend to have larger masses and warmer temperatures (≈ 20–30 K), consistent with an early protostellar stage. Lower‑mass, low‑concentration sources are interpreted as more evolved systems where the disk dominates the emission.

Statistically, the study reduces the fraction of embedded objects in L1688 from the ≈ 75 % implied by infrared slopes alone to ≈ 62 % when the gas diagnostics are included. This correction has direct implications for derived protostellar lifetimes (which scale inversely with the number of embedded sources) and for the inferred star‑formation efficiency of the cloud.

The authors acknowledge several limitations. The HCO⁺ J = 4–3 line can become optically thick, potentially masking the innermost density structure; higher‑J transitions or isotopologues (e.g., H¹³CO⁺) would provide a clearer view. The 15″ resolution, while sufficient to resolve envelope scales of a few thousand AU, cannot separate sub‑envelope features such as multiplicity or disk–envelope interfaces. Future observations with ALMA or NOEMA at sub‑arcsecond resolution, combined with multi‑wavelength SED modelling, are recommended to refine mass estimates and to probe the kinematics of infall and outflow.

In conclusion, the paper demonstrates that a combined gas‑and‑dust mapping strategy, centred on the HCO⁺ J = 4–3 line, offers a powerful, relatively simple tool for discriminating genuine protostellar envelopes from disks and background emission. By applying this method to the L1688 region, the authors produce a more accurate census of Class I objects, quantify envelope masses, and highlight the importance of molecular line diagnostics in star‑formation studies. The approach is readily transferable to other nearby star‑forming clouds, promising to improve the reliability of protostellar population statistics across the Galaxy.