2023ApJ...943...39L

Although episodic star formation (SF) has been suggested for nearby SF regions, a panoramic view of the recent episodic SF history in the solar neighborhood is still missing. By uniformly constraining the slope α of infrared spectral energy distributions (SEDs) of young stellar objects (YSOs) in the 13 largest Gould's Belt (GB) protoclusters surveyed by the Spitzer Space Telescope, we have constructed a cluster-averaged histogram of α representing the YSO evolution lifetime as a function of the α value. Complementary to the traditional SED classification scheme (0, I, F, II, and III) that is based on different α values, a staging scheme (A, B, C, D, and E) of SED evolution is proposed on the basis of the α statistical features that can be better matched to the physical stages of disk dissipation and giant planet formation. This has also allowed us to unravel the fluctuations of SF rate (SFR) in the 3 Myr history of these GB protoclusters. Diverse evolutionary patterns such as single peaks, double peaks, and ongoing acceleration of SFR are revealed. The SFR fluctuations are between 20% and 60% (∼40% on average) and no dependence on the average SFR or the number of SFR episodes is found. However, spatially close protoclusters tend to have similar SFR fluctuation trends, indicating that the driving force of the fluctuations should be at size scales beyond the typical cluster sizes of several parsecs.