Free Access
Issue
A&A
Volume 641, September 2020
Article Number A32
Number of page(s) 8
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/201936542
Published online 03 September 2020

© ESO 2020

1. Introduction

The study of galaxy clusters plays an important role in understanding cosmological structure formation and the astrophysics of galaxy evolution. Statistics concerning galaxy cluster size, mass, and redshift distribution provide constraints for cosmological models, while the properties of the galaxies and gas inside clusters give clues about galaxy evolution and the star formation history of the universe (e.g. Boylan-Kolchin et al. 2009). The progenitors of galaxy clusters, the so-called protoclusters, start off as over-dense regions and groups of galaxies at high redshift, which over time coalesce into the larger galaxy clusters we see today. Observing the early stages of cluster formation at higher redshifts has been challenging. Since protoclusters lack many of the observational properties of the massive virialised galaxy clusters of today, one of the best ways to find them is to identify galaxy over-densities at high redshift (Overzier 2016). Readily observable populations of galaxies include Lyman break galaxies (LBGs) and Lyα emitters (LAEs) (e.g. Overzier et al. 2008). These objects, which are compact galaxies that have strong emission in the Lyα line, are relatively easy to observe over a wide range of redshifts at z ∼ 2–6 (e.g. Taniguchi et al. 2005; Gronwall et al. 2007; Nilsson et al. 2009; Guaita et al. 2010; Cassata et al. 2015). LAEs are mainly star-forming, low-mass objects, and some may be the building blocks of Milky Way-like galaxies (Gawiser et al. 2007).

Giant Lyα-emitting nebulae, also known as Lyα “blobs” (LABs; Steidel et al. 2000; Matsuda et al. 2004) emit Lyα radiation on large scales (50–100 kpc) and have high Lyα luminosities of ∼1043 − 44 erg s−1. Moreover, it has been found that they are also apparent occasional signposts of LAE over-densities (e.g. Yang et al. 2010). The mechanism that powers the strong extended Lyα emission in these blobs is still poorly understood. Possible powering mechanisms include gravitational cooling radiation (e.g. Rosdahl & Blaizot 2012), the resonant scattering of Lyα photons produced by star formation (e.g. Cen & Zheng 2013), photo-ionising radiation from active galactic nuclei (AGNs; e.g. Yang et al. 2014), or a combination thereof (e.g. Reuland et al. 2003; Dey et al. 2005). Another potential source is shock heating from starburst-driven winds (Mori & Umemura 2006), although recent studies of the emission of non-resonant lines from eight Lyα blobs exclude models that require fast galactic winds driven by AGNs or supernovae (e.g. Prescott et al. 2015a). Regardless of the energy sources of Lyα blobs, the association of these blobs with compact LAE over-densities of ∼10–20 Mpc in size (e.g. Saito et al. 2015) suggests that LABs are good potential markers of large protoclusters. Furthermore, the number density and variance of Lyα blobs, as well as the 200–400 km s−1 relative velocities of their embedded galaxies, suggest that blobs themselves occupy individual group-like halos of ∼1013M (e.g. Prescott et al. 2015b). Thus, blobs may be sites of massive galaxy formation and trace significant components of the build-up of protoclusters.

These Lyα nebulae are usually found in the redshift range 2 <  z <  6 using different techniques, such as a narrow-band (NB) search or integral field spectroscopy (IFS). The narrow-band technique is the most successful method for detecting strong Lyα emission lines of galaxies, since it relies on a specific redshift interval as well as a selected low-sky background window. Using integral field spectroscopy is an optimal method as it does not require any spectroscopic follow-up and covers not only a large volume but also a wide range of redshifts.

To date, a few dozen of blobs have been reported. Bădescu et al. (2017) found three new z  ∼  2.3 blobs using Narrow-Band technique. Matsuda et al. (2004) and Arrigoni Battaia et al. (2019) found several tens of these objects around a redshift of z  ∼  3. Vanzella et al. (2017) and Lusso et al. (2019) found five LABs using the IFS technique with the VLT/Multi-Unit Spectroscopy Explorer (MUSE) around z  ∼  3.3. Borisova et al. (2016) detected one at z  ∼  3.5, and Saito et al. (2015) also detected one at z ∼ 4.1. The highest redshift LAB remains Himiko at z ∼ 6.6 discovered by Ouchi et al. (2009), also found using the narrow-band technique.

This paper is organised as follows: in Sect. 2, we first present the method and the data set we used. In Sect. 3, we explain the characteristics of the newly detected Lyα nebula. Throughout this study, we adopted the following cosmological parameters: H0 = 70 km s−1 Mpc−1, Ωm = 0.3, and ΩΛ = 0.7. All magnitudes are given using the AB system (Oke & Gunn 1983).

2. Method and data set

The ATACAMA project (PI: P. Hibon; Hibon et al., in prep.) aims to perform a systematic search for Lyα nebulae in a multi-wavelength broad-band imaging data set, instead of using the conventional narrow-band imaging technique. This strategy is similar to the one described in Prescott et al. (2012) who pioneered this method by confirming the detection of several LABs at 1.7 ≤ z ≤ 2.7 (Prescott et al. 2013). Based on this method, we designed a software (Thomas et al., in prep.) that consists of a GUI communicating with SExtractor (Bertin & Arnouts 1996) and optimised for LAB detection. We use a single MegaCam ultra deep pointing of 1 sq. degree of the Cosmic Evolution Survey (COSMOS, Scoville et al. 2007) field from the CFHT Legacy Deep Survey (CFHTLS Deep), reaching limiting magnitudes between 25 and 26.44 depending on the filter (80% completeness limit). This provides a high-quality and homogeneous optical data set, calibrated photometrically to 1.0% and astrometrically to 0.028 arcsec. This data set was reduced by Terapix dedicated to the processing of very large data flows from digital sky surveys. The transmission curves are available online1.

