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.

From the approximately $\sim$3,500 planetary nebulae (PNe) discovered in our Galaxy, only 14 are known to be members of the Galactic halo. Nevertheless, a systematic search for halo PNe has never been performed. In this study, we present new photometric diagnostic tools to identify compact PNe in the Galactic halo by making use of the novel 12-filter system projects, J-PLUS (Javalambre Photometric Local Universe Survey) and S-PLUS (Southern-Photometric Local Universe Survey). We reconstructed the IPHAS (Isaac Newton Telescope (INT) Photometric H${\alpha}$ Survey of the Northern Galactic Plane) diagnostic diagram and propose four new ones using i) the J-PLUS and S-PLUS synthetic photometry for a grid of photo-ionisation models of halo PNe, ii) several observed halo PNe, as well as iii) a number of other emission-line objects that resemble PNe. All colour-colour diagnostic diagrams are validated using two known halo PNe observed by J-PLUS during the scientific verification phase and the first data release (DR1) of S-PLUS and the DR1 of J-PLUS. By applying our criteria to the DR1s ($\sim$1,190 deg$^2$), we identified one PN candidate. However, optical follow-up spectroscopy proved it to be a H II region belonging to the UGC 5272 galaxy. Here, we also discuss the PN and two H II galaxies recovered by these selection criteria. Finally, the cross-matching with the most updated PNe catalogue (HASH) helped us to highlight the potential of these surveys, since we recover all the known PNe in the observed area. The tools here proposed to identify PNe and separate them from their emission-line contaminants proved to be very efficient thanks to the combination of many colours, even when applied -like in the present work- to an automatic photometric search that is limited to compact PNe.

We present a synthetic galaxy lightcone specially designed for narrow-band optical photometric surveys. To reduce time-discreteness effects, unlike previous works, we directly include the lightcone construction in the \texttt{L-Galaxies} semi-analytic model applied to the subhalo merger trees of the {\tt Millennium} simulation. Additionally, we add a model for the nebular emission in star-forming regions, which is crucial for correctly predicting the narrow/medium-band photometry of galaxies. Explicitly, we consider, individually for each galaxy, the contribution of 9 different lines: $\rm Ly{\alpha}$ (1216\AA), \Hb (4861\AA), \Ha (6563\AA), {\oii} (3727\AA, 3729\AA), {\oiii} (4959\AA, 5007\AA), $\rm [\ion{Ne}{III}]$ (3870\AA), {\oi} (6300\AA), $\rm [\ion{N}{II}]$ (6548\AA, 6583\AA), and $\rm [\ion{S}{II}]$ (6717\AA, 6731\AA). We validate our lightcone by comparing galaxy number counts, angular clustering, and \Ha, \Hb, {\oii} and {\oiiiFd} luminosity functions to a compilation of observations. As an application of our mock lightcones, we generate catalogues tailored for J-PLUS, a large optical galaxy survey featuring 5 broad and 7 medium band filters. We study the ability of the survey to correctly identify, with a simple \textit{three filter method}, a population of emission-line galaxies at various redshifts. We show that the $4000\AA$ break in the spectral energy distribution of galaxies can be misidentified as line emission. However, all significant excess (larger than 0.4 magnitudes) can be correctly and unambiguously attributed to emission line galaxies. Our catalogues are publicly released to facilitate their use in interpreting narrow-band surveys and for quantifying the impact of line emission in broad band photometry.

Our aim is to determine the distribution of stellar population parameters (extinction, age, metallicity, and star formation rate) of quiescent galaxies within the rest-frame stellar mass$-$colour and $UVJ$ colour$-$colour diagrams corrected for extinction up to $z\sim1$. These novel diagrams reduce the contamination in samples of quiescent galaxies owing to dust-reddened galaxies, and they provide useful constraints on stellar population parameters. We set constraints on the stellar population parameters of quiescent galaxies combining the ALHAMBRA multi-filter photo-spectra with our SED-fitting code MUFFIT, making use of composite stellar population models. The extinction obtained by MUFFIT allowed us to remove dusty star-forming (DSF) galaxies from the sample of red $UVJ$ galaxies. The distributions of stellar population parameters across these rest-frame diagrams are revealed after the dust correction and are fitted by the LOESS method to reduce uncertainty effects. Quiescent galaxy samples defined via classical $UVJ$ diagrams are typically contaminated by a $\sim20$% fraction of DSF galaxies. A significant part of the galaxies in the green valley are actually obscured star-forming galaxies ($\sim30-65$%). Consequently, the transition of galaxies from the blue cloud to the red sequence, and hence the related mechanisms for quenching, seems to be much more efficient and faster than previously reported. The rest-frame stellar mass$-$colour and $UVJ$ colour$-$colour diagrams are useful for constraining the age, metallicity, extinction, and star formation rate of quiescent galaxies by only their redshift, rest-frame colours, and/or stellar mass. Dust correction plays an important role in understanding how quiescent galaxies are distributed in these diagrams and is key to performing a pure selection of quiescent galaxies via intrinsic colours.

The next generation of galaxy surveys will allow us to test one of the most fundamental assumptions of the standard cosmology, i.e., that gravity is governed by the general theory of relativity (GR). In this paper we investigate the ability of the Javalambre Physics of the Accelerating Universe Astrophysical Survey (J-PAS) to constrain GR and its extensions. Based on the J-PAS information on clustering and gravitational lensing, we perform a Fisher matrix forecast on the effective Newton constant, $\mu$, and the gravitational slip parameter, $\eta$, whose deviations from unity would indicate a breakdown of GR. Similar analysis is also performed for the DESI and Euclid surveys and compared to J-PAS with two configurations providing different areas, namely an initial expectation with 4000 $\mathrm{deg}^2$ and the future best case scenario with 8500 $\mathrm{deg}^2$. We show that J-PAS will be able to measure the parameters $\mu$ and $\eta$ at a sensitivity of $2\% - 7\%$, and will provide the best constraints in the interval $z = 0.3 - 0.6$, thanks to the large number of ELGs detectable in that redshift range. We also discuss the constraining power of J-PAS for dark energy models with a time-dependent equation-of-state parameter of the type $w(a)=w_0+w_a(1-a)$, obtaining $\Delta w_0=0.058$ and $\Delta w_a=0.24$ for the absolute errors of the dark energy parameters.

We present the photometric calibration of the twelve optical passbands observed by the Javalambre Photometric Local Universe Survey (J-PLUS). The proposed calibration method has four steps: (i) definition of a high-quality set of calibration stars using Gaia information and available 3D dust maps; (ii) anchoring of the J-PLUS gri passbands to the Pan-STARRS photometric solution, accounting for the variation of the calibration with the position of the sources on the CCD; (iii) homogenization of the photometry in the other nine J-PLUS filters using the dust de-reddened instrumental stellar locus in (X - r) versus (g - i) colours, where X is the filter to calibrate. The zero point variation along the CCD in these filters was estimated with the distance to the stellar locus. Finally, (iv) the absolute colour calibration was obtained with the white dwarf locus. We performed a joint Bayesian modelling of eleven J-PLUS colour-colour diagrams using the theoretical white dwarf locus as reference. This provides the needed offsets to transform instrumental magnitudes to calibrated magnitudes outside the atmosphere. The uncertainty of the J-PLUS photometric calibration, estimated from duplicated objects observed in adjacent pointings and accounting for the absolute colour and flux calibration errors, are ~19 mmag in u, J0378 and J0395, ~11 mmag in J0410 and J0430, and ~8 mmag in g, J0515, r, J0660, i, J0861, and z. We present an optimized calibration method for the large area multi-filter J-PLUS project, reaching 1-2% accuracy within an area of 1 022 square degrees without the need for long observing calibration campaigns or constant atmospheric monitoring. The proposed method will be adapted for the photometric calibration of J-PAS, that will observe several thousand square degrees with 56 narrow optical filters.

