The near infrared SED of young star clusters in the FEAST galaxies: Missing ingredients at 1 - 5 $μ$m

Abstract: We present a combined HST and JWST 0.2 - to - 5 $\mu$m analysis of the spectral energy distributions (SEDs) of emerging young star clusters (eYSCs) in four nearby galaxies from the Feedback in Emerging extrAgalactic Star clusTers (FEAST) survey: M51, M83, NGC 628, and NGC 4449. These clusters, selected for their bright Pa$\alpha$ and 3.3 $\mu$m polycyclic aromatic hydrocarbon (PAH) emission, are still associated to their natal gas cloud and have been largely missed in previous HST optical campaigns. We modeled their SEDs using the CIGALE fitting code and identified: i) a systematic flux excess at 1.5 - 2.5 $\mu$m that is not accounted for by current stellar population models; ii) the preference for a set of dust model parameters that is not aligned with expectations from self-consistent analyses of star-forming regions, suggesting model shortcomings also in the 3 - 5 $\mu$m. The near-infrared (NIR) excess is most prominent in low-mass ($\leq 3000$ M$_\odot$) and young ($\leq 6$ Myr) clusters. Additionally, we see that the SED fitting analysis wrongly assigns ages $\geq 6$ Myr to a fraction of strong Pa$\alpha$ emitters with equivalent widths suggestive of significantly younger ages. A parallel analysis with the slug code suggests that stochastic initial mass function (IMF) sampling of pre-main-sequence stars combined with extinction might partially reduce the gap. We conclude that the inclusion of young stellar object SEDs, along with more realistic sampling of the cluster IMF, might be needed to fully account for the stellar population and dust properties of eYSCs.