We selected the COSMOS field due to the availability of a very wide range of wavelength bands in various data archives.

The ATACAMA tool processed the images through various steps: (1) the progressive elimination of bright, intermediate, and faint point-like objects in one specific band; (2) the production of wavelet planes allowing us to detect extended objects; (3) the creation of a photometric candidate catalogue with the spectral energy distribution in all available broad-band images; (4) the cross-correlation of the catalogue of candidates with the COSMOS2015 photometric redshift catalogue (Laigle et al. 2016), which contains photometric redshifts for over half a million sources in the CFHTLS Deep2 field; (5) the determination of the photometric redshift through spectral energy distribution (SED) fitting with LePhare (Arnouts & Ilbert 2011) in case of no matching object in the COSMOS15 catalogue; (6) the careful visual inspection of each object. These steps are performed by several team members and the final results are compared between the different parameter choices of each step.

3. The Lyα nebula

One object (LAB1) in particular caught our immediate interest, after we ran ATACAMA several times and independently of the parameter variations during the different tests.

3.1. Photometric information

As mentioned in the previous section, SExtractor is used to estimate the photometric information of this object. We chose to use the isophotal-corrected magnitudes as such an extended object does not have all its flux within neat boundaries. The double-image mode is used to derive the multi-band photometry of this LAB, with the g′-band image as the detection image and all other bands as photometric bands. Figure 1 shows the multi-wavelength 28″ × 28″ images of the LAB. We choose not to show the Herschel/SPIRE stamps at 250microns, 350microns and 500microns, as their resolution is not optimal to possibly find a detection (6″ pixel−1, 10″ pixel−1, 14″ pixel−1, respectively). Prescott et al. (2012) selected their highest 2 ≤ z ≤ 3 LAB candidate with the following colour: BW − R <  0.45 while searching for LABs in the BW band image. We verified the colour we obtained from its magnitudes. This LAB has g′−r′ = 0.45. The dropout has a magnitude difference of u′−r′ = 2.71. Our object is neither detected in J-, H- nor Ks-band images. van der Burg et al. (2010) led a search for Lyman-break Galaxies in the CFHTLS data. They defined the following criteria to detect u-dropout candidates:

1.0 < ( u g ) 1.0 < ( g r ) < 1.2 1.5 ( g r ) < ( u g ) 0.75 . $$ \begin{aligned}&1.0 < (u^{*}-g\prime )\\&-1.0 < (g\prime -r\prime ) < 1.2 \\&1.5(g\prime -r\prime ) < (u^{*}-g\prime )-0.75. \end{aligned} $$

thumbnail Fig. 1.

Thumbnails. North up, east left. Size of each stamp is 28 arcsec side-wise. We show here NUV and FUV from GALEX, the optical bands from CFHT/Megacam, the bands available from Subaru/Suprime-Cam, the NIR bands from VISTA/VIRCam and from CFHT/WIRCam, the four channels from Spitzer/IRAC and the 24 μm Spitzer/MIPS. The circles are 7.65 arcsec in radius and are encompassing the complete nebula.

The colours of our newly detected Lyα nebula respect their criteria and reinforce the detection and the high-redshift nature of this object.

3.2. Morphology of LAB1

This LAB is composed of one main clump and a halo, as can be seen in the g′- and r′-band (see Fig. 1). Its complete extent is 28 kpc  ×  59 kpc, which corresponds to 29.4 arcsec2, measured in the g′-band image. The main clump, also called the bright core, is located at RA = 149.97769, Dec = 1.9358. The halo or filamentary structure is orientated 45deg north-east of the bright core, located at RA = 149.9777428, Dec = 1.9357793, as seen and indicated in Fig. 2, where we show the contour of this LAB with the filamentary structure encircled in red (area of the red ellipse: 24.3 sq.arcsec) and the bright core in blue (area of the blue circle: 5.4 sq.arcsec). We also show in the same figure the r′-band image subtracted by the g′-band.

thumbnail Fig. 2.

Left: contour of LAB1 on the 15″ × 15″ r-band image. North is up, east is left. The identified bright core is encircled in blue. The filamentary structure is encircled in red. The second image shows the r − g bands subtraction. The third image is the Bw-band of object 59 of Prescott et al. (2012). Right: Bw-band of object 26 of Prescott et al. (2012). Both Bw-band stamps were extracted from Fig. 10 of Prescott et al. (2012).

The morphology of this Lyα nebula is very similar to the famous one found by Steidel et al. (2000, Fig. 6) (called Blob 1), as it can be resolved into a bright source and a diffuse gas halo. The bright core is not centred on the peak of the extended emission but it is located at 5.58 arcsec in the south-western direction (see Fig. 2). There is no established consensus on the relation morphology-powering source yet. Both radio-loud and radio-quiet LABs can show filamentary and/or bubble-like structures (Matsuda et al. 2004; Dijkstra & Loeb 2009; Geach et al. 2009). Prescott et al. (2012) divided their candidates in two morphological categories: GROUP and DIFFUSE. GROUP corresponds to a tight grouping of compact sources selected as a single large source, and DIFFUSE indicates spatially extended diffuse emission. When comparing our object (Figs. 1 and 2) to the complete sample of Prescott et al. (2012), we first remark that it is more similar in morphology with the DIFFUSE ones. From Fig. 10 of Prescott et al. (2012), 21 out of 85 blobs are diffuse. Comparing this sub-sample with LAB1 in more detail, it presents the same kind of structure as source numbers: 5, 10, 26, 59. Figure 2 shows the detection band for objects 26 and 59 on the side of our r′-band contoured object, to facilitate the visual morphology comparison.

