We present the results of a new study investigating the relationship between
observed Lyα equivalent width (Wλ(Lyα)) and the
metallicity of the ionizing stellar population (Z⋆) for a sample of
768 star-forming galaxies at 3≤z≤5 drawn from the VANDELS survey.
Dividing our sample into quartiles of rest-frame Wλ(Lyα)
across the range -58 \unicode{xC5} \lesssimWλ(Lyα)
\lesssim 110 \unicode{xC5} we determine Z⋆ from full spectral
fitting of composite far-ultraviolet (FUV) spectra and find a clear
anti-correlation between Wλ(Lyα) and Z⋆. Our results
indicate that Z⋆ decreases by a factor ≳3 between the lowest
Wλ(Lyα) quartile
(⟨Wλ(Lyα)\rangle=-18\unicode{xC5}) and the highest
Wλ(Lyα) quartile
(⟨Wλ(Lyα)\rangle=24\unicode{xC5}). Similarly,
galaxies typically defined as Lyman Alpha Emitters (LAEs;
Wλ(Lyα) >20\unicode{xC5}) are, on average, metal poor with
respect to the non-LAE galaxy population (Wλ(Lyα)
\leq20\unicode{xC5}) with Z⋆non−LAE≳2×Z⋆LAE. Finally, based on the best-fitting stellar models, we
estimate that the increasing strength of the stellar ionizing spectrum towards
lower Z⋆ is responsible for ≃15−25% of the observed variation
in Wλ(Lyα) across our sample, with the remaining contribution
(≃75−85%) being due to a decrease in the HI/dust covering fractions in
low Z⋆ galaxies.Comment: 10 pages, 6 figures, MNRAS accepte