Search

This is an Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited., Context. Galaxies with extreme emission lines (EELGs) may play a key role in the evolution of the Universe, as well as in our understanding of the star formation process itself. For this reason an accurate determination of their spatial density and fundamental properties in different epochs of the Universe will constitute a unique perspective towards a comprehensive picture of the interplay between star formation and mass assembly in galaxies. In addition to this, EELGs are also interesting in order to explain the reionization of the Universe, since their interstellar medium (ISM) could be leaking ionizing photons, and thus they could be low z, analogous of extreme galaxies at high z. Aims. This paper presents a method to obtain a census of EELGs over a large area of the sky by detecting galaxies with rest-frame equivalent widths ≥300 Å in the emission lines [O II]λλ3727,3729Å, [O III]λ5007Å, and Hα. For this, we aim to use the J-PAS survey, which will image an area of ≈8000 deg2 with 56 narrow band filters in the optical. As a pilot study, we present a methodology designed to select EELGs on the miniJPAS images, which use the same filter dataset as J-PAS, and thus will be exportable to this larger survey. Methods. We make use of the miniJPAS survey data, conceived as a proof of concept of J-PAS, and covering an area of ≈1 deg2. Objects were detected in the rSDSS images and selected by imposing a condition on the flux in a given narrow-band J-PAS filter with respect to the contiguous ones, which is analogous to requiring an observed equivalent width larger than 300 Å in a certain emission line within the filter bandwidth. The selected sources were then classified as galaxies or quasi-stellar objects (QSOs) after a comparison of their miniJPAS fluxes with those of a spectral database of objects known to present strong emission lines. This comparison also provided a redshift for each source, which turned out to be consistent with the spectroscopic redshifts when available (|Δz/(1 + zspec)| ≤ 0.01). Results. The selected candidates were found to show a compact appearance in the optical images, some of them even being classified as point-like sources according to their stellarity index. After discarding sources classified as QSOs, a total of 17 sources turned out to exhibit EW0 ≥ 300 Å in at least one emission line, thus constituting our final list of EELGs. Our counts are fairly consistent with those of other samples of EELGs in the literature, although there are some differences, which were expected due to biases resulting from different selection criteria. © J. Iglesias-Páramo et al. 2022., This work has been partially funded by projects PID2019-107408GB-C44 from the Spanish PNAYA, co-funded with FEDER, and grand P18-FR-2664, funded by Junta de Andalucía. We acknowledge financial support from the State Agency for Research of the Spanish MCIU through the “Center of Excellence Severo Ochoa” award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709). RGD and LADG acknowledge financial support from the State Agency for Research of the Spanish MCIU through the “Center of Excellence Severo Ochoa” award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709), and PID2019-109067-GB100. IM acknowledges financial support from the State Agency for Research of the Spanish MCIU through the PID2019-106027GB-C41. JCM acknowledges partial support from the Spanish Ministry of Science, Innovation and Universities (MCIU/AEI/FEDER, UE) through the grant PGC2018-097585-B-C22. SDP is grateful to the Fonds de Recherche du Québec – Nature et Technologies. LSJ acknowledges the support of CNPq (304819/2017-4) and FAPESP (2019/10923-5). JAFO acknowledges the financial support from the Spanish Ministry of Science and Innovation and the European Union – NextGenerationEU through the Recovery and Resilience Facility project ICTS-MRR-2021-03- CEFCA. Funding for the J-PAS Project has been provided by the Governments of España and Aragón though the Fondo de Inversión de Teruel, European FEDER funding and the MINECO and by the Brazilian agencies FINEP, FAPESP, FAPERJ and by the National Observatory of Brazil. Based on observations made with the JST/T250 telescope and PathFinder camera for the miniJPAS project at the Observatorio Astrofísico de Javalambre (OAJ), in Teruel, owned, managed, and operated by the Centro de Estudios de Física del Cosmos de Aragón (CEFCA). We acknowledge the OAJ Data Processing and Archiving Unit (UPAD) for reducing and calibrating the OAJ data used in this work. Funding for OAJ, UPAD, and CEFCA has been provided by the Governments of Spain and Aragón through the Fondo de Inversiones de Teruel; the Aragón Government through the Research Groups E96, E103, and E16_17R; the Spanish Ministry of Science, Innovation and Universities (MCIU/AEI/FEDER, UE) with grant PGC2018-097585-B-C21; the Spanish Ministry of Economy and Competitiveness (MINECO/FEDER, UE) under AYA2015-66211-C2-1-P, AYA2015-66211-C2-2, AYA2012-30789, and ICTS-2009-14; and European FEDER funding (FCDD10-4E-867, FCDD13-4E-2685). This research has made use of the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. Funding for SDSS-IV has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, and the U.S. Department of Energy Office of Science. The SDSS-IV web site is https://www.sdss.org/. This project has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 898633.