The morphologies in the sample of Matsuda et al. (2011) are very diverse, going from compact systems to filamentary, bubble-like, conical structures, including some with shapes similar to our candidate.

Borisova et al. (2016) and Arrigoni Battaia et al. (2019) Lyα nebulae with a very circular morphology, centred on the host quasi-stellar objects (QSOs). However, blob 2 of Borisova et al. (2016) is elongated, with the QSO on one of the extremity. It could also be seen as a more conical structure. From both samples, this is another very similar object to LAB1.

Cai et al. (2017) presented a blob at z = 2.3 with a filamentary structure possibly aligned with a large-scale structure found previously (Cai 2016).

The newly discovered z = 3.32 LAB from Marques-Chaves et al. (2019) shows an elongated morphology (and is extended over ≃11″, or ≃85 kpc at z = 3.33, within the 3σ detection limit.)

Reuland et al. (2003) reported the detection of Lyα nebulae centred on z >  3 radio galaxies, two at z ∼ 3.8, one at z ∼ 3.4. We focus our comparison on B2 0902+34 (z ∼ 3.4) which shows a diffuse morphology with an extension of 10″ × 8″ and some possible signs of large-scale filamentary emission (Figs. 7 and 8 of Reuland et al. 2003).

Vanzella et al. (2017) reported a z ∼ 3.3 LAB in the Hubble Ultra Deep Field (HUDF) with a size of ∼40 kpc × 80 kpc. Lusso et al. (2019) found three blobs at z ∼ 3.2: the most extended one has a 150 kpc and the smallest one 35 kpc. The measured size of our newly found object is therefore in agreement with the Lyα nebulae previously found in a similar redshift range.

3.3. Ancillary data

Based on the position of our blob candidate we looked for counterparts from different astronomical resources. We list below the result of this search from different facilities.

HST. The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) data set (Grogin et al. 2011; Koekemoer et al. 2011) covers 210 sq.arcmin in the following ACS (Advanced Camera for Surveys) and WFC3 (Wide Field Camera 3) filters: F606W, F814, F350LP, F125W and F160W. However, the CANDELS area does not overlap with our object.

CFHT. The CFHTLS also provides the WIRCam Deep Survey (WIRDS, Bielby et al. 2012), an extremely deep, high-quality (FWHM ∼ 0.6″) J, H, Ks imaging covering a total effective area of 2.1 sq.deg and reaching AB 50% completeness limits of ∼24.5.

Chandra. We verified the available merged images (Elvis et al. 2009) corresponding to the X-ray energy bands 0.5–2 keV (soft), 2 keV–7 keV (hard), and 0.5–7 keV (full) energy bands, reaching 2.2 × 10−16 cgs, 1.5−15 cgs and 8.9 × 10−16 cgs in depth, respectively. There is no detection associated with our object.

GALEX. We also searched for a possible counterpart in the near and far-UV images reaching a depth of 26 AB magnitude (Zamojski et al. 2007). We did not identify any detection in these bands.

Subaru. We checked the Subaru/Suprime-Cam mosaic images available in several filters (Taniguchi et al. 2015). We detected our object in the Subaru mosaic with the intermediate filter IA527, centred at 5262.3 Å, FWHM = 243.5 Å. As the photometric redshift of our Lyα nebula is 3.32 (see Sect. 3.4), we should expect to have a detection of the extended emission at λ ∼ 5225 Å, wavelength covered by this intermediate band filter. The Lyα nebula is detected in this image.

Spitzer. Looking at a possible detection in IRAC channels 3.6 and 4.5 microns, and in the MIPS image at 24 microns, we found a marginal detection in the 24 μm MIPS image (Sanders et al. 2007) at RA = 149.97831, Dec = 1.8353, at 2.86 arcsec south-east of the bright core of our nebula. Hence, we can assume that this could be a marginal detection of the bright core of this LAB. The 24 μm MIPS image has a 3σ detection limit of 54 μJy. Using the zero-point information given in the MIPS Instrument Handbook, we found a flux of 29.9 μJy ± 5.9 for the possible bright core detection.

Herschel. We looked for a possible detection in the SPIRE images reaching the following 5σ noise levels: 8.0 mJy, 6.6 mJy, 9.5 mJy in the SPIRE bands 250, 350 and 500 mJy, respectively (Oliver et al. 2012; Chapin et al. 2011; Levenson et al. 2010; Viero et al. 2013, 2014). The very low resolution of the SPIRE images (6″ pixel−1,10″ pixel−1 and 14″ pixel−1 for the 250, 350 and 500 mJy bands respectively) prevent us from concluding on the possible detection for any counterpart corresponding to the associated bright core, which could be the host galaxy.

We thus cannot rule out the presence of an associated active galactic nucleus, which could be the powering source of this extended Lyα emission.

VLA. We verified the available VLA images corresponding to the following frequencies: 327 MHz (P-band, 90 cm), 1.4 GHz ( = L-band, 20 cm) (Schinnerer et al. 2007; Bondi et al. 2008), and S-band at 2–4 GHz (Smolčić et al. 2017). These images reach the following 3σ detection limit: 0.12 mJy in P-band, 21 μJy in L-band and 3.75 mJy in S-band. There is no VLA detection associated either with the bright core or the filamentary component of our object.

