The ionizing properties of two bright Ly$\alpha$ emitters in the BDF reionized bubble at z=7

We investigate the ionizing properties of the pair of bright Ly$\alpha$ emitting galaxies BDF521 and BDF2195 at z=7.012 in order to constrain their contribution to the formation of the BDF"reionized bubble". We obtain constraints on UV emission lines (CIV$\lambda 1548$ doublet, HeII$\lambda 1640$, OIII]$\lambda 1660$ doublet, and CIII]$\lambda 1909$ doublet) from deep VLT-XSHOOTER observations and compare them to those available for other high-redshift objects, and to models with mixed stellar and AGN emission. We use this spectroscopic information together with the photometry available in the field to constrain the physical properties of the two objects using the spectro-photometric fitting code BEAGLE. We do not detect any significant emission at the expected position of UV lines, with 3$\sigma$ upper limits of EW$\lesssim$2-7 AA rest-frame. We find that the two objects have lower CIII] emission than expected on the basis of the correlation between the Ly$\alpha$ and CIII] EWs. The EW limits on CIV and HeII emission exclude pure AGN templates at $\sim2-3\sigma$ significance, and only models with a $\lesssim$40% AGN contribution are compatible with the observations. The two objects are found to be relatively young ($\sim$20-30 Myrs) and metal-poor ($\lesssim 0.3 Z_{\odot}$) with stellar masses of a few $10^9M_{\odot}$. Their production rate of hydrogen ionizing photons per intrinsic UV luminosity is log($\xi_{ion}^*$/Hz erg$^{-1}$)=25.02-25.26, consistent with values typically found in high-redshift galaxies, but more than twice lower than values measured in $z>$7 galaxies with strong CIII] and/or optical line emission ($\simeq$25.6-25.7). The two BDF emitters have no evidence of higher than average ionizing capabilities and are not capable of reionizing their surroundings by their own means under realistic assumptions on the escape fraction of ionizing photons. (abridged)


