We present a new prospective analysis of deep multi-band imaging with the James Webb Space Telescope (JWST). In this work, we investigate the recovery of high-redshift $5<z<12$ galaxies through extensive image simulations of accepted JWST programs such as CEERS in the EGS field and HUDF GTO. We introduce complete samples of $\sim300,000$ galaxies with stellar masses $\log(M_*/M_\odot)>6$ and redshifts $0<z<15$, as well as galactic stars, into realistic mock NIRCam, MIRI and HST images to properly describe the impact of source blending. We extract the photometry of the detected sources as in real images and estimate the physical properties of galaxies through spectral energy distribution fitting. We find that the photometric redshifts are primarily limited by the availability of blue-band and near-infrared medium-band imaging. The stellar masses and star-formation rates are recovered within $0.25$ and $0.3$ dex respectively, for galaxies with accurate photometric redshifts. Brown dwarfs contaminating the $z>5$ galaxy samples can be reduced to $<0.01$ arcmin$^{-2}$ with a limited impact on galaxy completeness. We investigate multiple high-redshift galaxy selection techniques and find the best compromise between completeness and purity at $5<z<10$ using the full redshift posterior probability distributions. In the EGS field, the galaxy completeness remains higher than $50\%$ for $m_\text{UV}<27.5$ sources at all redshifts, and the purity is maintained above $80$ and $60\%$ at $z\leq7$ and $10$ respectively. The faint-end slope of the galaxy UV luminosity function is recovered with a precision of $0.1-0.25$, and the cosmic star-formation rate density within $0.1$ dex. We argue in favor of additional observing programs covering larger areas to better constrain the bright end.

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 cold molecular gas in contemporary galaxies is structured in discrete cloud complexes. These giant molecular clouds (GMCs), with $10^4$-$10^7$ solar masses and radii of 5-100 parsecs, are the seeds of star formation. Highlighting the molecular gas structure at such small scales in distant galaxies is observationally challenging. Only a handful of molecular clouds were reported in two extreme submillimetre galaxies at high redshift. Here we search for GMCs in a typical Milky Way progenitor at z = 1.036. Using the Atacama Large Millimeter/submillimeter Array (ALMA), we mapped the CO(4-3) emission of this gravitationally lensed galaxy at high resolution, reading down to 30 parsecs, which is comparable to the resolution of CO observations of nearby galaxies. We identify 17 molecular clouds, characterized by masses, surface densities and supersonic turbulence all of which are 10-100 times higher than present-day analogues. These properties question the universality of GMCs and suggest that GMCs inherit their properties from ambient interstellar medium. The measured cloud gas masses are similar to the masses of stellar clumps seen in the galaxy in comparable numbers. This corroborates the formation of molecular clouds by fragmentation of distant turbulent galactic gas disks, which then turn into stellar clumps ubiquitously observed in galaxies at cosmic noon.

We present high-fidelity, 30 milliarcsecond (200-pc) resolution ALMA rest-frame 240 $\mu$m observations of cold dust emission in three typical main-sequence star-forming galaxies (SFGs) at $z \sim 3$ in the Hubble Ultra-Deep Field (HUDF). The cold dust is distributed within the smooth disk-like central regions of star formation $1 - 3$ kpc in diameter, despite their complex and disturbed rest-frame UV and optical morphologies. No dust substructures or clumps are seen down to $\simeq 1- 3$ $M_\odot$yr$^{-1}$ (1$\sigma$) per 200-pc beam. No dust emission is observed at the locations of UV-emitting clumps, which lie $\simeq 2-10$ kpc from the bulk of star formation. Clumpy substructures can contribute no more than $1-7$% of the total star formation in these galaxies (3$\sigma$ upper limits). The lack of star-forming substructures in our HUDF galaxies is to be contrasted with the multiple substructures characteristic of submillimeter-selected galaxies (SMGs) at the same cosmic epoch, particularly the far-IR-bright SMGs with similarly high-fidelity ALMA observations of Hodge et al. (2019). Individual star-forming substructures in these SMGs contain $\sim10-30$% of their total star formation. A substructure in these SMGs is often comparably bright in the far-infrared as (or in some cases brighter than) our typical SFGs, suggesting that these SMGs originate from a class of disruptive event involving multiple objects at the scale of our HUDF galaxies. The scale of the disruptive event found in our main-sequence SFGs, characterized by the lack of star-forming substructures at our resolution and sensitivity, could be less violent, e.g., gas-rich disk instability or minor mergers.

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.

A diverse range of dust attenuation laws is found in star-forming galaxies. In particular, Tress et al. (2018) studied the SHARDS survey to constrain the NUV bump strength (B) and the total-to selective ratio (Rv) of 1,753 star-forming galaxies in the GOODS-N field at 1.5<z<3. We revisit here this sample to assess the implications and possible causes of the correlation found between Rv and B. The UVJ bicolour plot and main sequence of star formation are scrutinised to look for clues into the observed trend. The standard boundary between quiescent and star-forming galaxies is preserved when taking into account the wide range of attenuation parameters. However, an additional degeneracy, regarding the effective attenuation law, is added to the standard loci of star-forming galaxies in the UVJ diagram. A simple phenomenological model with an age-dependent extinction (at fixed dust composition) is compatible with the observed trend between Rv and B, whereby the opacity decreases with the age of the populations, resulting in a weaker NUV bump when the overall attenuation is shallower (greyer). In addition, we compare the constraints obtained by the SHARDS sample with dust models from the literature, supporting a scenario where geometry could potentially drive the correlation between Rv and B

Stellar masses in future James Webb Space Telescope (JWST) deep blank-fields will be mainly derived fitting the spectral energy distribution with theoretical galaxy templates. We investigate the uncertainties and biases of the stellar masses derived by using the LePhare code for spectral energy distribution fitting and the Yggdrasil theoretical templates. We consider a sample of mock galaxies at z$=$7-10 with mock JWST observations with S/N$_{F150W}\geqslant$ 10. Our goal is to provide a list of statistical stellar mass corrections to include on the stellar mass derivation for different output galaxy properties and JWST filter combinations to correct for template degeneracies. Median statistical stellar mass corrections vary from -0.83 dex to 0.87 dex, while 25$\%$ (75$\%$) quartiles range from -0.83 dex (-0.67 dex) to 0.51 dex (0.88 dex), depending on filter combinations and galaxy models. The most challenging cases are galaxies with nebular emission lines, especially the ones that are wrongly identified as galaxies without, relative dust-free galaxies and galaxies with small metallicities (i.e. Z$=1/50Z_{\odot}$). The stellar mass estimation of galaxies correctly identified without emission lines is generally fine, except at z$=$10 when considering only the 8 NIRCam bands, which make the MIRI bands very valuable. We have tested our stellar mass corrections using the public JAGUAR galaxy catalogue, deriving that the average discrepancy in the recovered stellar mass distribution decreases by 20-50$\%$ at z$>$7 after the correction. We found that without the stellar-mass corrections the number of low-mass galaxies (M$^{*}<10^{7}M_{\odot}$) is overestimated, which can potentially lead to systematic errors in the calculation of the galaxy stellar mass function faint-end slope at high z.

Using the Spitzer Space Telescope and Herschel Space Observatory, we have conducted a survey of infrared galaxies in the field of the galaxy cluster Abell S1063 (AS1063) at $z=0.347$, which is one of the most massive clusters known and a target of the HST CLASH and Frontier-Field surveys. The Spitzer/MIPS 24 $\mu$m and Herschel/PACS & SPIRE images revealed that the core of AS1063 is surprisingly devoid of infrared sources, showing only a few detectable sources within the central r$\sim1^{\prime}$. There is, however, one particularly bright source (2.3 mJy at 24 $\mu$m; 106 mJy at 160 $\mu$m), which corresponds to a background galaxy at $z=0.61$. The modest magnification factor (4.0$\times$) implies that this galaxy is intrinsically IR-luminous (L$_{\rm IR}=3.1\times10^{11}\ \rm L_{\odot}$). What is particularly interesting about this galaxy is that HST optical/near-infrared images show a remarkably bright and large (1 kpc) clump at one edge of the disk. Our follow-up optical/near-infrared spectroscopy shows Balmer (H$\alpha$-H8) and forbidden emission from this clump ([OII] $\lambda$3727, [OIII] $\lambda\lambda$4959,5007, [NII] $\lambda\lambda$6548,6583), indicating that it is a HII region. The HII region appears to have formed in-situ, as kinematically it is part of a rotating disk, and there is no evidence of nearby interacting galaxies. With an extinction correction of A$_{\rm V}=1.5$ mag, the star formation rate of this giant HII region is $\sim$10 M$_{\odot}$ yr$^{-1}$, which is exceptionally large, even for high redshift HII regions. Such a large and luminous HII region is often seen at $z\sim2$ but quite rare in the nearby Universe.

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.

Galaxy mergers are expected to have a significant role in the mass assembly of galaxies in the early Universe, but there are very few observational constraints on the merger history of galaxies at $z>2$. We present the first study of galaxy major mergers (mass ratios $>$ 1:4) in mass-selected samples out to $z\approx6$. Using all five fields of the HST/CANDELS survey and a probabilistic pair count methodology that incorporates the full photometric redshift posteriors and corrections for stellar mass completeness, we measure galaxy pair-counts for projected separations between 5 and 30 kpc in stellar mass selected samples at $9.7 < \log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) < 10.3$ and $\log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) > 10.3$. We find that the major merger pair fraction rises with redshift to $z\approx6$ proportional to $(1+z)^{m}$, with $m = 0.8\pm0.2$ ($m = 1.8\pm0.2$) for $\log_{10}(\rm{M}_{*} / \rm{M}_{\odot}) > 10.3$ ($9.7 < \log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) < 10.3$). Investigating the pair fraction as a function of mass ratio between 1:20 and 1:1, we find no evidence for a strong evolution in the relative numbers of minor to major mergers out to $z<3$. Using evolving merger timescales we find that the merger rate per galaxy ($\mathcal{R}$) rises rapidly from $0.07\pm 0.01$ Gyr$^{-1}$ at $z < 1$ to $7.6\pm 2.7$ Gyr$^{-1}$ at $z = 6$ for galaxies at $\log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) > 10.3$. The corresponding co-moving major merger rate density remains roughly constant during this time, with rates of $\Gamma \approx 10^{-4}$ Gyr$^{-1}$ Mpc$^{-3}$. Based on the observed merger rates per galaxy, we infer specific mass accretion rates from major mergers that are comparable to the specific star-formation rates for the same mass galaxies at $z>3$ - observational evidence that mergers are as important a mechanism for building up mass at high redshift as in-situ star-formation.

We study the rest-frame optical mass-size relation of bulges and discs from z~2 to z~0 for a complete sample of massive galaxies in the CANDELS fields using 2 component S\'ersic models (Dimauro et a. 2018). Discs and star forming galaxies follow similar mass-size relations. The mass-size relation of bulges is less steep than the one of quiescent galaxies (best fit slope of ~ 0.7 for quiescent galaxies against ~ 0.4 for bulges). We find little dependence of the structural properties of massive bulges and discs with the global morphology of galaxies (disc vs. bulge dominated) and the star formation activity (star-forming vs. quiescent). This result suggests similar bulge formation mechanisms for most massive galaxies and also that the formation of the bulge component does not significantly affect the disc structure. Our findings pose a challenge to models envisioning multiple channels for massive bulge growth, such as disc instabilities and mergers.

We quantify the star formation (SF) in the inner cores ($\mathcal{R}$/$R_{200}$$\leq$0.3) of 24 massive galaxy clusters at 0.2$\lesssim$$z$$\lesssim$0.9 observed by the $Herschel$ Lensing Survey and the Cluster Lensing and Supernova survey with $Hubble$. These programmes, covering the rest-frame ultraviolet to far-infrared regimes, allow us to accurately characterize stellar mass-limited ($\mathcal{M}_{*}$$>$$10^{10}$ $M_{\odot}$) samples of star-forming cluster members (not)-detected in the mid- and/or far-infrared. We release the catalogues with the photometry, photometric redshifts, and physical properties of these samples. We also quantify the SF displayed by comparable field samples from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey. We find that in intermediate-$z$ cluster cores, the SF activity is suppressed with respect the field in terms of both the fraction ($\mathcal{F}$) of star-forming galaxies (SFG) and the rate at which they form stars ($\mathcal{SFR}$ and $s\mathcal{SFR} = \mathcal{SFR}/\mathcal{M}_{*}$). On average, the $\mathcal{F}$ of SFGs is a factor $\sim$$2$ smaller in cluster cores than in the field. Furthermore, SFGs present average $\mathcal{SFR}$ and $s\mathcal{SFR}$ typically $\sim$0.3 dex smaller in the clusters than in the field along the whole redshift range probed. Our results favour long time-scale quenching physical processes as the main driver of SF suppression in the inner cores of clusters since $z$$\sim$0.9, with shorter time-scale processes being very likely responsible for a fraction of the missing SFG population.