ALMA. We looked for corresponding radio sources in the ALMA Science Archive. We did not find any ALMA data taken in the surroundings of our LAB.

3.4. Photometric redshift

The COSMOS2015 catalogue (Laigle et al. 2016) contains precise photometric redshift for more than half a million objects over the whole COSMOS field. They derived them using the photometry of 30 bands, several tens of templates and several extinction laws. For a redshift in the range 3 <  z <  6, they reached a photometric redshift precision of 0.021. We verified any possible corresponding object for this LAB, and found no matching entry in the COSMOS2015 catalogue within a search radius of 2 arcsec. It is also absent from the COSMOS Master Spectroscopic Redshift Catalogue, which contains more than 65 000 sources (M. Salvato, priv. comm.). Hence, we need to estimate a photometric redshift.

We fitted its SED using the PHotometric Analysis for Redshift Estimations, also known as LePhare (Arnouts et al. 1999). The photometric redshift was then computed using LePhare (Arnouts & Ilbert 2011), which is based on a χ2 template-fitting method (Arnouts & Ilbert 2011). We used the galaxy templates optimised for the COSMOS survey (Ilbert et al. 2006) with the emission line contribution enabled.

Table 1 presents the 3σ magnitudes for the complete nebula, the bright core and the faint halo in all the available filters. We used the following bands to obtain necessary photometric information for the bright core, the faint core, and the complete nebula. We used aperture magnitudes obtained with SExtractor with a diameter aperture of 3 arcsec for the bright core, 15 pixels for the faint part, and 35 pixels for the whole nebula.

Table 1.

Aperture magnitudes of the complete nebula: 35 pixels aperture, 3σ detection, the bright core: 3 arcsec aperture, 3σ detection, and the faint halo: 15 pixels aperture, 3σ detection.

A 3σ detection threshold is used to determine the magnitudes for all of them.

For the bright core, different tests of filter combinations were realised with a minimum of 11 filters used each time. For some of these tests, we used the flux measured in the 24 μm MIPS image, but we chose not to include it in the figure as we have no certainty it really corresponds to the bright core. All the tests are consistent with a result of zphot = 3.32 ± 0.05 with χ2 = 25.9. Similarly to the bright core, we also realised several tests with different apertures (10-, 15- and 20-pixel diameters) and different filter combinations to estimate the photometric redshift of the faint core, taking advantage of the suit of data available in the COSMOS field and of the marginal MIPS24 μm detection. All the tests are consistent with a resdhift estimation of zphot = 3.28 ± 0.03 with χ2 = 20.88. The photometric redshift of the complete nebula was calculated with an aperture of 35 pixels in diameter. The same strategy was employed by using different combinations of filters. All confirmed a photometric redshift of zphot = 3.3 ± 0.05 with χ2 = 15.67. The resulting fits are shown in Fig. 3.

thumbnail Fig. 3.

Best-fit galaxy template in observed wavelengths from LePhare including the emission line contribution. Twelve filters were used to estimate these photometric redshift. The filled squares represent the ISOCOR magnitudes of our object, and the empty circles show the magnitudes from the best model. The downward facing triangles represent the detection limit used for the photometric redshift estimation. Top: for the complete nebula, magnitude of aperture 35 pixels. Centre: for the bright core, with 3 arcsec aperture, 3sigma detection. Bottom: for the faint halo with 15 pixels aperture. The three fits give a redshift of around 3.3 ± 0.05.

3.5. Lyα line

With the applied broad-band selection technique, estimating the Lyα equivalent width is a delicate task. Moreover, we know that all the flux coming from the g-band cannot be fully accounted to the Lyα line.

Based on the redshift estimation of the previous section, the Lyα line, at ∼5248 Å should be present in the intermediate Subaru filter IA527. Assuming that all the flux in this filter is associated with the Lyα emission, we compute that the Lyα flux is ∼1.5 × 10−15 erg s−1 cm−1. Obviously, this should be considered as an upper limit. We also tried to estimate the Lyα equivalent width using broad-band SED fitting. For this exercise, we used the SPARTAN software2 (Thomas, in prep.). SPARTAN is an SED fitting algorithm based on a χ2 statistic. It is able to use both photometry and spectroscopy, either separately or at the same time to infer the physical parameters of galaxies. For the purpose of the Lyα equivalent width estimation, we used Bruzual & Charlot (2003) models at low resolution with the Chabrier (2003) initial mass function. The metallicity can take values up to the solar metallicity (Z). We used a star formation history exponentially delayed with a timescale parameter from 0.1 to 2.0 Gyr. In terms of attenuation, we used the Calzetti’s dust attenuation (Calzetti et al. 2000) with E(B − V) ranging from 0.0 to 0.5, and we used the IGM prescription of Thomas et al. (2017) which made it possible to treat the IGM as a free parameter. SPARTAN includes the possibility to add emission lines. This is done following the method proposed in the paper by Schaerer & de Barros (2009). In this method both nebular continuum and emission lines are added converting the number of ionising photons to Hβ flux. Then, the other emission lines are added using the metallicity dependent line ratios given by Anders & Alvensleben (2003). While these line ratios are fixed during the fitting process, it is possible to test multiple possibilities by changing them manually. This allows us to compute the equivalent width directly on the theoretical model for different line ratios. We made the Lyα/Hβ ratio vary between 2 and 40 in steps of 2. At each step, we computed the equivalent width in the theoretical model and estimated the χ2.

