The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will soon start to scan thousands of square degrees of the northern extragalactic sky with a unique set of $56$ optical filters from a dedicated $2.55$m telescope, JST, at the Javalambre Astrophysical Observatory. Before the arrival of the final instrument (a 1.2 Gpixels, 4.2deg$^2$ field-of-view camera), the JST was equipped with an interim camera (JPAS-Pathfinder), composed of one CCD with a 0.3deg$^2$ field-of-view and resolution of 0.23 arcsec pixel$^{-1}$. To demonstrate the scientific potential of J-PAS, with the JPAS-Pathfinder camera we carried out a survey on the AEGIS field (along the Extended Groth Strip), dubbed miniJPAS. We observed a total of $\sim 1$ deg$^2$, with the $56$ J-PAS filters, which include $54$ narrow band (NB, $\rm{FWHM} \sim 145$Angstrom) and two broader filters extending to the UV and the near-infrared, complemented by the $u,g,r,i$ SDSS broad band (BB) filters. In this paper we present the miniJPAS data set, the details of the catalogues and data access, and illustrate the scientific potential of our multi-band data. The data surpass the target depths originally planned for J-PAS, reaching $\rm{mag}_{\rm {AB}}$ between $\sim 22$ and $23.5$ for the NB filters and up to $24$ for the BB filters ($5\sigma$ in a $3$~arcsec aperture). The miniJPAS primary catalogue contains more than $64,000$ sources extracted in the $r$ detection band with forced photometry in all other bands. We estimate the catalogue to be complete up to $r=23.6$ for point-like sources and up to $r=22.7$ for extended sources. Photometric redshifts reach subpercent precision for all sources up to $r=22.5$, and a precision of $\sim 0.3$% for about half of the sample. (Abridged)

We present the photometric determination of the bright-end (L_Lya>10^43.5 erg/s) of the Lya luminosity function (LF) within four redshifts windows in the interval 2.2<z<3.3. Our work is based on the Javalambre Photometric Local Universe Survey (J-PLUS) first data-release, which provides multiple narrow-band measurements over ~1000 deg^2, with limiting magnitude r~22. The analysis of high-z Lya-emitting sources over such a wide area is unprecedented, and allows to select a total of ~14,500 hyper-bright (L_Lya>10^43.3 erg/s) Lya-emitting candidates. We test our selection with two spectroscopic follow-up programs at the GTC telescope, confirming ~89% of the targets as line-emitting sources, with ~64% being genuine z~2.2 QSOs. We extend the 2.2<z<3.3 Lya LF for the first time above L_Lya~10^44 erg/s and down to densities of ~10^-8 Mpc^-3. Our results unveil with high detail the Schechter exponential-decay of the brightest-end of the Lya LF, complementing the power-law component of previous LF determinations at 43.3<Log_10(L_Lya / [erg/s])<44. We measure Phi^*=(3.33+-0.19)x10^-6, Log(L^*)=44.65+-0.65 and alpha=-1.35+-0.84 as an average over the redshifts we probe. These values are significantly different than the typical Schechter parameters measured for the Lya LF of high-z star-forming LAEs. This suggests that z>2 AGN/QSOs (likely dominant in our samples) are described by a structurally different LF than z>2 star-forming LAEs, namely with L^*_QSOs ~ 100 L^*_LAEs and Phi^*_QSOs ~ 10^-3 Phi^*_LAEs. Finally, our method identifies very efficiently as high-z line-emitters sources without previous spectroscopic confirmation, currently classified as stars (~2000 objects in each redshift bin, on average). Assuming a large predominance of Lya-emitting AGN/QSOs in our samples, this supports the scenario by which these are the most abundant class of z>2 Lya emitters at L_Lya>10^43.3 erg/s.

