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 present a WFC3 F160W ($H$-band) selected catalog in the CANDELS/GOODS-N field containing photometry from the ultraviolet (UV) to the far-infrared (IR), photometric redshifts and stellar parameters derived from the analysis of the multi-wavelength data. The catalog contains 35,445 sources over the 171 arcmin$^{2}$ of the CANDELS F160W mosaic. The 5$\sigma$ detection limits (within an aperture of radius 0\farcs17) of the mosaic range between $H=27.8$, 28.2 and 28.7 in the wide, intermediate and deep regions, that span approximately 50\%, 15\% and 35\% of the total area. The multi-wavelength photometry includes broad-band data from UV (U band from KPNO and LBC), optical (HST/ACS F435W, F606W, F775W, F814W, and F850LP), near-to-mid IR (HST/WFC3 F105W, F125W, F140W and F160W, Subaru/MOIRCS Ks, CFHT/Megacam K, and \spitzer/IRAC 3.6, 4.5, 5.8, 8.0 $\mu$m) and far IR (\spitzer/MIPS 24$\mu$m, HERSCHEL/PACS 100 and 160$\mu$m, SPIRE 250, 350 and 500$\mu$m) observations. In addition, the catalog also includes optical medium-band data (R$\sim50$) in 25 consecutive bands, $\lambda=500$ to 950~nm, from the SHARDS survey and WFC3 IR spectroscopic observations with the G102 and G141 grisms (R$\sim210$ and 130). The use of higher spectral resolution data to estimate photometric redshifts provides very high, and nearly uniform, precision from $z=0-2.5$. The comparison to 1,485 good quality spectroscopic redshifts up to $z\sim3$ yields $\Delta z$/(1+$z_{\rm spec}$)$=$0.0032 and an outlier fraction of $\eta=$4.3\%. In addition to the multi-band photometry, we release added-value catalogs with emission line fluxes, stellar masses, dust attenuations, UV- and IR- based star formation rates and rest-frame colors.

We present a sample of 33 Balmer Break Galaxies (BBGs) selected as HST/F160W dropouts in the deepest CANDELS/GOODS fields ($H\gtrsim27.3$~mag) but relatively bright in {\it Spitzer}/IRAC ($[3.6],[4.5]<24.5$~mag), implying red colors (median and quartiles: $\langle H-[3.6]\rangle=3.1^{3.4}_{2.8}$\,mag). Half of these BBGs are newly identified sources. Our BBGs are massive ($\langle \log(\rm{M}/\rm{M}_\odot)\rangle=10.8$) high redshift ($\langle z\rangle=4.8$) dusty ($\langle \rm{A(V)}\rangle=2.0$~mag) galaxies. The SEDs of half of our sample indicate that they are star-forming galaxies with typical specific SFRs 0.5-1.0~Gyr$^{-1}$, qualifying them as main sequence (MS) galaxies at $3<z<6$. One third of those SEDs indicates the presence of prominent emission lines (H$\beta$+$[OIII]$, H$\alpha$$+$[NII]) boosting the IRAC fluxes and red colors. Approximately 20\% of the BBGs are very dusty ($\rm{A(V)}\sim2.5$~mag) starbursts with strong mid-to-far infrared detections and extreme SFRs ($\rm{SFR}>10^{3}\,\rm{M}_\odot/yr$) that place them above the MS. The rest, 30\%, are post-starbursts or quiescent galaxies located $>2\sigma$ below the MS with mass-weighted ages older than 700~Myr. Only 2 of the 33 galaxies are X-ray detected AGN with optical/near-infrared SEDs dominated by stellar emission, but the presence of obscured AGN in the rest of sources cannot be discarded. Our sample accounts for 8\% of the total number density of $\log(\rm{M}/\rm{M}_\odot)>10$ galaxies at $z>3$, but it is a significant contributor (30\%) to the general population of red $\log(\rm{M}/\rm{M}_\odot)>11$ galaxies at $4<z<6$. Finally, our results point out that 1 of every 30 massive $\log(\rm{M}/\rm{M}_\odot)>11$ galaxies in the local Universe was assembled in the first 1.5~Gyr after the Big Bang, a fraction that is not reproduced by state-of-the-art galaxy formation simulations.

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.