(Work partially presented in Lumbreras-Calle et al. 2021) Extreme emission line galaxies (EELGs) constitute a population of unique systems undergoing very intense events of star formation, which are very rare in the local Universe. They nevertheless hold clues into the limits of star formation processes and can be considered analogs of the first galaxies in the Universe, given the physical properties they share. To find EELGs, we have used the unprecedented dataset of the second data release of the J-PLUS survey, with 2000 square degrees of the northern sky observed from the Observatorio Astronómico de Javalambre (Teruel) with a set of narrow, medium and broad-band filters. We select as emitters objects with excess of flux in mediumband filters, and remove contaminants using WISE and J-PLUS photometry. We have found 466 EELGs selected by their [OIII]5007 emission, with 410 of them previously unclassified as such. The derived [OIII]/Hα ratios show a wide range, more similar to previous simulations of reionization-era galaxies than Green Peas or Blueberries. The strong [OIII] and high [OIII]/[OII] ratios can be considered proxies for high ionization parameter and potential Lyman continuum leakage. Thanks to the use of mediumband filterrs, the volume density of the EELGs in this work is more than one order of magnitude higher than previous broadband-selected samples, and similar to magnitude-limited spectroscopic surveys. Additionally, we find 310 extreme Hα emitters at z~0.3, out of which 271 are newly discovered. Finally, we have obtained longslit spectroscopic data for a subsample these EELGs. Preliminary analysis of the spectrashowthe accuracy of the sample selection, and the low metallicity and strong ionizing spectrum of the sources.

We present the third data release of the Javalambre Local Universe Photometric Survey (J-PLUS) covering more than 3200 square degrees of the northern sky. J-PLUS is an ongoing photometric survey with 12 optical bands observing thousands of square degrees of the sky from the JAST80 telescope at the Observatorio Astrofísico de Javalambre (OAJ). T80Cam is a 2 deg² field-of-view (FoV) camera mounted on JAST80, and is equipped with a unique system of filters spanning the entire optical range (3500 − 10000 Å), optimally designed to extract the rest-frame spectral features that are key to both characterize stellar types and to deliver a low-resolution photo-spectrum for each observed object. With a typical depth (5σ in 3 arcsec aperture) of AB ∼ 20.7 mag per band, data containing about 15 million stars and 10 million galaxies at r < 21 mag are going to be released by the end of 2022. We also present announce some relevant scientific results obtained with J-PLUS., Missing links to J-PLUS e-posters added

We use a sample of 706 galaxies, selected as [O ii]λ3727 ([O ii]) emitters in the Survey for High-z Absorption Red and Dead Sources (SHARDS) on the CANDELS/GOODS-N field, to study the differential attenuation of the nebular emission with respect to the stellar continuum. The sample includes only galaxies with a counterpart in the infrared and log10(M∗/M⊙) > 9, over the redshift interval 0.3 ≲ z ≲ 1.5. Our methodology consists in the comparison of the star formation rates inferred from [O ii] and Hα emission lines with a robust quantification of the total star-forming activity (SFR TOT) that is independently estimated based on both infrared and ultraviolet (UV) luminosities. We obtain $f\, =\, E(B-V)_{\mathrm{stellar}}$/E(B - V)nebular = 0.69$^{0.71}_{0.69}$ and 0.55$^{0.56}_{0.53}$ for [O ii] and Hα, respectively. Our resulting f-factors display a significant positive correlation with the UV attenuation and shallower or not-significant trends with the stellar mass, the SFRTOT, the distance to the main sequence, and the redshift. Finally, our results favour an average nebular attenuation curve similar in shape to the typical dust curve of local starbursts., LR-M thanks Casiana Muñoz Tuñón for useful suggestions on this manuscript. LR-M, GR, and AF acknowledge the support from grant PRIN MIUR 2017-20173ML3WW_001. LR-M also acknowledges funding support from the Università degli Studi di Padova - Dipartimento di Fisica e Astronomia ‘G. Galilei’. PGP-G and LC acknowledge support from Spanish Ministerio de Ciencia, Innovación y Universidades through grant PGC2018-093499-B-I00. LC acknowledges also financial support from Comunidad de Madrid under Atracción de Talento grant 2018-T2/TIC-11612. AVG acknowledges support from the European Research Council through the Advanced Grant MIST (FP7/2017-2020, No. 742719). AP gratefully acknowledges financial support from STFC through grants ST/T000244/1 and ST/P000541/1. Analyses were performed in R 3.6.1 (R Core Team 2019).