The differences between the inherent stellar populations (SPs) of LAEs and LBGs are a key factor in understanding early galaxy formation and evolution. We have run a set of SP burst-like models for a sample of 1,558 sources at $3.4<z<6.8$ from the Survey for High-$z$ Absorption Red and Dead Sources (SHARDS) over the GOODS-N field. This work focuses on the differences between the three different observational subfamilies of our sample: LAE-LBGs, no-Ly$\alpha$ LBGs and pure LAEs. Single and double SP synthetic spectra were used to model the SEDs, adopting a Bayesian information criterion to analyse under which situations a second SP is required. We find that the sources are well modelled using a single SP in $\sim79\%$ of the cases. The best models suggest that pure LAEs are typically young low mass galaxies ($t\sim26^{+41}_{-25}$ Myr; $M_{\mathrm{star}}\sim5.6^{+12.0}_{-5.5}\times10^{8}\ M_{\odot}$), undergoing one of their first bursts of star formation. On the other hand, no-Ly$\alpha$ LBGs require older SPs ($t\sim71\pm12$ Myr), and they are substantially more massive ($M_{\mathrm{star}}\sim3.5\pm1.1\times10^{9}\ M_{\odot}$). LAE-LBGs appear as the subgroup that more frequently needs the addition of a second SP, representing an old and massive galaxy caught in a strong recent star-forming episode. The relative number of sources found from each subfamily at each $z$ supports an evolutionary scenario from pure LAEs and single SP LAE-LBGs to more massive LBGs. Stellar Mass Functions are also derived, finding an increase of $M^{*}$ with cosmic time and a possible steepening of the low mass slope from $z\sim6$ to $z\sim5$ with no significant change to $z\sim4$. Additionally, we have derived the SFR-$M_{\mathrm{star}}$ relation, finding a $\mathrm{SFR}\propto M_{\mathrm{star}}^{\beta}$ behaviour with negligible evolution from $z\sim4$ to $z\sim6$.

We characterize for the first time the torus properties of an ultra-hard X-ray (14-195 keV) volume-limited (DL<40 Mpc) sample of 24 Seyfert (Sy) galaxies (BCS40 sample). The sample was selected from the Swift/BAT nine month catalog. We use high angular resolution nuclear infrared (IR) photometry and N-band spectroscopy, the CLUMPY torus models and a Bayesian tool to characterize the properties of the nuclear dust. In the case of the Sy1s we estimate the accretion disk contribution to the subarcsecond resolution nuclear IR SEDs (~0.4'') which is, on average, 46+-28, 23+-13 and 11+-5% in the J-, H- and K-bands, respectively. This indicates that the accretion disk templates that assume a steep fall for longer wavelengths than 1 micron might underestimate its contribution to the near-IR emission. Using both optical (broad vs narrow lines) and X-ray (unabsorbed vs absorbed) classifications, we compare the global posterior distribution of the torus model parameters. We confirm that Sy2s have larger values of the torus covering factor (CT~0.95) than Sy1s (CT~0.65) in our volume-limited Seyfert sample. These findings are independent of whether we use an optical or X-ray classification. We find that the torus covering factor remains essentially constant within the errors in our luminosity range and there is no clear dependence with the Eddington ratio. Finally, we find tentative evidence that even an ultra hard X-ray selection is missing a significant fraction of highly absorbed type 2 sources with very high covering factor tori.

We present a sample of 33 Balmer Break Galaxies (BBGs) selected as HST/F160W dropouts in the deepest CANDELS/GOODS fields ($H\gtrsim27.3$~mag) but relatively bright in {\it Spitzer}/IRAC ($[3.6],[4.5]<24.5$~mag), implying red colors (median and quartiles: $\langle H-[3.6]\rangle=3.1^{3.4}_{2.8}$\,mag). Half of these BBGs are newly identified sources. Our BBGs are massive ($\langle \log(\rm{M}/\rm{M}_\odot)\rangle=10.8$) high redshift ($\langle z\rangle=4.8$) dusty ($\langle \rm{A(V)}\rangle=2.0$~mag) galaxies. The SEDs of half of our sample indicate that they are star-forming galaxies with typical specific SFRs 0.5-1.0~Gyr$^{-1}$, qualifying them as main sequence (MS) galaxies at $3<z<6$. One third of those SEDs indicates the presence of prominent emission lines (H$\beta$+$[OIII]$, H$\alpha$$+$[NII]) boosting the IRAC fluxes and red colors. Approximately 20\% of the BBGs are very dusty ($\rm{A(V)}\sim2.5$~mag) starbursts with strong mid-to-far infrared detections and extreme SFRs ($\rm{SFR}>10^{3}\,\rm{M}_\odot/yr$) that place them above the MS. The rest, 30\%, are post-starbursts or quiescent galaxies located $>2\sigma$ below the MS with mass-weighted ages older than 700~Myr. Only 2 of the 33 galaxies are X-ray detected AGN with optical/near-infrared SEDs dominated by stellar emission, but the presence of obscured AGN in the rest of sources cannot be discarded. Our sample accounts for 8\% of the total number density of $\log(\rm{M}/\rm{M}_\odot)>10$ galaxies at $z>3$, but it is a significant contributor (30\%) to the general population of red $\log(\rm{M}/\rm{M}_\odot)>11$ galaxies at $4<z<6$. Finally, our results point out that 1 of every 30 massive $\log(\rm{M}/\rm{M}_\odot)>11$ galaxies in the local Universe was assembled in the first 1.5~Gyr after the Big Bang, a fraction that is not reproduced by state-of-the-art galaxy formation simulations.