Introduction
The epoch of reionization (EoR) marked a major phase transition of the Universe, during which the intergalactic medium (IGM) became transparent to UV photons.Determining the physical processes involved in the reionization process, its timeline and topology represents the latest frontier in observational cosmology (Dayal & Ferrara 2018).
The first indication of an increased neutral hydrogen fraction (χ HI ) in the IGM at z 6 was obtained by observations of the Gunn-Peterson effect in distant quasars (Fan et al. 2002).However, a substantial step in our knowledge of the reionization timeline has been made possible only by constraints on the redshift evolution of the fraction of Lyman-break galaxies (LBGs) that have an appreciable Lyα emission line (e.g., Malhotra & Rhoads 2006;Stark et al. 2010;Fontana et al. 2010).In fact, the neutral hydrogen in the IGM scatters Lyα photons out of the line-ofsight, such that a lower line visibility indicates a higher IGM neutral fraction if the properties of the inter-stellar medium remain unchanged (e.g., Dayal et al. 2011;Dijkstra 2019).As of today, a decrease in the Lyα fraction at z 6 has been confirmed by many independent analyses and interpreted as indication of an increasing neutral hydrogen fraction in the IGM (e.g., Pentericci et al. 2011;Schenker et al. 2012;Ono et al. 2012;Mason et al. 2018;Pentericci et al. 2018;Mason et al. 2019).Together with measurements of the CMB Thomson optical depth (Planck Collaboration et al. 2020), of the clustering of Lyαemitters (LAEs, Hutter et al. 2015;Sobacchi & Mesinger 2015) and with measurements of ionized regions around the most distant QSOs (e.g., Bañados et al. 2018), the available observations point to a scenario in which the universe was still substantially neutral at z ∼10 (χ HI 90%) and was rapidly reionized in ∼500Myrs, ending the EoR at z ∼5.5-6 (e.g., Mitra et al. 2015;Greig & Mesinger 2017;Mitra et al. 2018;Hutter et al. 2021).Star-forming galaxies are currently considered as the most likely responsible of reionization (e.g., Robertson et al. 2015;Bouwens et al. 2015a;Dayal et al. 2020;Endsley et al. 2021b;Romanello et al. 2021), although our poor knowledge of galaxy physical properties does not allow us to rule out that active galactic nuclei (AGN) also contributed (Finkelstein et al. 2019;Giallongo et al. 2015).
The rate of ionizing photons escaping into the IGM from a given population is Ṅ = ρ UV ξ * ion f esc , where ρ UV is the total, dustcorrected UV luminosity density at 1500Å, ξ * ion is the ionizing photon production efficiency per unit UV luminosity, and f esc is the fraction of ionizing photons leaked into the surrounding environment.Currently, the UV luminosity density from starforming galaxies is well constrained by measurements of their UV luminosity function (LF) at z >6 down to L < L * (e.g., Bouwens et al. 2015b;Livermore et al. 2017;Ishigaki et al. 2018;Yue et al. 2018;Oesch et al. 2018).On the other hand, the escape fraction of ionizing photons can only be directly constrained at z 3-4 and on average is quite modest in LBGs with values lower than 5-10% (Boutsia et al. 2011;Marchi et al. 2017;Grazian et al. 2017;Steidel et al. 2018;Pahl et al. 2021).Only in few rare objects it reaches values of f esc 50% (de Barros et al. 2016;Vanzella et al. 2016;Naidu et al. 2017;Vanzella et al. 2018;Izotov et al. 2018a,b;Flury et al. 2022).As of today, no real constraints exist on a possible increase of f esc with redshift and we can only make assumptions on its value in the EoR.
The ionizing efficiency ξ * ion can be directly obtained by measuring, at the same time, the non-ionizing UV continuum and the Balmer emission lines, the latter yielding the emission rate of ionizing photons (e.g., Shivaei et al. 2018), or from the equivalent width of the [OIII]λ4959, 5007 doublet (Chevallard et al. 2018;Tang et al. 2019).At higher redshifts where optical emission lines are not directly observable through spectroscopy, constraints on ξ * ion were derived by the analysis of the rest-frame UV colors (e.g., Duncan & Conselice 2015) or from the contamination from strong emission lines to the mid-IR colors, finding a typical log(ξ * ion /Hz erg −1 ) 25.3 at z ∼4-5 (Bouwens et al. 2016;Lam et al. 2019), with indication of higher ionizing efficiencies in LAEs (Harikane et al. 2018;Sobral & Matthee 2019) and strong [OIII]λ4959, 5007 emitters (Tang et al. 2019).
For galaxies in the EoR, constraints on the ionizing capabilities can be derived using UV emission lines.These however are extremely faint and have been detected only in a handful of objects (e.g., Stark et al. 2015a,b;Laporte et al. 2017;Mainali et al. 2018;Hutchison et al. 2019;Topping et al. 2021).
Such features might be associated to AGN or metal poor stellar populations (Gutkin et al. 2016;Feltre et al. 2016;Nakajima et al. 2018b) as also shown at lower redshifts (e.g., Amorín et al. 2017;Calabrò et al. 2017;Nakajima et al. 2018a).The spectro-photometric analysis of a CIV emitter at z=7.045 allowed Stark et al. (2015b) to estimate a ionizing efficiency log(ξ * ion /Hz erg −1 ) 25.68 much higher than the typical value measured at lower redshifts, and comparable to values measured in rare, extreme emission line galaxies at z ∼3-4 (Nakajima et al. 2016).Similarly, spectroscopically confirmed z >7 galaxies with mid-IR colors suggestive of intense Hβ+[OIII] emission show a log(ξ * ion /Hz erg −1 ) 25.5 (Stark et al. 2017;Endsley et al. 2021a).These results point to an increased fraction of galaxies with high ionizing efficiency in the EoR.
A promising way to shed light on the sources of reionization is the investigation of regions where this process is more advanced and observations point to the presence of "reionized bubbles" (Castellano et al. 2016a;Higuchi et al. 2019;Tilvi et al. 2020;Jung et al. 2021;Leonova et al. 2021;Endsley & Stark 2021).The analysis of the Lyα fraction in independent lines-ofsight suggests that reionization was a spatially inhomogeneous process (Treu et al. 2012;Pentericci et al. 2014).A detailed spectroscopic investigation of the Bremer Deep Field (BDF) (Lehnert & Bremer 2003;Castellano et al. 2010), led to the first discovery of a likely reionized "bubble" at z ∼7 (Castellano et al. 2016a, hereafter C16).The BDF hosts three bright (L ∼ L * ) Lyα emitting galaxies with EW>50Å (Vanzella et al. 2011;Castellano et al. 2018, hereafter V11 and C18), and the density of z 7 faint LBGs in this area is 3-4 times higher than the average (C16).Two of the Lyα emitting galaxies (BDF521 and BDF2195) have exactly the same redshift and are at a projected physical separation of only 91.3kpc, the third one, BDF3299, being at 1.9 proper Mpc projected distance.
The high number density of LBGs in the BDF region is consistent with the clustered faint galaxies being key contributors to the local reionization history (C16, Rodríguez Espinosa et al. 2021), as expected in "inside-out" reionization scenarios (Choudhury et al. 2009;Dayal et al. 2009;Trebitsch et al. 2021).However, significant uncertainties remain due to the lack of Lyα detections in the faint companion galaxies, despite the significant observational efforts (C18).Lacking a spectroscopic confirmation of the faint LBGs in the field, it is fundamental to fully constrain the physical properties of the bright emitters and assess whether they can create the reionized region possibly thanks to hard-ionizing stellar populations or AGN.
In this paper we analyse deep VLT-XSHOOTER observations of the BDF521 and BDF2195 pair to put constraints on the ionizing budget of the two emitters from UV rest-frame metal emission lines, and to ascertain whether these L L * galaxies play a major role in the creation of the BDF bubble.
The paper is organised as follows: in Sect. 2 we present the observations and data reduction, in Sect. 3 we discuss the constraints on emission lines and compare the BDF galaxies to other high redshift sources.We place constraints on the contribution from AGN in Sect.4, and quantify their physical properties and their contribution to the creation of the BDF bubble in Sect. 5.The results are summarised in Sect.6.