We make use of SHARDS, an ultra-deep (<26.5AB) galaxy survey that provides optical photo-spectra at resolution R~50, via medium band filters (FWHM~150A). This dataset is combined with ancillary optical and NIR fluxes to constrain the dust attenuation law in the rest-frame NUV region of star-forming galaxies within the redshift window 1.5<z<3. We focus on the NUV bump strength (B) and the total-to-selective extinction ratio (Rv), targeting a sample of 1,753 galaxies. By comparing the data with a set of population synthesis models coupled to a parametric dust attenuation law, we constrain Rv and B, as well as the colour excess, E(B-V). We find a correlation between Rv and B, that can be interpreted either as a result of the grain size distribution, or a variation of the dust geometry among galaxies. According to the former, small dust grains are associated with a stronger NUV bump. The latter would lead to a range of clumpiness in the distribution of dust within the interstellar medium of star-forming galaxies. The observed wide range of NUV bump strengths can lead to a systematic in the interpretation of the UV slope ($\beta$) typically used to characterize the dust content. In this study we quantify these variations, concluding that the effects are $\Delta\beta$~0.4.

We present Star Formation Histories (SFHs) for a sample of 104 massive (stellar mass M > 10$^{10}$ M$_{\odot}$) quiescent galaxies (MQGs) at $z$=1.0-1.5 from the analysis of spectro-photometric data from the SHARDS and HST/WFC3 G102 and G141 surveys of the GOODS-N field, jointly with broad-band observations from ultraviolet (UV) to far-infrared (Far-IR). The sample is constructed on the basis of rest-frame UVJ colours and specific star formation rates (sSFR=SFR/Mass). The Spectral Energy Distributions (SEDs) of each galaxy are compared to models assuming a delayed exponentially declining SFH. A Monte Carlo algorithm characterizes the degeneracies, which we are able to break taking advantage of the SHARDS data resolution, by measuring indices such as MgUV and D4000. The population of MQGs shows a duality in their properties. The sample is dominated (85%) by galaxies with young mass-weighted ages, tM $<$ 2 Gyr, short star formation timescales, $\langle \tau \rangle$ $\sim$ 60-200 Myr, and masses log(M/M$_{\odot}$) $\sim$ 10.5. There is an older population (15%) with tM$=$2 - 4 Gyr, longer star formation timescales, $\langle \tau \rangle$$\sim$ 400 Myr, and larger masses, log(M/M$_{\odot}$) $\sim$ 10.7. The SFHs of our MQGs are consistent with the slope and the location of the Main Sequence (MS) of star-forming galaxies at $z$ $>$ 1.0, when our galaxies were 0.5-1.0 Gyr old. According to these SFHs, all the MQGs experienced a Luminous Infrared Galaxy (LIRG) phase that lasts for $\sim$ 500 Myr, and half of them an Ultra Luminous Infrared Galaxy (ULIRG) phase for $\sim$ 100 Myr. We find that the MQG population is almost assembled at $z$ $\sim$ 1, and continues evolving passively with few additions to the population.

In this paper, we present a comprehensive analysis of star-forming galaxies (SFGs) at intermediate redshifts (z~1). We combine the ultra-deep optical spectro-photometric data from the Survey for High-z Absorption Red and Dead Sources (SHARDS) with deep UV-to-FIR observations in the GOODS-N field. Exploiting two of the 25 SHARDS medium-band filters, F687W17 and F823W17, we select [OII] emission line galaxies at z~0.84 and z~1.23 and characterize their physical properties. Their rest-frame equivalent widths (EW$_{\mathrm{rf}}$([OII])), line fluxes, luminosities, star formation rates (SFRs) and dust attenuation properties are investigated. The evolution of the EW$_{\mathrm{rf}}$([OII]) closely follows the SFR density evolution of the universe, with a trend of EW$_{\mathrm{rf}}$([OII])$\propto$(1+z)$^3$ up to redshift z~1, followed by a possible flattening. The SF properties of the galaxies selected on the basis of their [OII] emission are compared with complementary samples of SFGs selected by their MIR and FIR emission, and also with a general mass-selected sample of galaxies at the same redshifts. We demonstrate observationally that the UVJ diagram (or, similarly, a cut in the specific SFR) is only partially able to distinguish the quiescent galaxies from the SFGs. The SFR-M$_*$ relation is investigated for the different samples, yelding a logarithmic slope ~1, in good agreement with previous results. The dust attenuations derived from different SFR indicators (UV(1600), UV(2800), [OII], IR) are compared and show clear trends with respect to both the stellar mass and total SFR, with more massive and highly star-forming galaxies being affected by stronger dust attenuation.