With the recent discovery of a dozen dusty star-forming galaxies and around 30 quasars at z>5 that are hyper-luminous in the infrared ($\mu$$L_{\rm IR}>10^{13}$ L$_{\odot}$, where $\mu$ is a lensing magnification factor), the possibility has opened up for SPICA, the proposed ESA M5 mid-/far-infrared mission, to extend its spectroscopic studies toward the epoch of reionization and beyond. In this paper, we examine the feasibility and scientific potential of such observations with SPICA's far-infrared spectrometer SAFARI, which will probe a spectral range (35-230 $\mu$m) that will be unexplored by ALMA and JWST. Our simulations show that SAFARI is capable of delivering good-quality spectra for hyper-luminous infrared galaxies (HyLIRGs) at z=5-10, allowing us to sample spectral features in the rest-frame mid-infrared and to investigate a host of key scientific issues, such as the relative importance of star formation versus AGN, the hardness of the radiation field, the level of chemical enrichment, and the properties of the molecular gas. From a broader perspective, SAFARI offers the potential to open up a new frontier in the study of the early Universe, providing access to uniquely powerful spectral features for probing first-generation objects, such as the key cooling lines of low-metallicity or metal-free forming galaxies (fine-structure and H2 lines) and emission features of solid compounds freshly synthesized by Population III supernovae. Ultimately, SAFARI's ability to explore the high-redshift Universe will be determined by the availability of sufficiently bright targets (whether intrinsically luminous or gravitationally lensed). With its launch expected around 2030, SPICA is ideally positioned to take full advantage of upcoming wide-field surveys such as LSST, SKA, Euclid, and WFIRST, which are likely to provide extraordinary targets for SAFARI.

We investigate the galaxy quenching process at intermediate redshift using a sample of $\sim4400$ galaxies with $M_{\ast} > 10^{9}M_{\odot}$ between redshift 0.5 and 1.0 in all five CANDELS fields. We divide this sample, using the integrated specific star formation rate (sSFR), into four sub-groups: star-forming galaxies (SFGs) above and below the ridge of the star-forming main sequence (SFMS), transition galaxies and quiescent galaxies. We study their $UVI$ ($U-V$ versus $V-I$) color gradients to infer their sSFR gradients out to twice effective radii. We show that on average both star-forming and transition galaxies at all masses are not fully quenched at any radii, whereas quiescent galaxies are fully quenched at all radii. We find that at low masses ($M_{\ast} = 10^{9}-10^{10}M_{\odot}$) SFGs both above and below the SFMS ridge generally have flat sSFR profiles, whereas the transition galaxies at the same masses generally have sSFRs that are more suppressed in their outskirts. In contrast, at high masses ($M_{\ast} > 10^{10.5}M_{\odot}$), SFGs above and below the SFMS ridge and transition galaxies generally have varying degrees of more centrally-suppressed sSFRs relative to their outskirts. These findings indicate that at $z\sim~0.5-1.0$ the main galaxy quenching mode depends on its already formed stellar mass, exhibiting a transition from "the outside-in" at $M_{\ast} \leq 10^{10}M_{\odot}$ to "the inside-out" at $M_{\ast} > 10^{10.5}M_{\odot}$. In other words, our findings support that internal processes dominate the quenching of massive galaxies, whereas external processes dominate the quenching of low-mass galaxies.

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.

Understanding how bulges grow in galaxies is critical step towards unveiling the link between galaxy morphology and star-formation. To do so, it is necessary to decompose large sample of galaxies at different epochs into their main components (bulges and disks). This is particularly challenging, especially at high redshifts, where galaxies are poorly resolved. This work presents a catalog of bulge-disk decompositions of the surface brightness profiles of ~17.600 H-band selected galaxies in the CANDELS fields (F160W<23, 0<z<2) in 4 to 7 filters covering a spectral range of 430-1600nm. This is the largest available catalog of this kind up to z = 2. By using a novel approach based on deep-learning to select the best model to fit, we manage to control systematics arising from wrong model selection and obtain less contaminated samples than previous works. We show that the derived structural properties are within $\sim10-20\%$ of random uncertainties. We then fit stellar population models to the decomposed SEDs (Spectral Energy Distribution) of bulges and disks and derive stellar masses (and stellar mass bulge-to-total ratios) as well as rest-frame colors (U,V,J) for bulges and disks separately. All data products are publicly released with this paper and through the web page https://lerma.obspm.fr/huertas/form_CANDELS and will be used for scientific analysis in forthcoming works.

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.

Type-III S0 galaxies present tight scaling relations between their surface brightness photometric and structural parameters. Several evolutionary models have been proposed for the formation of Type-III S0 galaxies but the observations of are usually limited to the local Universe. We study the evolution of the photometric and structural scaling relations found between the parameters of the surface brightness profiles in the rest-frame R-band of Type-III S0 galaxies with z and the possible differences between the rest-frame (B-R) colours of the inner and outer disc profiles. We make use of a sample of 14 Type-III E/S0--S0 galaxies at 0.2<z<0.6 to study if the correlations found in local Type-III S0 galaxies were present ~6 Gyr ago. We analyse the distribution of the surface brightness characteristic parameters as a function of the stellar mass and if there is a significant change with z. We also derive their rest-frame (B-R) colour profiles and we compare these results with the predictions from a grid of SSP models. We find that the inner and outer scale-lengths of Type-III S0 galaxies at 0.4<z<0.6 follow compatible trends and scaling relations with those observed in local S0 galaxies. We do not detect any significant differences between the location of Rbreak between z~0.6 and z=0 for a fixed stellar mass of the object, whereas the surface brightness at the break radius is ~1.5 mag arcsec-2 dimmer in the local Universe than at z~0.6 for a fixed stellar mass. In contrast to Type-II profiles, the Type-III surface brightness profiles of S0 galaxies present compatible Rbreak values and scaling relations during the last 6 Gyr. This result and the similarity of the colours of the inner and outer discs point to a highly scalable and stable formation process, probably more related to gravitational and dynamical processes than to the evolution of stellar populations (abridged).

We present results from the KMOS lensing survey-KLENS which is exploiting gravitational lensing to study the kinematics of 24 star forming galaxies at $1.4<z<3.5$ with a median mass of $\rm log(M_\star/M_\odot)=9.6$ and median star formation rate (SFR) of $\rm 7.5\,M_\odot\,yr^{-1}$. We find that 25% of these low-mass/low-SFR galaxies are rotation dominated, while the majority of our sample shows no velocity gradient. When combining our data with other surveys, we find that the fraction of rotation dominated galaxies increases with the stellar mass, and decreases for galaxies with a positive offset from the main sequence. We also investigate the evolution of the intrinsic velocity dispersion, $\sigma_0$, as a function of the redshift, $z$, and stellar mass, $\rm M_\star$, assuming galaxies in quasi-equilibrium (Toomre Q parameter equal to 1). From the $z-\sigma_0$ relation, we find that the redshift evolution of the velocity dispersion is mostly expected for massive galaxies ($\rm log(M_\star/M_\odot)>10$). We derive a $\rm M_\star-\sigma_0$ relation, using the Tully-Fisher relation, which highlights that a different evolution of the velocity dispersion is expected depending on the stellar mass, with lower velocity dispersions for lower masses, and an increase for higher masses, stronger at higher redshift. The observed velocity dispersions from this work and from comparison samples spanning $0<z<3.5$ appear to follow this relation, except at higher redshift ($z>2$), where we observe higher velocity dispersions for low masses ($\rm log(M_\star/M_\odot)\sim 9.6$) and lower velocity dispersions for high masses ($\rm log(M_\star/M_\odot)\sim 10.9$) than expected. This discrepancy could, for instance, suggest that galaxies at high-$z$ do not satisfy the stability criterion, or that the adopted parametrisation of the specific star formation rate and molecular properties fail at high redshift.

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.

Studying giant star-forming clumps in distant galaxies is important to understand galaxy formation and evolution. At present, however, observers and theorists have not reached a consensus on whether the observed "clumps" in distant galaxies are the same phenomenon that is seen in simulations. In this paper, as a step to establish a benchmark of direct comparisons between observations and theories, we publish a sample of clumps constructed to represent the commonly observed "clumps" in the literature. This sample contains 3193 clumps detected from 1270 galaxies at $0.5 \leq z < 3.0$. The clumps are detected from rest-frame UV images, as described in our previous paper. Their physical properties, e.g., rest-frame color, stellar mass (M*), star formation rate (SFR), age, and dust extinction, are measured by fitting the spectral energy distribution (SED) to synthetic stellar population models. We carefully test the procedures of measuring clump properties, especially the method of subtracting background fluxes from the diffuse component of galaxies. With our fiducial background subtraction, we find a radial clump U-V color variation, where clumps close to galactic centers are redder than those in outskirts. The slope of the color gradient (clump color as a function of their galactocentric distance scaled by the semi-major axis of galaxies) changes with redshift and M* of the host galaxies: at a fixed M*, the slope becomes steeper toward low redshift, and at a fixed redshift, it becomes slightly steeper with M*. Based on our SED-fitting, this observed color gradient can be explained by a combination of a negative age gradient, a negative E(B-V) gradient, and a positive specific star formation rate gradient of the clumps. We also find that the color gradients of clumps are steeper than those of intra-clump regions. [Abridged]

The SPICA mid and far-infrared telescope will address fundamental issues in our understanding of star formation and ISM physics in galaxies. A particular hallmark of SPICA is the outstanding sensitivity enabled by the cold telescope, optimized detectors, and wide instantaneous bandwidth throughout the mid- and far-infrared. The spectroscopic, imaging and polarimetric observations that SPICA will be able to collect will help in clarifying the complex physical mechanisms which underlie the baryon cycle of galaxies. In particular: (i) The access to a large suite of atomic and ionic fine-structure lines for large samples of galaxies will shed light on the origin of the observed spread in star formation rates within and between galaxies. (ii) Observations of HD rotational lines (out to $\sim$10 Mpc) and fine structure lines such as [CII] 158 $\mu$m (out to $\sim$100 Mpc) will clarify the main reservoirs of interstellar matter in galaxies, including phases where CO does not emit. (iii) Far-infrared spectroscopy of dust and ice features will address uncertainties in the mass and composition of dust in galaxies, and the contributions of supernovae to the interstellar dust budget will be quantified by photometry and monitoring of supernova remnants in nearby galaxies. (iv) Observations of far-infrared cooling lines such as [OI] 63 $\mu$m from star-forming molecular clouds in our Galaxy will evaluate the importance of shocks to dissipate turbulent energy. The paper concludes with requirements for the telescope and instruments, and recommendations for the observing strategy.

Giant stellar clumps are ubiquitous in high-redshift galaxies. They are thought to play an important role in the build-up of galactic bulges and as diagnostics of star formation feedback in galactic discs. Hubble Space Telescope (HST) blank field imaging surveys have estimated that these clumps have masses up to 10$^{9.5}$ M$_{\odot}$ and linear sizes larger than ~1 kpc. Recently, gravitational lensing has also been used to get higher spatial resolution. However, both recent lensed observations and models suggest that the clumps properties may be overestimated by the limited resolution of standard imaging techniques. A definitive proof of this observational bias is nevertheless still missing. Here we investigate directly the effect of resolution on clump properties by analysing multiple gravitationally-lensed images of the same galaxy at different spatial resolutions, down to 30 pc. We show that the typical mass and size of giant clumps, generally observed at $\sim$1 kpc resolution in high-redshift galaxies, are systematically overestimated. The high spatial resolution data, only enabled by strong gravitational lensing using currently available facilities, support smaller scales of clump formation by fragmentation of the galactic gas disk via gravitational instabilities.