The results are displayed in Fig. 4 where we display the evolution of the χ2 from the SPARTAN fit with respect to the equivalent width of Lyα (computed in the observed frame). It shows that there are two χ2 minima at EW(Lyα)∼–291 Å and EW(Lyα)∼–118 Å. This corresponds to line fluxes of 2.02 × 10−15 erg s−1 cm−1 and 1.38×10−15 erg s−1 cm−1, respectively. The difference in χ2 for these two solutions is less than 1. This is in good agreement with the estimation from the Subaru band presented above. In addition, these values are consistent with the LAB search of Matsuda et al. (2004) whose few candidates have a line flux ∼10−15 erg s−1 cm−1. The two nebulae, identified by Steidel et al. (2000) also have a typical Lyα flux of the order of ∼1.3 × 10−15 erg s−1 cm−1 (for a physical extent of ∼100 kpc). Lusso et al. (2019) estimated the Lyα flux for > 100 kpc blobs and found them to be at ∼0.2 − 0.8 × 10−15 erg s−1 cm−1. Finally, Marques-Chaves et al. (2019) found a Lyα flux of ∼1.8 × 10−15 erg s−1 cm−1 for the giant blob they found that extend to ∼100 kpc. Our estimations will have to be verified using spectroscopic data.

thumbnail Fig. 4.

Attempt to estimate the equivalent width of Lyα based on the variation of the line ratios during SED fitting. The two vertical dashed lines give the position of the two local minima. The y-axis has been reversed for clarity and we choose the emission line to be negative.

4. Environment and contaminants

In this section, we discuss the different aspects of the detection of our candidate.

Possible contaminants. There are several populations that may contaminate our Lyα nebula candidate sample. The morphological selection method tends to select any source of low surface-brightness extended emission that cannot be easily identified, or low surface-brightness nebulae within the Galaxy. In principle, lower redshift sources such as [OII]λλ3727,3729 emitters, are a contaminant population in Lyα-emitting galaxy surveys (e.g. Hibon et al. 2010, 2011, 2012). The [OII] emitters would have a redshift in the range [0.065;0.55] considering the transmission of the CFHT/MegaCam g′-band filter with a mean redshift of z ∼ 0.3 for λmean ∼ 4862.5 Å. The lack of detection in the broad-band COSMOS HST data very likely rules out all of the low-redshift options. Thanks to the result from LePhare, we can hence discard these usual contaminants in such a broad-band search.

Quasar. Using quasar templates during the photometric redshift estimation, we find a photometric redshift of z ∼ 3.3 ± 0.05. This reinforces the assumption that this Lyα nebula could be associated with a quasar at a similar redshift, as explained in Borisova et al. (2016) and Husemann et al. (2018) among others. Borisova et al. (2016) searched for giant Lyα nebulae around 12 of the brightest radio-quiet quasars. Husemann et al. (2018) studied the environment and the extended Lyα nebula of one specific luminous radio-quiet QSO. For such a redshift of z ∼ 3.3, we also expect to have the Lyman limit at λ  ∼ 3921 Å, which agrees with our SED (see Fig. 3).

Radio detection. The absence of any radio detection informs us that there is no radio core associated with the nebula. The marginal detection of the bright core in the MIPS 24 μm image with the absence of any radio detection could indicate the presence of a radio-quiet quasar associated with the bright core. However, with its very limited resolution, the MIPS image does not provide more information.

Sargent et al. (2010) differentiate star-forming galaxies and AGNs from their (u − K) colour: AGN sources should show (u − K) ≥ 2.42. The bright core of our LAB has a u − K = 1.53. Hence, following their criterion, it is likely that the bright core of LAB1 is a star-forming galaxy rather than an AGN. Using their definition of q24, obs (Eq. (2) of Sargent et al. 2010), determined as the ratio of observed MIPS 24 μm and VLA 1.4 GHz fluxes, and their Fig. 11 showing q24, obs in function of redshift for star-forming and AGN, we remark that at z ∼ 3.3, star-forming galaxies would have q24, obs <  1 and AGN q24, obs ∼ 1. Using the VLA detection limit for our bright core, we have an upper limit of its FIR/Radio ratio of q24, obs ∼ 1.42. Data with a better FIR resolution and deeper 1.4 GHz would help us to establish better constraints on the nature of LAB1’s bright core.

Clustering. From the COSMOS spectroscopic catalogue (M. Salvato, priv. comm.), there is no evidence of a cluster at a similar redshift or at another redshift in the vicinity of our object within a one arcmin radius. The redshift distribution in this same area is in the range of z = 0 to 1.5. No source is included in the COSMOS spectroscopic catalogue within 30 arcsec radius from the nebula, except three foreground objects all at z <  1: (09:59:545461; +01:56:36.395), (09:59:54.4717; +01:56:22.861) and (09:59:54.4967; +01:55:59.426). From the COSMOS2015 photometric redshift catalogue (Laigle et al. 2016), in a one arcmin radius around this LAB, there are 61 galaxies with an estimated photometric redshift in the range of 0.05 <  zphot <  3.68. Eight of these objects have zphot ≥ 3. Five of them are in a 2.4 arcmin radius of each other, at 36 arcsec from our object. These five objects have a photometric redshift in the range 3.035 <  zphot <  3.68.

From Ono et al. (2018), the only spectroscopically confirmed Lyα emitter at z <  4, and more precisely at z = 3.41 is 8.41 arcmin north-east of our object. We also verified the presence of Lyα emitter at z ∼ 3.23 selected in the Subaru/IA527 filter and presented in Sobral et al. (2018) and Santos et al. (2020). They created and studied a catalogue of 641 z ∼ 3.23 LAEs. Their four closest z ∼ 3.23 LAEs are at the following separations: 1.79 arcmin S-E, 1.82 arcmin N-W, 1.96 arcmin N-E and 2.56 arcmin S-W.

