Star Formation in the H II Region Sh 2-205: 3D Morphology and Kinematics from Young Stars and Molecular Gas

Abstract: Using Gaia astrometry of young stars combined with CO observations, we present the first systematic three-dimensional (3D) analysis of the structure, kinematics, and evolutionary history of the star-forming regions in the environs of the H II region Sh 2-205 (S205). S205 exhibits a complex morphology and coherent expansion on both global and subregional scales. We identify several O9-B1 stars and a 0.56 Myr old pulsar that are likely associated with the region. A momentum estimate suggests that feedback from these objects may account for the observed overall expansion. Trace-back analysis of the expansion, combined with color-magnitude diagram fitting for young star clusters, indicates at least two episodes of star formation. These results reveal a complex star-formation history of S205 and provide new insights into its 3D evolution.

• The paper combines Gaia DR3 astrometry for 218 young stars with MWISP 12CO and 13CO mapping to place Sh 2-205 at 1,070 ± 15 pc and reconstruct its three-dimensional structure.
• The paper finds global expansion dominated by radial motions averaging 2.9 ± 0.6 km s−1, with local cluster expansion of roughly 0.2–0.9 km s−1 and momentum plausibly supplied by massive-star feedback and a supernova.
• The paper identifies star-formation episodes around 11–12 and 4–5 million years ago, linking the younger episode to triggered formation and highlighting PSR B0355+54 as a possible recent supernova remnant from the region.

Motivation and scope

The H II region Sh 2-205 (S205), a $\sim$70–90 pc nebula in the Camelopardalis OB1 association, has long been studied in two dimensions, but its distance (previously estimated anywhere from 0.5 to 1.1 kpc), driving sources, and evolutionary history remained poorly constrained. This paper combines Gaia DR3 astrometry of young stars with $^{12}$CO and $^{13}$CO (1–0) data from the MWISP survey to construct the first systematic three-dimensional picture of S205's structure, kinematics, and star-formation history (2608.16179). The approach follows earlier Gaia-based 3D analyses of Orion and Scorpius-Centaurus, extending them to the environment of an H II region where prior 3D work was largely absent.

Data and sample construction

The stellar sample is built from seven published YSO catalogs (760 entries at Class II or earlier) cross-matched with Gaia DR3, combined with member stars of four young clusters from the Hunt & Reffert open cluster census. After quality cuts ($ruwe$, $ipd\_gof\_harmonic\_amplitude$, parallax-over-error $>3$), duplicate removal, and iterative outlier rejection, the final catalog contains 218 young stars. Their parallax distribution peaks near 0.9 mas, and the authors adopt a reference-center distance of $1070\pm15$ pc, consistent with extinction-based estimates of $\sim$1067–1114 pc. This resolves the historical distance ambiguity for the region.

Five major molecular clouds are identified with a DBSCAN-based decomposition of the MWISP data. A notable negative result concerns neighboring regions: Sh 2-209, NGC 1491, and IRAS 04000+5052 lie close to S205 on the sky but have CO velocities and distances inconsistent with it, so they are unlikely to be physically associated despite apparent proximity.

Global expansion in 3D

S205 is decomposed into six kinematically coherent subregions (A–F). Because only 16% of young stars have radial velocity measurements—and subregional stellar RVs show sensitivity to uncertainty thresholds—the authors adopt the $T_{\rm MB}$-weighted mean $^{13}$CO velocity as each subregion's RV, an approximation they flag as introducing relatively large uncertainties into the derived 3D motions.

Within this framework, all six subregions exhibit positive radial velocities relative to the reference center, yielding a mean radial expansion of $^{12}$0 km s$^{12}$1, compared to a mean tangential component of only $^{12}$21.4 km s$^{12}$3. Radial motions therefore dominate, providing the first 3D evidence for global expansion of S205. The result is robust to the choice of reference center: using subregions A or C (which host the oldest clusters) as alternative centers still yields positive $^{12}$4 for all other subregions, with subregion A emerging as the most likely expansion center. Unlike the approximately spherical shell found in Canis Major, the S205 subregions are asymmetrically distributed, spanning distances of $^{12}$515–77 pc from the center, so no single spherical shell describes its structure.

At smaller scales, internal kinematics of the clusters within each subregion reveal coherent expansion as well, with mean 2D expansion velocities of $^{12}$60.2–0.9 km s$^{12}$7—consistent with values reported for young stellar systems generally. S205 thus expands simultaneously on $^{12}$8100 pc global and $^{12}$910 pc local scales.

Driving sources and momentum budget

