We have obtained three-dimensional maps of the universe in $\sim200\times200\times80$ comoving Mpc$^3$ (cMpc$^3$) volumes each at $z=5.7$ and $6.6$ based on a spectroscopic sample of 179 galaxies that achieves $\gtrsim80$\% completeness down to the Ly$\alpha$ luminosity of $\log(L_{\rm Ly\alpha}/[\mathrm{erg\ s^{-1}}])=43.0$, based on our Keck and Gemini observations and the literature. The maps reveal filamentary large-scale structures and two remarkable overdensities made out of at least 44 and 12 galaxies at $z=5.692$ (z57OD) and $z=6.585$ (z66OD), respectively, making z66OD the most distant overdensity spectroscopically confirmed to date with $>10$ spectroscopically confirmed galaxies. We compare spatial distributions of submillimeter galaxies at $z\simeq 4-6$ with our $z=5.7$ galaxies forming the large-scale structures, and detect a $99.97\%$ signal of cross correlation, indicative of a clear coincidence of dusty star-forming galaxy and dust unobscured galaxy formation at this early epoch. The galaxies in z57OD and z66OD are actively forming stars with star formation rates (SFRs) $\gtrsim5$ times higher than the main sequence, and particularly the SFR density in z57OD is 10 times higher than the cosmic average at the redshift (a.k.a. the Madau-Lilly plot). Comparisons with numerical simulations suggest that z57OD and z66OD are protoclusters that are progenitors of the present-day clusters with halo masses of $\sim10^{14}\ \mathrm{M_\odot}$.

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

Far-infrared (FIR)--radio correlation is a well-established empirical connection between continuum radio and dust emission of star-forming galaxies, often used as a tool in determining star-formation rates. Here we expand the point made by Murphy (2013) that in the case of some interacting star-forming galaxies there is a non-thermal emission from the gas bridge in between them, which might cause a dispersion in this correlation. Galactic interactions and mergers have been known to give rise to tidal shocks and disrupt morphologies especially in the smaller of the interacting components. Here we point out that these shocks can also heat the gas and dust and will inevitably accelerate particles and result in a tidal cosmic-ray population in addition to standard galactic cosmic rays in the galaxy itself. This would result in a non-thermal emission not only from the gas bridges of interacting systems, but from interacting galaxies as a whole in general. Thus both tidal heating and additional non-thermal radiation will obviously affect the FIR-radio correlation of these systems, the only question is how much. In this scenario the FIR-radio correlation is not stable in interacting galaxies, but rather evolves as the interaction/merger progresses. To test this hypothesis and probe the possible impact of tidal cosmic ray population we have analyzed a sample of 43 infrared bright star-forming interacting galaxies at different merger stages. We have found that their FIR-radio correlation parameter and radio emission spectral index vary noticeably over different merger stages and behave as it would be expected from our tidal-shock scenario. Important implications of departure of interacting galaxies from the FIR-radio correlation are discussed.