We analyze the properties of a multiply-imaged Lyman-alpha (Lya) emitter at z=5.75 identified through SHARDS Frontier Fields intermediate-band imaging of the Hubble Frontier Fields (HFF) cluster Abell 370. The source, A370-L57, has low intrinsic luminosity (M_UV~-16.5), steep UV spectral index (\beta=-2.4+/-0.1), and extreme rest-frame equivalent width of Lya (EW(Lya)=420+180-120 \AA). Two different gravitational lens models predict high magnification (\mu~10--16) for the two detected counter-images, separated by 7", while a predicted third counter-image (\mu~3--4) is undetected. We find differences of ~50% in magnification between the two lens models, quantifying our current systematic uncertainties. Integral field spectroscopy of A370-L57 with MUSE shows a narrow (FWHM=204+/-10 km/s) and asymmetric Lya profile with an integrated luminosity L(Lya)~10^42 erg/s. The morphology in the HST bands comprises a compact clump (r_e<100 pc) that dominates the Lya and continuum emission and several fainter clumps at projected distances <1 kpc that coincide with an extension of the Lya emission in the SHARDS F823W17 and MUSE observations. The latter could be part of the same galaxy or an interacting companion. We find no evidence of contribution from AGN to the Lya emission. Fitting of the spectral energy distribution with stellar population models favors a very young (t<10 Myr), low mass (M*~10^6.5 Msun), and metal poor (Z<4x10^-3) stellar population. Its modest star formation rate (SFR~1.0 Msun/yr) implies high specific SFR (sSFR~2.5x10^-7 yr^-1) and SFR density (Sigma_SFR ~ 7-35 Msun/yr/kpc^2). The properties of A370-L57 make it a good representative of the population of galaxies responsible for cosmic reionization.

Imaging with the James Webb Space Telescope (JWST) will allow for observing the bulk of distant galaxies at the epoch of reionisation. The recovery of their properties, such as age, color excess E(B-V), specific star formation rate (sSFR) and stellar mass, will mostly rely on spectral energy distribution fitting, based on the data provided by JWST's two imager cameras, namely the Near Infrared Camera (NIRCam) and the Mid Infrared Imager (MIRI). In this work we analyze the effect of choosing different combinations of NIRCam and MIRI broad-band filters, from 0.6 {\mu}m to 7.7 {\mu}m, on the recovery of these galaxy properties. We performed our tests on a sample of 1542 simulated galaxies, with known input properties, at z=7-10. We found that, with only 8 NIRCam broad-bands, we can recover the galaxy age within 0.1 Gyr and the color excess within 0.06 mag for 70% of the galaxies. Besides, the stellar masses and sSFR are recovered within 0.2 and 0.3 dex, respectively, at z=7-9. Instead, at z=10, no NIRCam band traces purely the {\lambda}> 4000 {\AA} regime and the percentage of outliers in stellar mass (sSFR) increases by > 20% (> 90%), in comparison to z=9. The MIRI F560W and F770W bands are crucial to improve the stellar mass and the sSFR estimation at z=10. When nebular emission lines are present, deriving correct galaxy properties is challenging, at any redshift and with any band combination. In particular, the stellar mass is systematically overestimated in up to 0.3 dex on average with NIRCam data alone and including MIRI observations improves only marginally the estimation.

We investigate the environmental quenching of galaxies, especially those with stellar masses (M*)$<10^{9.5} M_\odot$, beyond the local universe. Essentially all local low-mass quenched galaxies (QGs) are believed to live close to massive central galaxies, which is a demonstration of environmental quenching. We use CANDELS data to test {\it whether or not} such a dwarf QG--massive central galaxy connection exists beyond the local universe. To this purpose, we only need a statistically representative, rather than a complete, sample of low-mass galaxies, which enables our study to $z\gtrsim1.5$. For each low-mass galaxy, we measure the projected distance ($d_{proj}$) to its nearest massive neighbor (M*$>10^{10.5} M_\odot$) within a redshift range. At a given redshift and M*, the environmental quenching effect is considered to be observed if the $d_{proj}$ distribution of QGs ($d_{proj}^Q$) is significantly skewed toward lower values than that of star-forming galaxies ($d_{proj}^{SF}$). For galaxies with $10^{8} M_\odot < M* < 10^{10} M_\odot$, such a difference between $d_{proj}^Q$ and $d_{proj}^{SF}$ is detected up to $z\sim1$. Also, about 10\% of the quenched galaxies in our sample are located between two and four virial radii ($R_{Vir}$) of the massive halos. The median projected distance from low-mass QGs to their massive neighbors, $d_{proj}^Q / R_{Vir}$, decreases with satellite M* at $M* \lesssim 10^{9.5} M_\odot$, but increases with satellite M* at $M* \gtrsim 10^{9.5} M_\odot$. This trend suggests a smooth, if any, transition of the quenching timescale around $M* \sim 10^{9.5} M_\odot$ at $0.5<z<1.0$.

The controversy about the origin of the structure of S0--E/S0 galaxies may be due to the difficulty of comparing surface brightness profiles with different depths, photometric corrections and PSF effects (almost always ignored). We aim to quantify the properties of Type-III (anti-truncated) discs in a sample of S0 galaxies at 0.2<z<0.6. In this paper, we present the sample selection and describe in detail the methods to robustly trace the structure in their outskirts and correct for PSF effects. We have selected and classified a sample of 150 quiescent galaxies at 0.2<z<0.6 in the GOODS-N field. We perform a quantitative structural analysis of 44 S0-E/S0 galaxies. We corrected their surface brightness profiles for PSF distortions and analysed the biases in the structural and photometric parameters when the PSF correction is not applied. Additionally, we have developed Elbow, an automatic statistical method to determine whether a possible break is significant - or not - and its type and made it publicly available. We found 14 anti-truncated S0-E/S0 galaxies in the range 0.2<z<0.6 (~30% of the final sample). This fraction is similar to the those reported in the local Universe. In our sample, ~25% of the Type-III breaks observed in PSF-uncorrected profiles are artifacts, and their profiles turn into a Type I after PSF correction. PSF effects also soften Type-II profiles. We found that the profiles of Type-I S0 and E/S0 galaxies of our sample are compatible with the inner profiles of the Type-III, in contrast with the outer profiles. We have obtained the first robust and reliable sample of 14 anti-truncated S0--E/S0 galaxies beyond the local Universe, in the range 0.2<z<0.6. PSF effects significantly affect the shape of the surface brightness profiles in galaxy discs even in the case of the narrow PSF of HST/ACS images, so future studies on the subject should make an effort to correct them.

We explore empirical constraints on the statistical relationship between the radial size of galaxies and the radius of their host dark matter halos from $z\sim 0.1$--3 using the GAMA and CANDELS surveys. We map dark matter halo mass to galaxy stellar mass using relationships from abundance matching, applied to the Bolshoi-Planck dissipationless N-body simulation. We define SRHR$\equiv r_e/R_h$ as the ratio of galaxy radius to halo virial radius, and SRHR$\lambda \equiv r_e/(\lambda R_h)$ as the ratio of galaxy radius to halo spin parameter times halo radius. At $z\sim 0.1$, we find an average value of SRHR $\simeq 0.018$ and SRHR$\lambda \simeq 0.5$ with very little dependence on stellar mass. SRHR and SRHR$\lambda$ have a weak dependence on cosmic time since $z\sim 3$. SRHR shows a mild decrease over cosmic time for low mass galaxies, but increases slightly or does not evolve for more massive galaxies. We find hints that at high redshift ($z\sim 2$--3), SRHR$\lambda$ is lower for more massive galaxies, while it shows no significant dependence on stellar mass at $z\lesssim 0.5$. We find that for both the GAMA and CANDELS samples, at all redshifts from $z\sim 0.1$--3, the observed conditional size distribution in stellar mass bins is remarkably similar to the conditional distribution of $\lambda R_h$. We discuss the physical interpretation and implications of these results.

We investigate the burstiness of star formation histories (SFHs) of galaxies at $0.4<z<1$ by using the ratio of star formation rates (SFRs) measured from H$\beta$ and FUV (1500 \AA) (H$\beta$--to--FUV ratio). Our sample contains 164 galaxies down to stellar mass (M*) of $10^{8.5} M_\odot$ in the CANDELS GOODS-N region, where Team Keck Redshift Survey DEIMOS spectroscopy and HST/WFC3 F275W images from CANDELS and Hubble Deep UV Legacy Survey are available. When the {\it ratio} of H$\beta$- and FUV-derived SFRs is measured, dust extinction correction is negligible (except for very dusty galaxies) with the Calzetti attenuation curve. The H$\beta$--to--FUV ratio of our sample increases with M* and SFR. The median ratio is $\sim$0.7 at M*$\sim10^{8.5} M_\odot$ (or SFR$\sim 0.5 M_\odot/yr$) and increases to $\sim$1 at M*$\sim10^{10} M_\odot$ (or SFR $\sim 10 M_\odot/yr$). At M*$<10^{9.5} M_\odot$, our median H$\beta$--to--FUV ratio is lower than that of local galaxies at the same M*, implying a redshift evolution. Bursty SFH on a timescale of a few tens of megayears on galactic scales provides a plausible explanation of our results, and the importance of the burstiness increases as M* decreases. Due to sample selection effects, our H$\beta$--to--FUV ratio may be an upper limit of the true value of a complete sample, which strengthens our conclusions. Other models, e.g., non-universal initial mass function or stochastic star formation on star cluster scales, are unable to plausibly explain our results.

The determination of galaxy redshifts in James Webb Space Telescope (JWST)'s blank-field surveys will mostly rely on photometric estimates, based on the data provided by JWST's Near-Infrared Camera (NIRCam) at 0.6-5.0 {\mu}m and Mid Infrared Instrument (MIRI) at {\lambda}>5.0 {\mu}m. In this work we analyse the impact of choosing different combinations of NIRCam and MIRI broad-band filters (F070W to F770W), as well as having ancillary data at {\lambda}<0.6 {\mu}m, on the derived photometric redshifts (zphot) of a total of 5921 real and simulated galaxies, with known input redshifts z=0-10. We found that observations at {\lambda}<0.6 {\mu}m are necessary to control the contamination of high-z samples by low-z interlopers. Adding MIRI (F560W and F770W) photometry to the NIRCam data mitigates the absence of ancillary observations at {\lambda}<0.6 {\mu}m and improves the redshift estimation. At z=7-10, accurate zphot can be obtained with the NIRCam broad bands alone when S/N>=10, but the zphot quality significantly degrades at S/N<=5. Adding MIRI photometry with one magnitude brighter depth than the NIRCam depth allows for a redshift recovery of 83-99%, depending on SED type, and its effect is particularly noteworthy for galaxies with nebular emission. The vast majority of NIRCam galaxies with [F150W]=29 AB mag at z=7-10 will be detected with MIRI at [F560W, F770W]<28 mag if these sources are at least mildly evolved or have spectra with emission lines boosting the mid-infrared fluxes.

We present spatially-resolved Atacama Large Millimeter/sub-millimeter Array (ALMA) 870 $\mu$m dust continuum maps of six massive, compact, dusty star-forming galaxies (SFGs) at $z\sim2.5$. These galaxies are selected for their small rest-frame optical sizes ($r_{\rm e, F160W}\sim1.6$ kpc) and high stellar-mass densities that suggest that they are direct progenitors of compact quiescent galaxies at $z\sim2$. The deep observations yield high far-infrared (FIR) luminosities of L$_{\rm IR}=10^{12.3-12.8}$ L$_{\odot}$ and star formation rates (SFRs) of SFR$=200-700$ M$_{\odot}$yr$^{-1}$, consistent with those of typical star-forming "main sequence" galaxies. The high-spatial resolution (FWHM$\sim$0.12"-0.18") ALMA and HST photometry are combined to construct deconvolved, mean radial profiles of their stellar mass and (UV+IR) SFR. We find that the dusty, nuclear IR-SFR overwhelmingly dominates the bolometric SFR up to $r\sim5$ kpc, by a factor of over 100$\times$ from the unobscured UV-SFR. Furthermore, the effective radius of the mean SFR profile ($r_{\rm e, SFR}\sim1$ kpc) is $\sim$30% smaller than that of the stellar mass profile. The implied structural evolution, if such nuclear starburst last for the estimated gas depletion time of $\Delta t=\pm100$ Myr, is a 4$\times$ increase of the stellar mass density within the central 1 kpc and a 1.6$\times$ decrease of the half-mass radius. This structural evolution fully supports dissipation-driven, formation scenarios in which strong nuclear starbursts transform larger, star-forming progenitors into compact quiescent galaxies.

