We present deep rest-frame UV spectroscopic observations using the Gran Telescopio Canarias of six gravitationally lensed Lya emitters (LAEs) at $2.36<z<2.82$ selected from the BELLS GALLERY survey. By taking the magnifications into account, we show that LAEs can be as luminous as L(Lya) = 30x10$^{42}$ erg s-1 and M(UV) = -23 (AB) without invoking an AGN component, in contrast with previous findings. We measure Lya rest-frame equivalent widths, EW(Lya), ranging from 16\AA to 50\AA and Lya escape fractions, fesc(Lya), from 10% to 40%. Large EW(Lya) and fesc(Lya) are found predominantly in LAEs showing weak low-ionization ISM absorption (EW < 1\AA) and narrow Lya profiles (< 300 km s-1 FWHM) with their peak close (< 80 km s-1) to their systemic redshifts, suggestive of less scatter from low HI column densities that favours the escape of Lya photons. We infer stellar metallicities of Z/Zsun ~ 0.2 in almost all LAEs by comparing the P-Cygni profiles of the wind lines NV1240\AA and CIV1549\AA with those from stellar synthesis models. We also find a trend between M(UV) and the velocity offset of ISM absorption lines, such as the most luminous LAEs experience stronger outflows. The most luminous LAEs show star formation rates up to 180 Msun yr-1, yet they appear relatively blue ($\beta$(UV) ~ -1.8 to -2.0) showing evidence of little dust attenuation (E(B-V) = 0.10-0.14). These luminous LAEs may be particular cases of young starburst galaxies that have had no time to form large amounts of dust. If so, they are ideal laboratories to study the early phase of massive star formation, stellar and dust mass growth, and chemical enrichment histories of starburst galaxies at high-z.

This work explores from a statistical point of view the rest-frame Far-ultraviolet (FUV) to Far-infrared (FIR) emission of a population of Lyman break galaxies (LBGs) at $z\sim3$ that cannot be individually detected from current FIR observations. We perform a stacking analysis over a sample of $\sim$17000 LBGs at redshift $2.5<z<3.5$ in the COSMOS field. The sample is binned as a function of UV luminosity ($L_{\mathrm{FUV}}$), UV continuum slope ($\beta_{\mathrm{UV}}$), and stellar mass (M$_{*}$), and then, stacked at optical ($BVriz$ bands), near-infrared ($YJHKs$ bands), IRAC (3.6, 4.5, 5.6 and 8.0 $\mu$m), MIPS (24$\mu$m), PACS (100 and 160~$\mu$m), SPIRE (250, 350, and 500~$\mu$m), and AzTEC (1.1mm) observations. We obtain thirty rest-frame FUV-to-FIR spectral energy distribution (SEDs) of LBGs at $z\sim3$, and analyse them with CIGALE SED-fitting analysis code. We are able to derive fully consistent physical parameters (M$_{*}$, $\beta_{\mathrm{UV}}$, $L_{\mathrm{FUV}}$, $L_{\mathrm{IR}}$, A$_{FUV}$, SFR, and slope of the dust attenuation law), and build a semi-empirical library of thirty rest-frame FUV-to-FIR stacked LBG SEDs as a function of $L_{\mathrm{FUV}}$, $\beta_{\mathrm{UV}}$, and M$_{*}$. We used the so-called IR-excess ($IRX \equiv L_{\mathrm{IR}} / L_{\mathrm{FUV}}$) to investigate the dust attenuation as a function of $\beta_{\mathrm{UV}}$ and M$_{*}$. Additionally, the SED-fitting analysis results provide a diversity of dust attenuation curve along the LBG sample, and their slope are well correlated with M$_{*}$. Stepper dust attenuations curves than Calzetti's are favoured in low stellar mass LBGs ($\log(M_{*} [M_{\odot}]) < 10.25$), while grayer ones are favoured in high stellar mass LBGs ($\log(M_{*} [M_{\odot}]) > 10.25$).