The XSHOOTER observations
The BDF521 and BDF2195 pair was observed by XSHOOTER in nodding mode with 11×0.9 arcsec slits and 900 and 865 seconds per single exposure in the NIR and VIS arms, respectively.Observations were acquired between 26 September 2019 and 8 June 2021.The total observing time was of 15 and 12 hours for BDF521 and BDF2195, respectively, corresponding to on-target exposure times of 12.0 (11.5) and 9.5 (9.1) hours in NIR (VIS).
Science frames were reduced with the official XSHOOTER pipeline (Modigliani et al. 2010), v. 3.3.5,using the associated raw calibrations from the ESO archive.Sky subtraction was performed using the XSHOOTER nodding strat-egy, combining frames acquired at two different positions in the sky for each OB execution.A fixed boxcar window was used to extract the target on the rectified 2D spectra, to ensure that the target was correctly localized along the slit (extractmethod=LOCALIZATION and localize-method=MANUAL in pipeline recipe xsh_scired_slit_nod).The position of the window with respect to the slit center and its half height (in arcsec) were fine-tuned for the two targets independently, based on the position of the Lyα line which is well detected in the 2D VIS median stacked spectra (and visible in each single OB); we used localize-slit-position=-1.0, localize-slit-hheight=0.5 for BDF521 and localize--slit-position=-3.2, localize-slit-hheight=0.6 for BDF2195 in pipeline recipe xsh_scired_slit_nod.
The 1D extracted spectra, one for each XSHOOTER arm and OB execution, were corrected to the barycentric reference frame in vacuum and combined with the data analysis package Astrocook (Cupani et al. 2020).The combined spectra were rebinned to 18 km s −1 in the NIR arm and to 11 km s −1 in the VIS arm, roughly corresponding to 3 pixels per resolution element, given the nominal resolution FWHMs of the instrument in the three arms for the adopted slit.Undersampled spectra with larger spectral bins (54 km s −1 in the NIR arm, 33 km s −1 in the VIS arm) were also produced to facilitate the identification of spectral features by visual inspection.
While the slits are not designed to include the total flux from the observed objects, following Lemaux et al. (2009) we expect slit losses to be small (∼10-15%) on the basis of the slit dimension, the compact size of the sources (∼0.15") and the median seeing (∼0.8") of the XSHOOTER observations.In order to quantitatively assess the significance of potential slit losses we smoothed the HST F125W image (C16) of the objects to the seeing of the spectroscopic observations and measured the fraction of the total flux observed within the slit size.We find corrections of ∼10% that are, however, uncertain due to effects of light contamination from other objects and local background subtraction of the images.Therefore, we decided not to apply correction factors in the following.

Constraints on UV rest-frame line emission
At the Lyα redshift estimated from FORS2 (z=7.008,C18), the VIS arm of the instrument covers the Lyα line, while the NIR arm covers the range of the CIV doublet (1548Å-1550Å), HeIIλ1640, the OIII] doublet (1661Å-1666Å), and the CIII] doublet (1907Å-1909Å).We first visually inspected the spectra to localize Lyα emission and obtain a precise positioning of the object within the slit.

Lyα emission
We detect Lyα line emission for both objects: the flux-weighted mean wavelength of the line corresponds to z=7.0121 and z=7.0124 for BDF521 and BDF2195, respectively.These redshifts are slightly higher than previously estimated from the peak of the Lyα emission in FORS2 spectra (z=7.008)but do confirm that the objects are approximately at the same redshift and, given the close separation on the sky (91.3kpc,C18), likely physically connected.We will use these new estimates in the present analysis.The total Lyα fluxes are 12.8 ± 1.0 and 21.8 ± 1.1 10 −18 erg/s/cm 2 for BDF521 and BDF2195, respectively.The widths in velocity space and the equivalent widths of the lines are consistent with previous analysis.We find EW(Lyα)=60±5Å,  Lenz & Ayres (1992).
The VIS arm covers also the expected range of the NVλ 1240 doublet but in a noisy, low efficiency, region of the spectrum that does not allow us to improve upon the limits previously obtained with FORS2 (C18).

