Taking advantage of HST CANDELS data, we analyze the lowest redshift (z<0.5) massive galaxies in order to disentangle their structural constituents and study possible faint non-axis-symmetric features. Due to the excellent HST spatial resolution for intermediate-z objects, they are hard to model by purely automatic parametric fitting algorithms. We performed careful single and double S\'ersic fits to their galaxy surface brightness profiles. We also compare the model color profiles with the observed ones and also derive multi-component global effective radii attempting to obtain a better interpretation of the mass-size relation. Additionally, we test the robustness of our measured structural parameters via simulations. We find that the S\'ersic index does not offer a good proxy for the visual morphological type for our sample of massive galaxies. Our derived multi-component effective radii give a better description of the size of our sample galaxies than those inferred from single S\'ersic models with GALFIT. Our galaxy population lays on the scatter of the local mass-size relation, indicating that these massive galaxies do not experience a significant growth in size since z~0.5. Interestingly the few outliers are late-type galaxies, indicating that spheroids must reach the local mass-size relation earlier. For most of our sample galaxies, both single and multi-component S\'ersic models with GALFIT show substantial systematic deviations from the observed SBPs in the outskirts. These residuals may be partly due to several factors, namely a non-optimal data reduction for low surface brightness features, the existence of prominent stellar haloes for massive galaxies and could also arise from conceptual shortcomings of parametric 2D image decomposition tools. They consequently propagate into galaxy color profiles.

The S\'ersic law (SL) offers a versatile functional form for the structural characterization of galaxies near and far. Whereas applying it to galaxies with a genuine SL luminosity distribution yields a robust determination of the S\'ersic exponent eta and effective surface brightness $\mu_{\rm eff}$, this is not necessarily the case for galaxies whose surface brightness profiles (SBPs) appreciably deviate from the SL (eg, early-type galaxies with a depleted core and nucleated dwarf ellipticals, or most late-type galaxies-LTGs). In this general case of "imperfect" SL profiles, the best-fitting solution may significantly depend on the radius (or surface brightness) interval fit and corrections for point spread function (PSF) convolution effects. Such uncertainties may then affect, in a non-easily predictable manner, automated structural studies of galaxies. We present a fitting concept (iFIT) that permits a robust determination of the equivalent SL model for the general case of galaxies with imperfect SL profiles. iFIT has been extensively tested on synthetic data with a S\'ersic index 0.3<${\eta}$<4.2 and an effective radius 1<$\rm{R}_{eff}$ (arcs)<20. Applied to non PSF-convolved data, iFIT can infer the S\'ersic exponent eta with an absolute error of <0.2 even for shallow SBPs. As for PSF-degraded data, iFIT can recover the input SL model parameters with a satisfactorily accuracy almost over the entire considered parameter space as long as FWHM(PSF)<$\rm{R}_{eff}$. Tests indicate that iFIT shows little sensitivity on PSF corrections and the SBP limiting surface brightness, and that subtraction of the best-fitting SL model in two different bands yields a good match to the observed radial color profile. The publicly available iFIT offers an efficient tool for the non-supervised structural characterization of large galaxy samples, as those expected to become available with Euclid and LSST.

We present an exploration of the expected detection of the earliest Active Galactic Nuclei (AGN) in the Universe from state-of-art galaxy formation and evolution semi-analytic models and hydro-dynamical simulations. We estimate the number and radiative characteristics of Super Massive Black Holes (SMBHs) at $z\geq 6$, a redshift range that will be intensively explored by the next generation of telescopes, in particular in the radio through the Square Kilometre Array (SKA) and at high energies with ESA's Athena X-ray Observatory. We find that Athena will be able to observe over 5000 AGN/$\rm deg^2$ at the Epoch of Re-ionization (EoR), $6\leq z \leq 10$. Similarly, for the same redshift range the models/simulations suggest that SKA will detect at least 400 AGN/$\rm deg^2$. Additionally, we stress the importance of the volume of the simulation box as well as the initial physical conditions of the models/simulations on their effect on the luminosity functions (LFs) and the creation of the most massive SMBHs that we currently observe at the EoR. Furthermore, following the evolution of the accretion mode of the SMBHs in each model/simulation, we show that, while the quasar dominates over the radio mode at the EoR, detection at radio wavelengths still reaches significant numbers even at the highest redshifts. Finally, we present the effect that the radiative efficiency has on the LFs by comparing results produced with a constant value for the radiative efficiency and more complex calculations based on the spin of each SMBH.