We present ALMA Band 6 $^{12}$CO(2--1) line and rest-frame 232GHz continuum observations of the nearby Compton-thick Seyfert galaxy NGC5643 with angular resolutions 0.11-0.26arcsec (9-21pc). The CO(2--1) integrated line map reveals emission from the nuclear and circumnuclear region with a two-arm nuclear spiral extending 10arcsec on each side. The circumnuclear CO(2--1) kinematics can be fitted with a rotating disk, although there are regions with large residual velocities and/or velocity dispersions. The CO(2--1) line profiles of these regions show two different velocity components. One is ascribed to the circular component and the other to the interaction of the AGN outflow, as traced by the [O III]5007AA emission, with molecular gas in the disk a few hundred parsecs from the AGN. On nuclear scales, we detected an inclined CO(2--1) disk (diameter 26 pc, FWHM) oriented almost in a north-south direction. The CO(2--1) nuclear kinematics can be fitted with a rotating disk which appears to be tilted with respect to the large scale disk. There are strong non-circular motions in the central 0.2-0.3 arcsec with velocities of up to 110km/s. In the absence of a nuclear bar, these motions could be explained as radial outflows in the nuclear disk. We estimate a total molecular gas mass for the nuclear disk of $M({\rm H}_2)=1.1\times 10^7\,M_\odot$ and an H$_2$ column density toward the location of the AGN of $N({\rm H}_2)\sim 5 \times 10^{23}\,{\rm cm}^{-2}$, for a standard CO-to-H$_2$ conversion factor. We interpret this nuclear molecular gas disk as the obscuring torus of NGC5643 as well as the collimating structure of the ionization cone.

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.

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.

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.

We use a semi-empirical model to reproduce the 0.1-10um spectral energy distribution (SED) of a sample of 85 luminous quasars. In the model, the continuum emission from the accretion disk as well as the nebular lines are represented by a single empirical template (disk), where differences in the optical spectral index are reproduced by varying the amount of extinction. The near- and mid-infrared emission of the AGN-heated dust is modelled as the combination of two black-bodies (dust). The model fitting shows that the disk and dust components are remarkably uniform among individual quasars, with differences in the observed SED largely accounted for by varying levels of obscuration in the disk as well as differences in the relative luminosity of the disk and dust components. By combining the disk-subtracted SEDs of the 85 quasars, we generate a template for the 1-10um emission of the AGN-heated dust. Additionally, we use a sample of local Seyfert 1 galaxies with full spectroscopic coverage in the 0.37um to 39um range to demonstrate a method for stitching together spectral segments obtained with different PSF and extraction apertures. We show that the disk and dust templates obtained from luminous quasars also reproduce the optical-to-mid-infrared spectra of local Seyfert 1s when the contribution from the host galaxy is properly subtracted.

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.