The nebula discovered by Marques-Chaves et al. (2019) and the B2 0902+34 nebula detected by Reuland et al. (2003) do not have a possible cluster or protocluster in their proximity. As mentioned previously, Borisova et al. (2016) searched especifically around bright QSOs, Matsuda et al. (2011) in the SSA22 protocluster, and Cai et al. (2017) in a Lyα absorption system. Their search and selection focused on LABs detected in specific kinds of surroundings. The absence of a cluster at similar redshift in the vicinity of our nebula candidate does not affect its legitimate detection. Its spectroscopic confirmation is essential to assessing the complete situation.

5. Conclusions

This paper presents the first result obtained from the ATACAMA project (PI: Hibon).

  1. We ran a broad-band search technique dedicated to searching for Lyα nebulae, presented in Prescott et al. (2012). After cross-correlating our candidates with the COSMOS2015 photometric redshift catalogue, we used the LePhare tool to determine the photometric redshift of the candidates absent from the COSMOS2015.

  2. We identified in the g′-band filter, from the CFHTLS COSMOS data set, one particular Lyα nebula covering an isophotal area of 29.4 arcsec2. Morphologically, we can identify a bright core and a faint core. Its redshift was photometrically estimated at z ∼ 3.3. We obtained a first estimation of the Lyα equivalent width and line flux. Both agree with the values from the Lyα nebulae found and spectroscopically observed by Matsuda et al. (2004).

  3. With VLT/XSHOOTER spectroscopy we will confirm the redshift of this exciting LAB, which will then be part of the high redshift Lyα blob sample discovered in a blind search. We are planning to use the following slit widths: 1.3″, 1.2″ and 1.2″ for the UVB, VIS, and NIR arms, respectively. A position angle of 100 deg will be ideal to observe both the bright core and part of the faint halo. We will set up the exposure times at 19 × 600 s for each arm. The study of such an object can significantly boost our understanding of this class of Lyα emitters and bring light to their still debated origin and formation. With a redshift of ∼3.3, the kinematics of the Lyα line will enhance our understanding of galaxies’ gas inflows and outflows and the escape of Lyα photons at the start of the cosmic noon, the epoch of the peak of star-forming activity.


Acknowledgments

We thank the reviewer for his/her thorough review and highly appreciate the comments and suggestions, which significantly contributed to improving the quality of the publication. Based on observations obtained with MegaPrime/MegaCam, a joint project of CFHT and CEA/IRFU, at the Canada-France-Hawaii Telescope (CFHT) which is operated by the National Research Council (NRC) of Canada, the Institut National des Science de l’Univers of the Centre National de la Recherche Scientifique (CNRS) of France, and the University of Hawaii. This work is based in part on data products produced at Terapix available at the Canadian Astronomy Data Centre as part of the Canada-France-Hawaii Telescope Legacy Survey, a collaborative project of NRC and CNRS. Based on observations obtained with WIRCam, a joint project of CFHT, Taiwan, Korea, Canada, France, at the Canada-France-Hawaii Telescope (CFHT) which is operated by the National Research Council (NRC) of Canada, the Institut National des Sciences de l’Univers of the Centre National de la Recherche Scientifique of France, and the University of Hawaii. This work is partly based on tools and data products produced by GAZPAR operated by CeSAM-LAM and IAP. The fits-file handling was done using the Thomas (2019) python module.