[abridged] We quantify the evolution of the stellar mass functions of star-forming and quiescent galaxies as a function of morphology from $z\sim 3$ to the present. Our sample consists of ~50,000 galaxies in the CANDELS fields ($\sim880$ $arcmin^2$), which we divide into four main morphological types, i.e. pure bulge dominated systems, pure spiral disk dominated, intermediate 2-component bulge+disk systems and irregular disturbed galaxies. Our main results are: Star-formation: At $z\sim 2$, 80\% of the stellar mass density of star-forming galaxies is in irregular systems. However, by $z\sim 0.5$, irregular objects only dominate at stellar masses below $10^9M\odot$. A majority of the star-forming irregulars present at $z\sim 2$ undergo a gradual transformation from disturbed to normal spiral disk morphologies by $z\sim 1$ without significant interruption to their star-formation. Rejuvenation after a quenching event does not seem to be common except perhaps for the most massive objects. Quenching: We confirm that galaxies reaching a stellar mass of $M_*\sim10^{10.8}M_\odot$ ($M^*$) tend to quench. Also, quenching implies the presence of a bulge: the abundance of massive red disks is negligible at all redshifts over 2~dex in stellar mass. However the dominant quenching mechanism evolves. At $z>2$, the SMF of quiescent galaxies above $M^*$ is dominated by compact spheroids. Quenching at this early epoch destroys the disk and produces a compact remnant unless the star-forming progenitors at even higher redshifts are significantly more dense. At $1<z<2$, the majority of newly quenched galaxies are disks with a significant central bulge. This suggests that mass-quenching at this epoch starts from the inner parts and preserves the disk. At $z<1$, the high mass end of the passive SMF is globally in place and the evolution mostly happens at stellar masses below $10^{10}M_\odot$.

We present a complete census of all Herschel-detected sources within the six massive lensing clusters of the HST Frontier Fields (HFF). We provide a robust legacy catalogue of 263 sources with Herschel fluxes, primarily based on imaging from the Herschel Lensing Survey (HLS) and PEP/HerMES Key Programmes. We optimally combine Herschel, Spitzer and WISE infrared (IR) photometry with data from HST, VLA and ground-based observatories, identifying counterparts to gain source redshifts. For each Herschel-detected source we also present magnification factor (mu), intrinsic IR luminosity and characteristic dust temperature, providing a comprehensive view of dust-obscured star formation within the HFF. We demonstrate the utility of our catalogues through an exploratory overview of the magnified population, including more than 20 background sub-LIRGs unreachable by Herschel without the assistance gravitational lensing.

We present the stellar mass ($M_{*}$)--gas-phase metallicity relation (MZR) and its scatter at intermediate redshifts ($0.5\leq z\leq0.7$) for 1381 field galaxies collected from deep spectroscopic surveys. The star formation rate (SFR) and color at a given $M_{*}$ of this magnitude-limited ($R\lesssim24$ AB) sample are representative of normal star-forming galaxies. For masses below $10^9 M_\odot$, our sample of 237 galaxies is $\sim$10 times larger than those in previous studies beyond the local universe. This huge gain in sample size enables superior constraints on the MZR and its scatter in the low-mass regime. We find a power-law MZR at $10^{8} M_\odot < M_{*} < 10^{11} M_\odot$: ${12+log(O/H) = (5.83\pm0.19) + (0.30\pm0.02)log(M_{*}/M_\odot)}$. Our MZR shows good agreement with others measured at similar redshifts in the literature in the intermediate and massive regimes, but is shallower than the extrapolation of the MZRs of others to masses below $10^{9} M_\odot$. The SFR dependence of the MZR in our sample is weaker than that found for local galaxies (known as the Fundamental Metallicity Relation). Compared to a variety of theoretical models, the slope of our MZR for low-mass galaxies agrees well with predictions incorporating supernova energy-driven winds. Being robust against currently uncertain metallicity calibrations, the scatter of the MZR serves as a powerful diagnostic of the stochastic history of gas accretion, gas recycling, and star formation of low-mass galaxies. Our major result is that the scatter of our MZR increases as $M_{*}$ decreases. Our result implies that either the scatter of the baryonic accretion rate or the scatter of the $M_{*}$--$M_{halo}$ relation increases as $M_{*}$ decreases. Moreover, our measures of scatter at $z=0.7$ appears consistent with that found for local galaxies.

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.

We introduce a new color-selection technique to identify high-redshift, massive galaxies that are systematically missed by Lyman-break selection. The new selection is based on the H_{160} and IRAC 4.5um bands, specifically H - [4.5] > 2.25 mag. These galaxies, dubbed "HIEROs", include two major populations that can be separated with an additional J - H color. The populations are massive and dusty star-forming galaxies at z > 3 (JH-blue) and extremely dusty galaxies at z < 3 (JH-red). The 350 arcmin^2 of the GOODS-N and GOODS-S fields with the deepest HST/WFC3 and IRAC data contain 285 HIEROs down to [4.5] < 24 mag. We focus here primarily on JH-blue (z > 3) HIEROs, which have a median photometric redshift z ~4.4 and stellar massM_{*}~10^{10.6} Msun, and are much fainter in the rest-frame UV than similarly massive Lyman-break galaxies (LBGs). Their star formation rates (SFRs) reaches ~240 Msun yr^{-1} leading to a specific SFR, sSFR ~4.2 Gyr^{-1}, suggesting that the sSFRs for massive galaxies continue to grow at z > 2 but at a lower growth rate than from z=0 to z=2. With a median half-light radius of 2 kpc, including ~20% as compact as quiescent galaxies at similar redshifts, JH-blue HIEROs represent perfect star-forming progenitors of the most massive (M_{*} > 10^{11.2} Msun) compact quiescent galaxies at z ~ 3 and have the right number density. HIEROs make up ~60% of all galaxies with M_{*} > 10^{10.5} Msun identified at z > 3 from their photometric redshifts. This is five times more than LBGs with nearly no overlap between the two populations. While HIEROs make up 15-25% of the total SFR density at z ~ 4-5, they completely dominate the SFR density taking place in M_{*} >10^{10.5} Msun galaxies, and are therefore crucial to understanding the very early phase of massive galaxy formation.

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.