Limits on UV metal lines
We inspected both 2D and 1D spectra to assess the presence of metal UV emission lines.We detect no continuum emission and find no significant features within 500 km/s from the position expected on the basis of the Lyα redshift.In general, no evident line emission is elsewhere found at the objects' positions in the NIR arm spectra.Considering the Lyα velocity offset range expected for bright z ∼7 galaxies (e.g.Maiolino et al. 2015;Mason et al. 2018) we estimate upper limits for the flux of all lines from the average flux uncertainty in a moving window of 100 km/s rest-frame across the range between -500 km/s and 0 km/s from the Lyα redshift, while masking sky line residuals.We derive corresponding EW limits after estimating the continuum flux at each position on the basis of the observed photometry (C16).Namely, we follow the standard assumption of approximating the UV continuum as a power law F λ ∝ λ β whose slope β can be measured from a regression across the observed bands (e.g.Castellano et al. 2012).The measured slopes are β =-2.25±0.25 and β =-1.66±0.16for BDF521 and BDF2195, respectively.We find rest-frame EW limits 2-7 Å at 3σ (Table 1).The limit on HeII emission for BDF521 is consistent with the estimate obtained by Cai et al. (2015) from narrow-band HST photometry.The spectra at the expected positions of the relevant lines are shown in Figs.A.1 and A.2.

Comparison with other high-redshift sources
In the last few years, several groups actively searched for UV emission lines in the most distant galaxies, as valuable diagnostics of the underlying radiation field (e.g., Mainali et al. 2018;Stark et al. 2017;Laporte et al. 2017;Endsley et al. 2021a) and in some case as redshift indicators, in the absence of Lyα (Jiang et al. 2021b).The most common line, besides Lyα, is the CIII] doublet which has been detected in several of the most distant galaxies including the most distant spectroscopically confirmed galaxy GN-z11 at z=10.957 (Oesch et al. 2016;Jiang  It is unclear whether the relation holds at z 6 where only a few detections are available.The reported detections we show in the left panel of Fig. 1 seem to agree with the lower redshift behavior, while our two BDF emitters do not follow the expected trend.Indeed, the galaxies have rather large Lyα EW and we would expect solid detections of the CIII] doublet at > 5σ from the average of lower redshift relations, albeit the BDF objects are consistent with the lower envelope of the measured distributions.As shown in the right panel of Fig. 1, all high-redshift objects with a secure CIII] detection are bright, L L * LBGs (i.e., M UV -20.5, e.g.Harikane et al. 2021), in most cases significantly brighter than the BDF emitters.Llerena et al. (2021) found that strong CIII] emitters at z ∼3 have strong radiation fields, low stellar metallicity but their gas is already partially enriched with C/O abundances in the range 35%-150% solar.In this respect, the lack of strong CIII] emission can be explained either by an extremely low C abundance in the gas phase, and/or by the objects being in a less active starformation phase (Berg et al. 2019).In the case of the BDF galaxies a low metallicity is likely the main factor affecting the low CIII] emission given the high Lyα EW suggesting active starformation.A low metallicity is also consistent with the lack of CII158µm detection in BDF521 (Maiolino et al. 2015).A strong starburst phase that photodissociates the molecular clouds could also lower CII158µm emission (Pallottini et al. 2019), but UV CIII] emission should be enhanced in such a situation if C is abundant in the ionized gas.In Fig. 2 we compare the limits we have obtained for the CIV emission on the two BDF galaxies to previous studies.In the EoR there are only two galaxies where a convincing CIV emission has been detected, the gravitationally lensed low-mass galaxies RXC J2248.7-4431-ID3(z=6.11,Mainali et al. 2017), and A1703-zd6 (z=7.045,Stark et al. 2015b).Several attempts to secure this line in small samples of galaxies at z=6-10 have failed to report detections (Laporte et al. 2017;Mainali et al. 2018;Shibuya et al. 2018), with typical EW limits around 3-10 Å comparable to the BDF galaxies.Although an AGN contribution to the Lyα flux cannot be completely ruled out, the CIV equivalent widths and f CIV / f Lyα ratios for the two BDF galaxies are low compared to most known Type II AGN at intermediate redshift (Alexandroff et al. 2013), and much lower than those for the two z >6 LBGs with CIV detections (left panel of Fig. 2).These results are consistent with the analysis of the Lyα/NV ratio from FORS2 spectroscopy (C18).The two BDF galaxies have Lyα/NV 8-10 which is higher than the values measured in suspected AGN at z 7 (e.g., Tilvi et al. 2016;Laporte et al. 2017;Mainali et al. 2018), albeit still consistent with values observed in AGN at lower redshifts (Humphrey et al. 2008;Hainline et al. 2011).

