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)

This work employs high-quality {\em Hubble Space Telescope} ({\em HST}) Advanced Camera for Surveys (ACS) F606W and F814W photometry to correct for the differential reddening affecting the colour-magnitude diagram (CMD) of the poorly-studied globular cluster (GC) Palomar\,2. Differential reddening is taken into account by assuming that morphological differences among CMDs extracted across the field of view of Palomar\,2 correspond essentially to shifts (quantified in terms of $\delta E(B-V)$) along the reddening vector due to a non-uniform dust distribution. The average reddening difference over all partial CMDs is $\overline{\delta E(B-V)}=0.24\pm0.08$, with the highest reaching $\delta E(B-V) =0.52$. The corrected CMD displays well-defined and relatively narrow evolutionary sequences, especially for the evolved stars, i.e. the red-giant, horizontal and asymptotic giant branches (RGB, HB and AGB, respectively). The average width of the upper main sequence and RGB profiles of the corrected CMD corresponds to 56\% of the original one. Parameters measured on this CMD show that Palomar\,2 is $\approx13.25$\,Gyr old, has the mass $M\sim1.4\times10^5\,m_\odot$ stored in stars, is affected by the foreground $E(B-V)\approx0.93$, is located at $d_\odot\approx26$\,Kpc from the Sun, and is characterized by the global metallicity $Z/Zo\approx0.03$, which corresponds to the range $-1.9\leq [Fe/H] \leq-1.6$ (for $0.0\leq[\alpha/Fe]\leq+0.4$), quite consistent with other outer halo GCs. Additional parameters are the absolute magnitude $M_V\approx-7.8$, and the core and half-light radii $r_C\approx2.6$\,pc and $R_{HL}\approx4.7$\,pc, respectively.

Despite numerous efforts, it is still unclear whether lenticular galaxies (S0s) evolve from spirals whose star formation was suppressed, or formed trough mergers or disk instabilities. In this paper we present a pilot study of 21 S0 galaxies in extreme environments (field and cluster), and compare their spatially-resolved kinematics and global stellar populations. Our aim is to identify whether there are different mechanisms that form S0s in different environments. Our results show that the kinematics of S0 galaxies in field and cluster are, indeed, different. Lenticulars in the cluster are more rotationally supported, suggesting that they are formed through processes that involve the rapid consumption or removal of gas (e.g. starvation, ram pressure stripping). In contrast, S0s in the field are more pressure supported, suggesting that minor mergers served mostly to shape their kinematic properties. These results are independent of total mass, luminosity, or disk-to-bulge ratio. On the other hand, the mass-weighted age, metallicity, and star formation time-scale of the galaxies correlate more with mass than with environment, in agreement with known relations from previous work such as the one between mass and metallicity. Overall, our results re-enforce the idea that there are multiple mechanisms that produce S0s, and that both mass $and$ environment play key roles. A larger sample is highly desirable to confirm or refute the results and the interpretation of this pilot study.

We present a star/galaxy classification for the Southern Photometric Local Universe Survey (S-PLUS), based on a Machine Learning approach: the Random Forest algorithm. We train the algorithm using the S-PLUS optical photometry up to $r$=21, matched to SDSS/DR13, and morphological parameters. The metric of importance is defined as the relative decrease of the initial accuracy when all correlations related to a certain feature is vanished. In general, the broad photometric bands presented higher importance when compared to narrow ones. The influence of the morphological parameters has been evaluated training the RF with and without the inclusion of morphological parameters, presenting accuracy values of 95.0\% and 88.1\%, respectively. Particularly, the morphological parameter {\rm FWHM/PSF} performed the highest importance over all features to distinguish between stars and galaxies, indicating that it is crucial to classify objects into stars and galaxies. We investigate the misclassification of stars and galaxies in the broad-band colour-colour diagram $(g-r)$ versus $(r-i)$. The morphology can notably improve the classification of objects at regions in the diagram where the misclassification was relatively high. Consequently, it provides cleaner samples for statistical studies. The expected contamination rate of red galaxies as a function of the redshift is estimated, providing corrections for red galaxy samples. The classification of QSOs as extragalactic objects is slightly better using photometric-only case. An extragalactic point-source catalogue is provided using the classification without any morphology feature (only the SED information) with additional constraints on photometric redshifts and {\rm FWHM/PSF} values.

The Southern Photometric Local Universe Survey (S-PLUS) is imaging ~9300 deg^2 of the celestial sphere in twelve optical bands using a dedicated 0.8 m robotic telescope, the T80-South, at the Cerro Tololo Inter-American Observatory, Chile. The telescope is equipped with a 9.2k by 9.2k e2v detector with 10 um pixels, resulting in a field-of-view of 2 deg^2 with a plate scale of 0.55"/pixel. The survey consists of four main subfields, which include two non-contiguous fields at high Galactic latitudes (8000 deg^2 at |b| > 30 deg) and two areas of the Galactic plane and bulge (for an additional 1300 deg^2). S-PLUS uses the Javalambre 12-band magnitude system, which includes the 5 u, g, r, i, z broad-band filters and 7 narrow-band filters centered on prominent stellar spectral features: the Balmer jump/[OII], Ca H+K, H-delta, G-band, Mg b triplet, H-alpha, and the Ca triplet. S-PLUS delivers accurate photometric redshifts (delta_z/(1+z) = 0.02 or better) for galaxies with r < 20 AB mag and redshift < 0.5, thus producing a 3D map of the local Universe over a volume of more than 1 (Gpc/h)^3. The final S-PLUS catalogue will also enable the study of star formation and stellar populations in and around the Milky Way and nearby galaxies, as well as searches for quasars, variable sources, and low-metallicity stars. In this paper we introduce the main characteristics of the survey, illustrated with science verification data highlighting the unique capabilities of S-PLUS. We also present the first public data release of ~336 deg^2 of the Stripe-82 area, which is available at http://datalab.noao.edu/splus.

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.