We present a catalog of visual like H-band morphologies of $\sim50.000$ galaxies ($H_{f160w}<24.5$) in the 5 CANDELS fields (GOODS-N, GOODS-S, UDS, EGS and COSMOS). Morphologies are estimated with Convolutional Neural Networks (ConvNets). The median redshift of the sample is $<z>\sim1.25$. The algorithm is trained on GOODS-S for which visual classifications are publicly available and then applied to the other 4 fields. Following the CANDELS main morphology classification scheme, our model retrieves the probabilities for each galaxy of having a spheroid, a disk, presenting an irregularity, being compact or point source and being unclassifiable. ConvNets are able to predict the fractions of votes given a galaxy image with zero bias and $\sim10\%$ scatter. The fraction of miss-classifications is less than $1\%$. Our classification scheme represents a major improvement with respect to CAS (Concentration-Asymmetry-Smoothness)-based methods, which hit a $20-30\%$ contamination limit at high z. The catalog is released with the present paper via the $\href{http://rainbowx.fis.ucm.es/Rainbow_navigator_public}{Rainbow\,database}$

Several authors have reported that the dynamical masses of massive compact galaxies ($M_\star \gtrsim 10^{11} \ \mathrm{M_\odot}$, $r_\mathrm{e} \sim 1 \ \mathrm{kpc}$), computed as $M_\mathrm{dyn} = 5.0 \ \sigma_\mathrm{e}^2 r_\mathrm{e} / G$, are lower than their stellar masses $M_\star$. In a previous study from our group, the discrepancy is interpreted as a breakdown of the assumption of homology that underlie the $M_\mathrm{dyn}$ determinations. Here, we present new spectroscopy of six redshift $z \approx 1.0$ massive compact ellipticals from the Extended Groth Strip, obtained with the 10.4 m Gran Telescopio Canarias. We obtain velocity dispersions in the range $161-340 \ \mathrm{km \ s^{-1}}$. As found by previous studies of massive compact galaxies, our velocity dispersions are lower than the virial expectation, and all of our galaxies show $M_\mathrm{dyn} < M_\star$ (assuming a Salpeter initial mass function). Adding data from the literature, we build a sample covering a range of stellar masses and compactness in a narrow redshift range $\mathit{z \approx 1.0}$. This allows us to exclude systematic effects on the data and evolutionary effects on the galaxy population, which could have affected previous studies. We confirm that mass discrepancy scales with galaxy compactness. We use the stellar mass plane ($M_\star$, $\sigma_\mathrm{e}$, $r_\mathrm{e}$) populated by our sample to constrain a generic evolution mechanism. We find that the simulations of the growth of massive ellipticals due to mergers agree with our constraints and discard the assumption of homology.

The Star Formation Rate (SFR) is one of the main parameters used to analyze the evolution of galaxies through time. The need for recovering the light reprocessed by dust commonly requires the use of low spatial resolution far-infrared data. Recombination-line luminosities provide an alternative, although uncertain dust-extinction corrections based on narrow-band imaging or long-slit spectroscopy have traditionally posed a limit to their applicability. Integral Field Spectroscopy (IFS) is clearly the way to overcome such limitation. We obtain integrated H{\alpha}, ultraviolet (UV) and infrared (IR)-based SFR measurements for 272 galaxies from the CALIFA survey at 0.005 < z < 0.03 using single-band and hybrid tracers. We provide updated calibrations, both global and split by properties (including stellar mass and morphological type), referred to H{\alpha}. The extinction-corrected H{\alpha} luminosity agrees with the updated hybrid SFR estimators based on either UV or H{\alpha} plus IR luminosity over the full range of SFRs (0.03-20 M$_{\odot}$ yr$^{-1}$). The coefficient that weights the amount of energy produced by newly-born stars that is reprocessed by dust on the hybrid tracers, a$_{IR}$, shows a large dispersion. However, it does not became increasingly small at high attenuations, as expected if significant highly-obscured H$\alpha$ emission would be missed. Lenticulars, early-type spirals and type-2 AGN host galaxies show smaller coefficients due to the contribution of optical photons and AGN to dust heating. In the Local Universe the H{\alpha} luminosity derived from IFS observations can be used to measure SFR, at least in statistically-significant, optically-selected galaxy samples. The analysis of the SFR calibrations by galaxies properties could be potentially used by other works to study the impact of different selection criteria in the SFR values derived.

[abridged] We quantify the morphological evolution of z~0 massive galaxies ($M*/M_\odot\sim10^{11}$) from z~3 in the 5 CANDELS fields. The progenitors are selected using abundance matching techniques to account for the mass growth. The morphologies strongly evolve from z~3. At z<1, the population matches the massive end of the Hubble sequence, with 30% of spheroids, 50% of galaxies with equally dominant disk and bulge components and 20% of disks. At z~2-3 there is a majority of irregular systems (~60-70%) with still 30% of spheroids. We then analyze the SFRs, gas fractions and structural properties for the different morphologies independently. Our results suggest two distinct channels for the growth of bulges in massive galaxies. Around 30-40% were already bulges at z~2.5, with low average SFRs and gas-fractions (10-15%), high Sersic indices (n>3-4) and small effective radii ($R_e$~1 kpc) pointing towards an early formation through gas-rich mergers or VDI. Between z~ 2.5 and z~0, they rapidly increase their size by a factor of ~4-5, become all passive but their global morphology remains unaltered. The structural evolution is independent of the gas fractions, suggesting that it is driven by ex-situ events. The remaining 60% experience a gradual morphological transformation, from clumpy disks to more regular bulge+disks systems, essentially happening at z>1. It results in the growth of a significant bulge component (n~3) for 2/3 of the systems possibly through the migration of clumps while the remaining 1/3 keeps a rather small bulge (n~1.5-2). The transition phase between disturbed and relaxed systems and the emergence of the bulge is correlated with a decrease of the star formation activity and the gas fractions. The growth of the effective radii scales roughly with $H(z)^{-1}$ and it is therefore consistent with the expected growth of disks in galaxy haloes.

We present the analysis of the integrated spectral energy distribution (SED) from the ultraviolet (UV) to the far-infrared and H$\alpha$ of a sample of 29 local systems and individual galaxies with infrared (IR) luminosities between 10^11 Lsun and 10^11.8 Lsun. We have combined new narrow-band H$\alpha$+[NII] and broad-band g, r optical imaging taken with the Nordic Optical Telescope (NOT), with archival GALEX, 2MASS, Spitzer, and Herschel data. The SEDs (photometry and integrated H$\alpha$ flux) have been fitted with a modified version of the MAGPHYS code using stellar population synthesis models for the UV-near-IR range and thermal emission models for the IR emission taking into account the energy balance between the absorbed and re-emitted radiation. From the SED fits we derive the star-formation histories (SFH) of these galaxies. For nearly half of them the star-formation rate appears to be approximately constant during the last few Gyrs. In the other half, the current star-formation rate seems to be enhanced by a factor of 3-20 with respect to that occured ~1 Gyr ago. Objects with constant SFH tend to be more massive than starbursts and they are compatible with the expected properties of a main-sequence (M-S) galaxy. Likewise, the derived SFHs show that all our objects were M-S galaxies ~1 Gyr ago with stellar masses between 10^10.1 and 10^11.5 Msun. We also derived from our fits the average extinction (A_v=0.6-3 mag) and the polycyclic aromatic hydrocarbons (PAH) luminosity to L(IR) ratio (0.03-0.16). We combined the A_v with the total IR and H$\alpha$ luminosities into a diagram which can be used to identify objects with rapidly changing (increasing or decreasing) SFR during the last 100 Myr.

We present a study of galaxies showing mid-infrared variability in the deepest Spitzer/MIPS 24 $\mu$m surveys in the GOODS-South field. We divide the dataset in epochs and subepochs to study the long-term (months-years) and the short-term (days) variability. We use a $\chi^2$-statistics method to select AGN candidates with a probability $\leq$ 1% that the observed variability is due to statistical errors alone. We find 39 (1.7% of the parent sample) sources that show long-term variability and 55 (2.2% of the parent sample) showing short-term variability. We compare our candidates with AGN selected in the X-ray and radio bands, and AGN candidates selected by their IR emission. Approximately, 50% of the MIPS 24 $\mu$m variable sources would be identified as AGN with these other methods. Therefore, MIPS 24 $\mu$m variability is a new method to identify AGN candidates, possibly dust obscured and low luminosity AGN that might be missed by other methods. However, the contribution of the MIPS 24 $\mu$m variable identified AGN to the general AGN population is small ($\leq$ 13%) in GOODS-South.

We explore the stellar initial mass function (IMF) of a sample of 49 massive quiescent galaxies (MQGs) at 0.9$<$z$<$1.5. We base our analysis on intermediate resolution spectro-photometric data in the GOODS-N field taken in the near-infrared and optical with the HST/WFC3 G141 grism and the Survey for High-z Absorption Red and Dead Sources (SHARDS). To constrain the slope of the IMF, we have measured the TiO$_2$ spectral feature, whose strength depends strongly on the content of low-mass stars, as well as on stellar age. Using ultraviolet to near-infrared individual and stacked spectral energy distributions, we have independently estimated the stellar ages of our galaxies. Knowing the age of the stellar population, we interpret the strong differences in the TiO$_2$ feature as an IMF variation. In particular, for the heaviest z$\sim$1 MQGs (M$>$10$^{11}$Msun) we find an average age of 1.7$\pm$0.3 Gyr and a bottom-heavy IMF ($\Gamma_b$=3.2$\pm$0.2). Lighter MQGs (2$\times$10$^{10}$$<$M$<$10$^{11}$ Msun) at the same redshift are younger on average (1.0$\pm$0.2 Gyr) and present a shallower IMF slope ($\Gamma_b=2.7^{+0.3}_{-0.4}$). Our results are in good agreement with the findings about the IMF slope in early-type galaxies of similar mass in the present-day Universe. This suggests that the IMF, a key characteristic of the stellar populations in galaxies, is bottom-heavier for more massive galaxies and has remained unchanged in the last $\sim$8 Gyr.

Lenticular galaxies (S0s) are more likely to host antitruncated (Type-III) stellar discs than galaxies of later Hubble types. Major mergers are popularly considered too violent mechanisms to form these breaks. We have investigated whether major mergers can result into S0-like remnants with realistic antitruncated stellar discs or not. We have analysed 67 relaxed S0 and E/S0 remnants resulting from dissipative N-body simulations of major mergers from the GalMer database. We have simulated realistic R-band surface brightness profiles of the remnants to identify those with antitruncated stellar discs. Their inner and outer discs and the breaks are quantitatively characterized to compare with real data. Nearly 70% of our S0-like remnants are antitruncated, meaning that major mergers that result in S0s have a high probability of producing Type-III stellar discs. Our remnants lie on top of the extrapolations of the observational trends (towards brighter magnitudes and higher break radii) in several photometric diagrams. In scale-free photometric diagrams, simulations and observations overlap and the remnants reproduce the observational trends, so the physical mechanism after antitruncations is highly scalable. We have found novel photometric scaling relations between the characteristic parameters of the antitruncations in real S0s, which are also reproduced by our simulations. The trends in all the photometric planes can be derived from three basic scaling relations that real and simulated Type-III S0s fulfill: h_i \prop R_brkIII, h_o \prop R_brkIII, and mu_brkIII \prop R_brkIII, where h_i and h_o are the scalelenghts of the inner and outer discs, and mu_brkIII and R_brkIII are the surface brightness and radius of the breaks. Mayor mergers provide a feasible mechanism to form realistic antitruncated S0 galaxies (abridged).

We present a study of galaxies showing mid-infrared variability in data taken in the deepest Spitzer/MIPS 24 $\mu$m surveys in the GOODS-South field. We divide the dataset in epochs and subepochs to study the long-term (months-years) and the short-term (days) variability. We use a $\chi^2$-statistics method to select AGN candidates with a probability $\leq$ 1% that the observed variability is due to statistical errors alone. We find 39 (1.7% of the parent sample) sources that show long-term variability and 55 (2.2% of the parent sample) showing short-term variability. That is, 0.03 sources $\times$ arcmin$^{-2}$ for both, long-term and short-term variable sources. After removing the expected number of false positives inherent to the method, the estimated percentages are 1.0% and 1.4% of the parent sample for the long-term and short-term respectively. We compare our candidates with AGN selected in the X-ray and radio bands, and AGN candidates selected by their IR emission. Approximately, 50% of the MIPS 24 $\mu$m variable sources would be identified as AGN with these other methods. Therefore, MIPS 24 $\mu$m variability is a new method to identify AGN candidates, possibly dust obscured and low luminosity AGN, that might be missed by other methods. However, the contribution of the MIPS 24 $\mu$m variable identified AGN to the general AGN population is small ($\leq$ 13%) in GOODS-South.

(Abridged) We probe the merging channel of massive galaxies over the z=0.3-1.3 redshift window by studying close pairs in a sample of 238 galaxies with stellar mass >1E11Msun, from the deep (m<26.5AB, 3 sigma) SHARDS survey. SHARDS provides medium band photometry equivalent to low-resolution optical spectra (R~50), allowing us to obtain extremely accurate photometric redshifts (|Dz|/(1+z)~0.55%) and to improve the constraints on the age distribution of the stellar populations. A strong correlation is found between the age difference of central and satellite galaxy and stellar mass ratio, from negligible age differences in major mergers to age differences ~4 Gyr for 1:100 minor mergers. However, this correlation is simply a reflection of the mass-age trend in the general population. The dominant contributor to the growth of massive galaxies corresponds to mass ratios mu=Msat/Mcen>0.3, followed by a decrease in the fractional mass growth rate linearly proportional to log mu, at least down to mu~0.01, suggesting a decreasing role of mergers involving low-mass satellites, especially if dynamical friction timescales are taken into account. A simple model results in an upper limit for the average mass growth rate of massive galaxies of DM/M/Dt~ 0.08+-0.02 per Gyr, over the z<1 range, with a ~70% fractional contribution from (major) mergers with mu>0.3. The majority of the stellar mass contributed by mergers does not introduce significantly younger populations, in agreement with the small radial age gradients observed in present-day early-type galaxies.

We analyse the stellar populations in the host galaxies of 53 X-ray selected optically dull active galactic nuclei (AGN) at 0.34<z<1.07 with ultra-deep (m=26.5) optical medium-band (R~50) photometry from the Survey for High-z Absorption Red and Dead Sources (SHARDS). The spectral resolution of SHARDS allows us to consistently measure the strength of the 4000 AA break, Dn(4000), a reliable age indicator for stellar populations. We confirm that most X-ray selected moderate-luminosity AGN (L_X<10^44 erg/s) are hosted by massive galaxies (typically M*>10^10.5 M_sun) and that the observed fraction of galaxies hosting an AGN increases with the stellar mass. A careful selection of random control samples of inactive galaxies allows us to remove the stellar mass and redshift dependencies of the AGN fraction to explore trends with several stellar age indicators. We find no significant differences in the distribution of the rest-frame U-V colour for AGN hosts and inactive galaxies, in agreement with previous results. However, we find significantly shallower 4000 AA breaks in AGN hosts, indicative of younger stellar populations. With the help of a model-independent determination of the extinction, we obtain extinction-corrected U-V colours and light-weighted average stellar ages. We find that AGN hosts have younger stellar populations and higher extinction compared to inactive galaxies with the same stellar mass and at the same redshift. We find a highly significant excess of AGN hosts with Dn(4000)~1.4 and light weighted average stellar ages of 300-500 Myr, as well as a deficit of AGN in intrinsic red galaxies. We interpret failure in recognising these trends in previous studies as a consequence of the balancing effect in observed colours of the age-extinction degeneracy.

We combine multiwavelength data in the AEGIS-XD and C-COSMOS surveys to measure the typical dark matter halo mass of X-ray selected AGN [Lx(2-10keV)>1e42 erg/s] in comparison with far-infrared selected star-forming galaxies detected in the Herschel/PEP survey (PACS Evolutionary Probe; Lir>1e11 solar) and quiescent systems at z~1. We develop a novel method to measure the clustering of extragalactic populations that uses photometric redshift Probability Distribution Functions in addition to any spectroscopy. This is advantageous in that all sources in the sample are used in the clustering analysis, not just the subset with secure spectroscopy. The method works best for large samples. The loss of accuracy because of the lack of spectroscopy is balanced by increasing the number of sources used to measure the clustering. We find that X-ray AGN, far-infrared selected star-forming galaxies and passive systems in the redshift interval 0.6<z<1.4 are found in halos of similar mass, $\log M_{DMH}/(M_{\odot}\,h^{-1})\approx13.0$. We argue that this is because the galaxies in all three samples (AGN, star-forming, passive) have similar stellar mass distributions, approximated by the J-band luminosity. Therefore all galaxies that can potentially host X-ray AGN, because they have stellar masses in the appropriate range, live in dark matter haloes of $\log M_{DMH}/(M_{\odot}\,h^{-1})\approx13.0$ independent of their star-formation rates. This suggests that the stellar mass of X-ray AGN hosts is driving the observed clustering properties of this population. We also speculate that trends between AGN properties (e.g. luminosity, level of obscuration) and large scale environment may be related to differences in the stellar mass of the host galaxies.

Chandra data in the COSMOS, AEGIS-XD and 4Ms CDFS are combined with optical/near-IR photometry to determine the rest-frame U-V vs V-J colours of X-ray AGN hosts at mean redshifts 0.40 and 0.85. This combination of colours (UVJ) provides an efficient means of separating quiescent from star-forming, including dust reddened, galaxies. Morphological information emphasises differences between AGN split by their UVJ colours. AGN in quiescent galaxies are dominated by spheroids, while star-forming hosts are split between bulges and disks. The UVJ diagram of AGN hosts is then used to set limits on the accretion density associated with evolved and star-forming systems. Most of the black hole growth since z~1 is associated with star-forming hosts. Nevertheless, ~15-20% of the X-ray luminosity density since z~1, is taking place in the quiescent region of the UVJ diagram. For the z~0.40 subsample, there is tentative evidence (2sigma significance), that AGN split by their UVJ colours differ in Eddington ratio. AGN in star-forming hosts dominate at high Eddington ratios, while AGN in quiescent hosts become increasingly important as a fraction of the total population toward low Eddington ratios. At higher redshift, z~0.8, such differences are significant at the 2sigma level only at Eddington ratios >1e-3. These findings are consistent with scenarios in which diverse accretion modes are responsible for the build-up of SMBHs at the centres of galaxies. We compare our results with the GALFORM semi-analytic model, which postulates two black hole fuelling modes, the first linked to star-formation and the second occuring in passive galaxies. GALFORM predicts a larger fraction of black hole growth in quiescent galaxies at z<1, compared to the data. Relaxing the strong assumption of the model that passive AGN hosts have zero star-formation rate could reconcile this disagreement.

In the course of our 870um APEX/LABOCA follow up of the Herschel Lensing Survey we have detected a source in AS1063 (RXC J2248.7-4431), that has no counterparts in any of the Herschel PACS/SPIRE bands, it is a Herschel 'drop-out' with S_870/S_500>0.5. The 870um emission is extended and centered on the brightest cluster galaxy suggesting either a multiply imaged background source or substructure in the Sunyaev-Zel'dovich (SZ) increment due to inhomogeneities in the hot cluster gas of this merging cluster. We discuss both interpretations with emphasis on the putative lensed source. Based on the observed properties and on our lens model we find that this source could be the first SMG with a moderate far infrared luminosity (L_FIR<10^12 L_sol) detected so far at z>4. In deep HST observations we identified a multiply imaged z~6 source and we measured its spectroscopic redshift z=6.107 with VLT/FORS. This source could be associated with the putative SMG but it is most likely offset spatially by 10-30kpc and they could be interacting galaxies. With a FIR luminosity in the range [5-15]x10^{11} L_sol corresponding to a star formation rate in the range [80-260]M_sol/yr, this SMG would be more representative than the extreme starbursts usually detected at z>4. With a total magnification of ~25 it would open a unique window to the 'normal' dusty galaxies at the end of the epoch of reionization.

Accretion onto SMBH is believed to occur mostly in obscured AGN. Such objects are proving rather elusive in surveys of distant galaxies, including those at X-ray energies. Our main goal is to determine whether the revised IRAC criteria of Donley et al. (2012) (objects with an IR power-law spectral shape), are effective at selecting X-ray type-2 AGN. We present the results from the X-ray spectral analysis of 147 AGN selected by cross-correlating the highest spectral quality ultra-deep XMM-Newton and the Spitzer/IRAC catalogues in the CDF-S. Consequently it is biased towards sources with high S/N X-ray spectra. In order to measure the amount of intrinsic absorption in these sources, we adopt a simple X-ray spectral model that includes a power-law modified by intrinsic absorption and a possible soft X-ray component. We find 21/147 sources to be heavily absorbed but the uncertainties in their obscuring column densities do not allow us to confirm their Compton-Thick nature without resorting to additional criteria. Although IR power-law galaxies are less numerous in our sample than IR non-power-law galaxies (60 versus 87 respectively), we find that the fraction of absorbed (N_{H} > 10^{22} cm^{-2}) AGN is significantly higher (at about 3 sigma level) for IR-power-law sources ($\sim$2/3) than for those sources that do not meet this IR selection criteria ($\sim$1/2). This behaviour is particularly notable at low luminosities, but it appears to be present, although with a marginal significance, at all luminosities. We therefore conclude that the IR power-law method is efficient in finding X-ray-absorbed sources. We would then expect that the long-sought dominant population of absorbed AGN is abundant among IR power-law spectral shape sources not detected in X-rays.

We report on results from the analysis of a stellar mass-selected (log M*>9.0) sample of 1644 galaxies at 0.65<z<1.1 with ultra-deep (m<26.5) optical medium-band (R~50) photometry from the Survey for High-z Absorption Red and Dead Sources (SHARDS). The spectral resolution of SHARDS allows us to consistently measure the strength of the 4000 Angstrom spectral break [Dn(4000), an excellent age indicator for the stellar populations of quiescent galaxies] for all galaxies at z~0.9 down to log M*9. The Dn(4000) index cannot be resolved from broad-band photometry, and measurements from optical spectroscopic surveys are typically limited to galaxies at least x10 more massive. When combined with the rest-frame U-V colour, Dn(4000) provides a powerful diagnostic of the extinction affecting the stellar population that is relatively insensitive to degeneracies with age, metallicity or star formation history. We use this novel approach to estimate de-reddened colours and light-weighted stellar ages for individual sources. We explore the relationships linking stellar mass, (U-V), and Dn(4000) for the sources in the sample, and compare them to those found in local galaxies. The main results are: a) both Dn(4000) and (U-V) correlate with M*. The dispersion in Dn(4000) values at a given M* increases with M*, while the dispersion for (U-V) decreases due to the higher average extinction prevalent in massive star-forming galaxies. b) for massive galaxies, we find a smooth transition between the blue cloud and red sequence in the intrinsic U-V colour, in contrast with other recent results. c) at a fixed stellar age, we find a positive correlation between extinction and stellar mass. d) the fraction of sources with declining or halted star formation increases steeply with the stellar mass, from ~5% at log M*~9.0-9.5 to ~80% at log M*>11, in agreement with downsizing scenarios.