We present new calculations of the CAT3D clumpy torus models, which now include a more physical dust sublimation model as well as AGN anisotropic emission. These new models allow graphite grains to persist at temperatures higher than the silicate dust sublimation temperature. This produces stronger near-infrared emission and bluer mid-infrared (MIR) spectral slopes. We make a statistical comparison of the CAT3D model MIR predictions with a compilation of sub-arcsecond resolution ground-based MIR spectroscopy of 52 nearby Seyfert galaxies (median distance of 36 Mpc) and 10 quasars. We focus on the AGN MIR spectral index $\alpha_{MIR}$ and the strength of the 9.7 $\mu$m silicate feature $S_{Sil}$. As with other clumpy torus models, the new CAT3D models do not reproduce the Seyfert galaxies with deep silicate absorption ($S_{Sil}<-1$). Excluding those, we conclude that the new CAT3D models are in better agreement with the observed $\alpha_{MIR}$ and $S_{Sil}$ of Seyfert galaxies and quasars. We find that Seyfert 2 are reproduced with models with low photon escape probabilities, while the quasars and the Seyfert 1-1.5 require generally models with higher photon escape probabilities. Quasars and Seyfert 1-1.5 tend to show steeper radial cloud distributions and fewer clouds along an equatorial line-of-sight than Seyfert 2. Introducing AGN anisotropic emission besides the more physical dust sublimation models alleviates the problem of requiring inverted radial cloud distributions (i.e., more clouds towards the outer parts of the torus) to explain the MIR spectral indices of type 2 Seyferts.

Dedicated searches generally find a decreasing fraction of obscured Active Galactic Nuclei (AGN) with increasing AGN luminosity. This has often been interpreted as evidence for a decrease of the covering factor of the AGN torus with increasing luminosity, the so-called receding torus models. Using a complete flux-limited X-ray selected sample of 199 AGN, from the Bright Ultra-hard XMM-Newton Survey, we determine the intrinsic fraction of optical type-2 AGN at 0.05$\leq$z$\leq$1 as a function of rest-frame 2-10 keV X-ray luminosity from 10$^{42}$ to 10$^{45}$ erg/s. We use the distributions of covering factors of AGN tori derived from CLUMPY torus models. Since these distributions combined over the total AGN population need to match the intrinsic type-2 AGN fraction, we reveal a population of X-ray undetected objects with high-covering factor tori, which are increasingly numerous at higher AGN luminosities. When these "missing" objects are included, we find that Compton-thick AGN account at most for 37$_{-10}^{+9}$% of the total population. The intrinsic type-2 AGN fraction is 58$\pm$4% and has a weak, non-significant (less than 2$\sigma$) luminosity dependence. This contradicts the results generally reported by AGN surveys, and the expectations from receding torus models. Our findings imply that the majority of luminous rapidly-accreting supermassive black holes at z<1 reside in highly-obscured nuclear environments but most of them are so deeply embedded that they have so far escaped detection in X-rays in <10 keV wide-area surveys.

Several authors have claimed that the less luminous active galactic nuclei (AGN) are not capable of sustaining the dusty torus structure. Thus, a gradual re-sizing of the torus is expected when the AGN luminosity decreases. Our aim is to confront mid-infrared observations of local AGN of different luminosities with this scenario. We decomposed about ~100 IRS/Spitzer spectra of LLAGN and powerful Seyferts in order to decontaminate the torus component from other contributors. We have used the affinity propagation (AP) method to cluster the data into five groups within the sample according to torus contribution to the 5-15 um range (Ctorus) and bolometric luminosity. The AP groups show a progressively higher torus contribution and an increase of the bolometric luminosity, from Group 1 (Ctorus~ 0% and logLbol ~ 41) and up to Group 5 (Ctorus ~80% and log(Lbol) ~44). We have fitted the average spectra of each of the AP groups to clumpy models. The torus is no longer present in Group 1, supporting the disappearance at low-luminosities. We were able to fit the average spectra for the torus component in Groups 3 (Ctorus~ 40% and log(Lbol)~ 42.6), 4 (Ctorus~ 60% and log(Lbol)~ 43.7), and 5 to Clumpy torus models. We did not find a good fitting to Clumpy torus models for Group 2 (Ctorus~ 18% and log(Lbol)~ 42). This might suggest a different configuration and/or composition of the clouds for Group 2, which is consistent with a different gas content seen in Groups 1, 2, and 3, according to the detections of H2 molecular lines. Groups 3, 4, and 5 show a trend to decrease of the width of the torus (which yields to a likely decrease of the geometrical covering factor), although we cannot confirm it with the present data. Finally, Groups 3, 4, and 5 show an increase on the outer radius of the torus for higher luminosities, consistent with a re-sizing of the torus according to the AGN luminosity.