The James Webb Space Telescope will provide deep imaging and spectroscopy for sources at redshifts above 6, covering the Epoch of Reionization (EoR, 6 < z < 10). The Mid-IR instrument (MIRI) integral field spectrograph (MRS) will be the only instrument on board JWST able to observe the brightest optical emission lines H$\alpha$ and [OIII]0.5007$\mu$m at redshifts above 7 and 9, respectively. This paper presents a study of the H$\alpha$ fluxes predicted by FIRSTLIGHT cosmological simulations for galaxies at redshifts of 6.5 to 10.5, and its detectability with MIRI. Deep (40 ks) spectroscopic integrations with MRS will be able to detect (S/N > 5) EoR sources at redshifts above 7 with intrinsic star formation rates of more than 2 M$_{\odot}$ yr$^{-1}$, and stellar masses above 4-9 $\times$ 10$^7$ M$_{\odot}$. In addition, the paper presents realistic MRS simulated observations of the expected (rest-frame) optical and near-infrared spectra for some spectroscopically confirmed EoR sources detected by ALMA as [OIII]88$\mu$m emitters. The MRS simulated spectra cover a wide range of low metallicities from about 0.2 to 0.02Z$_{\odot}$, and different [OIII]88$\mu$m/[OIII]0.5007$\mu$m line ratios. The simulated 10ks MRS spectra show S/N in the range of 5 to 90 for H$\beta$, [OIII]0.4959,0.5007$\mu$m, H$\alpha$ and HeI1.083$\mu$m emission lines of MACS1149-JD1 at z = 9.11, independent of metallicity. In addition, deep 40 ks simulated spectra of the luminous merger candidate B14-65666 at z=7.15 shows the MRS capabilities of detecting, or putting strong upper limits, on the [NII]0.6584$\mu$m, [SII]0.6717,0.6731$\mu$m, and [SIII]0.9069,0.9532$\mu$m emission lines. In summary, MRS will enable the detailed study of key physical properties like internal extinction, instantaneous star formation, hardness of the ionising continuum, and metallicity, in bright (intrinsic or lensed) EoR sources.

$\mathit{Herschel}$ extragalactic surveys offer a unique opportunity to efficiently select a significant number of rare and massive dusty objects, and thus gain insight into the prodigious star-forming activity that takes place in the very distant Universe. To search for $z\geq4$ dusty star-forming galaxies, in this work we consider red SPIRE objects with fluxes rising from 250 $\mu$m to $500\:\mu$m (so-called "500 $\mu$m-risers"). We aim to implement a novel method to obtain a statistical sample of "500 $\mu$m-risers" and fully evaluate our selection inspecting different models of galaxy evolution. We consider one of the largest and deepest ${\it Herschel}$ surveys, the Herschel Virgo Cluster Survey. We develop a novel selection algorithm which links the source extraction and spectral energy distribution fitting. We select 133 "500 $\mu$m-risers" over 55 deg$^{2}$, imposing the criteria: $S_{500}>S_{350}>S_{250}$, $S_{250}>13.2$ mJy and $S_{500}>$30 mJy. Differential number counts are in a fairly good agreement with models, displaying better match than other existing samples. In order to interpret the statistical properties of selected sources, which has been proven as a very challenging task due the complexity of observed artefacts, we make end-to-end simulations including physical clustering and lensing. The estimated fraction of strongly lensed sources is $24^{+6}_{-5}\%$ based on models. We present the faintest known statistical sample of "500 $\mu$m-risers" and show that noise and strong lensing have crucial impact on measured counts and redshift distribution of selected sources. We estimate the flux-corrected star formation rate density at $4<z<5$ with the "500 $\mu$m-risers" and found it close to the total value measured in far-infrared. It indicates that colour selection is not a limiting effect to search for the most massive, dusty $z>4$ sources.

In this work we investigate the far-UV to NIR shape of the dust attenuation curve of a sample of IR selected dust obscured (U)LIRGs at z$\sim$2. The spectral energy distributions (SEDs) are fitted with CIGALE, a physically-motivated spectral synthesis model based on energy balance. Its flexibility allows us to test a wide range of different analytical prescriptions for the dust attenuation curve, including the well-known Calzetti and Charlot & Fall curves, and modified versions of them. The attenuation curves computed under the assumption of our reference double power-law model are in very good agreement with those derived, in previous works, with radiative transfer (RT) SED fitting. We investigate the position of our galaxies in the IRX-$\beta$ diagram and find this to be consistent with grayer slopes, on average, in the UV. We also find evidence for a flattening of the attenuation curve in the NIR with respect to more classical Calzetti-like recipes. This larger NIR attenuation yields larger derived stellar masses from SED fitting, by a median factor of $\sim$ 1.4 and up to a factor $\sim$10 for the most extreme cases. The star formation rate appears instead to be more dependent on the total amount of attenuation in the galaxy. Our analysis highlights the need for a flexible attenuation curve when reproducing the physical properties of a large variety of objects.