We use a sample of 706 galaxies, selected as [OII]$\lambda$3727 ([OII]) emitters in the Survey for High-$z$ Absorption Red and Dead Sources (SHARDS) on the CANDELS/GOODS-N field, to study the differential attenuation of the nebular emission with respect to the stellar continuum. The sample includes only galaxies with a counterpart in the infrared and $\mathrm{log}_{10}(M_{*}/\mathrm{M}_{\odot})$ $>$ 9, over the redshift interval 0.3 $\lesssim$ $z$ $\lesssim$ 1.5. Our methodology consists in the comparison of the star formation rates inferred from [OII] and H$\alpha$ emission lines with a robust quantification of the total star-forming activity (${SFR}_{\mathrm{TOT}}$) that is independently estimated based on both infrared and ultraviolet (UV) luminosities. We obtain $f$$=$$E(B-V)_{\mathrm{stellar}}$/$E(B-V)_{\mathrm{nebular}}$ $=$ 0.69$^{0.71}_{0.69}$ and 0.55$^{0.56}_{0.53}$ for [OII] and H$\alpha$, respectively. Our resulting $f$-factors display a significant positive correlation with the UV attenuation and shallower or not-significant trends with the stellar mass, the $SFR_{\mathrm{TOT}}$, the distance to the main sequence, and the redshift. Finally, our results favour an average nebular attenuation curve similar in shape to the typical dust curve of local starbursts., Comment: 23 pages, 10 figures, 3 tables; accepted for publication in MNRAS

In this series of papers, we employ several machine learning (ML) methods to classify the point-like sources from the miniJPAS catalogue, and identify quasar candidates. Since no representative sample of spectroscopically confirmed sources exists at present to train these ML algorithms, we rely on mock catalogues. In this first paper we develop a pipeline to compute synthetic photometry of quasars, galaxies and stars using spectra of objects targeted as quasars in the Sloan Digital Sky Survey. To match the same depths and signal-to-noise ratio distributions in all bands expected for miniJPAS point sources in the range $17.5\leq r, 20 pages, 18 figures, submitted to MNRAS