In this paper we present a thorough discussion about the photometric redshift (photo-z) performance of the Southern Photometric Local Universe Survey (S-PLUS). This survey combines a 7 narrow + 5 broad passband filter system, with a typical photometric-depth of r$\sim$21 AB. For this exercise, we utilize the Data Release 1 (DR1), corresponding to 336 deg$^{2}$ from the Stripe-82 region. We rely on the \texttt{BPZ2} code to compute our estimates, using a new library of SED models, which includes additional templates for quiescent galaxies. When compared to a spectroscopic redshift control sample of $\sim$100k galaxies, we find a precision of $\sigma_{z}<$0.8\%, $<$2.0\% or $<$3.0\% for galaxies with magnitudes r$<$17, $<$19 and $<$21, respectively. A precision of 0.6\% is attained for galaxies with the highest \texttt{Odds} values. These estimates have a negligible bias and a fraction of catastrophic outliers inferior to 1\%. We identify a redshift window (i.e., 0.26$<z<$0.32) where our estimates double their precision, due to the simultaneous detection of two emission-lines in two distinct narrow-bands; representing a window opportunity to conduct statistical studies such as luminosity functions. We forecast a total of $\sim$2M, $\sim$16M and $\sim$32M galaxies in the S-PLUS survey with a photo-z precision of $\sigma_{z}<$1.0\%, $<$2.0\% and $<$2.5\% after observing 8000 $deg^{2}$. We also derive redshift Probability Density Functions, proving their reliability encoding redshift uncertainties and their potential recovering the $n(z)$ of galaxies at $z<0.4$, with an unprecedented precision for a photometric survey in the southern hemisphere.

We analyze a set of optical-to-near-infrared long-slit nuclear spectra of 16 infrared-luminous spiral galaxies. All of the studied sources present H$_2$ emission, which reflects the star-forming nature of our sample, and they clearly display H I emission lines in the optical. Their continua contain many strong stellar absorption lines, with the most common features due to Ca I, Ca II, Fe I, Na I, Mg I, in addition to prominent absorption bands of TiO, VO, ZrO, CN and CO. We report a homogeneous set of equivalent width (EW) measurements for 45 indices, from optical to NIR species for the 16 star-forming galaxies as well as for 19 early type galaxies where we collected the data from the literature. This selected set of emission and absorption-feature measurements can be used to test predictions of the forthcoming generations of stellar population models. We find correlations among the different absorption features and propose here correlations between optical and NIR indices, as well as among different NIR indices, and compare them with model predictions. While for the optical absorption features the models consistently agree with the observations,the NIR indices are much harder to interpret. For early-type spirals the measurements agree roughly with the models, while for star-forming objects they fail to predict the strengths of these indices.

We present the results of a systematic search and characterisation of jellyfish galaxy candidates in the multi-cluster system A901/2, at z ~ 0.165, as part of the OMEGA survey. By visual inspecting ACS/HST F606W images looking for morphological signatures of ram-pressure stripping events in Halpha-emitting galaxies, we identify a total of 70 candidates. Out of these, 53 are clearly star-forming galaxies and 5 are highly probable AGN hosts, the classification of the remaining galaxies is more uncertain. They have late-type and irregular morphologies and most of them are part of the blue cloud with only 4 being previously classified as dusty reds. The AGN activity is not prominent in the sample and, of the few cases of galaxies hosting AGN, such activity does not seem to be correlated to the gas stripping phenomenon. Our jellyfish galaxy candidates do not have a preferential pattern of motion within the multi-cluster system, although the most compelling cases appear to inhabit the inner regions of the most massive sub-cluster centres. The sSFR of these galaxies indicates that their star formation activity is enhanced, in contrast with what is observed for the rest of the star-forming galaxy population in the system. Half of the sample is forming stars at a higher rate than the main-sequence for field galaxies and this behaviour is more evident for the most compelling candidates. For some galaxies, the spatially resolved Halpha emission appears to be as disturbed and extended as their continuum counterparts. Our findings point towards a scenario where the ram pressure stripping is triggering a period of intense and extended star formation throughout the galaxy while it is also disturbing the morphology. This is the largest sample of jellyfish galaxy candidates found in a single system suggesting that cluster mergers might be the ideal environment for studying ram pressure stripping effects.

The A901/2 system is a rare case of galaxy cluster interaction, in which two galaxy clusters and two smaller groups are seen in route of collision with each other simultaneously. Within each of the four substructures, several galaxies with features indicative of jellyfish morphologies have been observed. In this paper, we propose a hydrodynamic model for the merger as a whole, compatible with its diffuse X-ray emission, and correlate the gas properties in this model with the locations of the jellyfish galaxy candidates in the real system. We find that jellyfish galaxies seem to be preferentially located near a boundary inside each subcluster where diffuse gas moving along with the subcluster and diffuse gas from the remainder of the system meet. The velocity change in those boundaries is such that a factor of up to $\sim$1000 increase in the ram pressure takes place within a few hundred kpc, which could trigger the high rate of gas loss necessary for a jellyfish morphology to emerge. A theoretical treatment of ram pressure stripping in the environment of galaxy cluster mergers has not been presented in the literature so far; we propose that this could be a common scenario for the formation of jellyfish morphologies in such systems.

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.