SHARDS is an ongoing ESO/GTC large program that is obtaining ultra-deep photometry of the GOODS-North field in 24 medium-band filters (reaching m=26.5 AB in all bands) in the 500-950 nm range with GTC/OSIRIS. It is designed to study the properties of high-z massive galaxies, but it can also provide very valuable information about the population of AGN at intermediate redshifts (z~0.5-2). Here we present preliminary results on a study of the stellar populations and star formation activity in the host galaxies of X-ray selected optically faint AGN at 0.6<z<1.1. We demonstrate that the SHARDS photometry provides a reliable measurement of the 4000 \AA{} break (Dn(4000), an indicator for the age of stellar populations) down to m=26.5 AB. We confirm that most X-ray selected AGN are hosted by massive galaxies (typically M*>10^10.5 Mo) and that the fraction of galaxies hosting an AGN increases with the stellar mass. AGN-hosts have restframe U-V and Dn(4000) comparable to those of non-active galaxies of the same mass, that is, they do not appear to have on average younger stellar populations, unlike in the Local Universe. Nevertheless, z~1 AGN hosts show an excess of IR emission at {\lambda}>3{\mu}m compared to non-AGN galaxies which might indicate increased star formation rates. In addition, the frequency of AGN in massive galaxies is about twice higher for the ones with young stellar populations compared to the older ones.

Using a compilation of 379 massive (stellar mass M > 10^{11} M_Sun) spheroid-like galaxies from the near-infrared Palomar/DEEP-2 survey, we have probed, up to z~1, whether the presence of companions depends on the size of the host galaxies. We have explored the presence of companions with mass ratios down to 1:10 and 1:100, with respect to the central massive galaxy, and within a projected distance of 30, 50 and 100 kpc of these objects. We find evidence for these companions being equally distributed around both compact and extended massive spheroids. This finding suggests that, at least since z~1, the merger activity in these objects is rather homogeneous across the whole population and its merger history is not affected for the size of the host galaxy. Our result could indicate that both compact and extended massive spheroid-like galaxies are growing in size at the same rate.

We demonstrate estimating the total infrared luminosity, LIR, and star formation rates (SFRs) of star-forming galaxies at redshift 0 < z < 2.8 from single-band 24 micron observations, using local spectral energy distributions (SED) templates without introducing additional free parameters. Our method is based on characterizing the SEDs of galaxies as a function of their LIR surface density, which is motivated by the indications that the majority of IR luminous star-forming galaxies at 1 < z < 3 have extended star-forming regions, in contrast to the strongly nuclear concentrated, merger-induced starbursts in local luminous and ultraluminous IR galaxies. We validate our procedure for estimating LIR by comparing the resulting LIRs with those measured from far-IR observations at 0 < z < 2.8. AGNs were excluded using X-ray and 3.6-8.0 micron observations, which are generally available in deep cosmological survey fields. The Gaussian fits to the distribution of the discrepancies between the LIR measurements from single-band 24 micron and Herschel observations have sigma < 0.1 dex, with ~10% of objects disagreeing by more than 0.2 dex. Since the 24 micron estimates are based on SEDs for extended galaxies, this agreement suggests that ~90% of IR galaxies at high z are indeed much more physically extended than local counterparts of similar LIR, consistent with recent independent studies of the fractions of galaxies forming stars in the main-sequence and starburst modes, respectively. Because we have not introduced empirical corrections to enhance these estimates, in principle, our method should be applicable to lower luminosity galaxies. This will enable use of the 21 micron band of the Mid-Infrared Instrument (MIRI) on board the JWST to provide an extremely sensitive tracer of obscured SFR in individual star-forming galaxies across the peak of the cosmic star formation history.

Nearby star-forming regions are ideal laboratories to study high-energy emission of different stellar populations, from very massive stars to brown dwarfs. NGC 2023 is a reflection nebula situated on the south of the Flame Nebula (NGC 2024) and at the edge of the H ii region IC 434, which also contains the Horsehead Nebula (Barnard 33). NGC 2023, NGC 2024, Barnard 33 and the surroundings of the O-type supergiant star {\zeta} Ori constitute the south part of the Orion B molecular complex. In this work, we present a comprehensive study of X-ray emitters in the region of NGC 2023 and its ?surroundings. We combine optical and infrared data to determine physical properties (mass, temperature, luminosity, presence of accretion disks) of the stars detected in an XMM-Newton observation. This study has allowed us to analyze spectral energy distribution of these stars for the first time and determine their evolutionary stage. Properties of the X-ray emitting plasma of these stars are compared to those found in other nearby star-forming regions. The results indicate that the stars that are being formed in this region have characteristics, in terms of physical properties and luminosity function, similar to those found in the Taurus-Auriga molecular complex.

Using far-infrared imaging from the "Herschel Lensing Survey", we derive dust properties of spectroscopically-confirmed cluster member galaxies within two massive systems at z~0.3: the merging Bullet Cluster and the more relaxed MS2137.3-2353. Most star-forming cluster sources (~90%) have characteristic dust temperatures similar to local field galaxies of comparable infrared (IR) luminosity (T_dust ~ 30K). Several sub-LIRG (L_IR < 10^11 L_sun) Bullet Cluster members are much warmer (T_dust > 37K) with far-infrared spectral energy distribution (SED) shapes resembling LIRG-type local templates. X-ray and mid-infrared data suggest that obscured active galactic nuclei do not contribute significantly to the infrared flux of these "warm dust" galaxies. Sources of comparable IR-luminosity and dust temperature are not observed in the relaxed cluster MS2137, although the significance is too low to speculate on an origin involving recent cluster merging. "Warm dust" galaxies are, however, statistically rarer in field samples (> 3sigma), indicating that the responsible mechanism may relate to the dense environment. The spatial distribution of these sources is similar to the whole far-infrared bright population, i.e. preferentially located in the cluster periphery, although the galaxy hosts tend towards lower stellar masses (M_* < 10^10 M_sun). We propose dust stripping and heating processes which could be responsible for the unusually warm characteristic dust temperatures. A normal star-forming galaxy would need 30-50% of its dust removed (preferentially stripped from the outer reaches, where dust is typically cooler) to recover a SED similar to a "warm dust" galaxy. These progenitors would not require a higher IR-luminosity or dust mass than the currently observed normal star-forming population.

We present the analysis of the full bi-dimensional optical spectral cube of the nearby spiral galaxy NGC 5668, observed with the PPAK IFU at the Calar Alto observatory 3.5m telescope. We make use of broad-band imaging to provide further constraints on the evolutionary history of the galaxy. This dataset will allow us to improve our understanding of the mechanisms that drive the evolution of disks. We investigated the properties of 62 H II regions and concentric rings in NGC 5668 and derived maps in ionized-gas attenuation and chemical (oxygen) abundances. We find that, while inwards of r\,$\sim\,36",\sim$\,4.4kpc\,$\sim$\,0.36\,$(\frac {D_{25}}{2})$ the derived O/H ratio follows the radial gradient typical of spiral galaxies, the abundance gradient beyond r$\sim36"$ flattens out. The analysis of the multi-wavelength surface brightness profiles of NGC 5668 is performed by fitting these profiles with those predicted by chemo-spectrophotometric evolutionary models of galaxy disks. From this, we infer a spin and circular velocity of $\lambda$=0.053 and v$_{c}$=167\,km\,s$^{-1}$, respectively. The metallicity gradient and rotation curve predicted by this best-fitting galaxy model nicely match the values derived from the IFU observations, especially within r\,$\sim36\arcsec$. The same is true for the colors despite of some small offsets and a reddening in the bluest colors beyond that radius. On the other hand, deviations of some of these properties in the outer disk indicate that a secondary mechanism, possibly gas transfer induced by the presence of a young bar, must have played a role in shaping the recent chemical and star formation histories of NGC 5668.