Aims. Our goal is to estimate the star formation main sequence (SFMS) and the star formation rate density (SFRD) at z = 0.017 (d . 75 Mpc) using the Javalambre Photometric Local Universe Survey (J-PLUS) first data release, that probes 897.4 deg2 with twelve optical bands. Methods. We extract the Ha emission flux of 805 local galaxies from the J-PLUS filter J0660, being the continuum level estimated with the other eleven J-PLUS bands, and the dust attenuation and nitrogen contamination corrected with empirical relations. Stellar masses (M), Ha luminosities (LHa), and star formation rates (SFRs) were estimated by accounting for parameters covariances. Our sample comprises 689 blue galaxies and 67 red galaxies, classified in the (u - g) vs. (g - z) color–color diagram, plus 49 AGN. Results. The SFMS is explored at log M 8 and it is clearly defined by the blue galaxies, with the red galaxies located below them. The SFMS is described as log SFR = 0.83 log M -8.44. We find a good agreement with previous estimations of the SFMS, especially those based on integral field spectroscopy. The Ha luminosity function of the AGN-free sample is well described by a Schechter function with log L * Ha = 41.34, log f * = -2.43, and a = -1.25. Our measurements provide a lower characteristic luminosity than several previous studies in the literature. Conclusions. The derived star formation rate density at d . 75 Mpc is log SFR = -2.10 ± 0.11, with red galaxies accounting for 15% of the SFRD. Our value is lower than previous estimations at similar redshift, and provides a local reference for evolutionary studies regarding the star formation history of the Universe © 2021 EDP Sciences. All rights reserved., Based on observations made with the JAST/T80 telescope for J-PLUS project at the Observatorio Astrofísico de Javalambre in Teruel, a Spanish Infraestructura Cientifico-Técnica Singular (ICTS) owned, managed and operated by the Centro de Estudios de Física del Cosmos de Aragón (CEFCA). Data has been processed and provided by CEFCA’s Unit of Processing and Archiving Data (UPAD). Funding for the J-PLUS Project has been provided by the Governments of Spain and Aragón through the Fondo de Inversiones de Teruel; the Aragón Government through the Research Groups E96, E103, and E16_17R; the Spanish Ministry of Science, Innovation and Universities (MCIU/AEI/FEDER, UE) with grants PGC2018-097585-B-C21 and PGC2018-097585-B-C22; the Spanish Ministry of Economy and Competitiveness (MINECO) under AYA2015-66211-C2-1-P, AYA2015-66211-C2-2, AYA2012-30789, and ICTS-2009-14; and European FEDER funding (FCDD10- 4E-867, FCDD13-4E-2685). The Brazilian agencies FINEP, FAPESP, and the National Observatory of Brazil have also contributed to this project. J.M.V. acknowledges financial support from the State Agency for Research of the Spanish MCIU through the “Center of Excellence Severo Ochoa” award to the Instituto de Astrofísica de Andalucía (SEV 2017-0709), and grant AYA2016-79724-C4-4P. R.A.D. acknowledges support from the Conselho Nacional de Desenvolvimento Científico e Tecnolágico –CNPq through BP grant 308105/2018-4, and the Financiadora de Estudos e Projetos – FINEP grants REF. 1217/13 – 01.13.0279.00 and REF 0859/10 – 01.10.0663.00 for hardware support for the J-PLUS project through the National Observatory of Brazil. We thank the help of Prof. G. Gavazzi in the analysis and understanding of Hα3 data. This research has made use of the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This research made use of Astropy, a community-developed core Python package for Astron omy (Astropy Collaboration 2013), and Matplotlib, a 2D graphics package used for Python for publication-quality image generation across user interfaces and operating systems (Hunter 2007). Some of the data sets discussed in the manuscript were recovered with WebPlotDigitizer (https://automeris.io/WebPlotDigitizer/index.html), a free software developed to facilitate easy and accurate data extraction from a variety of plot types.

The differences between the inherent stellar populations (SPs) of LAEs and LBGs are a key factor in understanding early galaxy formation and evolution. We have run a set of SP burst-like models for a sample of 1,558 sources at $3.4, Comment: Accepted for publication in MNRAS. 18 pages, 23 figures

We have undertaken a comprehensive search for both Lyman Alpha Emitters (LAEs) and Lyman Break Galaxies (LBGs) in the SHARDS Survey of the GOODS-N field. SHARDS is a deep imaging survey, made with the 10.4 m Gran Telescopio Canarias (GTC), employing 25 medium band filters in the range from 500 to 941 nm. This is the first time that both LAEs and LBGs are surveyed simultaneously in a systematic way in a large field. We draw a sample of 1558 sources; 528 of them are LAEs. Most of the sources (1434) show rest-frame UV continua. A minority of them (124) are pure LAEs with virtually no continuum detected in SHARDS. We study these sources from $z\sim3.35$ up to $z\sim6.8$, well into the epoch of reionization. Note that surveys done with just one or two narrow band filters lack the possibility to spot the rest-frame UV continuum present in most of our LAEs. We derive redshifts, Star Formation Rates (SFRs), Ly$\alpha$ Equivalent Widths (EWs) and Luminosity Functions (LFs). Grouping within our sample is also studied, finding 92 pairs or small groups of galaxies at the same redshift separated by less than 60 comoving kpc. In addition, we relate 87 and 55 UV-selected objects with two known overdensities at $z=4.05$ and $z=5.198$, respectively. Finally, we show that surveys made with broad band filters are prone to introduce many unwanted sources ($\sim20$% interlopers), which means that previous studies may be overestimating the calculated LFs, specially at the faint end., Comment: Accepted for publication in MNRAS. 18 pages, 19 figures

The morphological evolution of star-forming galaxies provides important clues to understand their physical properties, as well as the triggering and quenching mechanisms of star formation. We aim at connecting morphology and star-formation properties of low-mass galaxies (median stellar mass $\sim$ 10$^{8.5}$ M$_{\odot}$) at low redshift ($z, Comment: Accepted for publication in Astronomy & Astrophysics. 13 pages, 11 figures

The physical processes driving the evolution of star formation (SF) in galaxies over cosmic time still present many open questions. Recent galaxy surveys allow now to study these processes in great detail at intermediate redshift. In this work, we build a complete sample of star-forming galaxies and analyze their properties, reaching systems with low stellar masses and low star formation rates (SFRs) at intermediate-to-low redshift. We use data from the SHARDS multiband survey in the GOODS-North field. Its depth (up to magnitude $\langle m_{3\sigma}\rangle\sim26.5$) and its spectro-photometric resolution ($R\sim50$) provides us with an ideal dataset to search for emission line galaxies (ELGs). We develop a new algorithm to identify low-redshift ($z$, Comment: Accepted for publication in Astronomy & Astrophysics. 17 pages, 17 figures

Context. Anti-truncated lenticular galaxies (Type-III S0s) present tight scaling relations between their surface brightness photometric and structural parameters. Although several evolutionary models have been proposed for the formation of these structures, the observations of Type-III S0 galaxies are usually limited to the local Universe. Aims. We aim to compare the properties of Type-III discs in a sample of S0 galaxies at 0.2 < z < 0.6with those of the local Universe. In this paper, we study the evolution of the photometric and structural scaling relations measured in the rest-frame R-band with z and the possible differences between the rest-frame (B - R) colours of the inner and outer disc profiles. Methods. We make use of a sample of 14 Type-III E/ S0-S0 galaxies at 0.2 < z < 0 : 6 from the GOODS-N field identified and characterised in a previous paper. We study whether or not the correlations found in local Type-III S0 galaxies were present similar to 6 Gyr ago. We analyse the distribution of the surface brightness characteristic parameters (R-break, mu(break), h(i), h(o), mu(0;) i and mu(0; o)) as a function of the stellar mass and look to see if there is a significant change with z. We also derive their rest-frame (B - R) colour profiles. Finally, we compare these results with the predictions from a grid of SSP models. Results. We find that the inner and outer scale-lengths of Type-III S0 galaxies at 0.4 < z < 0.6follow compatible trends and scaling relations with those observed in local S0 galaxies as a function of the break radius, R-break. We do not detect any significant differences between the location of R-break between z similar to 0.6 and z similar to 0 for a fixed stellar mass of the object, whereas the surface brightness at the break radius mu(break) is similar to 1.5 mag arcsec 2 dimmer in the local Universe than at z similar to 0.6 for a fixed stellar mass. We find no significant differences in the (B - R) colour between the inner and outer profiles of the Type-III S0 galaxies at 0.2 < z < 0.6. Conclusions. In contrast to Type-II (down-bending) profiles, the anti-truncated surface brightness profiles of S0 galaxies present compatible R-break values and scaling relations during the last 6 Gyr. This result and the similarity of the colours of the inner and outer discs point to a highly scalable and stable formation process, probably more related to gravitational and dynamical processes than to the evolution of stellar populations.

International audience

Free access at publisher's website https://www.aanda.org/articles/aa/pdf/2019/01/aa34312-18.pdf., he Hubble Ultra Deep field (HUDF) is the deepest region ever observed with the Hubble Space Telescope. With the main objective of unveiling the nature of galaxies up to z ∼ 7 - 8, the observing and reduction strategy have focused on the properties of small and unresolved objects, rather than the outskirts of the largest objects, which are usually over-subtracted. Aims. We aim to create a new set of WFC3 IR mosaics of the HUDF using novel techniques to preserve the properties of the low surface brightness regions. Methods. We created ABYSS: A pipeline that optimises the estimate and modelling of low-level systematic effects to obtain a robust background subtraction. We have improved four key points in the reduction: 1) creation of new absolute sky flat fields, 2) extended persistence models, 3) dedicated sky background subtraction and 4) robust co-adding. Results. The new mosaics successfully recover the low surface brightness structure removed on the previous HUDF published reductions. The amount of light recovered with a mean surface brightness dimmer than μ = 26 mag arcsec -2 is equivalent to a m = 19 mag source when compared to the XDF and a m = 20 mag compared to the HUDF12. Conclusions. We present a set of techniques to reduce ultra-deep images (μ > 32.5 mag arcsec -2 , 3σ in 10 × 10 arcsec boxes), that successfully allow us to detect the low surface brightness structure of extended sources on ultra deep surveys. The developed procedures are applicable to HST, JWST, Euclid and many other space and ground-based observatories. © ESO 2019., The authors would like to thank the XDF and HUDF12 teams for their extraordinary work in which this contribution is based on, truly standing on the shoulders of giants. This work would not have been possible without the kind assistance of all the members of the STScI Help Desk. We specially thank Knox Long by its work on persistence effects and its extraordinary support during this work. A.B. acknowledges support from the International Space Science Institute (ISSI) as part of the Exploring the Ultra-Low Surface Brightness Universe International team. Some/all of the data presented in this paper were obtained from the Mikulski Archive for Space Telescopes (MAST). STScI is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. I.T. and C.M.L acknowledges the support from the SUNDIAL EU Network and from the European Union’s Horizon 2020 research and innovation programme under Marie Sklodowska-Curie grant agreement No 721463 and the SUNDIAL ITN network. C.M.L. acknowledges financial support from the Spanish Ministry of Economy and Competitiveness (MINECO) under grant number AYA2016-76219-P. N.C acknowledged support from the Spanish Programa Nacional de Astronimía y Astrofísica under grant AYA2016-75808-R. C.G.G. acknowledges support from the European Research Council (ERC) Consolidator Grant funding scheme (project ConTExt, grant number 648179). Support for MAST for non-HST data is provided by the NASA Office of Space Science via grant NNX09AF08G and by other grants and contracts. This work was partly done using GNU Astronomy Utilities Gnuastro version 0.5. Gnuastro is a generic package for astronomical data manipulation and analysis which was primarily created and developed for research funded by the Monbukagakusho (Japanese government) scholarship and ERC advanced grant 339659-MUSICOS. Support for this work was provided by the Spanish Ministe-rio de Economía y Competitividad (MINECO; grant AYA 2016-77237-C3-I-P). This research made use of Astropy (http://www.astropy.org) a community-developed core Python package for Astronomy (Astropy Collaboration 2013; Price-Whelan et al. 2018).

Type-III S0 galaxies present tight scaling relations between their surface brightness photometric and structural parameters. Several evolutionary models have been proposed for the formation of Type-III S0 galaxies but the observations of are usually limited to the local Universe. We study the evolution of the photometric and structural scaling relations found between the parameters of the surface brightness profiles in the rest-frame R-band of Type-III S0 galaxies with z and the possible differences between the rest-frame (B-R) colours of the inner and outer disc profiles. We make use of a sample of 14 Type-III E/S0--S0 galaxies at 0.2, Comment: Accepted for publishing in Astronomy & Astrophysics. 28 pages, 27 figures. Language editor corrections applied

We have undertaken a comprehensive search for both Lyman Alpha Emitters (LAEs) and Lyman Break Galaxies (LBGs) in the SHARDS Survey of the GOODS-N field. SHARDS is a deep imaging survey, made with the 10.4 m Gran Telescopio Canarias (GTC), employing 25 medium band filters in the range from 500 to 941 nm. This is the first time that both LAEs and LBGs are surveyed simultaneously in a systematic way in a large field. We draw a sample of 1558 sources; 528 of them are LAEs. Most of the sources (1434) show rest-frame UV continua. A minority of them (124) are pure LAEs with virtually no continuum detected in SHARDS. We study these sources from z ? 3.35 up to z ? 6.8, well into the epoch of reionization. Note that surveys done with just one or two narrow band filters lack the possibility to spot the rest-frame UV continuum present in most of our LAEs. We derive redshifts, Star Formation Rates (SFRs), Ly? Equivalent Widths (EWs) and Luminosity Functions (LFs). Grouping within our sample is also studied, finding 92 pairs or groups of galaxies at the same redshift separated by less than 60 comoving kpc. In addition, we relate 87 and 55 UV-selected objects with two known overdensities at z = 4.05 and z = 5.198, respectively. Finally, we show that surveys made with broad band filters are prone to introduce many unwanted sources (?20% interlopers), which means that previous studies may be overestimating the calculated LFs, specially at the faint end.