We present subarcsecond resolution mid-infrared (MIR) images obtained with 8-10 m-class ground-based telescopes of a complete volume-limited (DL<40 Mpc) sample of 24 Seyfert galaxies selected from the Swift/BAT nine month catalog. We use those MIR images to study the nuclear and circumnuclear emission of the galaxies. Using different methods to classify the MIR morphologies on scales of ~400 pc, we find that the majority of the galaxies (75-83%) are extended or possibly extended and 17-25% are point-like. This extended emission is compact and it has low surface brightness compared with the nuclear emission, and it represents, on average, ~30% of the total MIR emission of the galaxies in the sample. We find that the galaxies whose circumnuclear MIR emission is dominated by star formation show more extended emission (650+-700 pc) than AGN-dominated systems (300+-100 pc). In general, the galaxies with point-like MIR morphologies are face-on or moderately inclined (b/a~0.4-1.0), and we do not find significant differences between the morphologies of Sy1 and Sy2. We used the nuclear and circumnuclear fluxes to investigate their correlation with different AGN and SF activity indicators. We find that the nuclear MIR emission (the inner ~70 pc) is strongly correlated with the X-ray emission (the harder the X-rays the better the correlation) and with the [O IV] lambda 25.89 micron emission line, indicating that it is AGN-dominated. We find the same results, although with more scatter, for the circumnuclear emission, which indicates that the AGN dominates the MIR emission in the inner ~400 pc of the galaxies, with some contribution from star formation.

We analyse a sample of 85 luminous (log(nuLnu(3\um)/erg s-1)>45.5) quasars with restframe ~2-11 \um spectroscopy from AKARI and Spitzer. Their high luminosity allows a direct determination of the near-infrared quasar spectrum free from host galaxy emission. A semi-empirical model consisting of a single template for the accretion disk and two blackbodies for the dust emission successfully reproduces the 0.1-10 \um spectral energy distributions (SEDs). Excess emission at 1-2 \um over the best-fitting model suggests that hotter dust is necessary in addition to the ~1200 K blackbody and the disk to reproduce the entire near-infrared spectrum. Variation in the extinction affecting the disk and in the relative strength of the disk and dust components accounts for the diversity of individual SEDs. Quasars with higher dust-to-disk luminosity ratios show slightly redder infrared continua and less prominent silicate emission. We find no luminosity dependence in the shape of the average infrared quasar spectrum. We generate a new quasar template that covers the restframe range 0.1-11 \um, and separate templates for the disk and dust components. Comparison with other infrared quasar composites suggests that previous ones are less reliable in the 2-4 \um range. Our template is the first one to provide a detailed view of the infrared emission on both sides of the 4 \um bump.

We present far-infrared (FIR) $70-500\,\mu$m imaging observations obtained with Herschel/PACS and SPIRE of 33 nearby (median distance of 30 Mpc) Seyfert galaxies from the Revised Shapley-Ames (RSA) catalogue. We obtain the FIR nuclear ($r=1\,$kpc and $r=2\,$kpc) and integrated spectral energy distributions (SEDs). We estimate the unresolved nuclear emission at 70 $\mu$m and we fit the nuclear and integrated FIR SEDs with a grey body model. We find that the integrated FIR emission of the RSA Seyferts in our sample is dominated by emission from the host galaxy, with dust properties similar to those of normal galaxies (non AGN). We use four criteria to select galaxies whose nuclear $70\,\mu$m emission has a significant AGN contribution: (1) elevated 70/160 $\mu$m flux ratios, (2)spatially resolved, high dust temperature gradient, (3) $70\,\mu$m excess emission with respect to the fit of the FIR SEDs with a grey body, and (4) excess of nuclear SFR obtained from $70\,\mu$m over SFR from mid-infrared indicators. 16 galaxies (48 per cent of the initial sample) satisfy at least one of these conditions, whereas 10 satisfy half or more. After careful examination of these, we select six bona fide candidates (18 per cent of the initial sample) and estimate that $\sim 40-70$ per cent of their nuclear ($r=1-2\,$kpc) $70\,\mu$m emission is contributed by dust heated by the AGN.