The total radial momentum of the expanding region is estimated at %%%%2$\sim$2%%%%1 $^{13}$2 km s$^{13}$3, combining gas masses from $^{13}$4CO (with $^{13}$5 cm$^{13}$6(K km/s)$^{13}$7; total gas mass $^{13}$8$^{13}$9 $ruwe$0) and stellar masses from mass-function fits. Given an exposed-area fraction of roughly 4% under a simplified spherical-shell assumption, powering this expansion would require feedback from more than $ruwe$110 massive stars or about three supernovae. The authors note this momentum estimate is a lower limit, since a 10% star formation efficiency implies an initial gas reservoir about one-third larger than observed.

Seven candidate massive stars are identified: six O9–B1 stars within the region's boundaries (HD 24431, HD 23675, HD 23800, ALS 7793, HD 25348, HD 24094) plus one runaway candidate, HD 22253 (B0.5 III), whose trace-back trajectory reaches the S205 edge $ruwe$26.5 Myr ago. Four of these had been proposed previously as ionizing sources; HD 23800, HD 25348, and HD 22253 are new candidates. Distance uncertainties leave their true membership unconfirmed.

A particularly suggestive finding concerns PSR B0355+54. Its spin-down age is 0.56 Myr, its VLBI parallax distance ($ruwe$3 kpc) matches S205, FAST scintillation measurements point to a comparable scattering-screen distance, and its proper motion traces back to subregion A. Taken together, these properties strongly suggest the pulsar is the remnant of a massive star formed in S205—a supernova event that likely accelerated the region's recent evolution. Cometary cloud morphologies in subregions D and F, together with outward stellar motions there, further indicate feedback-shaped gas dynamics analogous to structures seen in Orion and Corona Australis.

Star formation history

Two independent trace-back methods—minimizing the sum of pairwise subregion distances and minimizing the mean distance from the reference center—both place the most compact configuration at $ruwe$4 and $ruwe$5 Myr ago, when the subregions were confined within $ruwe$627 pc, about half their current extent. CMD fitting with PARSEC isochrones yields ages of $ruwe$711 Myr for the clusters in subregions A and C, independently validating the trace-back onset. Subregions B and D host younger clusters aged 4–5 Myr, marking a second episode; notably, the age of subregion B matches the $ruwe$84 Myr expansion age of bubble LBN 148.11$ruwe$90.45, supporting a collect-and-collapse triggering scenario. Subregions E and F contain too few stars for reliable fitting and are assigned indicative ages of $ipd\_gof\_harmonic\_amplitude$03 Myr based on typical Class II YSO ages—an assumption rather than a measurement.

The massive-star census corroborates this two-episode picture: the O9 star HD 24431 must belong to the younger episode (O stars rarely survive past 5–6 Myr), while the B0.5–B1 stars are broadly consistent with the earlier epoch, though a more recent origin cannot be excluded. The resulting evolutionary scenario is: an early episode $ipd\_gof\_harmonic\_amplitude$111–12 Myr ago drove global expansion; feedback then triggered a second episode $ipd\_gof\_harmonic\_amplitude$24–5 Myr ago; and a supernova $ipd\_gof\_harmonic\_amplitude$30.56 Myr ago further shaped the region.

Limitations and open questions

Several caveats bear directly on the quantitative results. The substitution of gas RVs for stellar RVs makes the 3D velocities preliminary; only 25 stars (7% of the sample) have RV uncertainties below 5 km s$ipd\_gof\_harmonic\_amplitude$4, and in some subregions the stellar and gas RVs deviate or depend sensitively on the adopted threshold. The momentum analysis assumes uniform injection into a spherical shell with constant radial momentum, neglecting non-uniform motions, unknown initial surface areas, and deviations from symmetry. Cluster ages for E and F rest on assumed rather than fitted parameters, binary effects are ignored in the mass function, and the trace-back analyses assume constant velocities over timescales where Galactic potential effects are small but not zero. Whether the candidate massive stars truly formed in S205, and whether PSR B0355+54's association is physical rather than coincidental, remain to be confirmed with additional high-precision stellar RVs and deeper spectroscopy.

Conclusion

By uniting Gaia DR3 astrometry of 218 young stars with MWISP CO mapping, this work establishes S205 at $ipd\_gof\_harmonic\_amplitude$5 pc, demonstrates coherent multi-scale expansion driven plausibly by a combination of massive-star feedback and a recent supernova, and reconstructs a two-episode star-formation history beginning 11–12 Myr ago. The study provides a template for 3D analyses of H II region environments, while making clear that firmer kinematic conclusions await more complete stellar radial velocity coverage.