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<z<6.8$ from the Survey for High-$z$ Absorption Red and Dead Sources (SHARDS) over the GOODS-N field. This work focuses on the differences between the three different observational subfamilies of our sample: LAE-LBGs, no-Ly$\alpha$ LBGs and pure LAEs. Single and double SP synthetic spectra were used to model the SEDs, adopting a Bayesian information criterion to analyse under which situations a second SP is required. We find that the sources are well modelled using a single SP in $\sim79\%$ of the cases. The best models suggest that pure LAEs are typically young low mass galaxies ($t\sim26^{+41}_{-25}$ Myr; $M_{\mathrm{star}}\sim5.6^{+12.0}_{-5.5}\times10^{8}\ M_{\odot}$), undergoing one of their first bursts of star formation. On the other hand, no-Ly$\alpha$ LBGs require older SPs ($t\sim71\pm12$ Myr), and they are substantially more massive ($M_{\mathrm{star}}\sim3.5\pm1.1\times10^{9}\ M_{\odot}$). LAE-LBGs appear as the subgroup that more frequently needs the addition of a second SP, representing an old and massive galaxy caught in a strong recent star-forming episode. The relative number of sources found from each subfamily at each $z$ supports an evolutionary scenario from pure LAEs and single SP LAE-LBGs to more massive LBGs. Stellar Mass Functions are also derived, finding an increase of $M^{*}$ with cosmic time and a possible steepening of the low mass slope from $z\sim6$ to $z\sim5$ with no significant change to $z\sim4$. Additionally, we have derived the SFR-$M_{\mathrm{star}}$ relation, finding a $\mathrm{SFR}\propto M_{\mathrm{star}}^{\beta}$ behaviour with negligible evolution from $z\sim4$ to $z\sim6$.

We explore deep rest-frame UV to FIR data in the COSMOS field to measure the individual spectral energy distributions (SED) of the ~4000 SC4K (Sobral et al. 2018) Lyman-alpha (Lya) emitters (LAEs) at z~2-6. We find typical stellar masses of 10$^{9.3\pm0.6}$ M$_{\odot}$ and star formation rates (SFR) of SFR$_{SED}=4.4^{+10.5}_{-2.4}$ M$_{\odot}$/yr and SFR$_{Lya}=5.9^{+6.3}_{-2.6}$ M$_{\odot}$/yr, combined with very blue UV slopes of beta=-2.1$^{+0.5}_{-0.4}$, but with significant variations within the population. M$_{UV}$ and beta are correlated in a similar way to UV-selected sources, but LAEs are consistently bluer. This suggests that LAEs are the youngest and/or most dust-poor subset of the UV-selected population. We also study the Lya rest-frame equivalent width (EW$_0$) and find 45 "extreme" LAEs with EW$_0>240$ A (3 $\sigma$), implying a low number density of $(7\pm1)\times10^{-7}$ Mpc$^{-3}$. Overall, we measure little to no evolution of the Lya EW$_0$ and scale length parameter ($w_0$) which are consistently high (EW$_0=140^{+280}_{-70}$ A, $w_0=129^{+11}_{-11}$ A) from z~6 to z~2 and below. However, $w_0$ is anti-correlated with M$_{UV}$ and stellar mass. Our results imply that sources selected as LAEs have a high Lya escape fraction (f$_{esc, Lya}$) irrespective of cosmic time, but f$_{esc, Lya}$ is still higher for UV-fainter and lower mass LAEs. The least massive LAEs ($<10^{9.5}$ M$_{\odot}$) are typically located above the star formation "Main Sequence" (MS), but the offset from the MS decreases towards z~6 and towards $10^{10}$ M$_{\odot}$. Our results imply a lack of evolution in the properties of LAEs across time and reveals the increasing overlap in properties of LAEs and UV-continuum selected galaxies as typical star-forming galaxies at high redshift effectively become LAEs.

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.

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$<0.36) ELGs by detecting the [OIII]5007 and $H\alpha$ emission lines simultaneously. We fit the spectral energy distribution (SED) of the selected sample, using a model with two single stellar populations. We find 160 star-forming galaxies with equivalent widths (EWs) as low as 12 {\AA}, with median values for the sample of $\sim$ 35 {\AA} in [OIII]5007 and $\sim$ 56 {\AA} in $H\alpha$, respectively. Results from the SED fitting show a young stellar population with low median metallicity (36% of the solar value) and extinction ($A_V \sim$ 0.37), with median galaxy stellar mass $\sim$ 10$^{8.5}$ M$_{\odot}$. Gas-phase metallicities measured from available spectra are also low. ELGs in our sample present bluer colors in the UVJ plane than the median color-selected star-forming galaxy in SHARDS. We suggest a new (V-J) color criterion to separate ELGs from non-ELGs in blue galaxy samples. In addition, several galaxies present high densities of O-type stars. Robust fits to the full SEDs can only be obtained including an old stellar population, suggesting the young component is built up by a recent burst of SF in an otherwise old galaxy.