We study the effect of inclination on the apparent brightness of star-forming galaxies in spectral passbands that are commonly used as star-formation indicators. As diagnostics we use mass-to-light ratios in three passbands: the UV continuum at 280 nm, the H$\alpha$ emission line, and the FIR 24$\mu$-band. We include a study of inclination trends in the IR/UV ratio ("IRX") and the IR/H$\alpha$ ratio. Our sample comprises a few hundred galaxies from the region around the clusters Abell 901/902 with deep data and inclinations measured from outer disks in Hubble Space Telescope images. As a novelty, the H$\alpha$- and separately the NII-emission are measured by tunable-filter imaging and encompass galaxies in their entirety. At galaxy stellar masses above log $M_*/M_\odot > 10$ we find trends in the UV and H$\alpha$ mass-to-light ratio that suggest an inclination-induced attenuation from face-on to edge-on of $\sim 1$ mag and $\sim 0.7$ mag in UV and H$\alpha$, respectively, implying that star-formation rates of edge-on galaxies would be underestimated by $\sim 2.5\times$ in UV and $\sim 2\times$ in H$\alpha$. We find the luminosities in UV and H$\alpha$ to be well correlated, but the optical depth of diffuse dust that causes inclination dependence appears to be lower for stars emitting at 280 nm than for gas clouds emitting Balmer lines. For galaxies with log $M_*/M_\odot < 9.7$, we find no measurable effect at $>0.1$ mag. The absence of an inclination dependence at 24$\mu$ confirms that the average galaxy is optically thin in the FIR.

In the hierarchical two-phase formation scenario, the extended halos of early type galaxies (ETGs) are expected to have different physical properties from those of the galaxies' central regions. This work aims at characterizing the kinematic properties of ETG halos using planetary nebulae (PNe) as tracers, which allow us to overcome the limitations of absorption line spectroscopy of continuum at low surface brightness. We present two-dimensional velocity and velocity dispersion fields for 33 ETGs, including both fast (FRs) and slow rotators (SRs), making this the largest kinematic survey to-date of extragalactic PNe. The velocity fields extend out to a median 5.6 effective radii (Re), with a range [3Re-13Re] for the PN.S ETGs. We complemented the PN kinematics with absorption line data from the literature. We find that ETGs typically show a kinematic transition between inner regions and halo. Estimated transition radii in units of Re anticorrelate with stellar mass. SRs have increased but still modest rotational support at large radii, while most of the FRs show a decrease in rotation, due to the fading of the stellar disk in the outer, more slowly rotating spheroid. 30% of the FRs are dominated by rotation also at large radii. Most ETGs have flat or slightly falling halo velocity dispersion profiles, but 15% of the sample have steeply falling profiles. 40% of the FRs show kinematic twists, misalignments, or rotation along two axes, indicating that they turn from oblate at the center to triaxial in the halo, consistently with the distribution of isophote twists from extended photometry. ETGs have more diverse kinematic properties in their halos than in the central regions. The observed kinematic transition to the halo and its dependence on stellar mass is consistent with LambdaCDM simulations and supports a two-phase formation scenario. [abridged]

Recovering the origins of lenticular galaxies can shed light on the understanding of galaxy formation and evolution, since they present properties that can be found in both elliptical and spiral galaxies. In this work we study the kinematics of the globular cluster (GC) systems of three lenticular galaxies located in low density environments (NGC 2768, NGC 3115 and NGC 7457), and compare them with the kinematics of planetary nebulae (PNe). The PNe and GC data come from the Planetary Nebulae Spectrograph and the SLUGGS Surveys. Through photometric spheroid-disc decomposition and PNe kinematics we find the probability for a given GC to belong to either the spheroid or the disc of its host galaxy or be rejected from the model. We find that there is no correlation between the components that the GCs are likely to belong to and their colours. Particularly, for NGC 2768 we find that its red GCs display rotation preferentially at inner radii (Re<1). In the case of the GC system of NGC 3115 we find a group of GCs with similar kinematics that are not likely to belong to either its spheroid nor disc. For NGC 7457 we find that 70% of its GCs are likely to belong to the disc. Overall, our results suggest that these galaxies assembled into S0s through different evolutionary paths. Mergers seem to have been very important for NGC 2768 and NGC 3115 while NGC 7457 is more likely to have experienced secular evolution.

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

We aim to use multi-band imaging from the Phase-3 Verification Data of the J-PLUS survey to derive accurate photometric redshifts (photo-z) and look for potential new members in the surroundings of the nearby galaxy clusters A2589 (z=0.0414) & A2593 (z=0.0440), using redshift probability distribution functions. The ultimate goal is to demonstrate the usefulness of a 12-band filter system in the study of large-scale structure in the local universe. We present an optimized pipeline for the estimation of photo-z in clusters of galaxies. We tested our photo-z with a sample of 296 spectroscopically confirmed cluster members with a magnitude of <r>= 16.6 and redshift <z>=0.041. The combination of seven narrow and five broadband filters with a typical photometric-depth of r<21.5 provides dz/(1+z)=0.01 photo-z estimates. A precision of dz/(1+z)=0.005 is obtained for the 177 galaxies brighter than magnitude r<17. To foresee the precision beyond the spectroscopic sample, we designed a set of simulations in which real cluster galaxies are modeled and reinjected inside the images at different signal-to-noise. A precision of dz/(1+z)=0.02 and dz/(1+z)=0.03 is expected at <r>= 18-22, respectively. Complementarily, we used SDSS/DR12 data to derive photo-z estimates for the same galaxy sample, demonstrating that the wavelength-resolution of the J-PLUS can double the precision achieved by SDSS for galaxies with a high S/N. We find as much as 170 new candidates across the entire field. The spatial distribution of these galaxies may suggest an overlap between the systems with no evidence of a clear filamentary structure connecting the clusters. These preliminary results show the potential of J-PLUS data to revisit membership of groups and clusters from nearby galaxies, important for the determination of luminosity and mass functions and environmental studies at the intermediate and low-mass regime.

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.

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.

