Stellar associations powering H II regions – I. Defining an evolutionary sequence

Connecting the gas in H II regions to the underlying source of the ionizing radiation can help us constrain the physical processes of stellar feedback and how H II regions evolve over time. With PHANGS–MUSE, we detect nearly 24 000 H II regions across 19 galaxies and measure the physical properties of the ionized gas (e.g. metallicity, ionization parameter, and density). We use catalogues of multiscale stellar associations from PHANGS–HST to obtain constraints on the age of the ionizing sources. We construct a matched catalogue of 4177 H II regions that are clearly linked to a single ionizing association. A weak anticorrelation is observed between the association ages and the H α equivalent width EW(H α), the H α/FUV flux ratio, and the ionization parameter, log q. As all three are expected to decrease as the stellar population ages, this could indicate that we observe an evolutionary sequence. This interpretation is further supported by correlations between all three properties. Interpreting these as evolutionary tracers, we find younger nebulae to be more attenuated by dust and closer to giant molecular clouds, in line with recent models of feedback-regulated star formation. We also observe strong correlations with the local metallicity variations and all three proposed age tracers, suggestive of star formation preferentially occurring in locations of locally enhanced metallicity. Overall, EW(H α) and log q show the most consistent trends and appear to be most reliable tracers for the age of an H II region.