Accretion of minor satellites has been postulated as the most likely mechanism to explain the significant size evolution of the massive galaxies over cosmic time. Using a sample of 629 massive (Mstar~10^11 Msun) galaxies from the near-infrared Palomar/DEEP-2 survey, we explore which fraction of these objects has satellites with 0.01 Msat < Mcentral < 1 (1:100) up to z=1 and which fraction has satellites with 0.1 Msat < Mcentral < 1 (1:10) up to z=2 within a projected radial distance of 100 kpc. We find that the fraction of massive galaxies with satellites, after the background correction, remains basically constant and close to ~30% for satellites with a mass ratio down to 1:100 up to z=1, and ~15% for satellites with a 1:10 mass ratio up to z=2. The family of spheroid-like massive galaxies presents a 2-3 times larger fraction of objects with satellites than the group of disk-like massive galaxies. A crude estimation of the number of 1:3 mergers a massive spheroid-like galaxy experiences since z~2 is around 2. For a disk-like galaxy this number decreases to ~1.

SHARDS, an ESO/GTC Large Program, is an ultra-deep (26.5 mag) spectro-photometric survey with GTC/OSIRIS designed to select and study massive passively evolving galaxies at z=1.0-2.3 in the GOODS-N field using a set of 24 medium-band filters (FWHM~17 nm) covering the 500-950 nm spectral range. Our observing strategy has been planned to detect, for z>1 sources, the prominent Mg absorption feature (at rest-frame ~280 nm), a distinctive, necessary, and sufficient feature of evolved stellar populations (older than 0.5 Gyr). These observations are being used to: (1) derive for the first time an unbiased sample of high-z quiescent galaxies, which extends to fainter magnitudes the samples selected with color techniques and spectroscopic surveys; (2) derive accurate ages and stellar masses based on robust measurements of spectral features such as the Mg(UV) or D(4000) indices; (3) measure their redshift with an accuracy Delta(z)/(1+z)<0.02; and (4) study emission-line galaxies (starbursts and AGN) up to very high redshifts. The well-sampled optical SEDs provided by SHARDS for all sources in the GOODS-N field are a valuable complement for current and future surveys carried out with other telescopes (e.g., Spitzer, HST, and Herschel).

We present far-infrared (FIR) analysis of 68 Brightest Cluster Galaxies (BCGs) at 0.08 < z < 1.0. Deriving total infrared luminosities directly from Spitzer and Herschel photometry spanning the peak of the dust component (24-500um), we calculate the obscured star formation rate (SFR). 22(+6.2,-5.3)% of the BCGs are detected in the far-infrared, with SFR= 1-150 M_sun/yr. The infrared luminosity is highly correlated with cluster X-ray gas cooling times for cool-core clusters (gas cooling time <1 Gyr), strongly suggesting that the star formation in these BCGs is influenced by the cluster-scale cooling process. The occurrence of the molecular gas tracing Ha emission is also correlated with obscured star formation. For all but the most luminous BCGs (L_TIR > 2x10^11 L_sun), only a small (<0.4 mag) reddening correction is required for SFR(Ha) to agree with SFR_FIR. The relatively low Ha extinction (dust obscuration), compared to values reported for the general star-forming population, lends further weight to an alternate (external) origin for the cold gas. Finally, we use a stacking analysis of non-cool-core clusters to show that the majority of the fuel for star formation in the FIR-bright BCGs is unlikely to originate form normal stellar mass loss.

We present measurements of the mean mid-infrared-to-submillimeter flux densities of massive (M\ast \approx 2 \times 10^11 Msun) galaxies at redshifts 1.7 < z < 2.9, obtained by stacking positions of known objects taken from the GOODS NICMOS Survey (GNS) catalog on maps: at 24 {\mu}m (Spitzer/MIPS); 70, 100, and 160{\mu}m (Herschel/PACS); 250, 350, 500{\mu}m (BLAST); and 870{\mu}m (LABOCA). A modified blackbody spectrum fit to the stacked flux densities indicates a median [interquartile] star-formation rate of SFR = 63 [48, 81] Msun yr^-1 . We note that not properly accounting for correlations between bands when fitting stacked data can significantly bias the result. The galaxies are divided into two groups, disk-like and spheroid-like, according to their Sersic indices, n. We find evidence that most of the star formation is occurring in n \leq 2 (disk-like) galaxies, with median [interquartile] SFR = 122 [100,150] Msun yr^-1, while there are indications that the n > 2 (spheroid-like) population may be forming stars at a median [interquartile] SFR = 14 [9,20] Msun yr^-1, if at all. Finally, we show that star formation is a plausible mechanism for size evolution in this population as a whole, but find only marginal evidence that it is what drives the expansion of the spheroid-like galaxies.

We present the results of a comparison between the optical morphologies of a complete sample of 46 southern 2Jy radio galaxies at intermediate redshifts (0.05<z<0.7) and those of two control samples of quiescent early-type galaxies: 55 ellipticals at redshifts z<0.01 from the Observations of Bright Ellipticals at Yale (OBEY) survey, and 107 early-type galaxies at redshifts 0.2<z<0.7 in the Extended Groth Strip (EGS). Based on these comparisons, we discuss the role of galaxy interactions in the triggering of powerful radio galaxies (PRGs). We find that a significant fraction of quiescent ellipticals at low and intermediate redshifts show evidence for disturbed morphologies at relatively high surface brightness levels, which are likely the result of past or on-going galaxy interactions. However, the morphological features detected in the galaxy hosts of the PRGs (e.g. tidal tails, shells, bridges, etc.) are up to 2 magnitudes brighter than those present in their quiescent counterparts. Indeed, if we consider the same surface brightness limits, the fraction of disturbed morphologies is considerably smaller in the quiescent population (53% at z<0.2 and 48% at 0.2<z<0.7) than in the PRGs (93% at z<0.2 and 95% at 0.2<z<0.7 considering strong-line radio galaxies only). This supports a scenario in which PRGs represent a fleeting active phase of a subset of the elliptical galaxies that have recently undergone mergers/interactions. However, we demonstrate that only a small proportion (<20%) of disturbed early-type galaxies are capable of hosting powerful radio sources.

In this work we analyze the physical properties of a sample of 153 star forming galaxies at z~0.84, selected by their H-alpha flux with a NB filter. B-band luminosities of the objects are higher than those of local star forming galaxies. Most of the galaxies are located in the blue cloud, though some objects are detected in the green valley and in the red sequence. After the extinction correction is applied virtually all these red galaxies move to the blue sequence, unveiling their dusty nature. A check on the extinction law reveals that the typical extinction law for local starbursts is well suited for our sample but with E(B-V)_stars=0.55 E(B-V)_gas. We compare star formation rates (SFR) measured with different tracers (H-alpha, UV and IR) finding that they agree within a factor of three after extinction correction. We find a correlation between the ratios SFR_FUV/SFR_H-alpha, SFR_IR/SFR_H-alpha and the EW(H-alpha) (i.e. weighted age) which accounts for part of the scatter. We obtain stellar mass estimations fitting templates to multi-wavelength photometry. The typical stellar mass of a galaxy within our sample is ~10^10 Msun. The SFR is correlated with stellar mass and the specific star formation rate (sSFR) decreases with it, indicating that massive galaxies are less affected by star formation processes than less massive ones. This result is consistent with the downsizing scenario. To quantify this downsizing we estimated the quenching mass M_Q for our sample at z~0.84, finding that it declines from M_Q ~10^12 Msun to M_Q ~8x10^10 Msun at the local Universe.

The Sunyaev-Zel'dovich (SZ) effect is a spectral distortion of the cosmic microwave background as observed through the hot plasma in galaxy clusters. This distortion is a decrement in the CMB intensity for lambda > 1.3 mm, an increment at shorter wavelengths, and small again by lambda ~250 um. As part of the Herschel Lensing Survey (HLS) we have mapped 1E0657-56 (the Bullet cluster) with SPIRE with bands centered at 250, 350 and 500 um and have detected the SZ effect at the two longest wavelengths. The measured SZ effect increment central intensities are Delta I_{0} = 0.097 +- 0.019 MJy sr^{-1} at 350 um and Delta I_{0} = 0.268 +- 0.031 MJy sr^{-1} at 500 um, consistent with the SZ effect spectrum derived from previous measurements at 2 mm. No other diffuse emission is detected. The presence of the finite temperature SZ effect correction is preferred by the SPIRE data at a significance of 2.1 sigma, opening the possibility that the relativistic SZ effect correction can be constrained by SPIRE in a sample of clusters. The results presented here have important ramifications for both sub-mm measurements of galaxy clusters and blank field surveys with SPIRE.

We present preliminary results about the detection of high redshift (U)LIRGs in the Bullet cluster field by the PACS and SPIRE instruments within the Herschel Lensing Survey (HLS) Program. We describe in detail a photometric procedure designed to recover robust fluxes and deblend faint Herschel sources near the confusion noise. The method is based on the use of the positions of Spitzer/MIPS 24 um sources as priors. Our catalogs are able to reliably (5 sigma) recover galaxies with fluxes above 6 and 10 mJy in the PACS 100 and 160 um channels, respectively, and 12 to 18 mJy in the SPIRE bands. We also obtain spectral energy distributions covering the optical through the far-infrared/millimeter spectral ranges of all the Herschel detected sources, and analyze them to obtain independent estimations of the photometric redshift based on either stellar population or dust emission models. We exemplify the potential of the combined use of Spitzer position priors plus independent optical and IR photometric redshifts to robustly assign optical/NIR counterparts to the sources detected by Herschel and other (sub-)mm instruments.

We use deep, five band (100-500um) data from the Herschel Lensing Survey (HLS) to fully constrain the obscured star formation rate, SFR_FIR, of galaxies in the Bullet cluster (z=0.296), and a smaller background system (z=0.35) in the same field. Herschel detects 23 Bullet cluster members with a total SFR_FIR = 144 +/- 14 M_sun yr^-1. On average, the background system contains brighter far-infrared (FIR) galaxies, with ~50% higher SFR_FIR (21 galaxies; 207 +/- 9 M_sun yr^-1). SFRs extrapolated from 24um flux via recent templates (SFR_24) agree well with SFR_FIR for ~60% of the cluster galaxies. In the remaining ~40%, SFR_24 underestimates SFR_FIR due to a significant excess in observed S_100/S_24 (rest frame S_75/S_18) compared to templates of the same FIR luminosity.

The Herschel Lensing Survey (HLS) takes advantage of gravitational lensing by massive galaxy clusters to sample a population of high-redshift galaxies which are too faint to be detected above the confusion limit of current far-infrared/submillimeter telescopes. Measurements from 100-500 micron bracket the peaks of the far-infrared spectral energy distributions of these galaxies, characterizing their infrared luminosities and star formation rates. We introduce initial results from our science demonstration phase observations, directed toward the Bullet cluster (1E0657-56). By combining our observations with LABOCA 870 micron and AzTEC 1.1 mm data we fully constrain the spectral energy distributions of 19 MIPS 24 micron selected galaxies which are located behind the cluster. We find that their colors are best fit using templates based on local galaxies with systematically lower infrared luminosities.This suggests that our sources are not like local ultra-luminous infrared galaxies in which vigorous star formation is contained in a compact highly dust-obscured region. Instead, they appear to be scaled up versions of lower luminosity local galaxies with star formation occurring on larger physical scales.

The Herschel Lensing Survey (HLS) will conduct deep PACS and SPIRE imaging of ~40 massive clusters of galaxies. The strong gravitational lensing power of these clusters will enable us to penetrate through the confusion noise, which sets the ultimate limit on our ability to probe the Universe with Herschel. Here, we present an overview of our survey and a summary of the major results from our Science Demonstration Phase (SDP) observations of the Bullet Cluster (z=0.297). The SDP data are rich, allowing us to study not only the background high-redshift galaxies (e.g., strongly lensed and distorted galaxies at z=2.8 and 3.2) but also the properties of cluster-member galaxies. Our preliminary analysis shows a great diversity of far-infrared/submillimeter spectral energy distributions (SEDs), indicating that we have much to learn with Herschel about the properties of galaxy SEDs. We have also detected the Sunyaev-Zel'dovich (SZ) effect increment with the SPIRE data. The success of this SDP program demonstrates the great potential of the Herschel Lensing Survey to produce exciting results in a variety of science areas.