Limits on the contribution from AGN emission
The analysis of the limits on CIV emission presented in Fig. 2, and the existing constraints on the Lyα/NV ratio (C18), already suggest that it is unlikely that the BDF emitters are dominated by AGN emission.We perform here a more quantitative analysis by comparing the limits on CIV and HeII emission against narrow line AGN models by Feltre et al. (2016).We consider models with -5.0 ≤ log(U s ) ≤ -1.0 and interstellar gas metallicity Z = 0.025, 0.1, 0.4 Z , dust-to-metal mass ratio ξ d = 0.3, and hydrogen gas density n H =100 cm −3 .We restrict our analysis to models generated with incident spectra having UV spectral index α = −1.7,but we verified that results do not significantly change when adopting spectra with index α = −1.4.We first re-normalize the incident spectra to a chosen fraction of the observed flux of each object in the F125W filter (C16), i.e., the observable non-ionizing UV flux at ∼ 1500Å, and then compute the EW of the CIV and HeII lines on the basis of the predicted line flux for each model and of the observed flux at the relevant wavelength.We note that using the models in this way allows us to constrain scenarios in which the AGN is contributing to the UV continuum emission, which is the most important test to constrain the ionizing output of our targets.While the models do not include the broad line region emission, this is not a concern in our case given the narrow FWHM of the Lyα lines observed in the BDF galaxies.We show in Fig. 3 the position on the CIV(EW) versus HeII(EW) plane of AGN models with a 40%, 60% and 100% contribution to the observed UV emission highlighting the regions allowed at 1σ, 2σ and 3σ by the XSHOOTER observations.Models with pure AGN (100%) are always excluded at 2 − 3σ regardless of log(U s ) and metallicity.The more stringent constraints on BDF2195 allow us also to exclude a 60% AGN contribution at ∼3σ.Models with 40% contribution are more compatible with the observations, particularly for BDF521, for high log(U s ) and intermediate/high metallicities.
The narrow line models of Feltre et al. (2016) can be further used to investigate the case of an obscured AGN providing no contribution to the UV continuum emission.In such a case, we must consider the fractional contribution of the AGN to the observed Lyα emission.Since the radiative transfer effects of Lyα are severe, we therefore investigate what fraction of Lyα can escape for given fractional contribution of AGN without invalidating the measured line limits.For a similar parameter space of the models investigated above, we find that we can only rule out 100% narrow-line region contribution to Lyα if only 20% (15%) escapes for BDF2195 (BDF521).For lower escape fraction of Lyα, the line contribution can have 100% contribution from a Type-II AGN without the other rest-frame UV lines being detected.