References

  1. Anders, P. V., & Alvensleben, U. 2003, A&A, 401, 1063 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Arnouts, S., & Ilbert, O. 2011, LePHARE: Photometric Analysis for Redshift Estimate (Astrophysics Source Code Library) [Google Scholar]
  3. Arnouts, S., Cristiani, S., Moscardini, L., et al. 1999, MNRAS, 310, 540 [NASA ADS] [CrossRef] [Google Scholar]
  4. Arrigoni Battaia, F., Hennawi, J. F., Prochaska, J. X., et al. 2019, MNRAS, 482, 3162 [NASA ADS] [CrossRef] [Google Scholar]
  5. Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Bielby, R., Hudelot, P., McCracken, H. J., et al. 2012, A&A, 545, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Bondi, M., Ciliegi, P., Schinnerer, E., et al. 2008, ApJ, 681, 1129 [NASA ADS] [CrossRef] [Google Scholar]
  8. Borisova, E., Cantalupo, S., Lilly, S. J., et al. 2016, ApJ, 831, 39 [NASA ADS] [CrossRef] [Google Scholar]
  9. Boylan-Kolchin, M., Springel, V., White, S. D. M., Jenkins, A., & Lemson, G. 2009, MNRAS, 398, 1150 [NASA ADS] [CrossRef] [Google Scholar]
  10. Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000 [NASA ADS] [CrossRef] [Google Scholar]
  11. Bădescu, T., Yang, Y., Bertoldi, F., et al. 2017, ApJ, 845, 172 [CrossRef] [Google Scholar]
  12. Cai, Z. 2016, Imaging a Massive Galaxy Overdensity at z= 2.3: The Morphology-Density Relation at High Redshift, HST Proposal [Google Scholar]
  13. Cai, Z., Fan, X., Yang, Y., et al. 2017, ApJ, 837, 71 [NASA ADS] [CrossRef] [Google Scholar]
  14. Calzetti, D., Armus, L., Bohlin, R. C., et al. 2000, ApJ, 533, 682 [NASA ADS] [CrossRef] [Google Scholar]
  15. Cassata, P., Tasca, L. A. M., Le Fèvre, O., et al. 2015, A&A, 573, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Cen, R., & Zheng, Z. 2013, ApJ, 775, 112 [NASA ADS] [CrossRef] [Google Scholar]
  17. Chabrier, G. 2003, PASP, 115, 763 [Google Scholar]
  18. Chapin, E. L., Chapman, S. C., Coppin, K. E., et al. 2011, MNRAS, 411, 505 [NASA ADS] [CrossRef] [Google Scholar]
  19. Dey, A., Bian, C., Soifer, B. T., et al. 2005, ApJ, 629, 654 [NASA ADS] [CrossRef] [Google Scholar]
  20. Dijkstra, M., & Loeb, A. 2009, MNRAS, 400, 1109 [NASA ADS] [CrossRef] [Google Scholar]
  21. Elvis, M., Civano, F., Vignali, C., et al. 2009, ApJS, 184, 158 [NASA ADS] [CrossRef] [Google Scholar]
  22. Gawiser, E., Francke, H., Lai, K., et al. 2007, ApJ, 671, 278 [NASA ADS] [CrossRef] [Google Scholar]
  23. Geach, J. E., Alexander, D. M., Lehmer, B. D., et al. 2009, ApJ, 700, 1 [NASA ADS] [CrossRef] [Google Scholar]
  24. Grogin, N. A., Kocevski, D. D., Faber, S. M., et al. 2011, ApJS, 197, 35 [NASA ADS] [CrossRef] [Google Scholar]
  25. Gronwall, C., Ciardullo, R., Hickey, T., et al. 2007, ApJ, 667, 79 [NASA ADS] [CrossRef] [Google Scholar]
  26. Guaita, L., Gawiser, E., Padilla, N., et al. 2010, ApJ, 714, 255 [NASA ADS] [CrossRef] [Google Scholar]
  27. Hibon, P., Cuby, J. G., Willis, J., et al. 2010, A&A, 515, A97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Hibon, P., Malhotra, S., Rhoads, J., & Willott, C. 2011, ApJ, 741, 101 [NASA ADS] [CrossRef] [Google Scholar]
  29. Hibon, P., Kashikawa, N., Willott, C., Iye, M., & Shibuya, T. 2012, ApJ, 744, 89 [NASA ADS] [CrossRef] [Google Scholar]
  30. Husemann, B., Worseck, G., Arrigoni Battaia, F., & Shanks, T. 2018, A&A, 610, L7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Ilbert, O., Arnouts, S., McCracken, H. J., et al. 2006, A&A, 457, 841 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Koekemoer, A. M., Faber, S. M., Ferguson, H. C., et al. 2011, ApJS, 197, 36 [NASA ADS] [CrossRef] [Google Scholar]
  33. Laigle, C., McCracken, H. J., Ilbert, O., et al. 2016, ApJS, 224, 24 [NASA ADS] [CrossRef] [Google Scholar]
  34. Levenson, L., Marsden, G., Zemcov, M., et al. 2010, MNRAS, 409, 83 [NASA ADS] [CrossRef] [Google Scholar]
  35. Lusso, E., Fumagalli, M., Fossati, M., et al. 2019, MNRAS, 485, L62 [NASA ADS] [Google Scholar]
  36. Marques-Chaves, R., Pérez-Fournon, I., Villar-Martín, M., et al. 2019, A&A, 629, A23 [CrossRef] [EDP Sciences] [Google Scholar]
  37. Matsuda, Y., Yamada, T., Hayashino, T., et al. 2004, AJ, 128, 569 [NASA ADS] [CrossRef] [Google Scholar]
  38. Matsuda, Y., Yamada, T., Hayashino, T., et al. 2011, MNRAS, 410, L13 [NASA ADS] [CrossRef] [Google Scholar]
  39. Mori, M., & Umemura, M. 2006, Nature, 440, 644 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  40. Nilsson, K. K., Tapken, C., Møller, P., et al. 2009, A&A, 498, 13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Oke, J. B., & Gunn, J. E. 1983, ApJ, 266, 713 [NASA ADS] [CrossRef] [Google Scholar]
  42. Oliver, S. J., Bock, J., Altieri, B., et al. 2012, MNRAS, 424, 1614 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Ono, Y., Ouchi, M., Harikane, Y., et al. 2018, PASJ, 70, S10 [NASA ADS] [CrossRef] [Google Scholar]
  44. Ouchi, M., Ono, Y., Egami, E., et al. 2009, ApJ, 696, 1164 [NASA ADS] [CrossRef] [Google Scholar]
  45. Overzier, R. A. 2016, A&ARv, 24, 14 [NASA ADS] [CrossRef] [Google Scholar]
  46. Overzier, R. A., Bouwens, R. J., Cross, N. J. G., et al. 2008, ApJ, 673, 143 [NASA ADS] [CrossRef] [Google Scholar]
  47. Prescott, M. K. M., Dey, A., & Jannuzi, B. T. 2012, ApJ, 748, 125 [NASA ADS] [CrossRef] [Google Scholar]
  48. Prescott, M. K. M., Dey, A., & Jannuzi, B. T. 2013, ApJ, 762, 38 [NASA ADS] [CrossRef] [Google Scholar]
  49. Prescott, M. K. M., Martin, C. L., & Dey, A. 2015a, ApJ, 799, 62 [NASA ADS] [CrossRef] [Google Scholar]
  50. Prescott, M. K. M., Momcheva, I., Brammer, G. B., Fynbo, J. P. U., & Moller, P. 2015b, ApJ, 802, 32 [NASA ADS] [CrossRef] [Google Scholar]
  51. Reuland, M., van Breugel, W., Röttgering, H., et al. 2003, ApJ, 592, 755 [NASA ADS] [CrossRef] [Google Scholar]
  52. Rosdahl, J., & Blaizot, J. 2012, MNRAS, 423, 344 [NASA ADS] [CrossRef] [Google Scholar]
  53. Saito, T., Matsuda, Y., Lacey, C. G., et al. 2015, MNRAS, 447, 3069 [NASA ADS] [CrossRef] [Google Scholar]
  54. Sanders, D. B., Salvato, M., Aussel, H., et al. 2007, ApJS, 172, 86 [Google Scholar]
  55. Santos, S., Sobral, D., Matthee, J., et al. 2020, MNRAS, 493, 141 [CrossRef] [Google Scholar]
  56. Sargent, M. T., Schinnerer, E., Murphy, E., et al. 2010, ApJS, 186, 341 [Google Scholar]
  57. Schaerer, D., & de Barros, S. 2009, A&A, 502, 423 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Schinnerer, E., Smolčić, V., Carilli, C. L., et al. 2007, ApJS, 172, 46 [NASA ADS] [CrossRef] [Google Scholar]
  59. Scoville, N., Aussel, H., Brusa, M., et al. 2007, ApJS, 172, 1 [NASA ADS] [CrossRef] [Google Scholar]
  60. Smolčić, V., Novak, M., Bondi, M., et al. 2017, A&A, 602, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Sobral, D., Santos, S., Matthee, J., et al. 2018, MNRAS, 476, 4725 [NASA ADS] [CrossRef] [Google Scholar]
  62. Steidel, C. C., Adelberger, K. L., Shapley, A. E., et al. 2000, ApJ, 532, 170 [NASA ADS] [CrossRef] [Google Scholar]
  63. Taniguchi, Y., Ajiki, M., Nagao, T., et al. 2005, PASJ, 57, 165 [NASA ADS] [Google Scholar]
  64. Taniguchi, Y., Kajisawa, M., Kobayashi, M. A. R., et al. 2015, PASJ, 67, 104 [NASA ADS] [Google Scholar]
  65. Thomas, R. 2019, J. Open Source Soft., 4, 1249 [NASA ADS] [CrossRef] [Google Scholar]
  66. Thomas, R., Le Fèvre, O., Le Brun, V., et al. 2017, A&A, 597, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. van der Burg, R. F. J., Hildebrandt, H., & Erben, T. 2010, A&A, 523, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. Vanzella, E., Balestra, I., Gronke, M., et al. 2017, MNRAS, 465, 3803 [NASA ADS] [CrossRef] [Google Scholar]
  69. Viero, M. P., Wang, L., Zemcov, M., et al. 2013, ApJ, 772, 77 [NASA ADS] [CrossRef] [Google Scholar]
  70. Viero, M. P., Asboth, V., Roseboom, I. G., et al. 2014, ApJS, 210, 22 [NASA ADS] [CrossRef] [Google Scholar]
  71. Yang, Y., Zabludoff, A., Eisenstein, D., & Dave, R. 2010, ApJ, 719, 1654 [NASA ADS] [CrossRef] [Google Scholar]
  72. Yang, Y., Walter, F., Decarli, R., et al. 2014, ApJ, 784, 171 [CrossRef] [Google Scholar]
  73. Zamojski, M. A., Schiminovich, D., Rich, R. M., et al. 2007, ApJS, 172, 468 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1.