We conduct a comprehensive projected phase-space analysis of the A901/2 multi-cluster system at $z\sim0.165$. Aggregating redshifts from spectroscopy, tunable-filter imaging, and prism techniques, we assemble a sample of 856 cluster galaxies reaching $10^{8.5}M_\odot$ in stellar mass. We look for variations in cluster galaxy properties between virialised and non-virialised regions of projected phase-space (PPS). Our main conclusions point to relatively gentle environmental effects, expressed mainly on galaxy gas reservoirs. (1) Stacking the four subclusters in A901/2, we find galaxies in the virialised region are more massive, redder, and have marginally higher S\`ersic indices, but their half-light radii and Hubble types are not significantly different. (2) After accounting for trends in stellar mass, there is a remaining change in rest-frame colour across PPS. Primarily, the colour difference is due to an absence in the virialised region of galaxies with rest-frame $B-V<0.7$ and moderate-to-high ($M_\star>10^{9.85}M_\odot$) stellar mass. (3) There is an infalling population of lower-mass ($M_\star\leq10^{9.85}M_\odot$), relatively blue ($B-V<0.7$) elliptical or spheroidal galaxies that is strikingly absent in the virialised region. (4) The number of bona-fide star-forming and AGN galaxies in the PPS regions is strongly dictated by stellar mass. However, there remains a reduced fraction of star-forming galaxies in the centres of the clusters at fixed stellar mass, consistent with the star-formation-density relation in galaxy clusters. (5) There is no change in specific H$\alpha$-derived star-formation rates of star-forming galaxies at fixed mass across the cluster environment. This suggests that preprocessing of galaxies during infall plays a prominent role in quenching star formation.

We present a study of the star formation and AGN activity for galaxies in the Abell 901/2 multi-cluster system at z~0.167 as part of the OMEGA survey. Using Tuneable Filter data obtained with the OSIRIS instrument at the GTC we produce spectra covering the Halpha and [N II] spectral lines for more than 400 galaxies. Using optical emission-line diagnostics, we identify a significant number of galaxies hosting AGN, which tend to have high masses and a broad range of morphologies. Moreover, within the environmental densities probed by our study, we find no environmental dependence on the fraction of galaxies hosting AGN. The analysis of the integrated Halpha emission shows that the specific star formation rates (SSFRs) of a majority of the cluster galaxies are below the field values for a given stellar mass. We interpret this result as evidence for a slow decrease in the star formation activity of star-forming galaxies as they fall into higher-density regions, contrary to some previous studies which suggested a rapid truncation of star formation. We find that most of the intermediate- and high-mass spiral galaxies go through a phase in which their star formation is suppressed but still retain significant star-formation activity. During this phase, these galaxies tend to retain their spiral morphology while their colours become redder. The presence of this type of galaxies in high density regions indicates that the physical mechanism responsible for suppressing star-formation affects mainly the gas component of the galaxies, suggesting that ram-pressure stripping or starvation are potentially responsible.

We present new HST optical and near-infrared (NIR) photometry of the rich globular cluster (GC) system of NGC 4874, the cD galaxy in the core of the Coma cluster (Abell 1656). NGC 4874 was observed with the HST Advanced Camera for Surveys in the F475W (g) and F814W (I) passbands and the Wide Field Camera 3 IR Channel in F160W (H). The GCs in this field exhibit a bimodal optical color distribution with more than half of the GCs falling on the red side at g-I > 1. Bimodality is also present, though less conspicuously, in the optical-NIR I-H color. Consistent with past work, we find evidence for nonlinearity in the g-I versus I-H color-color relation. Our results thus underscore the need for understanding the detailed form of the color-metallicity relations in interpreting observational data on GC bimodality. We also find a very strong color-magnitude trend, or "blue tilt," for the blue component of the optical color distribution of the NGC 4874 GC system. A similarly strong trend is present for the overall mean I-H color as a function of magnitude; for M_814 < -10 mag, these trends imply a steep mass-metallicity scaling with $Z\propto M_{\rm GC}^{1.4\pm0.4}$, but the scaling is not a simple power law and becomes much weaker at lower masses. As in other similar systems, the spatial distribution of the blue GCs is more extended than that of the red GCs, partly because of blue GCs associated with surrounding cluster galaxies. In addition, the center of the GC system is displaced by 4+/-1 kpc towards the southwest from the luminosity center of NGC 4874, in the direction of NGC 4872. Finally, we remark on a dwarf elliptical galaxy with a noticeably asymmetrical GC distribution. Interestingly, this dwarf has a velocity of nearly -3000 km/s with respect to NGC 4874; we suggest it is on its first infall into the cluster core and is undergoing stripping of its GC system by the cluster potential.

Globular clusters (GCs) can be considered discrete, long-lived, dynamical tracers that retain crucial information about the assembly history of their parent galaxy. In this paper, we present a new catalogue of GC velocities and colours for the lenticular galaxy NGC 1023, we study their kinematics and spatial distribution, in comparison with the underlying stellar kinematics and surface brightness profile, and we test a new method for studying GC properties. Specifically, we decompose the galaxy light into its spheroid (assumed to represent the bulge + halo components) and disk components and use it to assign to each GC a probability of belonging to one of the two components. Then we model the galaxy kinematics, assuming a disk and spheroidal component, using planetary nebulae (PNe) and integrated stellar light. We use this kinematic model and the probability previously obtained from the photometry to recalculate for each GC its likelihood of being associated with the disk, the spheroid, or neither. We find that the reddest GCs are likely to be associated with the disk, as found for faint fuzzies in this same galaxy, suggesting that the disk of this S0 galaxy originated at z ~ 2. The majority of blue GCs are found likely to be associated with the spheroidal (hot) component. The method also allows us to identify objects that are unlikely to be in equilibrium with the system. In NGC1023 some of the rejected GCs form a substructure in phase space that is connected with NGC 1023 companion galaxy.

We describe techniques concerning wavelength calibration and sky subtraction to maximise the scientific utility of data from tunable filter instruments. While we specifically address data from the Optical System for Imaging and low Resolution Integrated Spectroscopy instrument (OSIRIS) on the 10.4~m Gran Telescopio Canarias telescope, our discussion is generalisable to data from other tunable filter instruments. A key aspect of our methodology is a coordinate transformation to polar coordinates, which simplifies matters when the tunable filter data is circularly symmetric around the optical centre. First, we present a method for rectifying inaccuracies in the wavelength calibration using OH sky emission rings. Using this technique, we improve the absolute wavelength calibration from an accuracy of 5 Angstroms to 1 Angstrom, equivalent to ~7% of our instrumental resolution, for 95% of our data. Then, we discuss a new way to estimate the background sky emission by median filtering in polar coordinates. This method suppresses contributions to the sky background from the outer envelopes of distant galaxies, maximising the fluxes of sources measured in the corresponding sky-subtracted images. We demonstrate for data tuned to a central wavelength of 7615~$\rm\AA$ that galaxy fluxes in the new sky-subtracted image are ~37% higher, versus a sky-subtracted image from existing methods for OSIRIS tunable filter data.

In this paper, the third in a series illustrating the power of generalized linear models (GLMs) for the astronomical community, we elucidate the potential of the class of GLMs which handles count data. The size of a galaxy's globular cluster population $N_{\rm GC}$ is a prolonged puzzle in the astronomical literature. It falls in the category of count data analysis, yet it is usually modelled as if it were a continuous response variable. We have developed a Bayesian negative binomial regression model to study the connection between $N_{\rm GC}$ and the following galaxy properties: central black hole mass, dynamical bulge mass, bulge velocity dispersion, and absolute visual magnitude. The methodology introduced herein naturally accounts for heteroscedasticity, intrinsic scatter, errors in measurements in both axes (either discrete or continuous), and allows modelling the population of globular clusters on their natural scale as a non-negative integer variable. Prediction intervals of 99% around the trend for expected $N_{\rm GC}$comfortably envelope the data, notably including the Milky Way, which has hitherto been considered a problematic outlier. Finally, we demonstrate how random intercept models can incorporate information of each particular galaxy morphological type. Bayesian variable selection methodology allows for automatically identifying galaxy types with different productions of GCs, suggesting that on average S0 galaxies have a GC population 35% smaller than other types with similar brightness.

Age, metallicity and spatial distribution of globular clusters (GCs) provide a powerful tool to reconstruct major star-formation episodes in galaxies. IKN is a faint dwarf spheroidal (dSph) in the M81 group of galaxies. It contains five old GCs, which makes it the galaxy with the highest known specific frequency (SN=126). We estimate the photometric age, metallicity and spatial distribution of the poorly studied IKN GCs. We search SDSS for GC candidates beyond the HST field of view, which covers half of IKN. To break the age-metallicity degeneracy in the V-I colour we use WHT/LIRIS Ks-band photometry and derive photometric ages and metallicities by comparison with SSP models in the V,I,Ks colour space. IKN GCs' VIKs colours are consistent with old ages ($\geq\!8$ Gyr) and a metallicity distribution with a higher mean than typical for such a dSph ([Fe/H$]\!\simeq\!-1.4_{-0.2}^{+0.6}$ dex). Their photometric masses range ($0.5 <{\cal M_{\rm GC}}<4\times10^5M_\odot$) implies a high mass ratio between GCs and field stars, of $10.6\%$. Mixture model analysis of the RGB field stars' metallicity suggests that 72\% of the stars may have formed together with the GCs. Using the most massive GC-SFR relation we calculate a SFR of $\sim\!10M_\odot/$yr during its formation epoch. We note that the more massive GCs are closer to the galaxy photometric centre. IKN GCs also appear spatially aligned along a line close to the IKN major-axis and nearly orthogonal to the plane of spatial distribution of galaxies in the M81 group. We identify one new IKN GC candidate based on colour and PSF analysis of the SDSS data. The evidence towards i) broad and high metallicity distribution of the field IKN RGB stars and its GCs, ii) high fraction and iii), spatial alignment of IKN GCs, supports a scenario for tidally triggered complex IKN's SFH in the context of interactions with galaxies in the M81 group.

We present an overview of and first results from the OMEGA survey: the OSIRIS Mapping of Emission-line Galaxies in the multi-cluster system A901/2. The ultimate goal of this project is to study star formation and AGN activity across a broad range of environments at a single redshift. Using the tuneable-filter mode of the OSIRIS instrument on GTC, we target Halpha and [NII] emission lines over a ~0.5 X 0.5 deg2 region containing the z~0.167 multi-cluster system A901/2. In this paper we describe the design of the survey, the observations and the data analysis techniques developed. We then present early results from two OSIRIS pointings centred on the cores of the A901a and A902 clusters. AGN and star-forming (SF) objects are identified using the [NII]/Halpha vs. W_Halpha (WHAN) diagnostic diagram. The AGN hosts are brighter, more massive, and possess earlier-type morphologies than SF galaxies. Both populations tend to be located towards the outskirts of the high density regions we study. The typical Halpha luminosity of these sources is significantly lower than that of field galaxies at similar redshifts, but greater than that found for A1689, a rich cluster at z~0.2. The Halpha luminosities of our objects translate into star-formation rates (SFRs) between ~0.02 and 6 Msun/yr. Comparing the relationship between stellar mass and Halpha-derived SFR with that found in the field indicates a suppression of star formation in the cores of the clusters. These findings agree with previous investigations of this multi-cluster structure, based on other star formation indicators, and demonstrate the power of tuneable filters for this kind of study.

The cores of clusters at 0 $\lesssim$ z $\lesssim$ 1 are dominated by quiescent early-type galaxies, whereas the field is dominated by star-forming late-type ones. Galaxy properties, notably the star formation (SF) ability, are altered as they fall into overdense regions. The critical issues to understand this evolution are how the truncation of SF is connected to the morphological transformation and the responsible physical mechanism. The GaLAxy Cluster Evolution Survey (GLACE) is conducting a study on the variation of galaxy properties (SF, AGN, morphology) as a function of environment in a representative sample of clusters. A deep survey of emission line galaxies (ELG) is being performed, mapping a set of optical lines ([OII], [OIII], H$\beta$ and H$\alpha$/[NII]) in several clusters at z $\sim$ 0.40, 0.63 and 0.86. Using the Tunable Filters (TF) of OSIRIS/GTC, GLACE applies the technique of TF tomography: for each line, a set of images at different wavelengths are taken through the TF, to cover a rest frame velocity range of several thousands km/s. The first GLACE results target the H$\alpha$/[NII] lines in the cluster ZwCl 0024.0+1652 at z = 0.395 covering $\sim$ 2 $\times$ r$_{vir}$. We discuss the techniques devised to process the TF tomography observations to generate the catalogue of H$\alpha$ emitters of 174 unique cluster sources down to a SFR below 1 M$_{\odot}$/yr. The AGN population is discriminated using different diagnostics and found to be $\sim$ 37% of the ELG population. The median SFR is 1.4 M$_{\odot}$/yr. We have studied the spatial distribution of ELG, confirming the existence of two components in the redshift space. Finally, we have exploited the outstanding spectral resolution of the TF to estimate the cluster mass from ELG dynamics, finding M$_{200}$ = 4.1 $\times$ 10$^{14}$ M$_{\odot} h^{-1}$, in agreement with previous weak-lensing estimates.

This work presents the first results from an ESO Large Programme carried out using the OSIRIS instrument on the 10m GTC telescope (La Palma). We have observed a large sample of galaxies in the region of the Abell 901/902 system (z ~ 0.165) which has been extensively studied as part of the STAGES project. We have obtained spectrally and spatially resolved H-alpha and [NII] emission maps for a very large sample of galaxies covering a broad range of environments. The new data are combined with extensive multi-wavelength observations which include HST, COMBO-17, Spitzer, Galex and XMM imaging to study star formation and AGN activity as a function of environment and galaxy properties such as luminosity, mass and morphology. The ultimate goal is to understand, in detail, the effect of the environment on star formation and AGN activity.

We combined new near-infrared VLT/HAWK-I data of the globular clusters (GCs) in the isolated edge-on S0 galaxy NGC 3115 with optical and spectroscopic ones taken from the literature, with the aim of analyzing the multiband GC color distributions. A recent study from the SLUGGS survey has shown that the GCs in this galaxy follow a bimodal distribution of Ca II triplet indices. Thus, NGC 3115 presents a critical example of a GC system with multiple, distinct, metallicity subpopulations, and this may argue against the "projection" scenario, which posits that the ubiquitous color bimodality mainly results from nonlinearities in the color-metallicity relations. Using optical, NIR, and spectroscopic data, we found strong and consistent evidence of index bimodality, which independently confirms the metallicity bimodality in NGC 3115 GCs. At the same time, we also found evidence for some color-color nonlinearity. Taken in the broader context of previous studies, the multicolor consistency of the GC bimodality in NGC 3115 suggests that in cases where GC systems exhibit clear differences between their optical and optical-NIR color distributions (as in some giant ellipticals), the apparent inconsistencies most likely result from nonlinearities in the color-metallicity relations.

Faint fuzzies are metal-rich apparently-old star clusters with unusually large radii (7-15 pc), found mostly in S0 galaxies, whose source remain obscure. To identify their origins, we compare planetary nebulae and neutral hydrogen with faint fuzzy positions and line-of-sight velocities in NGC1023. In this way, we rule out scenarios in which these objects are associated with an on-going merger or with a spheroid population in NGC1023. Their kinematics are indistinguishable from the stellar disk population in this galaxy, and we conclude that faint fuzzies are most likely just remnant open clusters. Their observed association with S0s then simply reflects the difficulty of identifying such objects in later-type disk galaxies.

The stellar kinematics of the spheroids and discs of S0 galaxies contain clues to their formation histories. Unfortunately, it is difficult to disentangle the two components and to recover their stellar kinematics in the faint outer parts of the galaxies using conventional absorption line spectroscopy. This paper therefore presents the stellar kinematics of six S0 galaxies derived from observations of planetary nebulae (PNe), obtained using the Planetary Nebula Spectrograph. To separate the kinematics of the two components, we use a maximum-likelihood method that combines the discrete kinematic data with a photometric component decomposition. The results of this analysis reveal that: the discs of S0 galaxies are rotationally supported; however, the amount of random motion in these discs is systematically higher than in comparable spiral galaxies; and the S0s lie around one magnitude below the Tully--Fisher relation for spiral galaxies, while their spheroids lie nearly one magnitude above the Faber--Jackson relation for ellipticals. All of these findings are consistent with a scenario in which spirals are converted into S0s through a process of mild harassment or "pestering," with their discs somewhat heated and their spheroid somewhat enhanced by the conversion process. In such a scenario, one might expect the properties of S0s to depend on environment. We do not see such an effect in this fairly small sample, although any differences would be diluted by the fact that the current location does not necessarily reflect the environment in which the transformation occurred. Similar observations of larger samples probing a broader range of environments, coupled with more detailed modelling of the transformation process to match the wide range of parameters that we have shown can now be measured, should take us from these first steps to the definitive answer as to how S0 galaxies form.

The origins of S0 galaxies remain obscure, with various mechanisms proposed for their formation, likely depending on environment. These mechanisms would imprint different signatures in the galaxies' stellar kinematics out to large radii, offering a method for distinguishing between them. We aim to study a sample of six S0 galaxies from a range of environments, and use planetary nebulae (PNe) as tracers of their stellar populations out to very large radii, to determine their kinematics in order to understand their origins. Using a special-purpose instrument, the Planetary Nebula Spectrograph, we observe and extract PNe catalogues for these six systems*. We show that the PNe have the same spatial distribution as the starlight, that the numbers of them are consistent with what would be expected in a comparable old stellar population in elliptical galaxies, and that their kinematics join smoothly onto those derived at smaller radii from conventional spectroscopy. The high-quality kinematic observations presented here form an excellent set for studying the detailed kinematics of S0 galaxies, in order to unravel their formation histories. We find that PNe are good tracers of stellar kinematics in these systems. We show that the recovered kinematics are largely dominated by rotational motion, although with significant random velocities in most cases.

We present new spectroscopic observations of a sample of globular clusters (GCs) in the Virgo Cluster elliptical NGC4365. The great majority of the objects analysed here, overlap on purpose with objects presented in two previous studies that came to contradictory conclusions. We measure Lick indices for the GCs and infer ages, abundances and alpha-enhancement values by comparing to population synthesis models. Our results do not support the existence of numerous intermediate-age GCs in this galaxy as suggested by some of the past studies. Instead, we find ages, consistent with the range ~8-14 Gyrs (model dependent and with large uncertainties). We also find that the newly observed GCs have [Fe/H]\sim-1.3 to 0.3 dex and the majority are consistent with [alpha/Fe]\sim0.3 dex, independent of metallicity. By comparing the metal sensitive index [MgFe] with optical colours, we find that our data is consistent with 14Gyr models. However for optical/near-infrared colours we find that the same models would have to be redder by ~0.2 mag in order to be consistent with our data. This offset is not found to be consistent with an age effect as it is even larger for younger models.

We present recent polarimetric images of the highly variable star R CrB using ExPo and archival WFPC2 images from the HST. We observed R CrB during its current dramatic minimum where it decreased more than 9 mag due to the formation of an obscuring dust cloud. Since the dust cloud is only in the line-of-sight, it mimics a coronograph allowing the imaging of the star's circumstellar environment. Our polarimetric observations surprisingly show another scattering dust cloud at approximately 1.3" or 2000 AU from the star. We find that to obtain a decrease in the stellar light of 9 mag and with 30% of the light being reemitted at infrared wavelengths (from R CrB's SED) the grains in R CrB's circumstellar environment must have a very low albedo of approximately 0.07%. We show that the properties of the dust clouds formed around R CrB are best fitted using a combination of two distinct populations of grains size. The first are the extremely small 5 nm grains, formed in the low density continuous wind, and the second population of large grains (~0.14 {\mu}m) which are found in the ejected dust clouds. The observed scattering cloud, not only contains such large grains, but is exceptionally massive compared to the average cloud.

The interpretation that bimodal colour distributions of globular clusters (GCs) reflect bimodal metallicity distributions has been challenged. Non-linearities in the colour to metallicity conversions caused by the horizontal branch (HB) stars may be responsible for transforming a unimodal metallicity distribution into a bimodal (optical) colour distribution. We study optical/near-infrared (NIR) colour distributions of the GC systems in 14 E/S0 galaxies. We test whether the bimodal feature, generally present in optical colour distributions, remains in the optical/NIR ones. The latter colour combination is a better metallicity proxy than the former. We use KMM and GMM tests to quantify the probability that different colour distributions are better described by a bimodal, as opposed to a unimodal distribution. We find that double-peaked colour distributions are more commonly seen in optical than in optical/NIR colours. For some of the galaxies where the optical (g-z) distribution is clearly bimodal, the (g-K) and (z-K) distributions are better described by a unimodal distribution. The two most cluster-rich galaxies in our sample, NGC4486 and NGC4649, show some interesting differences. The (g-K) distribution of NGC4649 is better described by a bimodal distribution, while this is true for the (g-K) distribution of NGC4486 GCs only if restricted to a brighter sub-sample with small K-band errors (< 0.05 mag). Formally, the K-band photometric errors cannot be responsible for blurring bimodal metallicity distributions to unimodal (g-K) colour distributions. However, simulations including the extra scatter in the colour-colour diagrams (not fully accounted for in the photometric errors) show that such scatter may contribute to the disappearance of bimodality in (g-K) for the full NGC4486 sample. For the less cluster-rich galaxies results are inconclusive due to poorer statistics. [Abridged]

Context: Some photometric studies of extragalactic globular cluster (GC) systems using the optical and near-infrared colour combination have suggested the presence of a large fraction of intermediate-age (2-8 Gyrs) GCs. Aims: We investigate the age distributions of GC systems in 14 E/S0 galaxies. Methods: We carry out a differential comparison of the (g-z) vs. (g-K) two-colour diagrams for GC systems in the different galaxies in order to see whether there are indications of age differences. We also compare the different GC systems with a few simple stellar population models. Results: No significant difference is detected in the mean ages of GCs among elliptical galaxies. S0 galaxies on the other hand, show evidence for younger GCs. Surprisingly, this appears to be driven by the more metal-poor clusters. The age distribution of GCs in NGC4365 seems to be similar to that of other large ellipticals (e.g. NGC4486, NGC4649). Padova SSPs with recently released isochrones for old ages (14 Gyrs) show less of an offset with respect to the photometry than previously published models. Conclusions: We suggest that E type galaxies assembled most of their GCs in a shorter and earlier period than S0 type galaxies. The latter galaxy type, seems to have a more extended period of GC formation/assembly.

Context: The combination of optical and near-infrared (NIR) colours has the potential to break the age/metallicity degeneracy and offers a better metallicity sensitivity than optical colours alone. Previous studies of extragalactic globular clusters (GCs) with this colour combination, however, have suffered from small samples or have been restricted to a few galaxies. Aims: We compile a homogeneous and representative sample of GC systems with multi-band photometry to be used in subsequent papers where ages and metallicity distributions will be studied. Methods: We acquired deep K-band images of 14 bright nearby early-type galaxies. The images were obtained with the LIRIS near-infrared spectrograph and imager at the William Herschel Telescope (WHT) and combined with optical ACS g and z images from the Hubble Space Telescope public archive. Results: For the first time GC photometry of 14 galaxies are observed and reduced homogeneously in this wavelength regime. We achieved a limiting magnitude of K~20-21. For the majority of the galaxies we detect about 70 GCs each. NGC4486 and NGC4649, the cluster-richest galaxies in the sample contain 301 and 167 GCs, respectively. We present tables containing coordinates, photometry and sizes of the GCs available.

Optical/near-infrared observations for 14 globular cluster (GC) systems in early- type galaxies are presented. We investigate the recent claims (Yoon, Yi & Lee 2006) of colour bimodality in GC systems being an artefact of the non linear colour - metallicity transformation driven by the horizontal branch morphology. Taking the advantage of the fact that the combination of optical and near-infrared colours can in principle break the age/metallicity degeneracy we also analyse age distributions in these systems.

NGC 1380 is a lenticular galaxy located near the centre of the Fornax Cluster northeast of NGC 1399. The globular cluster system of this galaxy was previously studied only from the ground. Recent studies of similar early-type galaxies, specially lenticular ones, reveal the existence of star clusters that apparently break up the traditional open/globular cluster dichotomy. With higher quality photometry from HST/WFPC2 we study the star clusters in NGC 1380, measuring their magnitudes, colours, sizes and projected distances from the centre of the galaxy. We used deep archival HST/WFPC2 in the B and V bands. We built colour magnitude diagrams from which we selected a sample of cluster candidates. We also analysed their colour distribution and measured their sizes. Based on their location in the luminosity-size diagram we estimated probabilities of them being typical globular clusters as those found in the Galaxy. A total of about 570 cluster candidates were found down to V=26.5. We measured sizes for approximately 200 of them. The observed colour distribution has three apparent peaks. Likewise for the size distribution. We identified the smaller population as being mainly typical globular clusters, while the more extended objects have small probabilities of being such objects. Different correlations between absolute magnitudes, sizes, colours and location were inferred for these cluster sub-populations. Most extended clusters (Reff > 4 pc) share similar properties to the diffuse star clusters reported to inhabit luminous early-type galaxies in the Virgo galaxy cluster such as being of low surface brightness and fainter than MV ~ -8. We also report on a small group of (Reff ~ 10 pc), -8< MV < -6, red clusters located near the centre of NGC 1380, which may be interpreted as faint fuzzies.

NGC 5846 is a giant elliptical galaxy with a previously well studied globular cluster system (GCS), known to have a bimodal colour distribution with a remarkably high red fraction. Here we revisit the central galaxy regions searching for new globular cluster (GC) candidates, and measuring, magnitudes, colours and sizes for them. We also search for their X-ray counterparts. We use archival Hubble Space Telescope WFPC2 images, from which we modelled and subtracted the host light distribution and increased the available sample of GCs. We performed photometry on the central objects, and measured sizes and equatorial coordinates for the entire system known in this galaxy. We detect two dozen previously unknown GC candidates in the central regions. Reliable sizes are obtained for about 60 GCs; their typical effective radii are in the range 3-5 pc. The largest clusters are located in the central regions. We find 7 X-ray counterparts to globular clusters, most of them in the central region. They are among the most luminous X-ray sources in NGC 5846. They are also optically luminous, compact and belong to the red subpopulation.