The ionizing emission budget of the BDF emitters
From the above analysis it appears that the BDF pair lacks highionization features found in some other sources at a similar redshift, and instead resembles standard star-forming galaxies at lower redshifts with low or absent contribution from AGN.We thus constrain the physical parameters of the two galaxies, including the ionizing budget, by means of a spectro-photometric fit performed with the BEAGLE v0.24.5 tool (Chevallard & Charlot 2016) using the most recent version of the Bruzual & Charlot (2003) stellar population synthesis models (see Vidal-García et al. 2017, for details).Nebular emission is modeled self-consistently as described in Gutkin et al. (2016) by processing stellar emission with the photoionization code CLOUDY (c13.03, Ferland et al. 2013).The fit is performed by fitting the upper limits on the integrated lines plus continuum fluxes measured as described in Sect.3, together with the broad-band photometric measurements redward of the Lyman break in the F105W, F125W, F160W and HAWK-I Ks bands (C16, see also Castellano et al. 2010;Cai et al. 2015).Given that the Lyα is very difficult to model in this framework, we subtract its flux from the affected bands and do not include it in the nebular emission mod- eled by BEAGLE.We fix the redshift at the value determined from Lyα.Although there is usually a shift between Lyα and systemic redshift, this in practice has no effect on the constraints derived here.The BEAGLE SED-fitting runs are performed with a configuration similar to the one discussed by Stark et al. (2017).The templates are based on a Chabrier (2003) initial mass function and have metallicity in the range −2.2 ≤ log(Z/Z ) ≤ 0.25.The star formation histories (SFH) are parametrized as an exponentially delayed function (SFR(t) ∝ t•exp(-t/τ)), which is the most flexible parametric SFH allowed by the code, plus an ongoing constant burst of 10 Myr duration.We adopt uniform priors on the SFH exponential timescale (7.0 ≤ τ/log(yr) ≤ 10.5), stellar mass (7.0 ≤ log(M/M ) ≤ 12), and maximum stellar age (7.0 ≤ log(Age/yr) ≤ age of the universe).Attenuation by dust is treated following the Charlot & Fall (2000) model combined with the Chevallard et al. (2013) prescriptions for geometry and inclination effects, assuming an effective V-band optical depth in the range -3.0 ≤ log(τ V ) ≤ 0.7 with a fixed fraction µ = 0.4 arising from dust in the diffuse ISM.Interstellar metallicity Z IS M is assumed to be identical to the stellar one, the dust-to-metal  mass ratio and ionization parameter are left free in the ranges 0.1 ≤ ξ d ≤ 0.5 and -4.0 ≤ log(U s ) ≤ -1.0, respectively.We exploit the observed Lyα flux to impose a prior on the ongoing SFR.Following Kennicutt (1998), the total line fluxes correspond to SFR(Lyα)=6.7±0.5 M /yr and SFR(Lyα)=11.4±0.6 M /yr for BDF521 and BDF2195, respectively.As we expect the SFR inferred from the detected Lyα to be a lower limit to the intrinsic SFR of BDF521 and BDF2195 depending on their Lyα escape fraction and attenuation by the IGM, we use the corresponding 5σ lower limits on SFR(Lyα), i.e. 4 M /yr and 8 M /yr, as conservative lower bounds for the uniform prior ranges on their current SFRs.
The best-fit parameters and relevant 68% confidence level uncertainties are shown in Table .2. The objects are found to be relatively young (mass-weighted age ∼20-30 Myrs) and metalpoor ( 0.3Z ) galaxies with stellar masses of a few 10 9 M .Both objects are actively star-forming (SFR∼15M /yr), but had experienced a higher star-formation activity in the past.We find a small degeneracy among the mass, extinction and age parameters, with older models being more massive and less extincted by dust.The results do not change significantly when adopting different parametrizations of the SFH, nor with updated treatment of dust within HII regions implemented in BEAGLE v0.27.1 described in Curtis-Lake et al. (2021).
The most important parameter to explore thanks to the combined spectroscopic and photometric information is the production rate of hydrogen ionizing photons per intrinsic (i.e.unattenuated) UV luminosity, ξ * ion .We find best-fit values for log(ξ * ion /Hz erg −1 ) of 25.26 and 25.02 for BDF521 and BDF2195, respectively.These values are consistent with the typical ionizing production efficiency of z 4 populations (Bouwens et al. 2015a(Bouwens et al. , 2016;;Lam et al. 2019), and significantly lower than the log(ξ * ion /Hz erg −1 ) 25.6-25.7 found in z >7 galaxies with CIII] detection and/or photometric evidence of strong optical emission lines (Stark et al. 2015b(Stark et al. , 2017)).The ionizing emission from the BDF emitters is lower, albeit consistent at ∼ 1σ, with the average value in z ∼5 LAEs (Harikane et al. 2018), and with the log(ξ * ion /Hz erg −1 )=25.43 estimated for the BDF galaxies by Rodríguez Espinosa et al. (2021) on the basis of the Lyα luminosity following Sobral & Matthee (2019).Most importantly, these results show that the two emitters are rather "standard" high-redshift star-forming galaxies, with no significant excess production of ionizing photons contributing to the local reionization history.
With the improved constraints coming from the new spectrophotometric fit described in the previous section, we can quantitatively explore whether the BDF pair is capable of re-ionizing its surroundings or not.The SFR and the age of a galaxy can be used to measure the size (R bubble ) of the resulting ionized region following, e.g., Shapiro & Giroux (1987) and Madau et al. (1999) by assuming an escape fraction of ionizing photons f esc , the hydrogen clumping factor C, and an average neutral hydrogen fraction χ HI at the onset of star-formation.As discussed in V11 and C18, the size R bubble must then be compared to the minimum size R min =1.1 Mpc enabling Lyα to be redshifted enough to reach the observers (Wyithe & Loeb 2005).The aforementioned R min is found under the assumption that the Lyα escapes from the galaxies at the systemic redshift.The size threshold decreases in presence of strong outflows, for example a 220km/s shift, which is the median value found by Mason et al. (2018) (see also Endsley et al. 2022) for galaxies in massive halos, results in R min ∼ 0.85 Mpc.
In practice, we compute the average SFR of the pair from all models within 95% c.l. from the best-fit as a function of time in bins of 5 Myr.We then convert this into R bubble by assuming C=2 and χ HI = 0.5 and let the escape fraction f esc vary from 10% to 80%.The assumption of a 50% neutral fraction is conservative as a higher χ HI would further decrease the size of the reionized region created by the sources.The results are shown in Fig. 4 as a function of the relevant maximum stellar age of the models (i.e., the onset of the star-formation episode): we find that in the absence of significant outflows only models with an extremely high f esc =80% and age larger than 200 Myr (formation redshift z >9) can create a large enough bubble.
When assuming a Lyα shift imprinted by 220km/s outflows, the two objects could grow a large enough bubble with an escape fraction f esc =40% on the same time frame.However, even in the presence of outflows, there are no spectro-photometric fitting solutions allowing the BDF pair to reionize the surrounding region with an escape fraction of f esc =10%.This is particularly rele-vant when considering that the physical properties of the BDF galaxies, i.e. moderate ξ * ion , and low EW of the UV metal lines, are typical for galaxies with low f esc (Naidu et al. 2022).Note that these results are considerably more stringent than those discussed in C18 based on the SED-fitting of the available photometry only.The C18 results were consistent with a combination of high SFRs and old ages which in turn allowed the formation of a large enough reionized bubble.Our current spectro-photometric analysis suggests that the two emitters are not solely responsible of the growth of the ionized bubble.In turn, this finding is consistent with the analysis of the BDF region presented by Rodríguez Espinosa et al. (2021) indicating that a dominant contribution is provided by the clustered, faint galaxies.