In a framework where galaxies mostly migrate on the colour-magnitude diagram from star-forming to quiescent, the green valley is considered a transitional galaxy stage. The details of the processes that drive galaxies from star-forming to passive systems still remain unknown. We developed a method that estimates empirically the star formation quenching times-scales of green valley galaxies, assuming an exponential decay model of the SFH and through a combination of narrow and broad bands from J-PLUS and GALEX. We correlate these quenching time-scales with the presence of bars. We find that the J-PLUS colours F0395-g and F0415-g are sensitive to different SFH, showing, a clear correlation with the Dn(4000) and H-delta,A spectral indices. We find that quenching time-scales obtained with our new approach are in agreement with those determined using spectral indices. We also find that galaxies with high bar probability tend to quench their star formation slowly. We conclude that: 1) J-PLUS filters can be used to measure quenching timescales in nearby green valley galaxies; and 2) the resulting star formation quenching time-scales are longer for barred green valley galaxies. Considering that the presence of a bar indicates that more violent processes (e.g., major mergers) are absent in host galaxies, we conclude that the presence of a bar can be used as a morphological signature for slow star formation quenching.

We aim at constraining the stellar population properties of quiescent galaxies. These properties reveal how these galaxies evolved and assembled since $z\sim1$ up to the present time. Combining the ALHAMBRA multi-filter photo-spectra with the SED-fitting code MUFFIT, we build a complete catalogue of quiescent galaxies via the dust-corrected stellar mass vs colour diagram. This catalogue includes stellar population properties, such as age, metallicity, extinction, stellar mass and photometric redshift, retrieved from the analysis of composited populations based on two independent sets of SSP models. We develop and apply a novel methodology to provide, for the first time, the analytic probability distribution functions (PDFs) of mass-weighted age, metallicity, and extinction of quiescent galaxies as a function of redshift and stellar mass. We adopt different star formation histories to discard potential systematics in the analysis. The number density of quiescent galaxies is found to increase since $z\sim1$, with a more substantial variation at lower mass. Quiescent galaxies feature extinction $A_V<0.6$, with median values in the range $A_V = 0.15\mathrm{-}0.3$. At increasing stellar mass, quiescent galaxies are older and more metal rich since $z\sim1$. A detailed analysis of the PDFs reveals that the evolution of quiescent galaxies is not compatible with passive evolution and a slight decrease is hinted at median metallicity $0.1\mathrm{-}0.2$~dex. The intrinsic dispersion of the age and metallicity PDFs show a dependence with stellar mass and/or redshift. These results are consistent with both sets of SSP models and the alternative SFH assumptions explored. Consequently, the quiescent population must undergo an evolutive pathway including mergers and/or remnants of star formation to reconcile the observed trends, where the `progenitor' bias should also be taken into account.

Ultracool dwarfs (UCDs) are objects with spectral types equal or later than M7. Most of them have been discovered using wide-field imaging surveys. The Virtual Observatory (VO) has proven to be of great utility to efficiently exploit these astronomical resources. We aim to validate a VO methodology designed to discover and characterize UCDs in the J-PLUS photometric survey. J-PLUS is a multiband survey carried out with the wide angle T80Cam optical camera mounted on the 0.83-m telescope JAST/T80 in the Observatorio Astrof\'isico de Javalambre. In this work we make use of the Internal Data Release (IDR) covering 528 deg$^2$. We complement J-PLUS photometry with other catalogues in the optical and IR using VOSA, a VO tool that estimates physical parameters from the spectral energy distribution fitting to collections of theoretical models. Objects identified as UCDs are distinguished from background M giants and highly reddened stars using parallaxes and proper motions from Gaia DR2. We identify 559 UCDs, ranging from i=16.2 to 22.4 mag, of which 187 are candidate UCDs not previously reported in the literature. This represents an increase in the number of known UCDs of about 50% in the studied region of the sky, particularly at the faint end of our sensitivity, which is interesting as reference for future wide and deep surveys such as Euclid. Three candidates constitute interesting targets for exoplanet surveys because of their proximity (<40 pc). We also analyze the kinematics of UCDs in our catalogue and find evidence that it is consistent with a Galactic thin-disk population, except for 6 objects that might be members of the thick disk. The results obtained validate the proposed methodology, which will be used in future J-PLUS and J-PAS releases. Considering the region of the sky covered by the IDR used, we foresee to discover 3,000-3,500 new UCDs at the end of the J-PLUS project.

This paper describes the characterization of the GALANTE photometric system, a seven intermediate- and narrow-band filter system with a wavelength coverage from 3000 $\r{A}$ to 9000 $\r{A}$ . We describe the photometric system presenting the full sensitivity curve as a product of the filter sensitivity, CCD, telescope mirror, and atmospheric transmission curves, as well as some first- and second-order moments of this sensitivity function. The GALANTE photometric system is composed of four filters from the J-PLUS photometric system, a twelve broad-to-narrow filter system, and three exclusive filters, specifically designed to measure the physical parameters of stars such as effective temperature $T_{\rm eff}$, $\log(g)$, metallicity, colour excess $E(4405-5495)$, and extinction type $R_{5495}$. Two libraries, the Next Generation Spectral Library (NGSL) and the one presented in Ma\'iz Apell\'aniz & Weiler (2018), have been used to determine the transformation equations between the Sloan Digital Sky Survey ($\textit{SDSS}$) $\textit{ugriz}$ photometry and the GALANTE photometric system. We will use this transformation to calibrate the zero points of GALANTE images. To this end, a preliminary photometric calibration of GALANTE has been made based on two different $\textit{griz}$ libraries ($\textit{SDSS}$ DR12 and ATLAS All-Sky Stellar Reference Catalog, hereinafter $\textit{RefCat2}$). A comparison between both zero points is performed leading us to the choice of $\textit{RefCat2}$ as the base catalogue for this calibration, and applied to a field in the Cyg OB2 association.

We present a new methodology for the estimation of stellar atmospheric parameters from narrow- and intermediate-band photometry of the Javalambre Photometric Local Universe Survey (J-PLUS), and propose a method for target pre-selection of low-metallicity stars for follow-up spectroscopic studies. Photometric metallicity estimates for stars in the globular cluster M15 are determined using this method. By development of a neural-network-based photometry pipeline, we aim to produce estimates of effective temperature, $T_{\rm eff}$, and metallicity, [Fe/H], for a large subset of stars in the J-PLUS footprint. The Stellar Photometric Index Network Explorer, SPHINX, is developed to produce estimates of $T_{\rm eff}$ and [Fe/H], after training on a combination of J-PLUS photometric inputs and synthetic magnitudes computed for medium-resolution (R ~ 2000) spectra of the Sloan Digital Sky Survey. This methodology is applied to J-PLUS photometry of the globular cluster M15. Effective temperature estimates made with J-PLUS Early Data Release photometry exhibit low scatter, \sigma($T_{\rm eff}$) = 91 K, over the temperature range 4500 < $T_{\rm eff}$ (K) < 8500. For stars from the J-PLUS First Data Release with 4500 < $T_{\rm eff}$ (K) < 6200, 85 $\pm$ 3% of stars known to have [Fe/H] <-2.0 are recovered by SPHINX. A mean metallicity of [Fe/H]=-2.32 $\pm$ 0.01, with a residual spread of 0.3 dex, is determined for M15 using J-PLUS photometry of 664 likely cluster members. We confirm the performance of SPHINX within the ranges specified, and verify its utility as a stand-alone tool for photometric estimation of effective temperature and metallicity, and for pre-selection of metal-poor spectroscopic targets.

The intracluster light (ICL) is a luminous component of galaxy clusters composed of stars that are gravitationally bound to the cluster potential but do not belong to the individual galaxies. Previous studies of the ICL have shown that its formation and evolution are intimately linked to the evolutionary stage of the cluster. Thus, the analysis of the ICL in the Coma cluster will give insights into the main processes driving the dynamics in this highly complex system. Using a recently developed technique, we measure the ICL fraction in Coma at several wavelengths, using the J-PLUS unique filter system. The combination of narrow- and broadband filters provides valuable information on the dynamical state of the cluster, the ICL stellar types, and the morphology of the diffuse light. We use the Chebyshev-Fourier Intracluster Light Estimator (CICLE) to disentangle the ICL from the light of the galaxies, and to robustly measure the ICL fraction in seven J-PLUS filters. We obtain the ICL fraction distribution of the Coma cluster at different optical wavelengths, which varies from $\sim 7\%-21\%$, showing the highest values in the narrowband filters J0395, J0410, and J0430. This ICL fraction excess is distinctive pattern recently observed in dynamically active clusters (mergers), indicating a higher amount of bluer stars in the ICL compared to the cluster galaxies. Both the high ICL fractions and the excess in the bluer filters are indicative of a merging state. The presence of younger/lower-metallicity stars the ICL suggests that the main mechanism of ICL formation for the Coma cluster is the stripping of the stars in the outskirts of infalling galaxies and, possibly, the disruption of dwarf galaxies during past/ongoing mergers.

Our goal is to characterise the dependence of the optical mass-to-light ratio on galaxy colour up to z = 1.5, expanding the redshift range explored in previous work. From the ALHAMBRA redshifts, stellar masses, and rest-frame luminosities provided by the MUFFIT code, we derive the mass-to-light ratio vs. colour relation (MLCR) both for quiescent and star-forming galaxies. The intrinsic relation and its physical dispersion are derived with a Bayesian inference model. The rest-frame i-band mass-to-light ratio of quiescent and star-forming galaxies presents a tight correlation with the rest-frame (g - i) colour up to z = 1.5. Such MLCR is linear for quiescent galaxies and quadratic for star-forming galaxies. The intrinsic dispersion in these relations is 0.02 dex for quiescent galaxies and 0.06 dex for star-forming ones. The derived MLCRs do not present a significant redshift evolution and are compatible with previous local results in the literature. Finally, these tight relations also hold for g- and r-band luminosities. The derived MLCRs in ALHAMBRA can be used to predict the mass-to-light ratio from a rest-frame optical colour up to z = 1.5. These tight correlations do not change with redshift, suggesting that galaxies have evolved along the derived relations during the last 9 Gyr.

The Javalambre Photometric Local Universe Survey (J-PLUS) provides wide field-of-view images in 12 narrow, intermediate and broad-band filters optimized for stellar photometry. Here we have applied J-PLUS data for the first time for the study of Galactic GCs using science verification data obtained for the very metal-poor GC M\,15. Our J-PLUS data provide low-resolution spectral energy distributions covering the near-UV to the near-IR, allowing us to search for MPs based on pseudo-spectral fitting diagnostics. J-PLUS CMDs are found to be particularly useful to search for splits in the sequences formed by the upper red giant branch (RGB) and asymptotic giant branch (AGB) stars. We interpret these split sequences as evidence for the presence of MPs. This demonstrates that the J-PLUS survey will have sufficient spatial coverage and spectral resolution to perform a large statistical study of GCs through multi-band photometry in the coming years.

In the present paper we aim to validate a methodology designed to extract the Halpha emission line flux from J-PLUS photometric data. J-PLUS is a multi narrow-band filter survey carried out with the 2 deg2 field of view T80Cam camera, mounted on the JAST/T80 telescope in the OAJ, Teruel, Spain. The information of the twelve J-PLUS bands, including the J0660 narrow-band filter located at rest-frame Halpha, is used over 42 deg2 to extract de-reddened and [NII] decontaminated Halpha emission line fluxes of 46 star-forming regions with previous SDSS and/or CALIFA spectroscopic information. The agreement of the inferred J-PLUS photometric Halpha fluxes and those obtained with spectroscopic data is remarkable, with a median comparison ratio R = 1.05 +- 0.25. This demonstrates that it is possible to retrieve reliable Halpha emission line fluxes from J-PLUS photometric data. With an expected area of thousands of square degrees upon completion, the J-PLUS dataset will allow the study of several star formation science cases in the nearby universe, as the spatially resolved star formation rate of nearby galaxies at z < 0.015, and how it is influenced by the environment, morphology or nuclear activity. As an illustrative example, the close pair of interacting galaxies NGC3994 and NGC3995 is analyzed, finding an enhancement of the star formation rate not only in the center, but also in outer parts of the disk of NGC3994.

The spatial variations of stellar population properties within a galaxy are intimately related to their formation process. Therefore, spatially resolved studies of galaxies are essential to uncover their formation and assembly. The Javalambre Photometric Local Universe Survey (J-PLUS) is a dedicated multi-filter designed to observed ~8500 deg2 using twelve narrow-, intermediate- and broad-band filters in the optical range. In this study, we test the potential of the multi-filter observation carried out with J-PLUS to investigate the properties of spatially-resolved nearby galaxies. We present detailed 2D maps of stellar population properties (age, metallicity, extinction, and stellar mass surface density) for two early-type galaxies observed in both, J-PLUS and CALIFA surveys: NGC 5473 and NGC 5485. Radial structures are also compared and luminosity- and mass-weighted profiles are derived. We use MUFFIT to process the J-PLUS observations, and two different techniques (STARLIGHT and STECKMAP) to analyze IFU CALIFA data. We demonstrate that this novel technique delivers radial stellar population gradients in good agreement with the IFU technique CALIFA/STECKMAP although comparison of the absolute values reveals the existence of intrinsic systematic differences. Radial stellar population gradients differ when CALIFA/STARLIGHT methodology is used. Age and metallicity radial profiles derived from J-PLUS/MUFFIT are very similar when luminosity- or mass-weighted properties are used, suggesting that the contribution of a younger component is small. Comparison between the three methodologies reveals some discrepancies suggesting that the specific characteristics of each method causes important differences. We conclude that the ages, metallicities and extinction derived for individual galaxies not only depend on the chosen models but also depend on the method used.

J-PLUS is an ongoing 12-band photometric optical survey, observing thousands of square degrees of the Northern hemisphere from the dedicated JAST/T80 telescope at the Observatorio Astrof\'isico de Javalambre. T80Cam is a 2 sq.deg field-of-view camera mounted on this 83cm-diameter telescope, and is equipped with a unique system of filters spanning the entire optical range. This filter system is a combination of broad, medium and narrow-band filters, optimally designed to extract the rest-frame spectral features (the 3700-4000\AA\ Balmer break region, H$\delta$, Ca H+K, the G-band, the Mgb and Ca triplets) that are key to both characterize stellar types and to deliver a low-resolution photo-spectrum for each pixel of the sky observed. With a typical depth of AB $\sim 21.25$ mag per band, this filter set thus allows for an indiscriminate and accurate characterization of the stellar population in our Galaxy, it provides an unprecedented 2D photo-spectral information for all resolved galaxies in the local universe, as well as accurate photo-z estimates ($\Delta\,z\sim 0.01-0.03$) for moderately bright (up to $r\sim 20$ mag) extragalactic sources. While some narrow band filters are designed for the study of particular emission features ([OII]/$\lambda$3727, H$\alpha$/$\lambda$6563) up to $z < 0.015$, they also provide well-defined windows for the analysis of other emission lines at higher redshifts. As a result, J-PLUS has the potential to contribute to a wide range of fields in Astrophysics, both in the nearby universe (Milky Way, 2D IFU-like studies, stellar populations of nearby and moderate redshift galaxies, clusters of galaxies) and at high redshifts (ELGs at $z\approx 0.77, 2.2$ and $4.4$, QSOs, etc). With this paper, we release $\sim 36$ sq.deg of J-PLUS data, containing about $1.5\times 10^5$ stars and $10^5$ galaxies at $r<21$ mag.

Our goal is to morphologically classify the sources identified in the images of the J-PLUS early data release (EDR) into compact (stars) or extended (galaxies) using a suited Bayesian classifier. J-PLUS sources exhibit two distinct populations in the r-band magnitude vs. concentration plane, corresponding to compact and extended sources. We modelled the two-population distribution with a skewed Gaussian for compact objects and a log-normal function for the extended ones. The derived model and the number density prior based on J-PLUS EDR data were used to estimate the Bayesian probability of a source to be star or galaxy. This procedure was applied pointing-by-pointing to account for varying observing conditions and sky position. Finally, we combined the morphological information from g, r, and i broad bands in order to improve the classification of low signal-to-noise sources. The derived probabilities are used to compute the pointing-by-pointing number counts of stars and galaxies. The former increases as we approach to the Milky Way disk, and the latter are similar across the probed area. The comparison with SDSS in the common regions is satisfactory up to r ~ 21, with consistent numbers of stars and galaxies, and consistent distributions in concentration and (g - i) colour spaces. We implement a morphological star/galaxy classifier based on PDF analysis, providing meaningful probabilities for J-PLUS sources to one magnitude deeper (r ~ 21) than a classical boolean classification. These probabilities are suited for the statistical study of 150k stars and 101k galaxies with 15 < r < 21 present in the 31.7 deg2 of the J-PLUS EDR. In a future version of the classifier, we will include J-PLUS colour information from twelve photometric bands.

Context. Knowing the exact shape of the UV luminosity function of high-redshift galaxies is important in order to understand the star formation history of the early universe. However, the uncertainties, especially at the faint and bright ends of the LFs, are still significant. Aims. In this paper, we study the UV luminosity function of redshift z = 2.5 - 4.5 galaxies in 2.38 deg^2 of ALHAMBRA data with I <= 24. Thanks to the large area covered by ALHAMBRA, we particularly constrain the bright end of the luminosity function. We also calculate the cosmic variance and the corresponding bias values for our sample and derive their host dark matter halo masses. Methods. We use a novel methodology based on redshift and magnitude probability distribution functions (PDFs). This methodology robustly takes into account the uncertainties due to redshift and magnitude errors, shot noise and cosmic variance, and models the luminosity function in two dimensions (z; M_UV ). Results. We find an excess of bright ~ M*_UV galaxies as compared to the studies based on broad-band photometric data. However, our results agree well with the luminosity function of the magnitude-selected spectroscopic VVDS data. We measure high bias values, b ~ 8 - 10, that are compatible with the previous measurements considering the redshifts and magnitudes of our galaxies and further reinforce the real high-redshift nature of our bright galaxies. Conclusions. We call into question the shape of the luminosity function at its bright end; is it a double power-law as suggested by the recent broad-band photometric studies or rather a brighter Schechter function, as suggested by our multi-filter analysis and the spectroscopic VVDS data.

We study the consistency of the physical properties of galaxies retrieved from SED-fitting as a function of spectral resolution and signal-to-noise ratio (SNR). Using a selection of physically motivated star formation histories, we set up a control sample of mock galaxy spectra representing observations of the local universe in high-resolution spectroscopy, and in 56 narrow-band and 5 broad-band photometry. We fit the SEDs at these spectral resolutions and compute their corresponding the stellar mass, the mass- and luminosity-weighted age and metallicity, and the dust extinction. We study the biases, correlations, and degeneracies affecting the retrieved parameters and explore the r\^ole of the spectral resolution and the SNR in regulating these degeneracies. We find that narrow-band photometry and spectroscopy yield similar trends in the physical properties derived, the former being considerably more precise. Using a galaxy sample from the SDSS, we compare more realistically the results obtained from high-resolution and narrow-band SEDs (synthesized from the same SDSS spectra) following the same spectral fitting procedures. We use results from the literature as a benchmark to our spectroscopic estimates and show that the prior PDFs, commonly adopted in parametric methods, may introduce biases not accounted for in a Bayesian framework. We conclude that narrow-band photometry yields the same trend in the age-metallicity relation in the literature, provided it is affected by the same biases as spectroscopy; albeit the precision achieved with the latter is generally twice as large as with the narrow-band, at SNR values typical of the different kinds of data.

We present a technique that permits the analysis of stellar population gradients in a relatively low cost way compared to IFU surveys analyzing a vastly larger samples as well as out to larger radii. We developed a technique to analyze unresolved stellar populations of spatially resolved galaxies based on photometric multi-filter surveys. We derived spatially resolved stellar population properties and radial gradients by applying a Centroidal Voronoi Tesselation and performing a multi-color photometry SED fitting. This technique has been applied to a sample of 29 massive (M$_{\star}$ > 10$^{10.5}$ M$_{\odot}$), early-type galaxies at $z$ < 0.3 from the ALHAMBRA survey. We produced detailed 2D maps of stellar population properties (age, metallicity and extinction). Radial structures have been studied and luminosity-weighted and mass-weighted gradients have been derived out to 2 - 3.5 R$_\mathrm{eff}$. We find the gradients of early-type galaxies to be on average flat in age ($\nabla$log Age$_\mathrm{L}$ = 0.02 $\pm$ 0.06 dex/R$_\mathrm{eff})$ and negative in metallicity ($\nabla$[Fe/H]$_\mathrm{L}$ = - 0.09 $\pm$ 0.06 dex/R$_\mathrm{eff}$). Overall, the extinction gradients are flat ($\nabla$A$_\mathrm{v}$ = - 0.03 $\pm$ 0.09 mag/R$_\mathrm{eff}$ ) with a wide spread. These results are in agreement with previous studies that used standard long-slit spectroscopy as well as with the most recent integral field unit (IFU) studies. According to recent simulations, these results are consistent with a scenario where early-type galaxies were formed through major mergers and where their final gradients are driven by the older ages and higher metallicity of the accreted systems. We demonstrate the scientific potential of multi-filter photometry to explore the spatially resolved stellar populations of local galaxies and confirm previous spectroscopic trends from a complementary technique.

Our goal is to study the evolution of the $B-$band luminosity function (LF) since $z=1$ using ALHAMBRA data. We used the photometric redshift and the $I-$band selection magnitude probability distribution functions (PDFs) of those ALHAMBRA galaxies with $I\leq24$ mag to compute the posterior LF. We statistically studied quiescent and star-forming galaxies using the template information encoded in the PDFs. The LF covariance matrix in redshift-magnitude-galaxy type space was computed, including the cosmic variance. That was estimated from the intrinsic dispersion of the LF measurements in the 48 ALHAMBRA sub-fields. The uncertainty due to the photometric redshift prior is also included in our analysis. We modelled the LF with a redshift-dependent Schechter function affected by the same selection effects than the data. The measured ALHAMBRA LF at $0.2\leq z<1$ and the evolving Schechter parameters both for quiescent and star-forming galaxies agree with previous results in the literature. The estimated redshift evolution of $M_B^* \propto Qz$ is $Q_{\rm SF}=-1.03\pm0.08$ and $Q_{\rm Q}=-0.80\pm0.08$, and of $\log \phi^* \propto Pz$ is $P_{\rm SF}=-0.01\pm0.03$ and $P_{\rm Q}=-0.41\pm0.05$. The measured faint-end slopes are $\alpha_{\rm SF}=-1.29\pm0.02$ and $\alpha_{\rm Q}=-0.53\pm0.04$. We find a significant population of faint quiescent galaxies, modelled by a second Schechter function with slope $\beta=-1.31\pm0.11$. We find a factor $2.55\pm0.14$ decrease in the luminosity density $j_B$ of star-forming galaxies, and a factor $1.25\pm0.16$ increase in the $j_B$ of quiescent ones since $z=1$, confirming the continuous build-up of the quiescent population with cosmic time. The contribution of the faint quiescent population to $j_B$ increases from 3% at $z=1$ to 6% at $z=0$. The developed methodology will be applied to future multi-filter surveys such as J-PAS.

The original ALHAMBRA catalogue contained over 400,000 galaxies selected using a synthetic F814W image, to the magnitude limit AB(F814W)$\approx$24.5. Given the photometric redshift depth of the ALHAMBRA multiband data (<z>=0.86) and the approximately $I$-band selection, there is a noticeable bias against red objects at moderate redshift. We avoid this bias by creating a new catalogue selected in the $K_s$ band. This newly obtained catalogue is certainly shallower in terms of apparent magnitude, but deeper in terms of redshift, with a significant population of red objects at $z>1$. We select objects using the $K_s$ band images, which reach an approximate AB magnitude limit $K_s \approx 22$. We generate masks and derive completeness functions to characterize the sample. We have tested the quality of the photometry and photometric redshifts using both internal and external checks. Our final catalogue includes $\approx 95,000$ sources down to $K_s \approx 22$, with a significant tail towards high redshift. We have checked that there is a large sample of objects with spectral energy distributions that correspond to that of massive, passively evolving galaxies at $z > 1$, reaching as far as $z \approx 2.5$. We have tested the possibility of combining our data with deep infrared observations at longer wavelengths, particularly Spitzer IRAC data.

We study the clustering of galaxies as a function of spectral type and redshift in the range $0.35 < z < 1.1$ using data from the Advanced Large Homogeneous Area Medium Band Redshift Astronomical (ALHAMBRA) survey. The data cover 2.381 deg$^2$ in 7 fields, after applying a detailed angular selection mask, with accurate photometric redshifts [$\sigma_z < 0.014(1+z)$] down to $I_{AB} < 24$. From this catalog we draw five fixed number density, redshift-limited bins. We estimate the clustering evolution for two different spectral populations selected using the ALHAMBRA-based photometric templates: quiescent and star-forming galaxies. For each sample, we measure the real-space clustering using the projected correlation function. Our calculations are performed over the range $[0.03,10.0] h^{-1}$ Mpc, allowing us to find a steeper trend for $r_p \lesssim 0.2 h^{-1}$ Mpc, which is especially clear for star-forming galaxies. Our analysis also shows a clear early differentiation in the clustering properties of both populations: star-forming galaxies show weaker clustering with evolution in the correlation length over the analysed redshift range, while quiescent galaxies show stronger clustering already at high redshifts, and no appreciable evolution. We also perform the bias calculation where similar segregation is found, but now it is among the quiescent galaxies where a growing evolution with redshift is clearer. These findings clearly corroborate the well known colour-density relation, confirming that quiescent galaxies are mainly located in dark matter halos that are more massive than those typically populated by star-forming galaxies.

The relative cosmic variance ($\sigma_v$) is a fundamental source of uncertainty in pencil-beam surveys and, as a particular case of count-in-cell statistics, can be used to estimate the bias between galaxies and their underlying dark-matter distribution. Our goal is to test the significance of the clustering information encoded in the $\sigma_v$ measured in the ALHAMBRA survey. We measure the cosmic variance of several galaxy populations selected with $B-$band luminosity at $0.35 \leq z < 1.05$ as the intrinsic dispersion in the number density distribution derived from the 48 ALHAMBRA subfields. We compare the observational $\sigma_v$ with the cosmic variance of the dark matter expected from the theory, $\sigma_{v,{\rm dm}}$. This provides an estimation of the galaxy bias $b$. The galaxy bias from the cosmic variance is in excellent agreement with the bias estimated by two-point correlation function analysis in ALHAMBRA. This holds for different redshift bins, for red and blue subsamples, and for several $B-$band luminosity selections. We find that $b$ increases with the $B-$band luminosity and the redshift, as expected from previous work. Moreover, red galaxies have a larger bias than blue galaxies, with a relative bias of $b_{\rm rel} = 1.4 \pm 0.2$. Our results demonstrate that the cosmic variance measured in ALHAMBRA is due to the clustering of galaxies and can be used to characterise the $\sigma_v$ affecting pencil-beam surveys. In addition, it can also be used to estimate the galaxy bias $b$ from a method independent of correlation functions.

We present the main steps that will be taken to extract H$\alpha$ emission flux from Javalambre Photometric Local Universe Survey (J-PLUS) photometric data. For galaxies with $z\lesssim0.015$, the H$\alpha$+[NII] emission is covered by the J-PLUS narrow-band filter $F660$. We explore three different methods to extract the H$\alpha$ + [NII] flux from J-PLUS photometric data: a combination of a broad-band and a narrow-band filter ($r'$ and $F660$), two broad-band and a narrow-band one ($r'$, $i'$ and $F660$), and a SED-fitting based method using 8 photometric points. To test these methodologies, we simulated J-PLUS data from a sample of 7511 SDSS spectra with measured H$\alpha$ flux. Based on the same sample, we derive two empirical relations to correct the derived H$\alpha$+[NII] flux from dust extinction and [NII] contamination. We find that the only unbiased method is the SED fitting based one. The combination of two filters underestimates the measurements of the H$\alpha$ + [NII] flux by a 28%, while the three filters method by a 9%. We study the error budget of the SED-fitting based method and find that, in addition to the photometric error, our measurements have a systematic uncertainty of a 4.3%. Several sources contribute to this uncertainty: differences between our measurement procedure and the one used to derive the spectroscopic values, the use of simple stellar populations as templates, and the intrinsic errors of the spectra, which were not taken into account. Apart from that, the empirical corrections for dust extinction and [NII] contamination add an extra uncertainty of 14%. Given the J-PLUS photometric system, the best methodology to extract H$\alpha$ + [NII] flux is the SED-fitting based one. Using this method, we are able to recover reliable H$\alpha$ fluxes for thousands of nearby galaxies in a robust and homogeneous way.

We consistently analyse for the first time the impact of survey depth and spatial resolution on the most used morphological parameters for classifying galaxies through non-parametric methods: Abraham and Conselice-Bershady concentration indices, Gini, M20 moment of light, asymmetry, and smoothness. Three different non-local datasets are used, ALHAMBRA and SXDS (examples of deep ground-based surveys), and COSMOS (deep space-based survey). We used a sample of 3000 local, visually classified galaxies, measuring their morphological parameters at their real redshifts (z ~ 0). Then we simulated them to match the redshift and magnitude distributions of galaxies in the non-local surveys. The comparisons of the two sets allow to put constraints on the use of each parameter for morphological classification and evaluate the effectiveness of the commonly used morphological diagnostic diagrams. All analysed parameters suffer from biases related to spatial resolution and depth, the impact of the former being much stronger. When including asymmetry and smoothness in classification diagrams, the noise effects must be taken into account carefully, especially for ground-based surveys. M20 is significantly affected, changing both the shape and range of its distribution at all brightness levels.We suggest that diagnostic diagrams based on 2 - 3 parameters should be avoided when classifying galaxies in ground-based surveys, independently of their brightness; for COSMOS they should be avoided for galaxies fainter than F814 = 23.0. These results can be applied directly to surveys similar to ALHAMBRA, SXDS and COSMOS, and also can serve as an upper/lower limit for shallower/deeper ones.

We present a catalogue of 348 galaxy clusters and groups with $0.2<z<1.2$ selected in the 2.78 $deg^2$ ALHAMBRA Survey. The high precision of our photometric redshifts, close to $1\%$, and the wide spread of the seven ALHAMBRA pointings ensure that this catalogue has better mass sensitivity and is less affected by cosmic variance than comparable samples. The detection has been carried out with the Bayesian Cluster Finder (BCF), whose performance has been checked in ALHAMBRA-like light-cone mock catalogues. Great care has been taken to ensure that the observable properties of the mocks photometry accurately correspond to those of real catalogues. From our simulations, we expect to detect galaxy clusters and groups with both $70\%$ completeness and purity down to dark matter halo masses of $M_h\sim3\times10^{13}\rm M_{\odot}$ for $z<0.85$. Cluster redshifts are expected to be recovered with $\sim0.6\%$ precision for $z<1$. We also expect to measure cluster masses with $\sigma_{M_h|M^*_{CL}}\sim0.25-0.35\, dex$ precision down to $\sim3\times10^{13}\rm M_{\odot}$, masses which are $50\%$ smaller than those reached by similar work. We have compared these detections with previous optical, spectroscopic and X-rays work, finding an excellent agreement with the rates reported from the simulations. We have also explored the overall properties of these detections such as the presence of a colour-magnitude relation, the evolution of the photometric blue fraction and the clustering of these sources in the different ALHAMBRA fields. Despite the small numbers, we observe tentative evidence that, for a fixed stellar mass, the environment is playing a crucial role at lower redshifts (z$<$0.5).

We present MUFFIT, a new generic code optimized to retrieve the main stellar population parameters of galaxies in photometric multi-filter surveys, and we check its reliability and feasibility with real galaxy data from the ALHAMBRA survey. Making use of an error-weighted $\chi^2$-test, we compare the multi-filter fluxes of galaxies with the synthetic photometry of mixtures of two single stellar populations at different redshifts and extinctions, to provide through a Monte Carlo method the most likely range of stellar population parameters (mainly ages and metallicities), extinctions, redshifts, and stellar masses. To improve the diagnostic reliability, MUFFIT identifies and removes from the analysis those bands that are significantly affected by emission lines. We highlight that the retrieved age-metallicity locus for a sample of $z \le 0.22$ early-type galaxies in ALHAMBRA at different stellar mass bins are in very good agreement with the ones from SDSS spectroscopic diagnostics. Moreover, a one-to-one comparison between the redshifts, ages, metallicities, and stellar masses derived spectroscopically for SDSS and by MUFFIT for ALHAMBRA reveals good qualitative agreements in all the parameters. In addition, and using as input the results from photometric-redshift codes, MUFFIT improves the photometric-redshift accuracy by $\sim 10$-$20\%$, and it also detects nebular emissions in galaxies, providing physical information about their strengths. Our results show the potential of multi-filter galaxy data to conduct reliable stellar population studies with the appropiate analysis techniques, as MUFFIT.

Context. Most observational results on the high redshift restframe UV-bright galaxies are based on samples pinpointed using the so called dropout technique or Ly-alpha selection. However, the availability of multifilter data allows now replacing the dropout selections by direct methods based on photometric redshifts. In this paper we present the methodology to select and study the population of high redshift galaxies in the ALHAMBRA survey data. Aims. Our aim is to develop a less biased methodology than the traditional dropout technique to study the high redshift galaxies in ALHAMBRA and other multifilter data. Thanks to the wide area ALHAMBRA covers, we especially aim at contributing in the study of the brightest, less frequent, high redshift galaxies. Methods. The methodology is based on redshift probability distribution functions (zPDFs). It is shown how a clean galaxy sample can be obtained by selecting the galaxies with high integrated probability of being within a given redshift interval. However, reaching both a complete and clean sample with this method is challenging. Hence, a method to derive statistical properties by summing the zPDFs of all the galaxies in the redshift bin of interest is introduced. Results. Using this methodology we derive the galaxy rest frame UV number counts in five redshift bins centred at z=2.5, 3.0, 3.5, 4.0, and 4.5, being complete up to the limiting magnitude at m_UV(AB)=24. With the wide field ALHAMBRA data we especially contribute in the study of the brightest ends of these counts, sampling well the surface densities down to m_UV(AB)=21-22. Conclusions. We show that using the zPDFs it is easy to select a clean sample of high redshift galaxies. We also show that statistical analysis of the properties of galaxies is better done using a probabilistic approach, which takes into account both the incompleteness and contamination in a natural way.

Our goal is to develop and test a novel methodology to compute accurate close pair fractions with photometric redshifts. We improve the current methodologies to estimate the merger fraction f_m from photometric redshifts by (i) using the full probability distribution functions (PDFs) of the sources in redshift space, (ii) including the variation in the luminosity of the sources with z in both the selection of the samples and in the luminosity ratio constrain, and (iii) splitting individual PDFs into red and blue spectral templates to deal robustly with colour selections. We test the performance of our new methodology with the PDFs provided by the ALHAMBRA photometric survey. The merger fractions and rates from the ALHAMBRA survey are in excellent agreement with those from spectroscopic work, both for the general population and for red and blue galaxies. With the merger rate of bright (M_B <= -20 - 1.1z) galaxies evolving as (1+z)^n, the power-law index n is larger for blue galaxies (n = 2.7 +- 0.5) than for red galaxies (n = 1.3 +- 0.4), confirming previous results. Integrating the merger rate over cosmic time, we find that the average number of mergers per galaxy since z = 1 is N_m = 0.57 +- 0.05 for red galaxies and N_m = 0.26 +- 0.02 for blue galaxies. Our new methodology exploits statistically all the available information provided by photometric redshift codes and provides accurate measurements of the merger fraction by close pairs only using photometric redshifts. Current and future photometric surveys will benefit of this new methodology.

The ALHAMBRA (Advance Large Homogeneous Area Medium Band Redshift Astronomical) survey has observed 8 different regions of the sky, including sections of the COSMOS, DEEP2, ELAIS, GOODS-N, SDSS and Groth fields using a new photometric system with 20 contiguous ~ $300\AA$ filters covering the optical range, combining them with deep $JHKs$ imaging. The observations, carried out with the Calar Alto 3.5m telescope using the wide field (0.25 sq. deg FOV) optical camera LAICA and the NIR instrument Omega-2000, correspond to ~700hrs on-target science images. The photometric system was designed to maximize the effective depth of the survey in terms of accurate spectral-type and photo-zs estimation along with the capability of identification of relatively faint emission lines. Here we present multicolor photometry and photo-zs for ~438k galaxies, detected in synthetic F814W images, complete down to I~24.5 AB, taking into account realistic noise estimates, and correcting by PSF and aperture effects with the ColorPro software. The photometric ZP have been calibrated using stellar transformation equations and refined internally, using a new technique based on the highly robust photometric redshifts measured for emission line galaxies. We calculate photometric redshifts with the BPZ2 code, which includes new empirically calibrated templates and priors. Our photo-zs have a precision of $dz/(1+z_s)=1%$ for I<22.5 and 1.4% for 22.5<I<24.5. Precisions of less than 0.5% are reached for the brighter spectroscopic sample, showing the potential of medium-band photometric surveys. The global $P(z)$ shows a mean redshift <z>=0.56 for I<22.5 AB and <z>=0.86 for I<24.5 AB. The data presented here covers an effective area of 2.79 sq. deg, split into 14 strips of 58.5'x15.5' and represents ~32 hrs of on-target.

We study the clustering of galaxies as function of luminosity and redshift in the range $0.35 < z < 1.25$ using data from the Advanced Large Homogeneous Area Medium Band Redshift Astronomical (ALHAMBRA) survey. The ALHAMBRA data used in this work cover $2.38 \mathrm{deg}^2$ in 7 independent fields, after applying a detailed angular selection mask, with accurate photometric redshifts, $\sigma_z \lesssim 0.014 (1+z)$, down to $I_{\rm AB} < 24$. Given the depth of the survey, we select samples in $B$-band luminosity down to $L^{\rm th} \simeq 0.16 L^{*}$ at $z = 0.9$. We measure the real-space clustering using the projected correlation function, accounting for photometric redshifts uncertainties. We infer the galaxy bias, and study its evolution with luminosity. We study the effect of sample variance, and confirm earlier results that the COSMOS and ELAIS-N1 fields are dominated by the presence of large structures. For the intermediate and bright samples, $L^{\rm med} \gtrsim 0.6L^{*}$, we obtain a strong dependence of bias on luminosity, in agreement with previous results at similar redshift. We are able to extend this study to fainter luminosities, where we obtain an almost flat relation, similar to that observed at low redshift. Regarding the evolution of bias with redshift, our results suggest that the different galaxy populations studied reside in haloes covering a range in mass between $\log_{10}[M_{\rm h}/(h^{-1}\mathrm{M}_{\odot})] \gtrsim 11.5$ for samples with $L^{\rm med} \simeq 0.3 L^{*}$ and $\log_{10}[M_{\rm h}/(h^{-1}\mathrm{M}_{\odot})] \gtrsim 13.0$ for samples with $L^{\rm med} \simeq 2 L^{*}$, with typical occupation numbers in the range of $\sim 1 - 3$ galaxies per halo.

There are many ways to solve the challenging problem of making a high performance robotic observatory from scratch. The Observatorio Astrof\'isico de Javalambre (OAJ) is a new astronomical facility located at the Sierra de Javalambre (Teruel, Spain) whose primary role will be to conduct all-sky astronomical surveys. The OAJ control system has been designed under a global point of view including not only astronomical subsystems but also infrastructure and other facilities. Three main factors have been considered in the design of a global control system for the robotic OAJ: quality, reliability and efficiency. We propose CIA (Control Integrated Architecture) design and OEE (Overall Equipment Effectiveness) as a key performance indicator in order to improve operation processes, minimizing resources and obtain high cost reduction maintaining quality requirements. The OAJ subsystems considered for the control integrated architecture are the following: two wide-field telescopes and their instrumentation, active optics subsystems, facilities for sky quality monitoring (seeing, extinction, sky background, sky brightness, clouds distribution, meteorological station), domes and several infrastructure facilities such as water supply, glycol water, water treatment plant, air conditioning, compressed air, LN2 plant, illumination, surveillance, access control, fire suppression, electrical generators, electrical distribution, electrical consumption, communication network, Uninterruptible Power Supply and two main control rooms, one at the OAJ and other remotely located in Teruel at 40km from the observatory, connected through a microwave radio-link. Here we present the OAJ strategy in control design to achieve maximum quality efficiency for the observatory processes and operations, giving practical examples of our approach.

Our goal is to estimate empirically, for the first time, the cosmic variance that affects merger fraction studies based on close pairs. We compute the merger fraction from photometric redshift close pairs with 10h^-1 kpc <= rp <= 50h^-1 kpc and Dv <= 500 km/s, and measure it in the 48 sub-fields of the ALHAMBRA survey. We study the distribution of the measured merger fractions, that follow a log-normal function, and estimate the cosmic variance sigma_v as the intrinsic dispersion of the observed distribution. We develop a maximum likelihood estimator to measure a reliable sigma_v and avoid the dispersion due to the observational errors (including the Poisson shot noise term). The cosmic variance of the merger fraction depends mainly on (i) the number density of the populations under study, both for the principal (n_1) and the companion (n_2) galaxy in the close pair, and (ii) the probed cosmic volume V_c. We find a significant dependence on neither the search radius used to define close companions, the redshift, nor the physical selection (luminosity or stellar mass) of the samples. We provide a parametrisation of the cosmic variance with n_1, n_2, and V_c, sigma_v = 0.48 n_1^{-0.54} V_c^{-0.48} (n_2/n_1)^{-0.37}. Thanks to this prescription, future merger fraction studies based on close pairs could account properly for the cosmic variance on their results.

ALHAMBRA is a photometric survey designed to trace the cosmic evolution and cosmic variance. It covers a large area of ~ 4 sq. deg in 8 fields, where 7 fields overlap with other surveys, allowing to have complementary data in other wavelengths. All observations were carried out in 20 continuous, medium band (30 nm width) optical and 3 near-infrared (JHK) bands, providing the precise measurements of photometric redshifts. In addition, morphological classification of galaxies is crucial for any kind of galaxy formation and cosmic evolution studies, providing the information about star formation histories, their environment and interactions, internal perturbations, etc. We present a morphological classification of > 40,000 galaxies in the ALHAMBRA survey. We associate to every galaxy a probability to be early-type using the automated Bayesian code galSVM. Despite of the spatial resolution of the ALHAMBRA images (~ 1 arcsec), for 22,051 galaxies we obtained the contamination by other type of less than 10%. Of those, 1,640 and 10,322 galaxies are classified as early- (down to redshifts ~ 0.5) and late-type (down to redshifts ~ 1.0), respectively, with magnitudes F613W < 22.0. In addition, for magnitude range 22.0 < F613W < 23.0 we classified other 10,089 late-type galaxies with redshifts < 1.3. We show that the classified objects populate the expected regions in the colour-mass and colour-magnitude planes. The presented dataset is especially attractive given the homogeneous multi-wavelength coverage available in the ALHAMBRA fields, and is intended to be used in a variety of scientific applications. The low-contamination catalogue (< 10%) is made publicly available with the present paper.

We aim to illustrate the potentiality of the Advanced Large, Homogeneous Area, Medium-Band Redshift Astronomical (ALHAMBRA) survey to investigate the high redshift universe through the detection of quasi stellar objects (QSOs) at redshifts larger than 5. The search for z>5 QSOs candidates was done by fitting an extensive library of spectral energy distributions --including active and non-active galaxy templates as well as stars-- to the photometric database of the ALHAMBRA survey (composed of 20 optical medium-band plus the 3 broad-band JHKs filters). Our selection over ~1 square degree of ALHAMBRA data (~1/4 of the total area covered by the survey), combined with GTC/OSIRIS spectroscopy, has yielded the identification of an optically faint QSO at very high redshift (z = 5.41). The QSO has an absolute magnitude of ~-24 at the 1450{\AA} continuum, a bolometric luminosity of ~2x10^46 erg/s and an estimated black hole mass of ~10^8 Msolar. This QSO adds itself to a reduced number of known UV faint sources at these redshifts. The preliminary derived space density is compatible with the most recent determinations of the high-z QSO luminosity functions (QLF). This new detection shows how ALHAMBRA, as well as forthcoming well designed photometric surveys, can provide a wealth of information on the origin and early evolution of this kind of objects.

We characterize the ability of the ALHAMBRA survey to assign accurate photo-z's to BLAGN and QSOs based on their ALHAMBRA very-low-resolution optical-NIR spectroscopy. A sample of 170 spectroscopically identified BLAGN and QSOs have been used together with a library of templates (including SEDs from AGN, normal, starburst galaxies and stars) in order to fit the 23 photometric data points provided by ALHAMBRA in the optical and NIR (20 medium-band optical filters plus the standard JHKs). We find that the ALHAMBRA photometry is able to provide an accurate photo-z and spectral classification for ~88% of the spectroscopic sources over 2.5 deg^2 in different areas of the survey, all of them brighter than m678=23.5 (equivalent to r(SLOAN)~24.0). The derived photo-z accuracy is better than 1% and comparable to the most recent results in other cosmological fields. The fraction of outliers (~12%) is mainly caused by the larger photometric errors for the faintest sources and the intrinsic variability of the BLAGN/QSO population. A small fraction of outliers may have an incorrectly assigned spectroscopic redshift. The definition of the ALHAMBRA survey in terms of the number of filters, filter properties, area coverage and depth is able to provide photometric redshifts for BLAGN/QSOs with a precision similar to any previous survey that makes use of medium-band optical photometry. In agreement with previous literature results, our analysis also reveals that, in the 0<z<4 redshift interval, very accurate photo-z can be obtained without the use of near-IR broadband photometry at the expense of a slight increase of outliers. The NIR importance is expected to increase at higher redshifts (z>4). These results are relevant for the design of future optical follow-ups of surveys with a large fraction of BLAGN, as it is the case for X-rays or radio surveys.

We present the stellar velocity dispersion measurements for 5 Luminous Compact Galaxies (LCGs) at z=0.5-0.7. These galaxies are vigorously forming stars with average SFR $\sim$ 40 M$_{\odot}$/yr. We find that their velocity dispersions range from $\sim137\ \rm{km/s}$ to $260\ \rm{km/s}$, while their stellar masses range between $4\times 10^{9}$ and $10^{11}$ M$_{\odot}$. If these LCGs evolve passively after this major burst of star formation, their masses and velocity dispersions, as well as their evolved colours and luminosities are most consistent with the values characteristic of early-type spiral galaxies today.

We present a new cataloge of EROs from the Groth strip and study the relation between their morphology and mass. We find 102 EROs (F814W-K=>4, K<=21.0), over a survey area of 155 arcmin^2. The photometric data include U,B,F606W,F814W,J,K bands. Morphologies are based on a by eye classification and we distinguish between 3 basic classes: compact objects, targets with a disc and/or a bulge component and irregular or merger candidates. The majority of our targets has either a very compact morphology (33+-6%), or show more or less distinct disc components (41+-6%). 14+-4% are merger or irregulars and 7 objects could not be classified. We also study the dependence of structural parameters on morphological appearance. EROs that are either compact or show a distinct bulge component have smaller median effective radii (1.22+-0.14 kpc and 3.31+-0.53 kpc) than disc dominated (5.50+-0.51 kpc) or possible irregular galaxies or merger candidates (4.92+-0.14 kpc). The Sersic index changes from 2.30+-0.34 and 3.24+-0.55, to 1.03+-0.24 and 1.54+-0.40 respectively. Most the EROs in our sample have redshifts between z=1 and z=2; however, compact EROs in our sample are found at redshifts as low as z=0.4 and as high as z=2.8; the latter qualify as well as DRGs. Disc-like EROs are also found up to z=2.8; however those with a bulge-disc structure are only seen at z<1.5. For each of these EROs we determined the stellar mass and mean population age by fitting synthetic Bruzual (2007) spectra to the SED. Mass estimates were obtained by assuming an exponentially declining star formation rate. Total stellar masses are in the range 9.1<log(M/M_sun)<11.6. We cannot detect significant differences between the stellar mass distribution of the morphological classes. EROs with masses of log(M/M_sun)>11.0 dominantly show compact morphologies, but also include a significant number of sources with a disc morphology.

(Abridged) We measure the fraction of galaxies undergoing disk-disk major mergers (f_m) at intermediate redshifts (0.35 <= z < 0.85) by studying the asymmetry index A of galaxy images. Results are provided for B- and Ks-band absolute magnitude selected samples from the Groth strip in the GOYA photometric survey. Three sources of systematic error are carefully addressed: (i) we avoid morphological K-corrections, (ii) we measure asymmetries in artificially redshifted to z_d = 0.75 galaxies to lead with loss of morphological information with redshift, and (iii) we take into account the observational errors in z and A, that tend to overestimate the merger fraction, by maximum likelihood techniques. We find: (i) our data allow for a robust merger fraction to be provided for a single redshift bin centered at z=0.6. (ii) Merger fractions have low values: f_m = 0.045 for M_B <= -20 galaxies, and f_m = 0.031 for M_Ks <= -23.5 galaxies. And, (iii) failure to address the effects of the observational errors leads to overestimating f_m by factors of 10%-60%. Combining our results with those on literature, and parameterizing the merger fraction evolution as f_m(z) = f_m(0)(1+z)^m, we obtain that m = 2.9 +- 0.8, and f_m(0) = 0.012 +- 0.004$. Assuming a Ks-band mass-to-light ratio not varying with luminosity, we infer that the merger rate of galaxies with stellar mass M >= 3.5x10^10 M_Sun is R_m = 1.6x10^-4 Mpc^-3 Gyr^-1. When we compare with previous studies at similar redshifts, we find that the merger rate decreases when mass increases.

(ABRIDGED) We describe the first results of the ALHAMBRA survey which provides cosmic tomography of the evolution of the contents of the Universe over most of Cosmic history. Our approach employs 20 contiguous, equal-width, medium-band filters covering from 3500 to 9700 A, plus the JHKs bands, to observe an area of 4 sqdeg on the sky. The optical photometric system has been designed to maximize the number of objects with accurate classification by SED and redshift, and to be sensitive to relatively faint emission lines. The observations are being carried out with the Calar Alto 3.5m telescope using the cameras LAICA and O-2000. The first data confirm that we are reaching the expected magnitude limits of AB<~25 mag in the optical filters from the blue to 8300 A, and from AB=24.7 to 23.4 for the redder ones. The limit in the NIR is (Vega) K_s~20, H~21, J~22. We expect to obtain accurate redshift values, Delta z/(1+z) <~ 0.03 for about 5x10^5 galaxies with I<~25 (60% complete), and z_med=0.74. This accuracy, together with the homogeneity of the selection function, will allow for the study of the redshift evolution of the large scale structure, the galaxy population and its evolution with redshift, the identification of clusters of galaxies, and many other studies, without the need for any further follow-up. It will also provide targets for detailed studies with 10m-class telescopes. Given its area, spectral coverage and its depth, apart from those main goals, the ALHAMBRA-Survey will also produce valuable data for galactic studies.

We summarize the multiwavelength properties of X-ray sources detected in the 80 ks XMM-Newton observation of the Groth-Westphal Strip. We find 23 XMM-Newton sources within the WFPC2 fields. Ten spectroscopic redshifts are available from the DEEP and CFRS projects and 4 of these show broad Mg II emission (type 1 AGNs). Two of those without any broad lines, nevertheless, have [NeV] emission which is an unambiguous signature of AGN activity, one of which is a narrow-line Seyfert 1 and the other a type 2 AGN. We have made near-infrared (NIR) spectroscopic observations using the Subaru OHS/CISCO spectrometer for five of the X-ray sources for which we found no indication of an AGN activity in the optical spectrum. We have detected H-alpha+[NII] emission in four of them. A broad H-alpha component and/or a large [NII]/H-alpha ratio is seen, suggestive of AGN activity. Nineteen sources have been detected in the Ks band and four of these are extremely red objects (I814-Ks>4). The optical counterparts for the majority of the X-ray sources are bulge-dominated with colors consistent with evolving elliptical galaxies, with starburst/AGN contamination. Assuming that the known local relations among the bulge luminosity,central velocity dispersion, and the mass of the central blackhole hold at about z=1, the AGN bolometric luminosity to Eddington luminosity ratio ranges from 0.3% to 10%. (abridged)