In the second work of this series, we analyse the connection between the availability of gas and the position of a region with respect to the spatially resolved main sequence (MS) relation. Following the procedure presented in Paper I we obtain 500pc scales estimates of stellar mass and star formation rate surface densities ($\Sigma_{\star}$ and $\Sigma_{\rm{SFR}}$). Our sample consists of five face-on, grand design spiral galaxies located on the MS. Thanks to HI 21cm and $^{12}$CO(2-1) maps, we connect the gas surface densities and gas fractions to the observed star formation properties of each region. We find that the spatially resolved MS ($\sigma=0.23$ dex) is the combination of two relations: the Kennicutt-Schmidt law ($\sigma=0.19$ dex) and the molecular gas MS (MGMS, $\sigma=0.22$ dex); $\Sigma_{\star}$, $\Sigma_{\rm{SFR}}$ and the surface density of the molecular gas, $\Sigma_{\rm{H_2}}$, define a 3D relation as proposed by \citet{2019ApJ...884L..33L}. We find that $\Sigma_{\rm{H_2}}$ steadily increases along the MS relation, varies little towards higher $\Sigma_{\rm{SFR}}$ at fixed stellar surface densities (not enough to sustain the change in SFR), and it is almost constant perpendicular to the relation. The surface density of neutral gas ($\Sigma_{\rm{HI}}$) is constant along the MS, and increases in its upper envelop. $\Sigma_{\rm{SFR}}$ can be expressed as a function of $\Sigma_{\star}$ and $\Sigma_{\rm{HI}}$, following the Equation: $\log\Sigma_{\rm{SFR}}$ = 0.97$\log\Sigma_{\star}$ + 1.99$\log\Sigma_{\rm{HI}}$ - 11.11. Finally, we show that f$_{\rm{gas}}$ increases significantly towards the starburst region in the $\log\Sigma_{\star}$ - $\log\Sigma_{\rm{SFR}}$ plane, accompanied by a slight increase in SFE.

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 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.