PRODIGE VI -- Envelope to Disk with NOEMA: VI. The Missing Sulfur Problem

Determining the amount of sulfur in volatiles and refractories in the ISM remains one of the main problems in astrochemistry. The detection of H$_2$S ices, which are thought to be one of the main sulfur reservoirs, has not been achieved yet, and the only S-bearing species detected in the ices to date is OCS. PRODIGE large survey observations with NOEMA of several Class 0/I protostars in the Perseus Molecular Cloud provide a perfect opportunity to study the H$_2$S and OCS composition of the ices through the volatiles sublimated in the warm inner core (T$>$100K, n $\sim10^6$cm$^{-3}$) of these protostars. Our aim is to determine the H$_2$S/OCS ratio in the warm inner core of 24 protostars in order to study how it is affected by different factors during its evolution. We used the NOEMA millimeter observations from the PRODIGE program of H$_2$S, H$_2^{33}$S, OCS, OC$^{33}$S and OC$^{34}$S to estimate the H$_2$S and OCS column densities in the warm inner cores. We used SO and SO$_2$ data from the ALMA archive to give a rough estimate of the total sulfur abundance. We explore the chemistry of H$_2$S and OCS in the warm cores using chemical and dynamical simulations of the collapse of a dense core to form a protostar. The estimated H$_2$S/OCS ratio reveals a segregation of the sources into ``OCS-poor'' and ``OCS-rich'' protostars, where the OCS-poor protostars present higher H$_2$S/OCS ratios than the OCS-rich ones. Total sulfur abundance is always dominated by either H$_2$S or OCS, grows with evolution during the Class 0 phase up to $D_S<8$, and decreases again in the Class I. Simulations show that temperature changes in the pre-stellar phase and during the collapse can produce substantial differences in the H$_2$S and OCS abundances and in the H$_2$S/OCS ratio. Our analysis shows that the H$_2$S/OCS ratio is strongly influenced by the environment and the initial conditions of the cloud.