Oxygen, sulfur, and iron radial abundance gradients of classical Cepheids across the Galactic thin disk

Context. Classical Cepheids (CCs) are solid distance indicators and tracers of young stellar populations. Dating back to the beginning of the 20th century, they have been safely adopted to trace the rotation, kinematics, and chemical enrichment history of the Galactic thin disk. Aims. The main aim of this investigation is to provide iron, oxygen, and sulfur abundances for the largest and most homogeneous sample of Galactic CCs analyzed so far (1118 spectra of 356 objects). The current sample, containing 70 CCs for which spectroscopic metal abundances are provided for the first time, covers a wide range in galactocentric distances, pulsation modes, and pulsation periods. Methods. Optical high-resolution spectra with a high signal-to-noise ratio that were collected with different spectrographs were adopted to provide homogeneous estimates of the atmospheric parameters (effective temperature, surface gravity, and microturbulent velocity) that are required to determine the abundance. Individual distances were based either on trigonometric parallaxes by the Gaia Data Release 3 (Gaia DR3) or on distances based on near-infrared period-luminosity relations. Results. We found that iron and α-element radial gradients based on CCs display a well-defined change in the slope for galactocentric distances larger than ~12 kpc. We also found that logarithmic regressions account for the variation in [X/H] abundances from the inner to the outer disk. Radial gradients for the same elements, but based on open clusters covering a wide range in cluster ages, display similar trends. This means that the flattening in the outer disk is an intrinsic feature of the radial gradients because it is independent of age. Empirical evidence indicates that the S radial gradient is steeper than the Fe radial gradient. The difference in the slope is a factor of two in the linear fit (-0.081 vs. -0.041 dex kpc-1) and changes from -1.62 to -0.91 in the logarithmic distance. Moreover, we found that S (explosive nucleosynthesis) is underabundant on average when compared with O (hydrostatic nucleosynthesis). The difference becomes clearer in the metal-poor regime and for the [O/Fe] and [S/Fe] abundance ratios. We performed a detailed comparison with Galactic chemical evolution models and found that a constant star formation efficiency for galactocentric distances larger than 12 kpc accounts for the flattening observed in both iron and α-elements. To further constrain the impact of the predicted S yields for massive stars on radial gradients, we adopted a toy model and found that the flattening in the outermost regions requires a decrease of a factor of four in the current S predictions. Conclusions. CCs are solid beacons for tracing the recent chemical enrichment of young stellar populations. Sulfur photospheric abundances, when compared with other α-elements, have the key advantage of being a volatile element. Therefore, stellar S abundances can be directly compared with nebular sulfur abundances in external galaxies.