Summary and future prospects
We have presented the analysis of deep VLT-XSHOOTER observations of BDF521 and BDF2195, a pair of L ∼ L * , bright Lyα emitting galaxies at ∼90 kpc separation shining within the BDF bubble, a z ∼7 overdensity of both LBGs and LAEs showing all expected properties of a "reionized bubble" embedded in a half-neutral universe (V11, C16, C18, Rodríguez Espinosa et al. 2021).Our main findings can be summarised as follows: -The Lyα lines are detected at high significance in the VIS arm spectra at z=7.0121 and z-7.0124 for BDF521 and BDF2195, respectively.The Lyα properties (EW 60Å and FWHM 250-260km/s) are consistent with previous analysis based on lower resolution FORS2 spectra.-No significant emission is detected in the NIR arm spectra at the expected position of the CIVλ1548 doublet, HeIIλ1640, OIII]λ1660 doublet, and CIII]λ1909 doublet.We can place stringent 3σ upper limits on the EW of each line or doublet component in the range ∼2-7Å (Table 1).
-A comparison with other high-redshift sources from the literature shows that our two BDF emitters have lower CIII] emission than expected on the basis of the average correlation between the Lyα and CIII] EWs measured at lower redshifts, albeit they are consistent with being at the lower envelope of the measured distributions (Fig. 1).-Although an AGN contribution cannot be completely ruled out, the constraints on the CIV(EW) and the flux ratio f CIV / f Lyα (Fig. 2) show that the objects are only marginally consistent with them being AGN when compared to known Type II QSOs at lower redshifts.Consistently, when comparing limits on HeII(EW) and CIV(EW) to narrow line AGN models by Feltre et al. (2016) we can exclude pure AGN (100%) templates at ∼ 2 − 3σ.Instead, models with a 40% AGN contribution could be compatible with the observations (Fig. 3).-A spectro-photometric fit on the available data indicate that the two objects are relatively young (∼20-30 Myrs) and metal-poor ( 0.3Z ) with stellar masses of a few 10 9 M (Table 2).Most importantly, we find a production rate of hydrogen ionizing photons per intrinsic UV luminosity of log(ξ * ion /Hz erg −1 )=25.02-25.26(under the assumption of a purely stellar origin), which is consistent with the value typically found in high-redshift populations, but significantly lower than the extreme values ( 25.6-25.7)measured in z >7 galaxies with CIII] detection (Stark et al. 2015b) or with evidence of strong optical emission lines from their mid-IR photometry (Stark et al. 2017).
-The range of physical parameters allowed by the spectrophotometric fit does not allow the BDF pair to reionize their surroundings by themselves with a low or even a moderate escape fraction of ionizing photons, as they would require extremely high f esc =40% and age larger than 200 Myr (formation redshift z >9) (Fig. 4).
The results summarised above indicate that the two emitters are typical high-redshift sources with no peculiar ionizing capabilities, either due to their stellar populations or to a significant contribution from AGN emission.The low metal line emission and the moderate ξ * ion are likely explained by a combination of low metallicity and the galaxies being caught in a quiet phase of star-formation activity.In addition, the lack of a significant AGN contribution is not surprising in the light of theoretical predictions showing a dominant contribution from star-formation to the UV emission of bright high-redshift galaxies (Piana et al. 2022).Most importantly for our purposes, it is unlikely that they are solely responsible of the growth of the ionized bubble.Instead, a dominant contribution to the local reionization history can be provided by the overdensity of faint galaxies.
These findings allow us to draw a plan for future observations capable of fully constraining the properties of this remarkable region.As of today, the most important missing ingredient to confirm that the faint galaxies in the overdensity reionized the BDF region is their spectroscopic confirmation.Given the absence of strong Lyα this is beyond the capabilities of ground-based telescopes, while attempting detection of submillimeter lines with ALMA in redshift scan mode would be exceedingly time consuming.However, spectroscopic confirmation is well within reach of the James Webb Space Telescope.Indeed, JWST-NIRSpec in ∼1.5 hours of integration time could detect [OIII]λ4959, 5007 from all currently known candidates down to M UV −19, while also fully constraining ionization state and metallicity of the LAEs through the detection of [OII]λ3727 and Hβ.In addition, JWST can extend the photo-metric mapping of the BDF in terms of area, depth and wavelength coverage to confirm the presence of an overdensity of LBGs down to fainter magnitudes.For example, ∼1-2 hours of NIRCam integration per filter can extended the constraints on the local UV LF down to M UV ∼ −18.0, while covering the optical rest-frame wavelength range which is crucial to reject low redshift interlopers and constrain the spectral energy distributions.
Finally, a significant contribution to our understanding of the formation of reionized regions will be provided by objects beyond the limit that can be reached on blank fields such as the BDF, e.g.objects analogous to the z ∼6 clustered ultra-faint dwarfs observed by Vanzella et al. (2019) (see also Vanzella et al. 2021b).In this respect, JWST surveys on lensed fields (Treu et al. 2017;Willott et al. 2017;Vanzella et al. 2021a) will provide crucial information, such as an accurate characterization of the faintest end of the UV LF and of its possible turn-over (e.g., Castellano et al. 2016b;Yue et al. 2018), and of the physical properties of high-redshift ultra-faint dwarfs (e.g., Vanzella et al. 2017a,b).
With the sample of candidate reionized regions steadily growing, the advent of JWST will enable a systematic investigation of the connection between the physical properties of high-redshift galaxies and the ionization state of the surrounding IGM, eventually providing an answer on which are the sources of reionization.