(Abriged) Our goal here is to provide merger frequencies that encompass both major and minor mergers, derived from close pair statistics. We use B-band luminosity- and mass-limited samples from an Spitzer/IRAC-selected catalogue of GOODS-S. We present a new methodology for computing the number of close companions, Nc, when spectroscopic redshift information is partial. We select as close companions those galaxies separated by 6h^-1 kpc < rp < 21h^-1 kpc in the sky plane and with a difference Delta_v <= 500 km s^-1 in redshift space. We provide Nc for four different B-band-selected samples. Nc increases with luminosity, and its evolution with redshift is faster in more luminous samples. We provide Nc of M_star >= 10^10 M_Sun galaxies, finding that the number including minor companions (mass ratio >= 1/10) is roughly two times the number of major companions alone (mass ratio >= 1/3) in the range 0.2 <= z < 1.1. We compare the major merger rate derived by close pairs with the one computed by morphological criteria, finding that both approaches provide similar merger rates for field galaxies when the progenitor bias is taken into account. Finally, we estimate that the total (major+minor) merger rate is ~1.7 times the major merger rate. Only 30% to 50% of the M_star >= 10^10 M_Sun early-type (E/S0/Sa) galaxies that appear z=1 and z=0 may have undergone a major or a minor merger. Half of the red sequence growth since z=1 is therefore unrelated to mergers.

We study the evolution of galaxy structure since z ~ 1 to the present. From a GOODS-S multi-band catalog we define (blue) luminosity- and mass-weighted samples, limited by M_B <= -20 and M_star >= 10^10 M_Sun, comprising 1122 and 987 galaxies, respectively. We extract early-type (E/S0/Sa) and late-type (Sb-Irr) subsamples by their position in the concentration-asymmetry plane, in which galaxies exhibit a clear bimodality. We find that the early-type fraction, f_ET, rises with cosmic time, with a corresponding decrease in the late-type fraction, f_LT, in both luminosity- and mass-selected samples. However, the evolution of the comoving number density is very different: the decrease in the total number density of M_B <= -20 galaxies since z = 1 is due to the decrease in the late-type population, which accounts for ~75% of the total star-formation rate in the range under study, while the increase in the total number density of M_star >= 10^10 M_Sun galaxies in the same redshift range is due to the evolution of early types. This suggests that we need a structural transformation between late-type galaxies that form stars actively and early-type galaxies in which the stellar mass is located. Comparing the observed evolution with the gas-rich major merger rate in GOODS-S, we infer that only ~20% of the new early-type galaxies with M_star >= 10^10 M_Sun appeared since z ~ 1 can be explained by this kind of mergers, suggesting that minor mergers and secular processes may be the driving mechanisms of the structural evolution of intermediate-mass (M_star ~ 4x10^10 M_Sun) galaxies since z ~ 1.

We present ultraviolet through far-infrared surface brightness profiles for the 75 galaxies in the Spitzer Infrared Nearby Galaxies Survey (SINGS). The imagery used to measure the profiles includes GALEX UV data, optical images from KPNO, CTIO and SDSS, near-IR data from 2MASS, and mid- and far-infrared images from Spitzer. Along with the radial profiles, we also provide multi-wavelength asymptotic magnitudes and several non-parametric indicators of galaxy morphology: the concentration index (C_42), the asymmetry (A), the Gini coefficient (G) and the normalized second-order moment of the brightest 20% of the galaxy's flux (M_20). Our radial profiles show a wide range of morphologies and multiple components (bulges, exponential disks, inner and outer disk truncations, etc.) that vary not only from galaxy to galaxy but also with wavelength for a given object. In the optical and near-IR, the SINGS galaxies occupy the same regions in the C_42-A-G-M_20 parameter space as other normal galaxies in previous studies. However, they appear much less centrally concentrated, more asymmetric and with larger values of G when viewed in the UV (due to star-forming clumps scattered across the disk) and in the mid-IR (due to the emission of Polycyclic Aromatic Hydrocarbons at 8.0 microns and very hot dust at 24 microns).

We present a detailed analysis of the radial distribution of dust properties in the SINGS sample, performed on a set of UV, IR and HI surface brightness profiles, combined with published molecular gas profiles and metallicity gradients. The internal extinction, derived from the TIR-to-FUV luminosity ratio, decreases with radius, and is larger in Sb-Sbc galaxies. The TIR-to-FUV ratio correlates with the UV spectral slope beta, following a sequence shifted to redder UV colors with respect to that of starbursts. The star formation history (SFH) is identified as the main driver of this departure. We have also derived radial profiles of the total dust mass surface density, the fraction of the dust mass contributed by PAHs, the fraction of the dust mass heated by very intense starlight and the intensity of the radiation field heating the grains. The dust profiles are exponential, their radial scale-length being constant from Sb to Sd galaxies (only ~10% larger than the stellar scale-length). Many S0/a-Sab galaxies have central depressions in their dust radial distributions. The PAH abundance increases with metallicity for 12+\log(O/H)<9, and at larger metallicities the trend flattens and even reverses, with the SFH being a plausible underlying driver for this behavior. The dust-to-gas ratio is also well correlated with metallicity and therefore decreases with galactocentric radius.

Aims: We study the major merger fraction in a SPITZER/IRAC-selected catalogue in the GOODS-S field up to z ~ 1 for luminosity- and mass-limited samples. Methods: We select disc-disc merger remnants on the basis of morphological asymmetries, and address three main sources of systematic errors: (i) we explicitly apply morphological K-corrections, (ii) we measure asymmetries in galaxies artificially redshifted to z_d = 1.0 to deal with loss of morphological information with redshift, and (iii) we take into account the observational errors in z and A, which tend to overestimate the merger fraction, though use of maximum likelihood techniques. Results: We obtain morphological merger fractions (f_m) below 0.06 up to z ~ 1. Parameterizing the merger fraction evolution with redshift as f_m(z) = f_m(0) (1+z)^m, we find that m = 1.8 +/- 0.5 for M_B <= -20 galaxies, while m = 5.4 +/- 0.4 for M_star >= 10^10 M_Sun galaxies. When we translate our merger fractions to merger rates (R_m), their evolution, parameterized as R_m(z) = R_m(0) (1+z)^n, is quite similar in both cases: n = 3.3 +/- 0.8 and n = 3.5 +/- 0.4, respectively. Conclusions: Our results imply that only ~8% of today's M_star >= 10^10 M_Sun galaxies have undergone a disc-disc major merger since z ~ 1. In addition, ~21% of this mass galaxies at z ~ 1 have undergone one of these mergers since z ~ 1.5. This suggests that disc-disc major mergers are not the dominant process in the evolution of M_star >= 10^10 M_Sun galaxies since z ~ 1, but may be an important process at z > 1.

We investigate the causes of the different shape of the $K$-band number counts when compared to other bands, analyzing in detail the presence of a change in the slope around $K\sim17.5$. We present a near-infrared imaging survey, conducted at the 3.5m telescope of the Calar Alto Spanish-German Astronomical Center (CAHA), covering two separated fields centered on the HFDN and the Groth field, with a total combined area of $\sim0.27$deg$^{2}$ to a depth of $K\sim19$ ($3\sigma$,Vega). We derive luminosity functions from the observed $K$-band in the redshift range [0.25-1.25], that are combined with data from the references in multiple bands and redshifts, to build up the $K$-band number count distribution. We find that the overall shape of the number counts can be grouped into three regimes: the classic Euclidean slope regime ($d\log N/dm\sim0.6$) at bright magnitudes; a transition regime at intermediate magnitudes, dominated by $M^{\ast}$ galaxies at the redshift that maximizes the product $\phi^{\ast}\frac{dV_{c}}{d\Omega}$; and an $\alpha$ dominated regime at faint magnitudes, where the slope asymptotically approaches -0.4($\alpha$+1) controlled by post-$M^{\ast}$ galaxies. The slope of the $K$-band number counts presents an averaged decrement of $\sim50%$ in the range $15.5<K<18.5$ ($d\log N/dm\sim0.6-0.30$). The rate of change in the slope is highly sensitive to cosmic variance effects. The decreasing trend is the consequence of a prominent decrease of the characteristic density $\phi^{\ast}_{K,obs}$ ($\sim60%$ from $z=0.5$ to $z=1.5$) and an almost flat evolution of $M^{\ast}_{K,obs}$ (1$\sigma$ compatible with $M^{\ast}_{K,obs}=-22.89\pm0.25$ in the same redshift range).

We present the results of an Halpha near-infrared narrow-band survey searching for star-forming galaxies at redshift z=0.84. This work is an extension of our previous narrow-band studies in the optical at lower redshifts. After removal of stars and redshift interlopers (using spectroscopic and photometric redshifts), we build a complete sample of 165 Halpha emitters in the Extended Groth strip and GOODS-N fields with L(Halpha)>10^41 erg/s. We compute the Halpha luminosity function at z=0.84 after corrections for [NII] flux contamination, extinction, systematic errors, and incompleteness. Our sources present an average dust extinction of A(Halpha)=1.5 mag. Adopting Halpha as a surrogate for the instantaneous star formation rate (SFR), we measure a extinction-corrected SFR density of 0.17+-0.03 M_sun/yr/Mpc3. Combining this result to our prior measurements at z=0.02, 0.24, and 0.40, we derive an Halpha-based evolution of the SFR density proportional to (1+z)^beta with beta=3.8+-0.5. This evolution is consistent with that derived by other authors using different SFR tracers.

We present the stellar calibrator sample and the conversion from instrumental to physical units for the 24 micron channel of the Multiband Imaging Photometer for Spitzer (MIPS). The primary calibrators are A stars, and the calibration factor based on those stars is 4.54*10^{-2} MJy sr^{-1} (DN/s)^{-1}, with a nominal uncertainty of 2%. We discuss the data-reduction procedures required to attain this accuracy; without these procdures, the calibration factor obtained using the automated pipeline at the Spitzer Science Center is 1.6% +/- 0.6% lower. We extend this work to predict 24 micron flux densities for a sample of 238 stars which covers a larger range of flux densities and spectral types. We present a total of 348 measurements of 141 stars at 24 micron. This sample covers a factor of ~460 in 24 micron flux density, from 8.6 mJy up to 4.0 Jy. We show that the calibration is linear over that range with respect to target flux and background level. The calibration is based on observations made using 3-second exposures; a preliminary analysis shows that the calibration factor may be 1% and 2% lower for 10- and 30-second exposures, respectively. We also demonstrate that the calibration is very stable: over the course of the mission, repeated measurements of our routine calibrator, HD 159330, show a root-mean-square scatter of only 0.4%. Finally, we show that the point spread function (PSF) is well measured and allows us to calibrate extended sources accurately; Infrared Astronomy Satellite (IRAS) and MIPS measurements of a sample of nearby galaxies are identical within the uncertainties.

We investigate why half of X-ray--selected active galactic nuclei (AGN) in deep surveys lack signs of accretion in their optical spectra. The majority of these ``optically--dull'' AGN are no more than 6 times fainter than their host galaxies in rest-frame R-band; as such, AGN lines are unlikely to be overwhelmed by stellar continuum in at least half the sample. We find that optically--dull AGN have the mid--infrared emission and L(x)/L(IR) ratios characteristic of local Seyferts, suggesting that the cause of optical dullness is not missing UV--optical continua. We compare the morphologies of 22 optically--dull and 9 optically--active AGN at 0.5<z<0.8, and find that optically--dull AGN show a wide range of axis ratio, but optically--active AGN have only very round axis ratios. We conclude that hard X-rays select AGN in host galaxies with a wide range of inclination angle, but only those AGN in the most face-on or spheroidal host galaxies show optical emission lines. Thus, extranuclear dust in the host galaxy plays an important role in hiding the emission lines of optically--dull AGN.

We stack Spitzer 24 micron images for ~7000 galaxies with 0.1<z<1 in the Chandra Deep Field South to probe the thermal dust emission in low-luminosity galaxies over this redshift range. Through stacking, we can detect mean 24 micron fluxes that are more than an order of magnitude below the individual detection limit. We find that the correlations for low and moderate luminosity galaxies between the average L_IR/L_UV and rest-frame B-band luminosity, and between the star formation rate (SFR) and L_IR/L_UV, are similar to those in the local Universe. This verifies that oft-used assumption in deep UV/optical surveys that the dust obscuration-SFR relation for galaxies with SFR < 20 solar mass per year varies little with epoch. We have used this relation to derive the cosmic IR luminosity density from z=1 to z=0.1. The results also demonstrate directly that little of the bolometric luminosity of the galaxy population arises from the faint end of the luminosity function, indicating a relatively flat faint-end slope of the IR luminosity function with a power law index of 1.2+-0.3.