We present the distributions of geometrical covering factors of active galactic nuclei (AGNs) dusty tori (f2) using an X-ray selected complete sample of 227 AGN drawn from the Bright Ultra-hard XMM-Newton Survey. The AGN have z from 0.05 to 1.7, 2-10 keV luminosities between 10^42 and 10^46 erg/s and Compton-thin X-ray absorption. Employing data from UKIDSS, 2MASS and the Wide-field Infrared Survey Explorer in a previous work we determined the rest-frame 1-20 microns continuum emission from the torus which we model here with the clumpy torus models of Nenkova et al. Optically classified type 1 and type 2 AGN are intrinsically different, with type 2 AGN having on average tori with higher f2 than type 1 AGN. Nevertheless, ~20 per cent of type 1 AGN have tori with large covering factors while ~23-28 per cent of type 2 AGN have tori with small covering factors. Low f2 are preferred at high AGN luminosities, as postulated by simple receding torus models, although for type 2 AGN the effect is certainly small. f2 increases with the X-ray column density, which implies that dust extinction an X-ray absorption takes place in material that shares an overall geometry and most likely belongs to the same structure, the putative torus. Based on our results, the viewing angle, AGN luminosity and also f2 determine the optical appearance of an AGN and control the shape of the rest-frame ~1-20 microns nuclear continuum emission. Thus, the torus geometrical covering factor is a key ingredient of unification schemes.

This is the first of a series of papers on the Infrared Database of Extragalactic Observables from Spitzer (IDEOS). In this work we describe the identification of optical counterparts of the infrared sources detected in Spitzer Infrared Spectrograph (IRS) observations, and the acquisition and validation of redshifts. The IDEOS sample includes all the spectra from the Cornell Atlas of Spitzer/IRS Sources (CASSIS) of galaxies beyond the Local Group. Optical counterparts were identified from correlation of the extraction coordinates with the NASA Extragalactic Database (NED). To confirm the optical association and validate NED redshifts, we measure redshifts with unprecedented accuracy on the IRS spectra ({\sigma}(dz/(1+z))=0.0011) by using an improved version of the maximum combined pseudo-likelihood method (MCPL). We perform a multi-stage verification of redshifts that considers alternate NED redshifts, the MCPL redshift, and visual inspection of the IRS spectrum. The statistics is as follows: the IDEOS sample contains 3361 galaxies at redshift 0<z<6.42 (mean: 0.48, median: 0.14). We confirm the default NED redshift for 2429 sources and identify 124 with incorrect NED redshifts. We obtain IRS-based redshifts for 568 IDEOS sources without optical spectroscopic redshifts, including 228 with no previous redshift measurements. We provide the entire IDEOS redshift catalog in machine-readable formats. The catalog condenses our compilation and verification effort, and includes our final evaluation on the most likely redshift for each source, its origin, and reliability estimates.

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 an analysis of the nuclear infrared (IR, 1.6 to 18 $\mu$m) emission of the ultraluminous IR galaxy UGC 5101 to derive the properties of its active galactic nucleus (AGN) and its obscuring material. We use new mid-IR high angular resolution ($0.3-0.5$ arcsec) imaging using the Si-2 filter ($\lambda_{C}=8.7\, \mu$m) and $7.5-13$ $\mu$m spectroscopy taken with CanariCam (CC) on the 10.4m Gran Telescopio CANARIAS. We also use archival HST/NICMOS and Subaru/COMICS imaging and Spitzer/IRS spectroscopy. We estimate the near- and mid-IR unresolved nuclear emission by modelling the imaging data with GALFIT. We decompose the Spitzer/IRS and CC spectra using a power-law component, which represents the emission due to dust heated by the AGN, and a starburst component, both affected by foreground extinction. We model the resulting unresolved near- and mid-IR, and the starburst subtracted CC spectrum with the CLUMPY torus models of Nenkova et al. The derived geometrical properties of the torus, including the large covering factor and the high foreground extinction needed to reproduce the deep $9.7\, \mu$m silicate feature, are consistent with the lack of strong AGN signatures in the optical. We derive an AGN bolometric luminosity $L_{bol}\sim1.9\times10^{45}\,$erg s$^{-1}$ that is in good agreement with other estimates in the literature.

