Probing the Stellar Populations and Star Formation History of Early-type Galaxies at 0 < z < 1.1 in the Rest-frame Ultraviolet

We measure the evolution of the rest-frame near-ultraviolet ( NUV ) − V colors for early-type galaxies in clusters at 0 < z < 1.1 using data from the Hyper Suprime-Cam Subaru Strategic Program, CFHT Large Area U -band Deep Survey, and local Sloan Digital Sky Survey clusters observed with Galaxy Evolution Explorer. Our results show that there is an excess in the ultraviolet spectrum in most quiescent galaxies ( compared to the expectations from models ﬁ tting their optical / infrared colors and spectra ) below z ∼ 0.6, beyond which the excess UV emission fades rapidly. This evolution of the UV color is only consistent with the presence of a highly evolved, hot horizontal branch subpopulation in these galaxies ( among the majority of cool and optically bright stars ) , comprising on average 10% of the total stellar mass and forming at z > 3. The blue UV colors of early-type galaxies at low – intermediate redshifts are likely driven by this subpopulation being enriched in helium up to ∼ 44%. At z > 0.8 ( when the extra UV component has not yet appeared ) the data allow us to constrain the star formation histories of galaxies by ﬁ tting models to the evolution of their UV colors: we ﬁ nd that the epoch at which the stellar populations formed lies in the range 3 < z form < 10 ( corresponding to 0.5 – 2.2 Gyr after the Big Bang ) with a star formation e -folding timescale of τ = 0.35 – 0.7 Gyr, suggesting that these galaxies formed the majority of stars at very high redshift, with a brief yet intense burst of star formation activity. The star formation history and chemical evolution of early-type galaxies resemble those of globular clusters, albeit on much larger scales.