Aperture magnitudes of the complete nebula: 35 pixels aperture, 3σ detection, the bright core: 3 arcsec aperture, 3σ detection, and the faint halo: 15 pixels aperture, 3σ detection.

All Figures

thumbnail Fig. 1.

Thumbnails. North up, east left. Size of each stamp is 28 arcsec side-wise. We show here NUV and FUV from GALEX, the optical bands from CFHT/Megacam, the bands available from Subaru/Suprime-Cam, the NIR bands from VISTA/VIRCam and from CFHT/WIRCam, the four channels from Spitzer/IRAC and the 24 μm Spitzer/MIPS. The circles are 7.65 arcsec in radius and are encompassing the complete nebula.

In the text
thumbnail Fig. 2.

Left: contour of LAB1 on the 15″ × 15″ r-band image. North is up, east is left. The identified bright core is encircled in blue. The filamentary structure is encircled in red. The second image shows the r − g bands subtraction. The third image is the Bw-band of object 59 of Prescott et al. (2012). Right: Bw-band of object 26 of Prescott et al. (2012). Both Bw-band stamps were extracted from Fig. 10 of Prescott et al. (2012).

In the text
thumbnail Fig. 3.

Best-fit galaxy template in observed wavelengths from LePhare including the emission line contribution. Twelve filters were used to estimate these photometric redshift. The filled squares represent the ISOCOR magnitudes of our object, and the empty circles show the magnitudes from the best model. The downward facing triangles represent the detection limit used for the photometric redshift estimation. Top: for the complete nebula, magnitude of aperture 35 pixels. Centre: for the bright core, with 3 arcsec aperture, 3sigma detection. Bottom: for the faint halo with 15 pixels aperture. The three fits give a redshift of around 3.3 ± 0.05.

In the text
thumbnail Fig. 4.

Attempt to estimate the equivalent width of Lyα based on the variation of the line ratios during SED fitting. The two vertical dashed lines give the position of the two local minima. The y-axis has been reversed for clarity and we choose the emission line to be negative.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.