We present a comprehensive study of the silicate features at 9.7 and 18 micron of a sample of almost 800 active galactic nuclei (AGN) with available spectra from the Spitzer InfraRed Spectrograph (IRS). We measure the strength of the silicate feature at 9.7 micron, S9.7, before and after subtracting the host galaxy emission from the IRS spectra. The numbers of type 1 and 2 AGN with the feature in emission increase by 20 and 50%, respectively, once the host galaxy is removed, while 35% of objects with the feature originally in absorption exhibit it in even deeper absorption. The peak of S9.7, lambda_peak, has a bimodal distribution when the feature is in emission, with about 65% of the cases showing lambda_peak > 10.2 micron. Silicates can appear in emission in objects with mid-infrared (MIR) luminosity spanning over six orders of magnitude. The derived distributions of the strength of the silicate features at 9.7 and 18 micron provide a solid test bed for modeling the dust distribution in AGN. Clumpiness is needed in order to produce absorption features in unobscured AGN and can also cause the silicates to be in absorption at 9.7 micron and in emission at 18 micron in type 1 sources. We find the `cosmic' silicates of Ossenkopf et al. to be more consistent with the observations than Draine's `astronomical' silicates. Finally, we discuss the possibility of a foreground absorber to explain the deep silicate absorption features in the MIR spectra of some type 2 AGN.

We present results on the spectral decomposition of 118 Spitzer Infrared Spectrograph (IRS) spectra from local active galactic nuclei (AGN) using a large set of Spitzer/IRS spectra as templates. The templates are themselves IRS spectra from extreme cases where a single physical component (stellar, interstellar, or AGN) completely dominates the integrated mid-infrared emission. We show that a linear combination of one template for each physical component reproduces the observed IRS spectra of AGN hosts with unprecedented fidelity for a template fitting method, with no need to model extinction separately. We use full probability distribution functions to estimate expectation values and uncertainties for observables, and find that the decomposition results are robust against degeneracies. Furthermore, we compare the AGN spectra derived from the spectral decomposition with sub-arcsecond resolution nuclear photometry and spectroscopy from ground-based observations. We find that the AGN component derived from the decomposition closely matches the nuclear spectrum, with a 1-sigma dispersion of 0.12 dex in luminosity and typical uncertainties of ~0.19 in the spectral index and ~0.1 in the silicate strength. We conclude that the emission from the host galaxy can be reliably removed from the IRS spectra of AGN. This allows for unbiased studies of the AGN emission in intermediate and high redshift galaxies -currently inaccesible to ground-based observations- with archival Spitzer/IRS data and in the future with the Mid-InfraRed Instrument of the James Webb Space Telescope. The decomposition code and templates are available at http://www.denebola.org/ahc/deblendIRS.

We have determined the relation between the AGN luminosities at rest-frame 6 {\mu}m associated to the dusty torus emission and at 2-10 keV energies using a complete, X-ray flux limited sample of 232 AGN drawn from the Bright Ultra-hard XMM-Newton Survey. The objects have intrinsic X-ray luminosities between 10^42 and 10^46 erg/s and redshifts from 0.05 to 2.8. The rest-frame 6 {\mu}m luminosities were computed using data from the Wide-Field Infrared Survey Explorer and are based on a spectral energy distribution decomposition into AGN and galaxy emission. The best-fit relationship for the full sample is consistent with being linear, L_6 {\mu}m $\propto$ L_2-10 keV^0.99$\pm$0.032, with intrinsic scatter, ~0.35 dex in log L_6 {\mu}m. The L_6 {\mu}m/L_2-10 keV luminosity ratio is largely independent on the line-of-sight X-ray absorption. Assuming a constant X-ray bolometric correction, the fraction of AGN bolometric luminosity reprocessed in the mid-IR decreases weakly, if at all, with the AGN luminosity, a finding at odds with simple receding torus models. Type 2 AGN have redder mid-IR continua at rest-frame wavelengths <12 {\mu}m and are overall ~1.3-2 times fainter at 6 {\mu}m than type 1 AGN at a given X-ray luminosity. Regardless of whether type 1 and type 2 AGN have the same or different nuclear dusty toroidal structures, our results imply that the AGN emission at rest-frame 6 {\mu}m is not isotropic due to self-absorption in the dusty torus, as predicted by AGN torus models. Thus, AGN surveys at rest-frame 6 {\mu}m are subject to modest dust obscuration biases.

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