ETGs may also show a degree of residual star formation (RSF) from young stellar populations (Vazdekis et al. 2016;Atlee et al. 2009;Salvador-Rusiñol et al. 2020;De Propris et al. 2022) especially for low mass galaxies (Ree et al. 2012) and lenticulars (Salim et al. 2012).Young stars (age ≲ 1 Gyr) contribute strongly to the rest-frame near ultraviolet.However, RSF and the short-wavelength tail of thermal emission from old stellar populations are not the only contributors to the UV spectra of ETGs in the local universe.It has long been known that ETGs exhibit an UV 'excess' (or 'upturn ' -Code & Welch 1979;O'Connell 1999) and that this is most likely produced by blue or extreme horizontal branch (HB) stars (Yi et al. 1997(Yi et al. , 1999)), that emit efficiently in the far and near-UV spectral range.Note that blue or extreme HB stars may also affect age determinations from spectroscopic indices, via their broad Balmer lines and high luminosities (De Propris 2000; Lee et al. 2000;Percival & Salaris 2011).This requires sampling a redshift range where blue HB stars have not yet appeared or modelling the contribution to the UV light from old stellar populations if we wish to constrain star formation histories and the presence of young stars in ETGs.
In this paper we derive rest-frame N U V photometry for galaxies in clusters at 0 < z < 1.1.As shown in our previous work (Ali et al. 2021) we can use N U V photometry to constrain the strength and evolution of the UV excess component and confirm the mechanisms that create a blue HB feature in otherwise old and metal rich stellar systems (where such stars should not exist).At redshifts beyond z = 0.8 where there is no doubt, from our data, that the blue HB stars have not yet appeared in the stellar population, we can use our rest-frame N U V photometry to set stringent limits to the age, star formation history and RSF of cluster ETGs.Furthermore, the rest-frame N U V is much less affected by the age-metallicity degeneracy than the optical (Worthey 1994) and multi-wavelength datasets that include the N U V maintain their sensitivity to age and τ (the e-folding time of the star formation histories) across a wide range of masses and redshifts.The lack of evolved UV-emitting stellar populations at z > 0.8 therefore makes our conclusions more robust than can be achieved by similar studies in the local universe (e.g.Kaviraj et al. 2007).
Data from this paper come from GALEX (Morrissey et al. 2007) N U V imaging of nearby clusters of galaxies  (Sawicki et al. 2019) and g band imaging of higher redshift clusters (z = 0.85 − 1.1) from the HyperSuprimeCam Subaru Strategic Program (Aihara et al. 2018(Aihara et al. , 2022)).These are described in section 2 of this paper.The analysis and results are presented in section 3 and discussed in section 4 and 5. Conclusions are summarized in section 6. Cosmological parameters are assumed to be H 0 = 67 km s −1 Mpc −1 , Ω m = 0.3 and Ω Λ = 0.7 (Planck Collaboration et al. 2020).All magnitudes quoted are in the AB system.

DATA AND PHOTOMETRY
In order to analyse the evolution of the UV colors of galaxies, we made use of data from a number of surveys -primarily the Subaru Hyper Suprime-Cam Subaru Strategic Program (HSC SSP) and the CFHT Large Area U-band Deep Survey (CLAUDS 1 ).The HSC SSP is a grizy (plus several narrow-band) imaging survey consisting of three layers -wide, deep and ultra-deep, covering 1400, 26 and 3.5 deg 2 and reaching r-band 5σ point source depths of roughly 26, 27 and 28 respectively.The wide layer targets the North field, the Spring Equatorial field and the Fall Equatorial field.The deep layer targets the XMM-LSS, E-COSMOS, ELAIS-N1 and DEEP2-3 fields each with four pointings (apart from XMM-LSS which has 3), while the ultra-deep layer adds an extra pointing to the XMM-LSS field and a fifth overlapping pointing to E-COSMOS, dubbed 'SXDS' and 'COSMOS' respectively.For the purposes of this paper, we made use of the data from the deep fields, which already have a slew of existing multi-wavelength imaging and spectroscopy datasets in various archives.CLAUDS is a partner survey to the HSC-SSP, adding u-band data using the CFHT Megacam u and u * filters (with central wavelengths at 3538 Å and 3743 Å respectively).Similar to the HSC-SSP, CLAUDS also consists of a deep and ultra-deep layer that target the exact same fields as the aforementioned survey.Though CLAUDS coverage largely overlaps with HSC-SSP, the deep and ultra-deep fields span somewhat smaller areas of 18.6 and 1.36 deg2 , reaching u-band 5σ (2" aperture) depths of 27.1 and 27.7 respectively.The ELAIS N1 and DEEP2-3 fields were imaged with the Megacam u filter, while the XMM-LSS was imaged with the u * filter, and observations of the E-COSMOS field include both filters, which were combined in separate u and u * stacks.The clusters used in the analysis and their proposed member galaxies are drawn from the CAMIRA cluster catalog, which uses the CAMIRA red sequence finding algorithm (see Oguri 2014;Oguri et al. 2018 for details).The most up to date catalogs for the deep/ultra-deep (DUD hereafter) layers 2 are derived from Data Release 4, which uses a similar setup to the cluster catalogs from Data Release 1 (Oguri et al. 2018).The catalogs contain the photometric redshift of each cluster and when available the spectroscopic redshift of the brightest cluster galaxy (BCG).Generally for clusters with confirmed spec-z, the photo-z is within ±0.1 of the spec-z (Fig. 1) and as such we use the photo-z for the rest of the clusters when necessary for our analysis given that we are only dealing with broadband data and do not require highly precise redshift measurements.Crucially, the catalogs also provide a list of proposed member galaxies for each cluster, which are the main targets of our study.
Photometric data was taken from the internal Data Release 4 (Data Release 3 is currently publicly available).The u-band (both MegaCam u and u * ) data was taken from the HSC-joint catalog that is part of Data Release 2. For each cluster galaxy in the DUD catalog we extract the ugrizy cmodel magnitudes (a measure of the total flux), magnitudes in fixed 1.0/2.0/3.0/4.0/5.7"diameter apertures and their associated errors.
For our analysis we have chosen all clusters and their associated member galaxies between z ∼ 0.4 − 0.65 and z ∼ 0.85 − 1.1 (i.e., where the CLAUDS u/u * bands and the HSC g band match the rest-frame N U V in each redshift interval respectively).In total we have used data of 378 galaxy clusters, consisting of approximately 8000 potential member ETGs.We split the clusters up into small redshift bins of width z = 0.05 and as in Ali et al. (2021) first determine the red sequence in a rest-frame optical color -in this case observed z − y.Then from the z − y red sequence galaxies, we select a further red sequence in a rest-frame u-band color -which are observed g − y at z = 0.4 − 0.65 and observed r − y at z = 0.85 − 1.1.The red sequence in z − y generally has a width of ±0.1 and the g − y/r − y ±0.15, though this changes slightly with redshift, as shown in Figs 2 and 3.This double red sequence selection method allows us to select almost entirely a passive sample of ETGs, rejecting even most residual star-forming galaxies.Furthermore, the dual red sequence selection helps to retain the most likely cluster members in the sample given that red sequence members are found to have a 80% likelihood of being cluster members (Annunziatella et al. 2014(Annunziatella et al. , 2016;;Rozo et al. 2015;De Propris et al. 2016).Finally, for the passive sample that remained after the color cuts, we compute the rest-frame N U V − optical color using the appropriate observed band (u/u * /g) that corresponds to the rest-frame NUV and the observed y-band data.The specific bands that were used for each redshift bin and their rest-frame wavelengths are shown in Table 1.
For nearby (z < 0.15) clusters we used galaxies in the sample of De Propris et al. ( 2021) that are confirmed spectroscopic members of each cluster (using public redshifts from the NASA Extragalactic Database) within r 200 of the brightest cluster galaxy.We have collated photometry from SDSS (Alam et al. 2015) with data from GALEX (Morrissey et al. 2007) by matching positions within 3".The galaxies were split into two redshift bins -z = 0 − 0.05 and 0.05 − 0.10, then a red sequence was selected for each bin in g − r, followed by u − r as with the HSC data described above in order to select as purely a quiescent sample as possible (Meng et al. 2023).The GALEX NUV data was then used to calculate the N U V − r color.We also carried out redvisual classification on all galaxies and rejected any non-ETGs.This sample will be the subject of a future paper analyzing environmental effects and RSF in nearby clusters.

RESULTS
For the red sequence galaxies selected from the optical and u-band colors, we calculate the N U V − optical colors at each redshift bin using the rest-frame UV band as shown in Table 1 and the HSC y-band as the optical band (the exact rest-frame wavelengths are given in the table).We use fixed aperture magnitudes in all bands to calculate the colors -2.0" diameter between z = 0.4 − 0.6 redshift bins, 1.5" between z = 0.6 − 0.7 and 1.0" between z = 0.85 − 1.15.These aperture sizes correspond roughly to a fixed physical diameter of 10 kpcs across redshift.Fixed apertures are particularly useful as the UV emission is often centrally concentrated in these systems and does not follow the same surface brightness profile as the optical, due to different stellar populations being responsible for the output at each wavelength region (Carter et al. 2011;Jeong et al. 2012).As such fixing the aperture to a diameter that covers most of the UV emission and using this aperture in all other bands allows for the calculation of a more accurate color that represents the sub-population responsible for the upturn in these galaxies.
The UV excess has classically and ideally been measured using the GALEX F U V band, previous studies have shown shown that the N U V is also a suitable proxy and is largely sensitive to the phenomenon (Schombert 2016;Ali et al. 2021).In this paper we have selected galaxies which have been observed in the rest-frame N U V centred between 2250 − 2500 Å.As seen in the analysis of Chavez & Bertone (2011), the emission at < 2500 Å in ETGs is largely dominated by the stellar sub-population responsible for the upturn, with the output from the standard population of stars that makes up the bulk of ETGs and dominates the optical wavelengths > 3000 Å being nearly negligible in comparison (see also detailed explanation in Ali et al. 2021).In Ali et al. (2018a) we showed that the observed SEDs of Coma ETGs were best fit by 2-component models -one component being a standard SSP model which represents the majority main sequence/red giant branch stars and a second blackbody component superimposed on top, representing the highly evolved blue horizontal branch stars that are believed to be the driver of the UV output.It could clearly be seen that the first component had little effect below 3000 Å, while the same was true for the second component above 3000 Å.This suggests that the UV-bright stellar population in ETGs is entirely distinct from the optically-bright stellar population, as is postulated in previous studies (Ali et al. 2018a,b).The N U V region is also much more sensitive to age than the optical or N IR.As such, in the absence of an upturn population, the N U V can be used to measure the age and star-formation history of the majority stellar population in ETGs at large, more efficiently than any other wavelength regime (Kaviraj et al. 2007).We will take advantage of this key feature of the N U V to explore the star-formation history of ETGs in our analysis of higher redshift galaxies.
In Fig. 4 we show the N U V − optical color for each individual redshift bin from z = 0.4 − 1.1.At z < 0.6 we detect objects in the UV at or above 5σ up to M * +2, but at z > 0.8 we reach only up to about M * +1.5, where M * is the characteristic magnitude in the Schechter luminosity function of red sequence galaxies in the observed y-band.Virtually all red sequence galaxies in the rest-frame optical and u bands are detected to a 5σ level or higher.Of these galaxies, > 80% (the percentage being higher at lower redshift bins) also have 5σ or higher detections in the rest-frame UV, down to ∼ M * +1.5.As such the vast majority of our optically selected galaxies also have strong UV detections.In the case of non-detections, which are mostly low luminosity (> M * ) UV-red galaxies -if we extend the detection limit down to 3σ, this includes most such galaxies and their inclusion does not change in any way the overall statistics and trends discussed henceforth.However, we exclude them from our analyses for consistency.
For each redshift bin we also include the color in the observed bands from the YEPS composite stellar popu-   lation (CSP) model (Chung et al. 2013) using an infall prescription and standard cosmological He abundance (see section 4 for further details).The models are for Z=Z ⊙ , Z=2Z ⊙ and Z=0.5Z ⊙ , with a redshift of formation (z f orm ) of 3 and a delta burst.Given that our cluster data roughly encompasses magnitudes of M * ±2, the metallicity range is derived from the Coma cluster for the same magnitude range (Thomas et al. 2005(Thomas et al. , 2010)).
For instance, a typical M * galaxy in Coma has an average metallicity of about Z=1.15Z ⊙ , while M * + 2 and M * −2 roughly correspond to Z=1.75Z ⊙ and Z=0.65Z ⊙ respectively.We show a slightly broader metallicity range of 0.5 − 2Z ⊙ in order to take into account the error in the M * r calculation (∼ 0.15 mags -Beijersbergen et al. 2002) and the spread in the mass-metallicity relation.As such, the models shown in Fig. 4 have a metallicity range and a formation history typical of ETGs in our luminosity range and provide a reference point of what the (upper and lower limits of) UV colors should be without any upturn component, i.e. the colors are driven purely by the age and metallicity of the overall stellar population in the galaxy.Below z ∼ 0.6 there is a large excess of galaxies that have colors bluer than even the Z=0.5Z ⊙ model, which for most ETGs in our magnitude range is already lower than their known metallicities from optical studies in the Coma cluster as noted above.Changing z f orm to 2 (a difference of about 1 Gyr from z f orm = 3) only changes the color by about ∼ 0.1 mags blueward, which still cannot account for most galaxies.Given that a large number of galaxies have colors bluer than even the bluest model suggests an external parameter at play that cannot simply be explained by an age/metallicity effect.Conversely, in all redshift bins above z > 0.8, the vast majority of galaxies are reasonably well encapsulated by the 0.5Z ⊙ < Z <2Z ⊙ models, with most galaxies no longer exhibiting a strong UV excess.Hence the UV colors at these redshifts are mostly driven by age/metallicity.In essence there is a clear evolution in the UV color of ETGs at around z = 0.6 − 0.8, which will be explored through the lens of several models in the next section in order to understand the fundamental physical properties of these galaxies.

THE ULTRAVIOLET UPTURN
As introduced above, ETGs are known to exhibit a significant excess of UV flux compared to the extrapolation of the spectral energy distribution of their optical and infrared stellar populations (e.g., see O'Connell 1999; Yi 2010 for a review).Recent RSF may provide some of the UV flux in a few cases (Yi et al. 2005;Boselli et al. 2005), but this explanation does not account for the majority of the galaxies where UV upturns have been observed (Ferguson et al. 1991;Brown et al. 1997Brown et al. , 2000;;Boselli et al. 2005;Han et al. 2007).The broad consensus is that hot, or extreme, horizontal branch (EHB) stars must be the key contributor to the UV fluxes of ETGs (Dorman et al. 1993(Dorman et al. , 1995;;Han et al. 2007;Lisker & Han 2008).Blue HB stars are directly observed in the bulge of M32 (Brown 2004) and identified as the source of the UV light, while Rosenfield et al. (2012) argue that such stars must also provide the bulk of the unresolved UV emission in the bulge of M31.
These blue HB stars may originate from a metal-poor stellar population (Park & Lee 1997), through normal stellar evolution channels, a metal-rich stellar population with extra mass loss (Yi et al. 1997), binaries stripped of their envelope during the first ascent of the red giant branch (Han et al. 2007;Hernández-Pérez & Bruzual 2014) or helium enriched stars (Tantalo et al. 1996;Chung et al. 2011Chung et al. , 2017)).The observed evolution of the UV color in Ali et al. (2018cAli et al. ( , 2021) ) shows that the blue HB stars appear in the population at z ≲ 0.7 and this is only consistent with a model where they originate from a population of stars enriched in helium at up to 40%, with a mass fraction of typically 10%, consistent with the model of Buzzoni & González-Lópezlira (2008).Similar stars likely produce the anomalously blue HBs observed in globular clusters in the Galaxy (e.g., Lee et al. 2005;Dalessandro et al. 2013;Marino et al. 2014).
In Fig. 5 we plot the evolution of the N U V − V color with redshift of all red sequence galaxies up to ∼ M * + 1.5, such that we consistently analyse cluster ETGs in the same mass range across redshift.The rest-frame NUV data at all redshifts is centred between ∼2300-2500 Å(as seen in Table 1), so k−corrections were not large for the UV data.In the optical even though we use the same observed y-band at all redshifts, in rest-frame the effective band changes significantly in the wavelength regime, spanning rest g to i from z = 0.4−1.1.Hence larger corrections were required, though this could be achieved with reasonable accuracy using standard models given that ETGs have consistent optical colors that evolve passively over cosmic time.In order to account for the slight difference in bandpass shapes (in the UV) and rest-frame wavelengths, we followed a similar k−correction prescription from Ali et al. (2021) to correct the U V −optical colors in each redshift bin to bring them in line with the rest-frame HSC g−y at z = 1 (rest-frame N U V − V ), such that the data across all redshifts could be directly compared.To do so we made use of a standard YEPS CSP model with Z=Z ⊙ , z f orm = 3, Y ini = 0.23 and age corresponding to the redshift of the galaxy bin being corrected.The model was then shifted to the redshift of the bin and z = 1, from which the k−correction to the observed color was calculated.Note that we used a fixed age as we do not wish to apply an evolutionary correction (given that we look to measure it in the data) and only apply a k−correction (to mostly account for the bandpass differences in the observed bands at different redshifts).Ultimately once the corrections are applied, it is possible to observe the evolution of the N U V − V color over redshift and fit models to the observations.Beside the main HSC SSP and CFHT CLAUDS data that is the focus of this paper, we have also included in the plot the data from our previous analyses of cluster galaxies from the HST Frontier Fields and CLASH surveys (Ali et al. 2018c(Ali et al. , 2021)), Abell 1689 (Ali et al. 2018b), the 2dF survey (Ali et al. 2019) and the SDSS/GALEX clusters as described The N U V − V is plotted in the form of box plots where the box represents the middle 50% of the color distribution of all galaxies in the redshift bin and the notch within the box represents the median, which also roughly corresponds to the typical M * galaxy in each bin.The whiskers then show the upper and lower 25% of the color distribution, with any outliers plotted as filled circles.The box plot has the benefit of allowing us to track both the overall evolution of galaxy colors with redshift (i.e. using the median), while also keeping an eye on the scatter in the color over redshift.Alongside the observed colors we have also included a number of YEPS He-enhanced spectrophotometric models (Chung et al. 2011(Chung et al. , 2017) ) for comparison.The metallicity distribution function derived from the infall chemical evolution model proposed by Kodama & Arimoto (1997) is employed in these composite stellar population (CSP) models for ETGs.The solid lines in black, dark gray and light gray represent models with the mean metallicities of Z=Z ⊙ , Z=2Z ⊙ and Z=0.5Z ⊙ -these act as the baseline model for the standard stellar population with no Helium enhancement as in Fig. 4. As previously noted, this metallicity range was chosen as it roughly aligns with that of Coma cluster galaxies in a similar magnitude range, though the sub-solar model shown here is slightly bluer than would be expected of the lower mass galaxies in each redshift bin on average, but can account for any outliers in the lower end of The plot shows that the difference between the observed and model color is large at low redshift (∼ 0.5 mags) but become nearly zero at z > 0.8, as would be expected when the upturn has faded.This evolution is not seen in the rest-frame B − V or U − V colors, where a simple passively evolving model fits the observed data reasonably well at all redshifts sampled by our data, albeit with some stochasticity.
the luminosity function.Also plotted in shaded regions are the colors of increasingly He-enhanced models, with Y ini 3 = 0.23, 0.28, 0.33, 0.38, 0.43 (and solar metallicity).The blue envelope of the shaded region indicates the color of the He-enhanced model when 100% of the population is He-enhanced (to the model's Y ini ) and the red envelope is simply the solar metallicity model with no He-enhancement, with the shaded regions representing between 0−100% of the population being enhanced.All models have a z f orm of 3, though it should be noted that increasing z f orm to 4 (or a similar reasonable value) only moves the color redward by approximately 0.1 mags (as shown in Ali et al. 2021) and hence the results would remain almost the same, with the strength of the upturn appearing slightly stronger in comparison to the CSP.The He-enhanced models initially trace the standard solar metallicity model at high redshift when most stars in the population are still on the main sequence, as the He-enhanced population only become extremely hot in late stages of stellar evolution.However stars with increased He abundance progress through the MS and 3 The helium abundance Y of a stellar population is related to the initial helium abundance Y ini and the metallicity Z through the following equation: Y = ∆Y /∆Z× Z + Y ini , where ∆Y /∆Z is the galactic helium enrichment parameter, assumed to be 2.0.
RGB phases at a faster rate than standard stars and reach the HB phase earlier.As such models with the highest He abundance start becoming bluer at higher redshifts as their blue HBs have already become populated, with those with the lowest He abundance only start becoming blue at much lower redshift.
In order to test whether the red sequence selection criterion (as detailed in section 2 and shown in Figs. 2 and 3) in the rest-frame optical and u bands had any impact on the evolutionary trend observed in the UV, we both doubled and halved the width of the red sequence in the observed z − y and g − y/r − y bands to re-run all of our analyses using the alternative selection criteria and galaxy samples.We find that at every redshift bin altering the selection criteria only changed the median U V −optical color of all galaxies within ±0.1 magnitude and hence did not have any clear impact on the overall UV color evolution as seen in Figs. 5 and 6.Significantly broadening the red sequence width does introduce some outliers in the U V − optical both in the red and blue end of the color distribution -likely dusty red galaxies in the former case and residual star-forming galaxies in the latter.
To visualise the evolution of the N U V − V color of galaxies in comparison to conventional galaxy evolution models, in Fig. 6 we subtract median N U V − V measured for each redshift bin (as seen in the box plots in Fig. 5) from the model N U V − V color of a solar metallicity model with z f orm = 3 and standard He abundance (shown as the solid black line in the top figure), which is then plotted against redshift.As mentioned previously, this model represents the evolution of the N U V − V of a standard M * ETG (with the models being calibrated in the optical/NIR), without any extreme HB stars.As such, if the evolution of the observed color matches this standard model, the difference in the model and observed colors should remain roughly constant with redshift.To demonstrate that this is indeed the case in the optical wavelengths (where these models are known to replicate the results well), we also plot the difference in the model and observed median colors for both rest-frame U -band and optical colors (the exact rest-frame wavelengths are shown in Table 1) for all redshift bins.At a glance it can be clearly seen that the difference between the model and U -band/optical colors remains roughly constant and nearly zero between z = 0 − 1.1, suggesting that the standard solar metallicity CSP model with z f orm = 3 can reasonably replicate the median color of quiescent cluster ETGs at all redshift bins in our sample.This is of course to be expected as most galaxy evolution models are designed to primarily interpret the evolution of the optical and NIR data.
However, when we look at the N U V −V color, we find a difference of about ∼ 0.5 magnitude between the model and observed colors particularly at z < 0.6, suggesting that the galaxies are far bluer in the UV than is expected from a passively evolving stellar population (with standard He abundance).At z > 0.6 we start to see the difference become smaller gradually until eventually at z > 0.8, the observed colors match reasonably closely the prediction from the model, as is the case in the Uband and optical colors, indicating that at higher redshift the N U V − V color behaves as is expected from a quiescent, non-starforming stellar population.This unusual form of evolution in the UV color can be physically explained by invoking a sub-population of He-enhanced blue HB stars driving the UV output amidst the majority red MS/RGB stars that are traditionally known to occupy ETGs.An evolutionary pattern where the UV color remains consistently blue from z = 0 ∼ 0.6, then declines and eventually reaches its reddest by z ∼ 0.8 and beyond requires the hot HB sub-population to be He-enhanced to Y = 0.44 or higher (Y ini = 0.40 as seen in Fig. 5), formed at z f orm ∼ 3 − 4 and making up roughly 10% of the overall population.From a physical standpoint the strong UV emission is explained by the presence of the highly evolved HB stars at low redshift, which naturally starts to disappear at intermediate redshifts (starting at z ∼ 0.6 and fully by z ∼ 0.8) as there simply has not been enough time at that epoch for these stars to have evolved through the MS/RGB (with them becoming sufficiently hot only in the HB phase).Given that the UV emission is a near ubiquitous feature in the vast majority of old stellar populations at low redshift, this poses further questions on how a sub-population of stars can become He-enhanced to nearly twice the cosmological level by at least a redshift of 3, which suggests that the chemical evolution must have happened in earlier generations at even higher redshifts.
The observed evolution in Fig. 5 is inconsistent with low metallicity HB stars (Park & Lee 1997) as these would only appear at z ∼ 0.3.Such stars, if present in the amounts needed to explain the UV flux, would also dramatically weaken the strong metal lines in the observed spectra of ETGs.The results are also inconsistent with high metallicity HB stars with extra mass loss (Yi et al. 1997) as these would also appear at much later times than observed here -note also that there is no evidence for a metallicity-dependent mass loss rate in the data of Miglio et al. (2012) and McDonald & Zijlstra (2015).Binary stars (Han et al. 2007;Hernández-Pérez & Bruzual 2014) would always be present in the population even at high redshift and therefore the match to the standard N U V colors at z ∼ 1 and the onset of the UV upturn at z < 0.8 also excludes this hypothesis.In the case of binary stars, no significant evolution would have been expected in UV colors at low-intermediate redshifts after its initial onset at very high redshift, which is quite clearly not the case in our observations.These observations strongly imply that a sub-population of He-rich stars is present within most ETGs at early epochs (see also Ali et al. 2018c for a direct comparison of observed UV colors of cluster ETGs with predictions from alternative models).
Helium-rich stars are believed to contribute to the anomalously blue horizontal branch morphology of several clusters in our Galaxy (e.g., Piotto et al. 2015), the Magellanic Clouds (Mucciarelli et al. 2009;Dalessandro et al. 2016), the Sagittarius dwarf (Carretta et al. 2014) and the Fornax dwarf (Larsen et al. 2014).A large number of these GCs are considerably metal-poor and hence have lower requisite degrees of He-enhancement (observed to be up to Y ∼ 0.35) in order to produce blue HB stars compared to more metal-rich systems.However in one of the most extreme cases, NGC 2808, more than half the population is He-rich, with around 10% having He abundance of ∼ 0.43, consistent with that inferred in ETGs.Metal-rich clusters in our Galaxy such as NGC 6388 and NGC 6441 (Yoon et al. 2008;Tailo et al. 2020) have a large fraction of stars with Y > 0.38 and F U V /N U V colors similar to those of ETGs.Metalrich clusters in M31 (Schiavon et al. 2013) and M87 (Peacock et al. 2017 -the closest comparable system to our galaxies) also have F U V /N U V colors typical of local ETGs.This led Goudfrooij (2018) to propose that the UV upturn population in ETGs originates from the disruption of large numbers of metal rich clusters in the early universe.
Given their small remaining envelopes, a fraction of low mass HB stars may not undergo the traditional AGB phase, thus refereed to as 'AGB-manqué' and are the hot UV-bright progeny of the most He-enriched HB stars (Lagioia et al. 2021;Carlos et al. 2023).Abundances of light elements in Galactic globular clusters show evidence of enrichment in N and Na and depletion in C and O (Milone et al. 2014) suggestive of hot bottom burning at temperatures of 70-100 MK, consistent with deep dredge-up in massive AGB stars (Ventura et al. 2001) or fast-rotating massive stars (Decressin et al. 2007) whose Hydrogen envelope is removed.One of the AGB-manqué stars in Carlos et al. (2023) is observed to be O-poor, a hallmark of the hot bottom burning process (see also Gratton et al. 2010).Such processes, at much higher metallicities, may produce the hot HB/EHB stars that contribute to the UV upturn in ETGs.

THE STAR FORMATION HISTORY OF ETGS
Given that the upturn has almost entirely faded by z ∼ 0.8 and the observed N U V − V color matches that of a standard solar metallicity CSP model -it can be reasonably postulated that in this redshift regime the UV colors are mostly driven by the standard age/metallicity of the overall stellar population in these galaxies.The average stellar masses that we are probing can be estimated to have roughly solar metallicity, particularly for the median galaxy in our redshift bins, as this corresponds roughly to the M * characteristic magnitude in the luminosity function.We can thus use this knowledge to estimate the star-formation history of ETGs at z > 0.8 using N U V − V , which is far more sensitive to age than optical/NIR colors, as it is likely being driven by the main sequence turnoff stars in the absence of young hot stars or old He-enhanced HB stars (Lotz & Ferguson 1999).We use the data of all ETGs in redshift bins above z > 0.8 (as in Fig. 5) in order to fit models of varying z f orm (redshift of formation) and τ (e-folding timescale in Gyrs) and find the most likely combination of the two parameters which explains the star-formation history of cluster ETGs, attempting to break the degeneracy in z f orm and τ .
For the purpose of this analysis we make use of the Flexible Stellar Population Synthesis (F SP S) models from Conroy et al. (2009); Conroy & Gunn (2010) -C09 hereafter.The results from a Z=Z ⊙ C09 CSP with the power of the metallicity distribution function (pmetals) set to the default value of 2 match reasonably with the YEPS CSP infall models (Y ini = 0.23) used before, with the benefit of being able to adjust the early starformation history of the CSPs.We generated a range of CSPs with 0 < τ < 2 Gyrs (in 0.025 Gyr intervals) and 2.5 < z f orm < 13 (with the highest redshift corresponding to the earliest observed galaxy - Robertson et al. 2023).We minimize least squares to find the best fitting value of τ for each z f orm that produces UV colors matching the median observed color of galaxies in the z > 0.8 redshift bins.Of particular importance is the color of the highest redshift bin (z = 1.1) given that if z f orm is too late or τ is too long then the model will not have had enough time to evolve and become red enough to match the color of the highest redshift galaxies, which effectively places a strong limit on how long the e-folding timescale can be and how late the ETG can start forming stars based on the highly age sensitive UV data.Fig. 7 shows a sample of the best fitting models with Z=Z ⊙ at z f orm = 3, 5 & 10 (roughly 1.5 Gyrs in range) as given by the dashed red line, with the Z=2Z ⊙ and Z=0.5Z ⊙ models given by the black and blue dashed lines for comparison.For z f orm = 3, 5 & 10, the best fit τ = 0.35, 0.55 & 0.675.Fig. 8 (black line) shows the best fit τ with increasing z f orm , where generally it is found that a higher formation redshift allows for and needs longer e-folding timescales to match the observed UV data.From our results we can reasonably infer that galaxies likely did not have z f orm < 3, as at z f orm = 2.5, τ needs to be 0.1 Gyr, indicating that stars would have to have formed almost instantaneously in order for the model colors to match the observed colors.It is also not possible to push z f orm much higher than 10 as the value of τ appears to plateau at this redshift (Fig. 8) and we would start to push past the epoch of reionization.As such, our results suggest that the average M * cluster ETG realistically formed between z f orm = 3 − 10 (corresponding to about 1.5 Gyr range in time), with e-folding timescales of τ = 0.35 − 0.7 Gyr, meaning that the vast majority of stars in these galaxies formed rapidly at high redshift and the galaxies became passive thereafter -consistent with some of the high redshift quiescent galaxies now observed with JWST (Nanayakkara et al. 2022;Antwi-Danso et al. 2023;Glazebrook et al. 2023;Looser et al. 2023;Strait et al. 2023).The results are also in general agreement with optical studies of massive quiescent galaxies at intermediate to high redshifts, where such galaxies are found to have started forming their stars at z f orm = 5 − 10 ( Tacchella et al. 2022) with timescales shorter than 1 Gyr, and the peak of starformation activity potentially occurring at z ∼ 2.3 − 3.5 (Onodera et al. 2015;D'Eugenio et al. 2021).However, our calculations of the star-formation history of ETGs can be further pinpointed with deep rest-frame UV data of cluster galaxies at z = 1.1 − 2, in the redshift range where large galaxy clusters are still observed (such as in the MaDCoWs catalog; Gonzalez et al. 2019).The HSC SSP survey is not deep enough to observe such galaxies, but surveys such as EUCLID may provide the requisite rest-frame UV data to further this analysis.
In addition, the presence of a He-rich population, if formed by the ejecta of stars in previous generations, and even assuming the maximum possible efficiency where 100% of the material is recycled into stars (which is not unfeasible for galaxies that may be considered as closed boxes over the likely formation period of a few hundred Myr), implies that a very large fraction of the stellar mass must also have been present and in situ at z > 3.

Can Residual Star Formation Explain the UV
Excess?
Studies have found that the star-formation rate (SFR) in galaxies shows a steady increase in lookback time from present day out to the cosmic noon around z ∼ 2, when the SFR in the Universe peaked (e.g., see compilation  in Yüksel et al. 2008).On average the cluster galaxies we observe follow the separate and opposite trend than would be expected from even residual star-formation (RSF) being the driver of the blue UV colors.Residual star formation (Vazdekis et al. 2016;Salvador-Rusiñol et al. 2020) would need to increase with decreasing redshift to explain the increasing UV excess compared to the passively evolving main stellar population, opposite to what is observed for the evolution of SFRs in galaxies with redshift (e.g.Finke et al. 2022).Kaviraj et al. (2008) show that in their study that RSF increases with increasing redshift, from ∼ 7% at z = 0.5 to ∼ 13% at z = 1.The intensity of (recent) star formation in the most active early-types has therefore halved between Figure 9. Plot showing the evolution of the median starformation rate at each redshift in Fig. 5 calculated using the C09 CSP models (Z=Z⊙, z f = 3) with an initial burst of star-formation followed by a constant SFR as given by the fraction of total mass converted to stars in the models.This assumes that the entire UV emission observed in galaxies is attributed to constant star-formation and not the upturn.In order for star-formation to reasonably replicate the evolution of the UV output in ETGs across redshift, it would have to decrease around z = 0.6 − 0.8, which would contradict the well established expectations of the SFR increasing at higher redshifts until the cosmic noon (Yüksel et al. 2008).
z ∼ 0.7 and present day.Ree et al. (2012) find that the number fraction of ETGs with RSF is negligible (< 5%) in their brightest (M r < −22) subgroup, but gradually increases toward the fainter galaxies (∼ 30%) at z < 0.12.It is therefore unlikely that our findings can be explained by fine-tuning RSF in bright ETGs, especially in cluster environments where galaxies are likely to be more strongly quenched.
As an experiment we use the C09 CSP models to calculate the average SFR in each redshift bin.The models have a fixed Z=Z ⊙ and z f = 3, with varying levels of constant SFRs from which the best fit to the median UV color in each redshift bin is derived, hence giving us the median SFR across redshift for roughly M * galaxies.The SFR given by C09 is defined as the fraction of mass converted to stars in a constant mode of star-formation.In Fig. 9 we plot the best fit SFR to the median color of each redshift bin in Fig. 5.As expected, the SFR in these galaxies is very low -less than 1% given that most galaxies in our sample are selected to be passive from their optical colors.However, it can also be clearly seen that given the rise in the UV output at lower redshift in ETGs, the SFR increases below z ∼ 0.6 up to 0.4 − 0.5%, whereas at z > 0.6 the SFR is nearly negligible at 0.1%.Altering the Z or z f orm would change the SFR values slightly but the overall trend remains intact.Such an increase is opposite to the trend in SFR seen in key studies based on the more thoroughly analysed multi-wavelength data where the SFR is found to gradually increase with increasing redshift.It would also be difficult to physically explain why there would be a sudden jump in the SFR particularly around z ∼ 0.6 from the perspective of stellar population evolution, keeping in mind the results shown here are median values derived from hundreds to thousands of cluster galaxies at each bin and shows the general trend in quiescent ETGs at large.As such, an older stellar sub-population that evolves off the MS/RGB phase between z = 0.6 − 0.8 and turns UV-bright in the HB phase can physically explain the results in a consistent and non-ad hoc manner.While the purpose of this study is not to derive highly accurate SFR values of cluster ETGs, this experiment does illustrate that while at any given redshift bin it may be possible to assume a small fixed SFR to account for the average UV color of ETGs, it is difficult to use star-formation to explain the overall evolution of the observed UV colors in ETGs over redshift.

Effect of Environment on the UV emission & SFH of ETGs
It is well understood that the environment has a significant impact on the evolution of the morphology and stellar populations of galaxies, as is seen by the morphology-density relation (e.g.Dressler et al. 1997) and the Butcher-Oemler effect (Butcher & Oemler 1978).Clusters are found to host mostly an earlytype galaxy population, particularly at the core where star-formation is quenched to a much stronger degree than galaxies in the outskirts, and even more so than in the field (see Boselli et al. 2016 and references therein).
Given this, we were motivated to test whether the cluster size in which a galaxy resides has any significant impact on its UV color or its star-formation history derived from therein.For this analysis we make use of the cluster size information from the HSC DUD cluster catalog, in which the estimated richness of each cluster and its constituent member galaxies are provided.Thus we make a cut of all clusters with a richness above and below 20, dubbing for the purpose of this specific analysis those with richness > 20 as 'clusters' and richness < 20 as 'groups'.As would be expected, at higher redshift, particularly at z > 1 the number of clusters declines significantly, with most galaxies being part of smaller groups that have 10-20 members.We re-create Fig. 5 -the evolution of the rest-frame N U V − V color with redshift in Fig. 10, except for cluster and group galaxies separately in order to identify any difference in the color evolution between the two galaxy populations.
In general both the cluster and group galaxies show the same trend in regards to the evolution in their UV color, in that there is an excess in the UV emission below z = 0.6, which largely disappears beyond z = 0.8.This suggests, as with previous studies in local galaxy clusters (Ali et al. 2019) and groups (Phillipps et al. 2020) that the UV upturn -driven by a sub-population of blue HB stars, is an intrinsic feature of all ETGs, or old stellar populations at large, which develops given sufficient time irrespective of the external environment.This would also explain the observation of blue HB stars in similarly old GCs in the Milky Way and elsewhere (e.g., the LMC and SMC, M31, M87) as discussed previously.If He-enhancement is the primary driver of blue colors, then such anomalous chemical enhancement likely arose in these extremely old stars through some universal mechanism in the early Universe.
Furthermore, using the evolution of the UV data at higher redshift, we wanted to test whether there is any clear difference in the star-formation history parameters, i.e. age and star-formation timescales between galaxies residing in larger clusters vs. smaller groups.As such we once again used the data of all ETGs above z = 0.8 in order to fit models of varying τ for z f orm between 2.5 − 13 and find the most likely combination of the two parameters which explains the star-formation history of ETGs in both cluster and group environments.The results are shown in Fig. 11 for a few select z f orm (similar to Fig. 7).For clusters with z f orm = 3−10, τ is between 0.425 − 0.725 Gyrs.For groups with z f orm = 3 − 10, τ is between 0.325 − 0.65 Gyrs.Fig. 8 shows the change in τ with z f orm between the two sets of galaxies for comparison.On average galaxies in both clusters and groups share similar formation histories, likely having finished  most of their initial star-formation activity in very short timescales of < 0.8 Gyrs and assembled at very high redshift -beyond z = 3.The cluster galaxies do appear to on average have slightly longer e-folding timescales compared to group galaxies but this difference is minute (∼ 0.1 Gyrs).It should also be noted that at increasing redshift (particularly z > 1), the number of galaxies within large cluster environments in the HSC SSP dataset is quite low -most ETGs at these redshifts are in smaller groups.In general there does not appear to be a significant divergence in star-formation history between group and cluster environments, though this finding may change with more data of large clusters at z > 1, e.g.such as those identified in the MaDCoWS survey (Gonzalez et al. 2019).It may then be of more significance to explore the UV emission and star-formation activity of field ETGs at z > 1 to those of ETGs as part of larger structures, where the environmental differences may be more apparent.

CONCLUSIONS
From an analysis of rest-frame N U V − V colors of numerous ETGs in hundreds of clusters at 0 < z < 1.1 from multiple surveys/datasets, we conclude: • ETGs show an UV excess that has an onset epoch at about z ∼ 0.6 − 0.8 and reaches colors as observed locally within less than 1 Gyr.This can best be explained by the presence of a subpopulation of stars (about 10% by mass but with a strong dependence on galaxy stellar mass) enriched in He at about the 44% level or above, comparable to the more extreme blue HB stars observed in the globular clusters of the Milky Way.
• At z > 0.8 this component is not yet present and the colors of galaxies may be used to set stringent limits to their star formation timescales.
Our data strongly constrain the majority of star formation to have occurred before z = 3 and with e-folding timescales shorter than 0.75 Gyrs in z < 1.2 ETGs.Galaxies that are part of smaller 'groups' and larger 'clusters' appear to have similar star-formation redshifts and timescales, with only marginal differences.
• Based on the presence of He-rich stars at z < 0.6, known evolutionary timescales and yields for chemical species, we also infer that the vast majority of stellar mass in ETGs must also have been present and assembled within galaxies by z = 3.We note here that the stars may have also formed at an even earlier redshift within lower mass systems, which then merged hierarchically into the galaxies observed here, though this scenario may be less certain given the general metal-rich nature of the galaxies.
• The UV excess appears to be an intrinsic feature in ETGs and is unaffected by environment.Furthermore, residual star-formation cannot simply explain the evolutionary trend observed in the UV data, reinforcing the need of a highly evolved population of blue HB stars to be the primary driver of the UV emission in quiescent ETGs.

Figure 1 .
Figure1.Spectrosopic vs photometric redshift of cluster galaxies for which spectroscopic data is readily available in the HSC SSP catalog.

Figure 2 .
Figure2.Observed z − y (rest-frame optical as noted in brackets) color-magnitude diagrams of cluster galaxies between z = 0.4 − 1.1 from the HSC SSP survey, separated into bins of 0.05 in redshift.The red sequence galaxies are denoted with the red filled circles within the dashed lines and have photometric uncertainties of < 0.05 magnitudes in their optical colors.The dashed lines show the selection region for quiescent cluster ETGs (as described in the text, roughly within ±0.1 mag. of the mean ridge line for the red sequence).

Figure 3 .
Figure 3. Observed g − y and r − y (rest-frame u-band as noted in brackets) color-magnitude diagrams of cluster galaxies between z = 0.4 − 1.1 from the HSC SSP survey, separated into bins of 0.05 in redshift.The red sequence galaxies are denoted with the red filled circles within the dashed lines and have photometric uncertainties of < 0.1 magnitudes in their optical colors.The dashed lines show the selection region for quiescent cluster ETGs (as described in the text, roughly within ±0.15 mag. of the mean ridge line for the red sequence).

Figure 4 .
Figure 4. Observed Megacam u/u * − y and g − y (rest-frame N U V − optical as noted in brackets) color-magnitude diagrams of cluster galaxies between z = 0.4 − 1.1, separated into bins of 0.05 in redshift.The red sequence galaxies are those selected from the optical colors in Fig. 2 and 3. Photometric uncertainties are < 0.15 magnitudes.

Figure 5 .
Figure 5. Evolution of the rest-frame N U V − V (observed g − y at z = 1) color over redshift/lookback time as given by the YEPS spectrophotometric (infall) models for a range of initial helium abundances -Yini = 0.28, 0.33, 0.38, 0.43 (shaded areas), with z f orm = 3 and metallicities as detailed in the figure legends.Also included is the evolution of the same color for infall models with Yini = 0.23 (i.e. standard cosmological He abundance with no upturn) for Z=Z⊙, 0.5Z⊙ and 2Z⊙ (solid lines).Plotted on top are box plots which show the rest-frame N U V − V colors of cluster galaxies in bins between z = 0 − 1.1, with different colors indicating the datasets from which the galaxies are derived, as labelled in the figure.Photometric uncertainties in color are < 0.2 magnitudes. in Section 2 (De Propris et al. 2021).The datasets are labelled with different colors in the figure.The N U V − V is plotted in the form of box plots where the box represents the middle 50% of the color distribution of all galaxies in the redshift bin and the notch within the box represents the median, which also roughly corresponds to the typical M * galaxy in each bin.The whiskers then show the upper and lower 25% of the color distribution, with any outliers plotted as filled circles.The box plot has the benefit of allowing us to track both the overall evolution of galaxy colors with redshift (i.e. using the median), while also keeping an eye on the scatter in the color over redshift.Alongside the observed colors we have also included a number of YEPS He-enhanced spectrophotometric models(Chung

Figure 6 .
Figure6.The difference between the median N U V − V color of each redshift bin (white notches in each box plot in the Fig.5) subtracted from the N U V −V color of the Yini = 0.23 solar metallicity infall model formed at z f orm = 3 (black line in Fig.5, i.e. standard passively evolving model with no upturn).The plot shows that the difference between the observed and model color is large at low redshift (∼ 0.5 mags) but become nearly zero at z > 0.8, as would be expected when the upturn has faded.This evolution is not seen in the rest-frame B − V or U − V colors, where a simple passively evolving model fits the observed data reasonably well at all redshifts sampled by our data, albeit with some stochasticity.

Figure 7 .
Figure7.Evolution of the rest-frame N U V − V (observed g − y at z = 1) color of cluster ETGs over redshift/lookback time as in Fig.5, but only for z > 0.8 where the upturn is expected to no longer be present.Also plotted for comparison are C09 CSP models for z f orm = 3, 5, 13 (covering roughly 1.5 Gyrs in range) and varying e-folding timescales (τ ) that best fit the data for each formation redshift.The red dashed lines show the closest fit to the median at all redshift bins in each case (the red envelope indicates the uncertainty due to photometric error), with the black and blue lines showing the Z=2Z⊙ and Z=0.5Z⊙ models respectively.The results suggest that ETGs have z f orm = 3 − 10 and τ = 0.35 − 0.7 Gyrs respectively.

Figure 8 .
Figure 8.The star-formation e-folding timescale (for solar metallicity C09 CSP models) that would be required at each formation redshift to fit the rest-frame N U V − V colors of cluster galaxies between z = 0.85 − 1.1.Several examples of the model fits to the data are shown in Fig. 7.

Table 1 .
Table giving details of the observed bands and rest-frame central wavelengths of HSC clusters, divided into redshift bins of 0.05.