Fig. 2 .
Fig. 2. Left: the position of BDF521 and BDF2195 on the f CIV / f Lyα versus Lyα(EW) plane compared to z 6 objects with detected CIV emission (Stark et al. 2015b; Mainali et al. 2017), and to Type II QSOs at intermediate redshift from Alexandroff et al. (2013).Right: the position of the BDF LAEs on the CIV(EW) versus M UV plane compared to z 6 sources from the literature including upper limits from Laporte et al. (2017), Mainali et al. (2018), and Shibuya et al. (2018).All upper limits are at 3σ, see legend for symbols and colors.

Fig. 3 .
Fig. 3.The position in the CIV(EW) versus HeII(EW) plane of models with 40% (crosses), 60% (triangles) and 100% (circles) contribution to the observed UV emission of the two BDF galaxies from narrow-line AGN (from Feltre et al. 2016) with different metallicities and ionization parameter.Symbols at increasing size are for models with Z=0.025, 0.1, 0.4 Z , whose color indicate the relevant ionization parameter log(U s ).The shaded regions highlight the parameter space allowed at 1σ, 2σ and 3σ by XSHOOTER observations of BDF521 (top) and BDF2195 (bottom).

a
Physical properties (mean values and 68% c.l. uncertainties) are from the BEAGLE spectro-photometric fit.bMass-weighted values.

Fig. 4 .
Fig. 4. The size (R bubble ) of the ionized bubble created by the combined emission of the BDF pair, as a function of the maximum stellar age of the stellar population for BEAGLE SED models within 95% c.l. from the best-fit.The cases for escape fraction f esc =0.1,0.4,0.8 are shown as squares, circles and stars, respectively, with color indicating the relevant, total, SFR of the system.The horizontal black line marks the minimum HII size R min =1.1 Mpc enabling Lyα to escape.The dark grey line indicates the minimum HII size R min =0.85Mpc in the case Lyα is affected by a 220km/s outflow.Vertical lines, from left to right, indicate formation redshifts z=9, 12 and 20.

Table 1 .
3-σ upper limits on the observed fluxes and the rest-frame EW of the UV emission lines of BDF521 and BDF2195

Table 2 .
Properties of the BDF pair a