In this work we explore the impact of the presence of an active galactic nucleus (AGN) on the mid- and far-infrared (IR) properties of galaxies as well as the effects of simultaneous AGN and starburst activity in these same galaxies. To do this we apply a multi-component, multi-band spectral synthesis technique to a sample of 250 micron selected galaxies of the Herschel Multi-tiered Extragalactic Survey, with IRS spectra available for all galaxies. Our results confirm that the inclusion of the IRS spectra plays a crucial role in the spectral analysis of galaxies with an AGN component improving the selection of the best-fit hot dust model (torus). We find a correlation between the obscured star formation rate (SFR) derived from the IR luminosity of the starburst component, SFR_IR and SFR_PAH, derived from the luminosity of the PAH features, L_PAH, with SFR_FIR taking higher values than SFR_PAH. The correlation is different for AGN- and starburst-dominated objects. The ratio of L_PAH to that of the starburst component, L_PAH/L_SB, is almost constant for AGN-dominated objects but decreases with increasing L_SB for starburst-dominated objects. SFR_FIR increases with the accretion luminosity, L_acc, with the increase less prominent for the very brightest, unobscured AGN-dominated sources. We find no correlation between the masses of the hot and cold dust components. We interpret this as a non-constant fraction of gas driven by the gravitational effects to the AGN while the starburst is ongoing. We also find no evidence of the AGN affecting the temperature of the cold dust component, though this conclusion is mostly based on objects with a non-dominant AGN component. We conclude that our findings do not provide evidence that the presence of AGN affects the star formation process in the host galaxy, but rather that the two phenomena occur simultaneously over a wide range of luminosities.

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.

We present a new SED-fitting based routine for redshift determination that is optimised for mid-infrared (MIR) low-resolution spectroscopy. Its flexible template scaling increases the sensitivity to slope changes and small scale features in the spectrum, while a new selection algorithm called Maximum Combined Pseudo-Likelihood (MCPL) provides increased accuracy and a lower number of outliers compared to the standard maximum-likelihood (ML) approach. Unlike ML, MCPL searches for local (instead of absolute) maxima of a 'pseudo-likelihood' (PL) function, and combines results obtained for all the templates in the library to weed out spurious redshift solutions. The capabilities of MCPL are demonstrated by comparing its results to those of regular ML and to the optical spectroscopic redshifts of a sample of 491 Spitzer/IRS spectra from sources at 0<z<3.7. MCPL achieves a redshift accuracy dz/(1+z)<0.005 for 78% of the galaxies in the sample compared to 68% for ML. The rate of outliers (dz/(1+z)>0.02) is 14% for MCPL and 22% for ML. chi^2 values for ML solutions are found to correlate with the SNR of the spectra, but not with redshift accuracy. By contrast, the peak value of the normalised combined PL (gamma) is found to provide a good indication on the reliability of the MCPL solution for individual sources. The accuracy and reliability of the redshifts depends strongly on the MIR SED. Sources with significant polycyclic aromatic hydrocarbon emission obtain much better results compared to sources dominated by AGN continuum. Nevertheless, for a given gamma the frequency of accurate solutions and outliers is largely independent on their SED type. This reliability indicator for MCPL solutions allows to select subsamples with highly reliable redshifts. In particular, a gamma>0.15 threshold retains 79% of the sources with dz/(1+z)<0.005 while reducing the outlier rate to 3.8% (abridged).

We present a panoramic atlas of Spitzer/IRS spectra of extragalactic sources collected from the recent literature, with value added measurements of their spectral features obtained in a homogeneous and concise manner. The atlas covers the full spectrum of the extragalactic universe and includes star forming galaxies, obscured and unobscured active galaxies, luminous and ultra-luminous infrared galaxies, and hybrid objects. Measured features such as the polycyclic aromatic hydrocarbons, the strength of the silicates in emission or absorption around 9.7 micron, rest-frame monochromatic luminosities or colours, combined with measurements derived from spectral decomposition, are used to establish diagnostics that allow for classification of sources, based on their infrared properties alone. Average templates of the various classes are also derived. The full atlas with the value added measurements and ancillary archival data are publicly available at http://www.denebola.org/atlas, with full references to the original data.