We present \ci\,(2--1) and multi-transition $^{12}$CO observations of a dusty star-forming galaxy, ACT\,J2029+0120, which we spectroscopically confirm to lie at $z$\,=\,2.64. We detect CO(3--2), CO(5--4), CO(7--6), CO(8--7), and \ci\,(2--1) at high significance, tentatively detect HCO$^{+}$(4--3), and place strong upper limits on the integrated strength of dense gas tracers (HCN(4--3) and CS(7--6)). Multi-transition CO observations and dense gas tracers can provide valuable constraints on the molecular gas content and excitation conditions in high-redshift galaxies. We therefore use this unique data set to construct a CO spectral line energy distribution (SLED) of the source, which is most consistent with that of a ULIRG/Seyfert or QSO host object in the taxonomy of the \textit{Herschel} Comprehensive ULIRG Emission Survey. We employ RADEX models to fit the peak of the CO SLED, inferring a temperature of T$\sim$117 K and $n_{\text{H}_2}\sim10^5$ cm$^{-3}$, most consistent with a ULIRG/QSO object and the presence of high density tracers. We also find that the velocity width of the \ci\ line is potentially larger than seen in all CO transitions for this object, and that the $L'_{\rm C\,I(2-1)}/L'_{\rm CO(3-2)}$ ratio is also larger than seen in other lensed and unlensed submillimeter galaxies and QSO hosts; if confirmed, this anomaly could be an effect of differential lensing of a shocked molecular outflow.

We explore from a statistical point of view the far-infrared (far-IR) and sub-millimeter (sub-mm) properties of a large sample of LBGs (22,000) at z~3 in the COSMOS field. The large number of galaxies allows us to split it in several bins as a function of UV luminosity, UV slope, and stellar mass to better sample their variety. We perform stacking analysis in PACS (100 and 160 um), SPIRE (250, 350 and 500 um) and AzTEC (1.1 mm) images. Our stacking procedure corrects the biases induced by galaxy clustering and incompleteness of our input catalogue in dense regions. We obtain the full IR spectral energy distributions (SED) of subsamples of LBGs and derive the mean IR luminosity as a function of UV luminosity, UV slope, and stellar mass. The average IRX is roughly constant over the UV luminosity range, with a mean of 7.9 (1.8 mag). However, it is correlated with UV slope, and stellar mass. We investigate using a statistically-controlled stacking analysis as a function of (stellar mass, UV slope) the dispersion of the IRX-UVslope and IRX-M* plane. Our results enable us to study the average relation between star-formation rate (SFR) and stellar mass, and we show that our LBG sample lies on the main sequence of star formation at z~3.

We report contributions to cosmic infrared background (CIB) intensities originating from known galaxies and their faint companions at submillimeter wavelengths. Using the publicly-available UltraVISTA catalog, and maps at 250, 350, and 500 {\mu}m from the \emph{Herschel} Multi-tiered Extragalactic Survey (HerMES), we perform a novel measurement that exploits the fact that uncatalogued sources may bias stacked flux densities --- particularly if the resolution of the image is poor --- and intentionally smooth the images before stacking and summing intensities. By smoothing the maps we are capturing the contribution of faint (undetected in K_S ~ 23.4) sources that are physically associated, or correlated, with the detected sources. We find that the cumulative CIB increases with increased smoothing, reaching 9.82 +- 0.78, 5.77 +- 0.43, and 2.32 +- 0.19$\, \rm nW m^{-2} sr^{-1}$ at 250, 350, and 500 {\mu}m at 300 arcsec FWHM. This corresponds to a fraction of the fiducial CIB of 0.94 +- 0.23, 1.07 +- 0.31, and 0.97 +- 0.26 at 250, 350, and 500 {\mu}m, where the uncertainties are dominated by those of the absolute CIB. We then propose, with a simple model combining parametric descriptions for stacked flux densities and stellar mass functions, that emission from galaxies with log(M/Msun) > 8.5 can account for the most of the measured total intensities, and argue against contributions from extended, diffuse emission. Finally, we discuss prospects for future survey instruments to improve the estimates of the absolute CIB levels, and observe any potentially remaining emission at z > 4.