Cosmological N-body simulations represent an excellent tool to study the formation and evolution of dark matter (DM) halos and the mechanisms that have originated the universal profile at the largest mass scales in the Universe. In particular, the combination of the velocity dispersion $\sigma_\mathrm{v}$ with the density $\rho$ can be used to define the pseudo-entropy $S(r)=\sigma_\mathrm{v}^2/\rho^{2/3}$, whose profile is well-described by a simple power-law $S\propto r^{\alpha}$. We analyze a set of cosmological hydrodynamical re-simulations of massive galaxy clusters and study the pseudo-entropy profiles as traced by different collisionless components in simulated galaxy clusters: DM, stars, and substructures. We analyze four sets of simulations, exploring different resolution and physics (N-body and full hydrodynamical simulations) to investigate convergence and the impact of baryons. We find that baryons significantly affect the inner region of pseudo-entropy profile as traced by substructures, while pseudo-entropy profiles as traced by DM particles are characterized by an almost universal behavior and can be used as a low-scatter ($\sim6-7\%$) mass-proxy. We compare observed and simulated pseudo-entropy profiles and find good agreement in both normalization and slope. We demonstrate, however, that the method used to derive observed pseudo-entropy profiles could introduce biases and underestimate the impact of mergers. Finally, we investigate the pseudo-entropy traced by the stars focusing our interest in the dynamical distinction between intracluster light (ICL) and the stars bound to the brightest cluster galaxy (BCG): the combination of these two pseudo-entropy profiles is well-described by a single power-law out to almost the entire cluster virial radius.

We present weak lensing shear catalogues from the fourth data release of the Kilo-Degree Survey, KiDS-1000, spanning 1006 square degrees of deep and high-resolution imaging. Our `gold-sample' of galaxies, with well calibrated photometric redshift distributions, consists of 21 million galaxies with an effective number density of $\sim 6.2$ galaxies per square arcminute. We quantify the accuracy of the spatial, temporal and flux-dependent point-spread function (PSF) model, verifying that the model meets our requirements to induce less than a $0.1\sigma$ change in the inferred cosmic shear constraints on the clustering cosmological parameter $S_8 = \sigma_8\sqrt{\Omega_{\rm m}/0.3}$. Through a series of two-point null-tests we validate the shear estimates, finding no evidence for significant non-lensing B-mode distortions in the data. PSF residuals are detected in the highest-redshift bins, originating from object selection and/or weight bias. The amplitude is however shown to be sufficiently low and within our stringent requirements. With a shear-ratio null-test we verify the expected redshift scaling of the galaxy-galaxy lensing signal around luminous red galaxies. We conclude that the joint KiDS-1000 shear and photometric redshift calibration is sufficiently robust for combined-probe gravitational lensing and spectroscopic clustering analyses.

The Probability Density Function (PDF) provides an estimate of the photometric redshift (zphot) prediction error. It is crucial for current and future sky surveys, characterized by strict requirements on the zphot precision, reliability and completeness. The present work stands on the assumption that properly defined rejection criteria, capable of identifying and rejecting potential outliers, can increase the precision of zphot estimates and of their cumulative PDF, without sacrificing much in terms of completeness of the sample. We provide a way to assess rejection through proper cuts on the shape descriptors of a PDF, such as the width and the height of the maximum PDF's peak. In this work we tested these rejection criteria to galaxies with photometry extracted from the Kilo Degree Survey (KiDS) ESO Data Release 4, proving that such approach could lead to significant improvements to the zphot quality: e.g., for the clipped sample showing the best trade-off between precision and completeness, we achieve a reduction in outliers fraction of $\simeq 75\%$ and an improvement of $\simeq 6\%$ for NMAD, with respect to the original data set, preserving the $\simeq 93\%$ of its content.

We present the methodology for a joint cosmological analysis of weak gravitational lensing from the fourth data release of the ESO Kilo-Degree Survey (KiDS-1000) and galaxy clustering from the partially overlapping BOSS and 2dFLenS surveys. Cross-correlations between galaxy positions and ellipticities are incorporated into the analysis, developing a hybrid model of non-linear scales that blends perturbative and non-perturbative approaches, and assessing signal contributions by astrophysical effects. All weak lensing signals are measured consistently via Fourier-space statistics that are insensitive to the survey mask and display low levels of mode mixing. The calibration of photometric redshift distributions and multiplicative gravitational shear bias is updated, and a more complete tally of residual calibration uncertainties is propagated into the likelihood. A dedicated suite of more than 20000 mocks is used to assess the performance of covariance models and to quantify the impact of survey geometry and spatial variations of survey depth on signals and their errors. The sampling distributions for the likelihood and the $\chi^2$ goodness-of-fit statistic are validated, with proposed changes to calculating the effective number of degrees of freedom. Standard weak lensing point estimates on $S_8=\sigma_8\,(\Omega_{\rm m}/0.3)^{1/2}$ derived from its marginal posterior are easily misinterpreted to be biased low, and an alternative estimator and associated credible interval are proposed. Known systematic effects pertaining to weak lensing modelling and inference are shown to bias $S_8$ by no more than 0.1 standard deviations, with the caveat that no conclusive validation data exist for models of intrinsic galaxy alignments. Compared to the previous KiDS analyses, $S_8$ constraints are expected to improve by 20% for weak lensing alone and by 29% for the joint analysis. [abridged]

Relic galaxies are thought to be the progenitors of high-redshift red nuggets that for some reason missed the channels of size growth and evolved passively and undisturbed since the first star formation burst (at $z>2$). These local ultracompact old galaxies are unique laboratories for studying the star formation processes at high redshift and thus the early stage of galaxy formation scenarios. Counterintuitively, theoretical and observational studies indicate that relics are more common in denser environments, where merging events predominate. To verify this scenario, we compared the number counts of a sample of ultracompact massive galaxies (UCMGs) selected within the third data release of the Kilo Degree Survey, that is, systems with sizes $R_{\rm e} < 1.5 \, \rm kpc$ and stellar masses $M_{\rm \star} > 8 \times 10^{10}\, \rm M_{\odot}$, with the number counts of galaxies with the same masses but normal sizes in field and cluster environments. Based on their optical and near-infrared colors, these UCMGs are likely to be mainly old, and hence representative of the relic population. We find that both UCMGs and normal-size galaxies are more abundant in clusters and their relative fraction depends only mildly on the global environment, with denser environments penalizing the survival of relics. Hence, UCMGs (and likely relics overall) are not special because of the environment effect on their nurture, but rather they are just a product of the stochasticity of the merging processes regardless of the global environment in which they live.

We present the study of the south-west group in the Fornax cluster centred on the brightest group galaxy (BGG) Fornax A, observed as part of the Fornax Deep Survey (FDS). This includes the analysis of the bright group members (mB < 16 mag) and the intra-group light (IGL). The main objective of this work is to investigate the assembly history of the Fornax A group and to compare its physical quantities as a function of the environment to that of the Fornax cluster core. For all galaxies, we extract the azimuthally averaged surface brightness profiles in three optical bands (g, r, i) by modelling the galaxy's isophotes. We derive their colour profiles and structural parameters in all respective bands. The long integration time and large covered area of the FDS allow us to also estimate the amount of IGL. The majority of galaxies in the Fornax A group are late-type galaxies (LTGs), spanning a range of stellar mass of $8 < log (M_* M_{\odot}) < 10.5$. Six out of nine LTGs show a Type III (up-bending) break in their light profiles, which is either suggestive of strangulation halting star-formation in their outskirts or their HI-richness causing enhanced star-formation in their outer-discs. The estimated luminosity of the IGL is $6 \pm 2 \times 10^{10} L_{\odot}$ in g-band, which corresponds to about 16% of the total light in the group. The Fornax A group appears to be in an early stage of assembly with respect to the cluster core. The environment of the Fornax A group is not as dense as that of the cluster core, with all galaxies except the BGG showing similar morphology, comparable colours and stellar masses, and Type III disc-breaks, without any clear trend of these properties with group-centric distances. The main contribution to the IGL is from the minor merging in the outskirts of the BGG NGC1316 and, probably, the disrupted dwarf galaxies close to the group centre.

We report new high-quality galaxy scale strong lens candidates found in the Kilo Degree Survey data release 4 using Machine Learning. We have developed a new Convolutional Neural Network (CNN) classifier to search for gravitational arcs, following the prescription by \cite{2019MNRAS.484.3879P} and using only $r-$band images. We have applied the CNN to two "predictive samples": a Luminous red galaxy (LRG) and a "bright galaxy" (BG) sample ($r<21$). We have found 286 new high probability candidates, 133 from the LRG sample and 153 from the BG sample. We have then ranked these candidates based on a value that combines the CNN likelihood to be a lens and the human score resulting from visual inspection (P-value) and we present here the highest 82 ranked candidates with P-values $\ge 0.5$. All these high-quality candidates have obvious arc or point-like features around the central red defector. Moreover, we define the best 26 objects, all with scores P-values $\ge 0.7$ as a "golden sample" of candidates. This sample is expected to contain very few false positives and thus it is suitable for follow-up observations. The new lens candidates come partially from the the more extended footprint adopted here with respect to the previous analyses, partially from a larger predictive sample (also including the BG sample). These results show that machine learning tools are very promising to find strong lenses in large surveys and more candidates that can be found by enlarging the predictive samples beyond the standard assumption of LRGs. In the future, we plan to apply our CNN to the data from next-generation surveys such as the Large Synoptic Survey Telescope, Euclid, and the Chinese Space Station Optical Survey.

Several semi-analytic models (SAMs) try to explain how galaxies form, evolve and interact inside the dark matter large-scale structure. These SAMs can be tested by comparing their predictions for galaxy-galaxy-galaxy-lensing (G3L), which is weak gravitational lensing around galaxy pairs, with observations. We evaluate the SAMs by Henriques et al. (2015; H15) and by Lagos et al. (2012; L12), implemented in the Millennium Run, by comparing their predictions for G3L to observations at smaller scales than previous studies and also for pairs of lens galaxies from different populations. We compare the G3L signal predicted by the SAMs to measurements in the overlap of the Galaxy And Mass Assembly survey (GAMA), the Kilo-Degree Survey (KiDS), and the VISTA Kilo-degree Infrared Galaxy survey (VIKING), splitting lens galaxies into two colour and five stellar-mass samples. Using an improved G3L estimator, we measure the three-point correlation of the matter distribution for mixed lens pairs with galaxies from different samples, and unmixed lens pairs with galaxies from the same sample. Predictions by the H15 SAM agree with the observations for all colour-selected and all but one stellar-mass-selected sample with 95% confidence. Deviations occur for lenses with stellar masses below $9.5h^{-2}\mathrm{M}_\odot$ at scales below $0.2h^{-1}\mathrm{Mpc}$. Predictions by the L12 SAM for stellar-mass selected samples and red galaxies are significantly higher than observed, while the predicted signal for blue galaxy pairs is too low. The L12 SAM predicts more pairs of small stellar-mass and red galaxies than the H15 SAM and the observations, as well as fewer pairs of blue galaxies. Likely explanations are different treatments of environmental effects by the SAMs and different models of the initial mass function. We conclude that G3L provides a stringent test for models of galaxy formation and evolution.

This paper continues the series of the Fornax Deep Survey (FDS). Following the previous studies dedicated to extended Fornax cluster members, we present the catalogs of compact stellar systems in the Fornax cluster as well as extended background sources and point-like sources. We derive ugri photometry of ~1.7 million sources over the $\sim$21 sq. degree area of FDS centered on NGC1399. For a wider area, of $\sim$27 sq. degs extending in the direction of NGC1316, we provide gri data for ~3.1 million sources. To improve the morphological characterization of sources we generate multi-band image stacks by coadding the best seeing gri-band single exposures with a cut at FWHM<=0.9 arcsec. We use the multi-band stacks as detection frames. The identification of compact sources is obtained from a combination of photometric and morphometric selection criteria taking as reference the properties of sources with well-defined classification from the literature. We present a preliminary analysis of globular cluster (GC) distributions in the Fornax area. The study confirms and extends further previous results. We observe the inter-galactic population of GCs, a population of mainly blue GCs centered on NGC1399, extends over $\sim$0.9Mpc, with an ellipticity $\sim$0.65. Several sub-structures extend over $\sim$0.5Mpc along various directions. Two of these structures do not cross any bright galaxy; one of them appears to be connected to NGC1404, a bright galaxy close to the cluster core and particularly poor of GCs. Using the gri catalogs we analyze the GC distribution over the extended FDS area, and do not find any obvious GC sub-structure bridging the two brightest cluster galaxies, NGC1316 and NGC1399. Although NGC1316 is twice brighter of NGC1399 in optical bands we estimate a factor of 3-4 richer GC population around NGC1399 compared to NGC1316, out to galactocentric distances of 40 arcmin

This work is based on deep multi-band (g, r, i) data from the Fornax Deep Survey with VST. We analyse the surface brightness profiles of the 19 bright ETGs inside the virial radius of the Fornax cluster. The main aim of this work is to identify signatures of accretion onto galaxies by studying the presence of outer stellar halos, and understand their nature and occurrence. Our analysis also provides a new and accurate estimate of the intra-cluster light inside the virial radius of Fornax. We performed multi-component fits to the azimuthally averaged surface brightness profiles available for all sample galaxies. This allows to quantify the relative weight of all components in the galaxy structure that contribute to the total light. In addition, we derived the average g-i colours in each component identified by the fit, as well as the azimuthally averaged g-i colour profiles, to correlate them with the stellar mass of each galaxy and the location inside the cluster. We find that in the most massive and reddest ETGs the fraction of light in, probably accreted, halos is much larger than in the other galaxies. Less-massive galaxies have an accreted mass fraction lower than 30%, bluer colours and reside in the low-density regions of the cluster. Inside the virial radius of the cluster, the total luminosity of the intra-cluster light, compared with the total luminosity of all cluster members, is about 34%. Inside the Fornax cluster there is a clear correlation between the amount of accreted material in the stellar halos of galaxies and the density of the environment in which those galaxies reside. By comparing this quantity with theoretical predictions and previous observational estimates, there is a clear indication that the driving factor for the accretion process is the total stellar mass of the galaxy, in agreement with the hierarchical accretion scenario.

Ultra-compact massive galaxies UCMGs, i.e. galaxies with stellar masses $M_{*} > 8 \times 10^{10} M_{\odot}$ and effective radii $R_{e} < 1.5$ kpc, are very rare systems, in particular at low and intermediate redshifts. Their origin as well as their number density across cosmic time are still under scrutiny, especially because of the paucity of spectroscopically confirmed samples. We have started a systematic census of UCMG candidates within the ESO Kilo Degree Survey, together with a large spectroscopic follow-up campaign to build the largest possible sample of confirmed UCMGs. This is the third paper of the series and the second based on the spectroscopic follow-up program. Here, we present photometrical and structural parameters of 33 new candidates at redshifts $0.15 \lesssim z \lesssim 0.5$ and confirm 19 of them as UCMGs, based on their nominal spectroscopically inferred $M_{*}$ and $R_{e}$. This corresponds to a success rate of $\sim 58\%$, nicely consistent with our previous findings. The addition of these 19 newly confirmed objects, allows us to fully assess the systematics on the system selection, and finally reduce the number density uncertainties. Moreover, putting together the results from our current and past observational campaigns and some literature data, we build the largest sample of UCMGs ever collected, comprising 92 spectroscopically confirmed objects at $0.1 \lesssim z \lesssim 0.5$. This number raises to 116, allowing for a $3\sigma$ tolerance on the $M_{*}$ and $R_{e}$ thresholds for the UCMG definition. For all these galaxies we have estimated the velocity dispersion values at the effective radii which have been used to derive a preliminary mass-velocity dispersion correlation.

The VST Optical Imaging of the CDFS and ES1 Fields (VOICE) Survey, in synergy with the SUDARE survey, is a deep optical $ugri$ imaging of the CDFS and ES1 fields using the VLT Survey Telescope (VST). The observations for the CDFS field comprise about 4.38 deg$^2$ down to $r\sim26$ mag. The total on-sky time spans over four years in this field, distributed over four adjacent sub-fields. In this paper, we use the multi-epoch $r$-band imaging data to measure the variability of the detected objects and search for transients. We perform careful astrometric and photometric calibrations and point spread function (PSF) modeling. A new method, referring to as differential running-average photometry, is proposed to measure the light curves of the detected objects. With the method, the difference of PSFs between different epochs can be reduced, and the background fluctuations are also suppressed. Detailed uncertainty analysis and detrending corrections on the light curves are performed. We visually inspect the light curves to select variable objects, and present some objects with interesting light curves. Further investigation of these objects in combination with multi-band data will be presented in our forthcoming paper.

Variability has proven to be a powerful tool to detect active galactic nuclei (AGN) in multi-epoch surveys. The new-generation facilities expected to become operational in the next few years will mark a new era in time-domain astronomy and their wide-field multi-epoch campaigns will favor extensive variability studies. We present our analysis of AGN variability in the second half of the VST survey of the Wide Chandra Deep Field South (W-CDF-S), performed in the r band and covering a 2 sq. deg. area. The analysis complements a previous work, in which the first half of the area was investigated. We provide a reliable catalog of variable AGN candidates, which will be critical targets in future variability studies. We selected a sample of optically variable sources and made use of infrared data from the Spitzer mission to validate their nature by means of color-based diagnostics. We obtain a sample of 782 AGN candidates among which 12 are classified as supernovae, 54 as stars, and 232 as AGN. We estimate a contamination $\lesssim 20\%$ and a completeness $\sim 38\%$ with respect to mid-infrared selected samples.

In the last years, Astroinformatics has become a well defined paradigm for many fields of Astronomy. In this work we demonstrate the potential of a multidisciplinary approach to identify globular clusters (GCs) in the Fornax cluster of galaxies taking advantage of multi-band photometry produced by the VLT Survey Telescope using automatic self-adaptive methodologies. The data analyzed in this work consist of deep, multi-band, partially overlapping images centered on the core of the Fornax cluster. In this work we use a Neural-Gas model, a pure clustering machine learning methodology, to approach the GC detection, while a novel feature selection method ($\Phi$LAB) is exploited to perform the parameter space analysis and optimization. We demonstrate that the use of an Astroinformatics based methodology is able to provide GC samples that are comparable, in terms of purity and completeness with those obtained using single band HST data (Brescia et al. 2012) and two approaches based respectively on a morpho-photometric (Cantiello et al. 2018b) and a PCA analysis (D'Abrusco et al. 2015) using the same data discussed in this work.

The formation of Low mass X-ray binaries (LMXB) is favored within dense stellar systems such as Globular Clusters (GCs). The connection between LMXB and Globular Clusters has been extensively studied in the literature, but these studies have always been restricted to the innermost regions of galaxies. We present a study of LMXB in GCs within the central 1.5 deg^2 of the Fornax cluster with the aim of confirming the existence of a population of LMXB in intra-cluster GCs and understand if their properties are related to the host GCs, to the environment or/and to different formation channels.

We analyse the mass density distribution in the centres of galaxies across five orders of magnitude in mass range. Using high-quality spiral galaxy rotation curves and infrared photometry from SPARC, we conduct a systematic study of their central dark matter fraction ($f_{\rm DM}$) and their mass density slope ($\alpha$), within their effective radius. We show that lower-mass spiral galaxies are more dark matter dominated and have more shallow mass density slopes when compared with more massive galaxies, which have density profiles closer to isothermal. Low-mass ($M_{\star} \lesssim 10^{10}\, \rm M_{\odot}$) gas-rich spirals span a wide range of \fdm\ values, but systematically lower than in gas-poor systems of similar mass. With increasing galaxy mass, the values of \fdm\ decrease and the density profiles steepen. In the most massive late-type gas-poor galaxies, a possible flattening of these trends is observed. When comparing these results to massive ($M_{\star} \gtrsim 10^{10}M_{\odot}$) elliptical galaxies from SPIDER and to dwarf ellipticals from SMACKED, these trends result to be inverted. Hence, the values of both $f_{\rm DM}$ and $\alpha$, as a function of $M_{\star}$, exhibit a U-shape trend. At a fixed stellar mass, the mass density profiles in dwarf ellipticals are steeper than in spirals. These trends can be understood by stellar feedback from a more prolonged star formation period in spirals, causing a transformation of the initial steep density cusp to a more shallow profile via differential feedback efficiency by supernovae, and by galaxy mergers or AGN feedback in higher-mass galaxies.

We present the curation and verification of a new combined optical and near infrared dataset for cosmology and astrophysics, derived from the combination of $ugri$-band imaging from the Kilo Degree Survey (KiDS) and $ZY\!J\!H\!K_{\rm s}$-band imaging from the VISTA Kilo degree Infrared Galaxy (VIKING) survey. This dataset is unrivaled in cosmological imaging surveys due to its combination of area ($458$ deg$^2$ before masking), depth ($r\le25$), and wavelength coverage ($ugriZY\!J\!H\!K_{\rm s}$). The combination of survey depth, area, and (most importantly) wavelength coverage allows significant reductions in systematic uncertainties (i.e. reductions of between 10 and 60\% in bias, outlier rate, and scatter) in photometric-to-spectroscopic redshift comparisons, compared to the optical-only case at photo-$z$ above $0.7$. The complementarity between our optical and NIR surveys means that over $80\%$ of our sources, across all photo-$z$, have significant detections (i.e. not upper limits) in our $8$ reddest bands. We derive photometry, photo-$z$, and stellar masses for all sources in the survey, and verify these data products against existing spectroscopic galaxy samples. We demonstrate the fidelity of our higher-level data products by constructing the survey stellar mass functions in 8 volume-complete redshift bins. We find that these photometrically derived mass functions provide excellent agreement with previous mass evolution studies derived using spectroscopic surveys. The primary data products presented in this paper are publicly available at http://kids.strw.leidenuniv.nl/.

Intrinsic galaxy alignments are a source of bias for weak lensing measurements as well as a tool for understanding galaxy formation and evolution. In this work, we measure the alignment of shapes of satellite galaxies, in galaxy groups, with respect to the brightest group galaxy (BGG), as well as alignments of the BGG shape with the satellite positions, using the highly complete Galaxy And Mass Assembly (GAMA) spectroscopic survey and deep imaging from the Kilo Degree Survey. We control systematic errors with dedicated image simulations and measure accurate shapes using the DEIMOS shape measurement method. We find a significant satellite radial alignment signal, which vanishes at large separations from the BGG. We do not identify any strong trends of the signal with galaxy absolute magnitude or group mass. The alignment signal is dominated by red satellites. We also find that the outer regions of galaxies are aligned more strongly than their inner regions, by varying the radial weight employed during the shape measurement process. This behaviour is evident for both red and blue satellites. BGGs are also found to be aligned with satellite positions, with this alignment being stronger when considering the innermost satellites, using red BGGs and the shape of the outer region of the BGG. Lastly, we measure the global intrinsic alignment signal in the GAMA sample for two different radial weight functions and find no significant difference.

We present the study of a magnitude limited sample (mB < 16.6 mag) of 13 late type galaxies (LTGs), observed inside the virial radius, Rvir 0.7 Mpc, of the Fornax cluster within the Fornax Deep Survey (FDS). The main objective is to use surface brightness (SB) profiles and g-i colour maps to obtain information on the internal structure of these galaxies and find signatures of the mechanisms that drive their evolution in high-density environment, which is inside the virial radius of the cluster. By modelling galaxy isophotes, we extract the azimuthally averaged surface brightness profiles in four optical bands. We also derive g-i colour profiles, and structural parameters like total magnitude and effective radius. For 10 of the galaxies in this sample, we observe a clear discontinuity in their SB profiles, derive their break radius (Br), and classify their disc-breaks into Type-II (down-bending) or Type-III (up-bending). We find that Type-II galaxies have bluer average (g-i) colour in their outer discs while Type-III galaxies are redder. Br increases with stellar mass and molecular gas mass while decreases with molecular gas-fractions. The inner and outer scale-lengths increase monotonically with absolute magnitude, as found in other works. In Fornax, galaxies with morphological type 5< T< 9 (~60 % of the sample) are located beyond the high-density, ETG-dominated regions, but no correlation found between T and the disc-break type. The main results of this work suggest that the disc-breaks of LTGs inside the virial radius of the Fornax cluster seem to have arisen through a variety of mechanisms, which is evident in their outer-disc colours and the absence of molecular gas beyond their break radius in some cases. This can result in a variety of stellar populations inside and outside the break radii.

The KiDS Strongly lensed QUAsar Detection project (KiDS-SQuaD) aims at finding as many previously undiscovered gravitational lensed quasars as possible in the Kilo Degree Survey. This is the second paper of this series where we present a new, automatic object classification method based on machine learning technique. The main goal of this paper is to build a catalogue of bright extragalactic objects (galaxies and quasars), from the KiDS Data Release 4, with a minimum stellar contamination, preserving the completeness as much as possible, to then apply morphological methods to select reliable gravitationally lensed (GL) quasar candidates. After testing some of the most used machine learning algorithms, decision trees based classifiers, we decided to use CatBoost, that was specifically trained with the aim of creating a sample of extragalactic sources as clean as possible from stars. We discuss the input data, define the training sample for the classifier, give quantitative estimates of its performances, and finally describe the validation results with Gaia DR2, AllWISE, and GAMA catalogues. We have built and make available to the scientific community the KiDS Bright EXtraGalactic Objects catalogue (KiDS-BEXGO), specifically created to find gravitational lenses. This is made of $\approx6$ millions of sources classified as quasars ($\approx 200\,000$) and galaxies ($\approx 5.7$M), up to $r<22^m$. From this catalog we selected 'Multiplets': close pairs of quasars or galaxies surrounded by at least one quasar, presenting the 12 most reliable gravitationally lensed quasar candidates, to demonstrate the potential of the catalogue, which will be further explored in a forthcoming paper. We compared our search to the previous one, presented in the first paper from this series, showing that employing a machine learning method decreases the stars-contaminators within the GL candidates.

The analysis of the variability of active galactic nuclei (AGNs) at different wavelengths and the study of possible correlations among different spectral windows are nowadays a major field of inquiry. Optical variability has been largely used to identify AGNs in multivisit surveys. The strength of a selection based on optical variability lies in the chance to analyze data from surveys of large sky areas by ground-based telescopes. However the effectiveness of optical variability selection, with respect to other multiwavelength techniques, has been poorly studied down to the depth expected from next generation surveys. Here we present the results of our r-band analysis of a sample of 299 optically variable AGN candidates in the VST survey of the COSMOS field, counting 54 visits spread over three observing seasons spanning > 3 yr. This dataset is > 3 times larger in size than the one presented in our previous analysis (De Cicco et al. 2015), and the observing baseline is ~8 times longer. We push towards deeper magnitudes (r(AB) ~23.5 mag) compared to past studies; we make wide use of ancillary multiwavelength catalogs in order to confirm the nature of our AGN candidates, and constrain the accuracy of the method based on spectroscopic and photometric diagnostics. We also perform tests aimed at assessing the relevance of dense sampling in view of future wide-field surveys. We demonstrate that the method allows the selection of high-purity (> 86%) samples. We take advantage of the longer observing baseline to achieve great improvement in the completeness of our sample with respect to X-ray and spectroscopically confirmed samples of AGNs (59%, vs. ~15% in our previous work), as well as in the completeness of unobscured and obscured AGNs. The effectiveness of the method confirms the importance to develop future, more refined techniques for the automated analysis of larger datasets.

Many ultra diffuse galaxies (UDGs) have now been identified in clusters of galaxies. However, the number of nearby UDGs suitable for detailed follow-up remain rare. Our aim is to begin to identify UDGs in the environments of nearby bright early-type galaxies from the VEGAS survey. Here we use a deep g band image of the NGC 5846 group, taken as part of the VEGAS survey, to search for UDGs. We found one object with properties of a UDG if it associated with the NGC 5846 group, which seems likely. The galaxy, we name NGC 5846$\_$UDG1, has an absolute magnitude of M$_g$ = -14.2, corresponding to a stellar mass of $\sim$10$^8$ M$_{\odot}$. It also reveals a system of compact sources which are likely globular clusters. Based on the number of globular clusters detected we estimate a halo mass that is greater than 8$\times$10$^{10}$ M$_{\odot}$ for UDG1.

The Fornax Deep Survey Dwarf galaxy Catalog (FDSDC) includes 564 dwarf galaxies in the Fornax cluster and the in-falling Fornax A subgroup. We use the FDSDC galaxies for statistical analysis of the structural and stellar population differences in the range of galactic environments within the Fornax cluster. We present the standard scaling relations for the dwarfs and analyze trends as a function of cluster-centric radius. We find a different behavior for the bright dwarfs (-18.5 mag < M$_r$ < -16 mag) as compared to the fainter ones (M$_r$ > -16 mag): While considering galaxies in the same magnitude-bins, we find that, while for fainter dwarfs the g'-r' color is redder for lower surface brightness objects (as expected from fading stellar populations), for brighter dwarfs the color is redder for the higher surface brightness and higher S\'ersic n objects. The trend of the bright dwarfs might be explained by those galaxies being affected by harassment and by slower quenching of star formation in their inner parts. As the fraction of early-type dwarfs with respect to late-types increases toward the central parts of the cluster, the color-surface brightness trends are also manifested in the cluster-centric trends, confirming that it is indeed the environment that changes the galaxies. We also estimate the strengths of the ram-pressure stripping, tidal disruption, and harassment in the Fornax cluster, and find that our observations are consistent with the theoretically expected ranges of galaxy properties where each of those mechanisms dominate. We furthermore find that the luminosity function, color-magnitude relation, and axis-ratio distribution of the dwarfs in the center of the Fornax cluster are similar to those in the center of the Virgo cluster.

We present a catalog of quasars selected from broad-band photometric ugri data of the Kilo-Degree Survey Data Release 3 (KiDS DR3). The QSOs are identified by the random forest (RF) supervised machine learning model, trained on SDSS DR14 spectroscopic data. We first cleaned the input KiDS data from entries with excessively noisy, missing or otherwise problematic measurements. Applying a feature importance analysis, we then tune the algorithm and identify in the KiDS multiband catalog the 17 most useful features for the classification, namely magnitudes, colors, magnitude ratios, and the stellarity index. We used the t-SNE algorithm to map the multi-dimensional photometric data onto 2D planes and compare the coverage of the training and inference sets. We limited the inference set to r<22 to avoid extrapolation beyond the feature space covered by training, as the SDSS spectroscopic sample is considerably shallower than KiDS. This gives 3.4 million objects in the final inference sample, from which the random forest identified 190,000 quasar candidates. Accuracy of 97%, purity of 91%, and completeness of 87%, as derived from a test set extracted from SDSS and not used in the training, are confirmed by comparison with external spectroscopic and photometric QSO catalogs overlapping with the KiDS footprint. The robustness of our results is strengthened by number counts of the quasar candidates in the r band, as well as by their mid-infrared colors available from WISE. An analysis of parallaxes and proper motions of our QSO candidates found also in Gaia DR2 suggests that a probability cut of p(QSO)>0.8 is optimal for purity, whereas p(QSO)>0.7 is preferable for better completeness. Our study presents the first comprehensive quasar selection from deep high-quality KiDS data and will serve as the basis for versatile studies of the QSO population detected by this survey.

The Kilo-Degree Survey (KiDS) is an ongoing optical wide-field imaging survey with the OmegaCAM camera at the VLT Survey Telescope, specifically designed for measuring weak gravitational lensing by galaxies and large-scale structure. When completed it will consist of 1350 square degrees imaged in four filters (ugri). Here we present the fourth public data release which more than doubles the area of sky covered by data release 3. We also include aperture-matched ZYJHKs photometry from our partner VIKING survey on the VISTA telescope in the photometry catalogue. We illustrate the data quality and describe the catalogue content. Two dedicated pipelines are used for the production of the optical data. The Astro-WISE information system is used for the production of co-added images in the four survey bands, while a separate reduction of the r-band images using the theli pipeline is used to provide a source catalogue suitable for the core weak lensing science case. All data have been re-reduced for this data release using the latest versions of the pipelines. The VIKING photometry is obtained as forced photometry on the theli sources, using a re-reduction of the VIKING data that starts from the VISTA pawprints. Modifications to the pipelines with respect to earlier releases are described in detail. The photometry is calibrated to the Gaia DR2 G band using stellar locus regression. In this data release a total of 1006 square-degree survey tiles with stacked ugri images are made available, accompanied by weight maps, masks, and single-band source lists. We also provide a multi-band catalogue based on r-band detections, including homogenized photometry and photometric redshifts, for the whole dataset. Mean limiting magnitudes (5 sigma in a 2" aperture) are 24.23, 25.12, 25.02, 23.68 in ugri, respectively, and the mean r-band seeing is 0.70".

The halo masses $M_{halo}$ of low surface brightness (LSB) galaxies are critical measurements for understanding their formation processes. One promising method to estimate a galaxy's $M_{halo}$ is to exploit the empirical scaling relation between $M_{halo}$ and the number of associated globular clusters ($N_{\mathrm{GC}}$). We use a Bayesian mixture model approach to measure $N_{\mathrm{GC}}$ for 175 LSB ($23\leq\left\langle \mu_{e,r} \right\rangle [\mathrm{mag\ arcsec}^{-2}]\leq 28$) galaxies in the Fornax cluster using the Fornax Deep Survey (FDS) data; this is the largest sample of low mass galaxies so-far analysed for this kind of study. The proximity of the Fornax cluster means that we can measure galaxies with much smaller physical sizes ($0.3\leq r_{e,r}\ [\mathrm{kpc}]\leq 9.5$) compared to previous studies of the GC systems of LSB galaxies, probing stellar masses down to $M_{*}\sim10^{5}\mathrm{M_{\odot}}$. The sample also includes \nudg\ ultra-diffuse galaxies (UDGs), with projected $r$-band half-light radii greater than 1.5 kpc. Our results are consistent with an extrapolation of the $M_{*}-M_{halo}$ relation predicted from abundance matching. In particular, our UDG measurements are consistent with dwarf sized halos, having typical masses between $10^{10}$ and $10^{11}\mathrm{M_{\odot}}$. Overall, our UDG sample is statistically indistinguishable from smaller LSB galaxies in the same magnitude range. We do not find any candidates likely to be as rich as some of those found in the Coma cluster. We suggest that environment might play a role in producing GC-rich LSB galaxies.

We present a new sample of galaxy-scale strong gravitational-lens candidates, selected from 904 square degrees of Data Release 4 of the Kilo-Degree Survey (KiDS), i.e., the "Lenses in the Kilo-Degree Survey" (LinKS) sample. We apply two Convolutional Neural Networks (ConvNets) to $\sim88\,000$ colour-magnitude selected luminous red galaxies yielding a list of 3500 strong-lens candidates. This list is further down-selected via human inspection. The resulting LinKS sample is composed of 1983 rank-ordered targets classified as "potential lens candidates" by at least one inspector. Of these, a high-grade subsample of 89 targets is identified with potential strong lenses by all inspectors. Additionally, we present a collection of another 200 strong lens candidates discovered serendipitously from various previous ConvNet runs. A straightforward application of our procedure to future Euclid or LSST data can select a sample of $\sim3000$ lens candidates with less than 10 per cent expected false positives and requiring minimal human intervention.

Quadruply lensed quasars are extremely rare objects, but incredibly powerful cosmological tools. Only few dozen are known in the whole sky. Here we present the spectroscopic confirmation of two new quadruplets WG0214-2105 and WG2100-4452 discovered by Agnello & Spiniello (2018) within the Dark Energy Survey (DES) public footprints. We have conducted spectroscopic follow-up of these systems with the Southern African Large Telescope as part of a program that aims at confirming the largest possible number of optically selected strong gravitational lensing systems in the Equatorial and Southern Hemisphere. For both systems, we present the spectra for the sources and deflectors that allowed us to estimate the source redshifts and unambiguously confirm their lensing nature. For the brighter deflector (WG2100-4452), we measure the stellar velocity dispersion from the spectrum. We also obtain photometry for both lenses, directly from DES multi-band images, isolating the lens galaxies from the quasar images. One of the quadruplets, WG0214-2105, was also observed by Pan-STARRS, allowing us to estimate the apparent brightness of each quasar image at two different epochs, and thus to find evidence for flux variability. This result could suggest a microlensing event for the faintest components, although intrinsic variability cannot be excluded with only two epochs. Finally, we present simple lens models for both quadruplets, obtaining Einstein radii, SIE velocity dispersions, ellipticities, and position angles of the lens systems, as well as time delay predictions assuming a concordance cosmological model.

We present the mass calibration for galaxy clusters detected with the AMICO code in KiDS DR3 data. The cluster sample comprises $\sim$ 7000 objects and covers the redshift range 0.1 < $z$ < 0.6. We perform a weak lensing stacked analysis by binning the clusters according to redshift and two different mass proxies provided by AMICO, namely the amplitude $A$ (measure of galaxy abundance through an optimal filter) and the richness $\lambda^*$ (sum of membership probabilities in a consistent radial and magnitude range across redshift). For each bin, we model the data as a truncated NFW profile plus a 2-halo term, taking into account uncertainties related to concentration and miscentring. From the retrieved estimates of the mean halo masses, we construct the $A$-$M_{200}$ and the $\lambda^*$-$M_{200}$ relations. The relations extend over more than one order of magnitude in mass, down to $M_{200} \sim 2 (5) \times 10^{13} M_\odot/h$ at $z$ = 0.2 (0.5), with small evolution in redshift. The logarithmic slope is $\sim 2.0$ for the $A$-mass relation, and $\sim 1.7$ for the $\lambda^*$-mass relation, consistent with previous estimations on mock catalogues and coherent with the different nature of the two observables.

Gravitationally lensed quasars are valuable, but extremely rare, probes of observational cosmology and extragalactic astrophysics. Progress in these fields has been limited just by the paucity of systems with good ancillary data. Here we present a first spectroscopic confirmation of lenses discovered in the Southern Sky from the DES and KiDS-DR3 footprints. % optical ground based surveys in the Southern Hemisphere. We have targeted 7 high-graded candidates, selected with new techniques, with NTT-EFOSC2, and confirmed 5 of them. We provide source spectroscopic redshifts, image separations, $gri$ photometry and first lens model parameters. The success rate of ~70% confirms our forecasts, based on the comparison between the number of candidate doubles and quadruplets in our searches over a ~5000 sq.deg footprint and theoretical predictions.

This paper is based on the multi-band (ugri) Fornax Deep Survey (FDS) with the VLT Survey Telescope (VST). We study bright early-type galaxies (m_B<15 mag) in the 9 square degrees around the core of the Fornax cluster, which covers the virial radius. The main goal of the present work is to provide the analysis of the light distribution for all galaxies out to unprecedented limits (in radius and surface brightness) and to release the main products resulting from this analysis in all FDS bands. From the isophote fit, we derive the azimuthally averaged surface brightness profiles, the position angle and ellipticity profiles as a function of the semi-major axis. In each band, we derive the total magnitudes, effective radii, integrated colors and stellar mass-to-light ratios. The long integration times, the arcsec-level angular resolution of OmegaCam@VST and the large covered area of FDS allow us to map the light and color distributions out to large galactocentric distances (up to about 10-15 R_e) and surface brightness levels beyond mu_r = 27 mag/arcsec^2 (mu_B > 28 mag/arcsec^2). Therefore, the new FDS data allow us to explore in great detail the morphology and structure of cluster galaxies out to the region of the stellar halo. The observations suggest that the Fornax cluster is not completely relaxed inside the virial radius. The bulk of the gravitational interactions between galaxies happens in the W-NW core region of the cluster, where most of the bright early-type galaxies are located and where the intra-cluster baryons (diffuse light and GCs) are found. We suggest that the W-NW sub-clump of galaxies results from an infalling group onto the cluster, which has modified the structure of the galaxy outskirts (making asymmetric stellar halos) and has produced the intra-cluster baryons (ICL and GCs), concentrated in this region of the cluster.

We present a weak gravitational lensing analysis of 815 square degree of $i$-band imaging from the Kilo-Degree Survey (KiDS-$i$-800). In contrast to the deep $r$-band observations, which take priority during excellent seeing conditions and form the primary KiDS dataset (KiDS-$r$-450), the complementary yet shallower KiDS-$i$-800 spans a wide range of observing conditions. The overlapping KiDS-$i$-800 and KiDS-$r$-450 imaging therefore provides a unique opportunity to assess the robustness of weak lensing measurements. In our analysis, we introduce two new `null' tests. The `nulled' two-point shear correlation function uses a matched catalogue to show that the calibrated KiDS-$i$-800 and KiDS-$r$-450 shear measurements agree at the level of $1 \pm 4$\%. We use five galaxy lens samples to determine a `nulled' galaxy-galaxy lensing signal from the full KiDS-$i$-800 and KiDS-$r$-450 surveys and find that the measurements agree to $7 \pm 5$\% when the KiDS-$i$-800 source redshift distribution is calibrated using either spectroscopic redshifts, or the 30-band photometric redshifts from the COSMOS survey.

Despite the high accuracy of photometric redshifts (zphot) derived using Machine Learning (ML) methods, the quantification of errors through reliable and accurate Probability Density Functions (PDFs) is still an open problem. First, because it is difficult to accurately assess the contribution from different sources of errors, namely internal to the method itself and from the photometric features defining the available parameter space. Second, because the problem of defining a robust statistical method, always able to quantify and qualify the PDF estimation validity, is still an open issue. We present a comparison among PDFs obtained using three different methods on the same data set: two ML techniques, METAPHOR (Machine-learning Estimation Tool for Accurate PHOtometric Redshifts) and ANNz2, plus the spectral energy distribution template fitting method, BPZ. The photometric data were extracted from the KiDS (Kilo Degree Survey) ESO Data Release 3, while the spectroscopy was obtained from the GAMA (Galaxy and Mass Assembly) Data Release 2. The statistical evaluation of both individual and stacked PDFs was done through quantitative and qualitative estimators, including a dummy PDF, useful to verify whether different statistical estimators can correctly assess PDF quality. We conclude that, in order to quantify the reliability and accuracy of any zphot PDF method, a combined set of statistical estimators is required.

We report the independent discovery and characterization of two K2 planets: K2-180b, a mini-Neptune-size planet in an 8.9-day orbit transiting a V = 12.6 mag, metal-poor ([Fe/H] =$-0.65\pm0.10$) K2V star in K2 campaign 5; K2-140b, a transiting hot Jupiter in a 6.6-day orbit around a V = 12.6 mag G6V ([Fe/H] = $+0.10\pm0.10$) star in K2 campaign 10. Our results are based on K2 time-series photometry combined with high-spatial resolution imaging and high-precision radial velocity measurements. We present the first mass measurement of K2-180b. K2-180b has a mass of $M_\mathrm{p}=11.3\pm1.9$ ${M_{\oplus}}$ and a radius of $R_\mathrm{p}=2.2\pm0.1$ ${R_{\oplus}}$ , yielding a mean density of $\rho_\mathrm{p}=5.6\pm1.9\,g\,cm^{-3}$, suggesting a rock composition. Given its radius, K2-180b is above the region of the so-called `planetary radius gap'. K2-180b is in addition not only one of the densest mini-Neptune-size planets, but also one of the few mini-Neptune-size planets known to transit a metal-poor star. We also constrain the planetary and orbital parameters of K2-140b and show that, given the currently available Doppler measurements, the eccentricity is consistent with zero, contrary to the results of a previous study.

The Fornax Deep Survey (FDS), an imaging survey in the u', g', r', and i'-bands, has a supreme resolution and image depth compared to the previous spatially complete Fornax Cluster Catalog (FCC). Our new data allows us to study the galaxies down to r'-band magnitude m$_{r'}\approx$21 mag (M$_{r'}\approx$-10.5 mag). These data provide an important legacy dataset to study the Fornax cluster. We aim to present the Fornax Deep Survey (FDS) dwarf galaxy catalog, focusing on explaining the data reduction and calibrations, assessing the quality of the data, and describing the methods used for defining the cluster memberships for the catalog objects. As a first step we used the SExtractor fine-tuned for dwarf galaxy detection, to find galaxies from the FDS data, covering a 26 deg$^2$ area of the main cluster, and the area around the Fornax A substructure. We made 2D-decompositions of the identified galaxies using GALFIT. We used color-magnitude, luminosity-radius and luminosity-concentration relations to separate the cluster galaxies from the background galaxies. We then divided the cluster galaxies into early- and late-type galaxies according to their morphology and gave first order morphological classifications. Our final catalog includes 14,095 galaxies. We classify 590 galaxies as being likely Fornax cluster galaxies, of which 564 are dwarfs (M$_{r'}$ > -18.5 mag) consisting our Fornax dwarf catalog. Of the cluster dwarfs we classify 470 as early-types, and 94 as late-type galaxies. Our final catalog reaches its 50% completeness limit at magnitude M$_{r'}$ = -10.5 mag and surface brightness $\bar{\mu}_{e,r'}$ = 26 mag arcsec-2, which is approximately three magnitudes deeper than the FCC. Based on previous works and comparison with a spectroscopically confirmed subsample, we estimate that our final Fornax dwarf galaxy catalog has < 10% contamination from the background objects.

We report the discovery of a candidate to quadrupole gravitationally lensed system KiDS0239-3211 based on the public data release 3 of the KiDS survey and machine learning techniques.

Convolutional Neural Networks (ConvNets) are one of the most promising methods for identifying strong gravitational lens candidates in survey data. We present two ConvNet lens-finders which we have trained with a dataset composed of real galaxies from the Kilo Degree Survey (KiDS) and simulated lensed sources. One ConvNet is trained with single \textit{r}-band galaxy images, hence basing the classification mostly on the morphology. While the other ConvNet is trained on \textit{g-r-i} composite images, relying mostly on colours and morphology. We have tested the ConvNet lens-finders on a sample of 21789 Luminous Red Galaxies (LRGs) selected from KiDS and we have analyzed and compared the results with our previous ConvNet lens-finder on the same sample. The new lens-finders achieve a higher accuracy and completeness in identifying gravitational lens candidates, especially the single-band ConvNet. Our analysis indicates that this is mainly due to improved simulations of the lensed sources. In particular, the single-band ConvNet can select a sample of lens candidates with $\sim40\%$ purity, retrieving 3 out of 4 of the confirmed gravitational lenses in the LRG sample. With this particular setup and limited human intervention, it will be possible to retrieve, in future surveys such as Euclid, a sample of lenses exceeding in size the total number of currently known gravitational lenses.

We present results from an ongoing investigation using the Kilo Degree Survey (KiDS) on the VLT Survey Telescope (VST) to provide a census of ultra-compact massive galaxies (UCMGs), defined as galaxies with stellar masses $M_{\rm \star} > 8 \times 10^{10} \rm M_{\odot}$ and effective radii $R_{\rm e} < 1.5\,\rm kpc$. UCMGs, which are expected to have undergone very few merger events, provide a unique view on the accretion history of the most massive galaxies in the Universe. Over an effective sky area of nearly 330 square degrees, we select UCMG candidates from KiDS multi-colour images, which provide high quality structural parameters, photometric redshifts and stellar masses. Our sample of $\sim 1000$ photometrically selected UCMGs at $z < 0.5$ represents the largest sample of UCMG candidates assembled to date over the largest sky area. In this paper we present the first effort to obtain their redshifts using different facilities, starting with first results for 28 candidates with redshifts $z < 0.5$, obtained at NTT and TNG telescopes. We confirmed, as bona fide UCMGs, 19 out of the 28 candidates with new redshifts. A further 46 UCMG candidates are confirmed with literature spectroscopic redshifts (35 at $z < 0.5$), bringing the final cumulative sample of spectroscopically-confirmed lower-z UCMGs to 54 galaxies, which is the largest sample at redshifts below $0.5$. We use these spectroscopic redshifts to quantify systematic errors in our photometric selection, and use these to correct our UCMG number counts. We finally compare the results to independent datasets and simulations.

Most of the galaxies in the Universe at present day are in groups, which are key to understanding the galaxy evolution. In this work we present a new deep mosaic of 1.2 x 1.0 square degrees of the group of galaxies centered on NGC 5018, acquired at the ESO VLT Survey Telescope. We use u, g, r images to analyse the structure of the group members and to estimate the intra-group light. Taking advantage of the deep and multiband photometry and of the large field of view of the VST telescope, we studied the structure of the galaxy members and the faint features into the intra-group space and we give an estimate of the intragroup diffuse light in the NGC 5018 group of galaxies. We found that ~ 41% of the total g-band luminosity of the group is in the form of intragroup light (IGL). The IGL has a (g - r) color consistent with those of other galaxies in the group, indicating that the stripping leading to the formation of IGL is ongoing. From the study of this group we can infer that there are at least two different interactions involving the group members: one between NGC 5018 and NGC 5022, which generates the tails and ring-like structures detected in the light, and another between NGC 5022 and MCG-03-34-013 that have produced the HI tail. A minor merging event also happened in the formation history of NGC 5018 that have perturbed the inner structure of this galaxy.

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]

We have obtained structural parameters of about 340,000 galaxies from the Kilo Degree Survey (KiDS) in 153 square degrees of data release 1, 2 and 3. We have performed a seeing convolved 2D single S\'ersic fit to the galaxy images in the 4 photometric bands (u, g, r, i) observed by KiDS, by selecting high signal-to-noise ratio (S/N > 50) systems in every bands. We have classified galaxies as spheroids and disc-dominated by combining their spectral energy distribution properties and their S\'ersic index. Using photometric redshifts derived from a machine learning technique, we have determined the evolution of the effective radius, \Re\ and stellar mass, \mst, versus redshift, for both mass complete samples of spheroids and disc-dominated galaxies up to z ~ 0.6. Our results show a significant evolution of the structural quantities at intermediate redshift for the massive spheroids ($\mbox{Log}\ M_*/M_\odot>11$, Chabrier IMF), while almost no evolution has found for less massive ones ($\mbox{Log}\ M_*/M_\odot < 11$). On the other hand, disc dominated systems show a milder evolution in the less massive systems ($\mbox{Log}\ M_*/M_\odot < 11$) and possibly no evolution of the more massive systems. These trends are generally consistent with predictions from hydrodynamical simulations and independent datasets out to redshift z ~ 0.6, although in some cases the scatter of the data is large to drive final conclusions. These results, based on 1/10 of the expected KiDS area, reinforce precedent finding based on smaller statistical samples and show the route toward more accurate results, expected with the the next survey releases.

New methods have been recently developed to search for strong gravitational lenses, in particular lensed quasars, in wide-field imaging surveys. Here, we compare the performance of three different, morphology- and photometry- based methods to find lens candidates over the Kilo-Degree Survey (KiDS) DR3 footprint (440 deg$^2$). The three methods are: i) a multiplet detection in KiDS-DR3 and/or Gaia-DR1, ii) direct modeling of KiDS cutouts and iii) positional offsets between different surveys (KiDS-vs-Gaia, Gaia-vs-2MASS), with purpose-built astrometric recalibrations. The first benchmark for the methods has been set by the recovery of known lenses. We are able to recover seven out of ten known lenses and pairs of quasars observed in the KiDS DR3 footprint, or eight out of ten with improved selection criteria and looser colour pre-selection. This success rate reflects the combination of all methods together, which, taken individually, performed significantly worse (four lenses each). One movelty of our analysis is that the comparison of the performances of the different methods has revealed the pros and cons of the approaches and, most of all, the complementarities. We finally provide a list of high-grade candidates found by one or more methods, awaiting spectroscopic follow-up for confirmation. Of these, KiDS 1042+0023 is to our knowledge the first confirmed lensed quasar from KiDS, exhibiting two quasar spectra at the same source redshift at either sides of a red galaxy, with uniform flux-ratio $f\approx1.25$ over the wavelength range $0.45\mu\mathrm{m}<\lambda<0.75\mu\mathrm{m}.$

The VST Optical Imaging of the CDFS and ES1 Fields (VOICE) Survey is a Guaranteed Time program carried out with the ESO/VST telescope to provide deep optical imaging over two 4 deg$^2$ patches of the sky centred on the CDFS and ES1 pointings. We present the cosmic shear measurement over the 4 deg$^2$ covering the CDFS region in the $r$-band using LensFit. Each of the four tiles of 1 deg$^2$ has more than one hundred exposures, of which more than 50 exposures passed a series of image quality selection criteria for weak lensing study. The $5\sigma$ limiting magnitude in $r$- band is 26.1 for point sources, which is $\sim$1 mag deeper than other weak lensing survey in the literature (e.g. the Kilo Degree Survey, KiDS, at VST). The photometric redshifts are estimated using the VOICE $u,g,r,i$ together with near-infrared VIDEO data $Y,J,H,K_s$. The mean redshift of the shear catalogue is 0.87, considering the shear weight. The effective galaxy number density is 16.35 gal/arcmin$^2$, which is nearly twice the one of KiDS. The performance of LensFit on such a deep dataset was calibrated using VOICE-like mock image simulations. Furthermore, we have analyzed the reliability of the shear catalogue by calculating the star-galaxy cross-correlations, the tomographic shear correlations of two redshift bins and the contaminations of the blended galaxies. As a further sanity check, we have constrained cosmological parameters by exploring the parameter space with Population Monte Carlo sampling. For a flat $\Lambda$CDM model we have obtained $\Sigma_8$ = $\sigma_8(\Omega_m/0.3)^{0.5}$ = $0.68^{+0.11}_{-0.15}$.

The VST Optical Imaging of the CDFS and ES1 Fields (VOICE) Survey is proposed to obtain deep optical $ugri$ imaging of the CDFS and ES1 fields using the VLT Survey Telescope (VST). At present, the observations for the CDFS field have been completed, and comprise in total about 4.9 deg$^2$ down to $r_\mathrm{AB}$$\sim$26 mag. In the companion paper by Fu et al. (2018), we present the weak lensing shear measurements for $r$-band images with seeing $\le$ 0.9 arcsec. In this paper, we perform image simulations to calibrate possible biases of the measured shear signals. Statistically, the properties of the simulated point spread function (PSF) and galaxies show good agreements with those of observations. The multiplicative bias is calibrated to reach an accuracy of $\sim$3.0%. We study the bias sensitivities to the undetected faint galaxies and to the neighboring galaxies. We find that undetected galaxies contribute to the multiplicative bias at the level of $\sim$0.3%. Further analysis shows that galaxies with lower signal-to-noise ratio (SNR) are impacted more significantly because the undetected galaxies skew the background noise distribution. For the neighboring galaxies, we find that although most have been rejected in the shape measurement procedure, about one third of them still remain in the final shear sample. They show a larger ellipticity dispersion and contribute to $\sim$0.2% of the multiplicative bias. Such a bias can be removed by further eliminating these neighboring galaxies. But the effective number density of the galaxies can be reduced considerably. Therefore efficient methods should be developed for future weak lensing deep surveys.

We present a machine-learning photometric redshift analysis of the Kilo-Degree Survey Data Release 3, using two neural-network based techniques: ANNz2 and MLPQNA. Despite limited coverage of spectroscopic training sets, these ML codes provide photo-zs of quality comparable to, if not better than, those from the BPZ code, at least up to zphot<0.9 and r<23.5. At the bright end of r<20, where very complete spectroscopic data overlapping with KiDS are available, the performance of the ML photo-zs clearly surpasses that of BPZ, currently the primary photo-z method for KiDS. Using the Galaxy And Mass Assembly (GAMA) spectroscopic survey as calibration, we furthermore study how photo-zs improve for bright sources when photometric parameters additional to magnitudes are included in the photo-z derivation, as well as when VIKING and WISE infrared bands are added. While the fiducial four-band ugri setup gives a photo-z bias $\delta z=-2e-4$ and scatter $\sigma_z<0.022$ at mean z = 0.23, combining magnitudes, colours, and galaxy sizes reduces the scatter by ~7% and the bias by an order of magnitude. Once the ugri and IR magnitudes are joined into 12-band photometry spanning up to 12 $\mu$, the scatter decreases by more than 10% over the fiducial case. Finally, using the 12 bands together with optical colours and linear sizes gives $\delta z<4e-5$ and $\sigma_z<0.019$. This paper also serves as a reference for two public photo-z catalogues accompanying KiDS DR3, both obtained using the ANNz2 code. The first one, of general purpose, includes all the 39 million KiDS sources with four-band ugri measurements in DR3. The second dataset, optimized for low-redshift studies such as galaxy-galaxy lensing, is limited to r<20, and provides photo-zs of much better quality than in the full-depth case thanks to incorporating optical magnitudes, colours, and sizes in the GAMA-calibrated photo-z derivation.

We aim to investigate the galaxy environment in GAMA Galaxy Groups Catalogue (G3C) using a volume-limited galaxy sample from the Kilo Degree Survey Data Release 3. The k-Nearest Neighbour technique is adapted to take into account the probability density functions (PDFs) of photometric redshifts in our calculations. This algorithm was tested on simulated KiDS tiles, showing its capability of recovering the relation between galaxy colour, luminosity and local environment. The characterization of the galaxy environment in G3C groups shows systematically steeper density contrasts for more massive groups. The red galaxy fraction gradients in these groups is evident for most of group mass bins. The density contrast of red galaxies is systematically higher at group centers when compared to blue galaxy ones. In addition, distinct group center definitions are used to show that our results are insensitive to center definitions. These results confirm the galaxy evolution scenario which environmental mechanisms are responsible for a slow quenching process as galaxies fall into groups and clusters, resulting in a smooth observed colour gradients in galaxy systems.

Using a sample of 472 local Universe (z<0.06) galaxies in the stellar mass range 10.25 < log M*/M_sun < 10.75, we explore the variation in galaxy structure as a function of morphology and galaxy colour. Our sample of galaxies is sub-divided into red, green and blue colour groups and into elliptical and non-elliptical (disk-type) morphologies. Using KiDS and VIKING derived postage stamp images, a group of eight volunteers visually classified bars, rings, morphological lenses, tidal streams, shells and signs of merger activity for all systems. We find a significant surplus of rings ($2.3\sigma$) and lenses ($2.9\sigma$) in disk-type galaxies as they transition across the green valley. Combined, this implies a joint ring/lens green valley surplus significance of $3.3\sigma$ relative to equivalent disk-types within either the blue cloud or the red sequence. We recover a bar fraction of ~44% which remains flat with colour, however, we find that the presence of a bar acts to modulate the incidence of rings and (to a lesser extent) lenses, with rings in barred disk-type galaxies more common by ~20-30 percentage points relative to their unbarred counterparts, regardless of colour. Additionally, green valley disk-type galaxies with a bar exhibit a significant $3.0\sigma$ surplus of lenses relative to their blue/red analogues. The existence of such structures rules out violent transformative events as the primary end-of-life evolutionary mechanism, with a more passive scenario the favoured candidate for the majority of galaxies rapidly transitioning across the green valley.

We present the largest and most spatially extended planetary nebulae (PNe) catalog ever obtained for the Fornax cluster. We measured velocities of 1452 PNe out to 200 kpc in the cluster core using a counter-dispersed slitless spectroscopic technique with data from FORS2 on the VLT. With such extended spatial coverage, we can study separately the stellar halos of some of the cluster main galaxies and the intracluster light. In this second paper of the Fornax Cluster VLT Spectroscopic Survey (FVSS), we identify and classify the emission-line sources, describe the method to select PNe and calculate their coordinates and velocities from the dispersed slitless images. From the PN 2D velocity map we identify stellar streams that are possibly tracing the gravitational interaction of NGC1399 with NGC1404 and NGC1387. We also present the velocity dispersion profile out to $\sim 200$ kpc radii, which shows signatures of a superposition of the bright central galaxy and the cluster potential, with the latter clearly dominating the regions outside R $\sim 1000$" ($\sim 100$ kpc).

We present the results of a wide spectroscopic survey aimed at detecting extragalactic globular clusters (GCs) in the core of the Fornax cluster. About 4500 low resolution spectra (from 4800 to 10000 \AA ) were observed in 25 VLT/VIMOS masks covering the central 1 deg$^{2}$ around the dominant galaxy NGC 1399 corresponding to $\sim$175 kpc galactocentric radius. We describe the methodology used for data reduction and data analysis. We found a total of 387 unique physical objects (372 GCs and 15 ultra compact dwarfs) in the field covered by our observations. Most of these objects are associated with NGC 1399, with only 10% likely belonging to other giant galaxies. The new VIMOS dataset is complementary to the many GC catalogues already present in the literature and it brings the total number of tracer particles around NGC 1399 to more than 1130 objects. With this comprehensive radial velocity sample we have found that the velocity dispersion of the GC population (equally for red and blue GC populations) shows a relatively sharp increase from low velocity dispersion ($\sim250$-$350$ kms$^{-1}$) to high velocity dispersion ($\sim300$-$400$ kms$^{-1}$) at projected radius of $\approx10$ arcmin ($\sim60$ kpc) from the galaxy centre. This suggests that at a projected radius of $\approx60$ kpc both blue and red GC populations begin to be governed by the dominating Fornax cluster potential, rather than by the central NGC 1399 galaxy potential. This kinematic evidence corroborates similar results found using surface brightness analysis and planetary nebulae kinematics.

In the core of the Fornax cluster, on the West side of NGC1399, we have detected a previously unknown region of intra-cluster light (ICL). It is made up by several faint ($\mu_r \simeq 28 - 29$~mag/arcsec$^2$) {\it patches} of diffuse light. The bulk of the ICL is located in between the three bright galaxies in the core, NGC1387, NGC1379 and NGC1381, at $10\leq R \leq40$~arcmin ($\sim 58 - 230$~kpc) from the central galaxy NGC~1399. We show that the ICL is the counterpart in the diffuse light of the known over-density in the population of blue globular clusters (GCs). The total g-band luminosity of the ICL is $L_g\simeq 8.3 \times 10^{9}$ $L_{\odot}$, which is $\sim5\%$ of the total luminosity of NGC1399. This is consistent with the fraction of the blue GCs in the same region of the cluster. The ICL has $g-r \sim 0.7$~mag, which is similar to the colors in the halo of the bright galaxies in the cluster core. The new findings were compared with theoretical predictions for the ICL formation and they support a scenario in which the intra-cluster population detected in the core of the Fornax cluster is build up by the tidal stripping of material (stars and GCs) from galaxy outskirts in a close passage with the cD. Moreover, the diffuse form of the ICL and its location close to the core of the cluster is expected in a dynamically evolved cluster as Fornax.

We analyze the globular cluster (GC) systems in two very different galaxies, N3115 and N1399. With the papers of this series, we aim at highlighting common and different properties in the GC systems in galaxies covering a wide range of parameter space. We compare the GCs in N3115 and N1399 as derived from the analysis of u, g, and i band images taken with the VST telescope as part of the VEGAS and FDS surveys. We selected GC candidates using as reference the morpho-photometric and color properties of confirmed GCs. The surface density of GCs in N3115 reveal a morphology similar to the profile of field stars; the same is true when blue and red GCs are taken separately. The GC maps for N1399 are richer in structure and confirm the existence of an intracluster GC component. We confirm the presence of a spatial offset in the N1399 GC centroid and find that the centroid of the GCs for N3115 coincides well with the galaxy center. Both GC systems show unambiguous color bimodality in g-i and u-i. The azimuthal average of the radial density profiles in both galaxies reveals a larger spatial extent for the total GCs population with respect to the galaxy surface brightness profile. As for the specific frequency of GCs, Sn, we find it is a factor of two higher in N1399 than for N3115. By inspecting the radial behavior of the specific frequency, Sn(<r), for the total, blue, and red GCs, we find notable similarities between the trends for red GCs in the two targets. In spite of extremely different host environments, the red GCs in both cases appear closely linked to the light distribution of field stars. Blue GCs extend to larger galactocentric scales than red GCs, marking a significant difference between the two galaxies: the blue/red GCs and field stellar components of N3115 appear well thermalized with each other and the blue GCs in N1399 appear to fade into an unrelaxed intracluster GC population.

We describe the observing strategy, data reduction tools and early results of a supernova (SN) search project, named SUDARE, conducted with the ESO VST telescope aimed at measuring the rate of the different types of SNe in the redshift range 0.2<z<0.8. The search was performed in two of the best-studied extragalactic fields, CDFS and COSMOS, for which a wealth of ancillary data are available in the literature or public archives. (abridged) We obtained a final sample of 117 SNe, most of which are SNIa (57%) and the remaining core collapse events of which 44% type II, 22% type IIn and 34% type Ib/c. In order to link the transients, we built a catalog of ~1.3x10^5 galaxies in the redshift range 0<z<1 with a limiting magnitude K_AB=23.5 mag. We measured the SN rate per unit volume for SN Ia and core collapse SNe in different bin of redshifts. The values are consistent with other measurements from the literature. The dispersion of the rate measurements for SNe Ia is comparable with the scatter of the theoretical tracks for single (SD) and double degenerate (DD) binary systems models, therefore the data do not allow to disentangle among the two different progenitor scenarios. However, we may notice that among the three tested models, SD and two flavours of DD, either with a steep (DDC) or a wide (DDW) delay time distribution, the SD gives a better fit across the whole redshift range whereas the DDC better matches the steep rise up to redshift ~1.2. The DDW appears instead less favoured. The core collapse SN rate is fully consistent, unlike recent claims, with the prediction based on recent estimates of the star formation history, and standard progenitor mass range.

Studies of low surface brightness (LSB) galaxies in nearby clusters have revealed a sub-population of extremely diffuse galaxies with central surface brightness $\mu_{0,g'}$ > 24 mag arcsec$^{-2}$ and effective radius between 1.5 kpc < R$_{e}$ < 10 kpc. The origin of these Ultra Diffuse Galaxies (UDGs) is still unclear, although several theories have been suggested. We exploit the deep g', r' and i'-band images of the Fornax Deep Survey (FDS), in order to identify LSB galaxies in the center of the Fornax cluster. We identified visually all extended structures having r'-band central surface brightness of $\mu_{0,r'}$ > 23 mag arcsec$^{-2}$. We classified the objects based on their appearance and performed 2D S\'ersic model fitting with GALFIT. We analyzed their distribution and orientations in the cluster, and studied their colors and compared the LSB galaxies in Fornax with those in other environments. Our sample consists of 205 galaxies of which 196 are LSB dwarfs (with R$_e$ < 1.5kpc) and nine are UDGs (R$_e$ > 1.5 kpc). We show that the UDGs have g'-r' colors similar to those of LSB dwarfs. The largest UDGs in our sample appear different from the other LSB galaxies, in that they are significantly more elongated and extended, whereas the smaller UDGs differ from the LSB dwarfs only by their effective radii. We do not find clear differences between the structural parameters of the UDGs in our sample and those of UDGs in other galaxy environments. We find that the dwarf LSB galaxies in our sample are less concentrated in the cluster center than the galaxies with higher surface brightness, and that their number density drops in the core of the cluster. Our findings are consistent with the small UDGs forming the tail of a continuous distribution of less extended LSB galaxies. However, the elongated and distorted shapes of the large UDGs could imply that they are tidally disturbed galaxies.

We study the dark matter (DM) assembly in the central regions of massive early-type galaxies up to $z\sim 0.65$. We use a sample of $\sim 3800$ massive ($\log M_{\rm \star}/M_{\rm \odot} > 11.2$) galaxies with photometry and structural parameters from 156 sq. deg. of the Kilo Degree Survey, and spectroscopic redshifts and velocity dispersions from SDSS. We obtain central total-to-stellar mass ratios, $M_{\rm dyn}/M_{\rm \star}$, and DM fractions, by determining dynamical masses, $M_{\rm dyn}$, from Jeans modelling of SDSS aperture velocity dispersions and stellar masses, $M_{\rm \star}$, from KiDS galaxy colours. We first show how the central DM fraction correlates with structural parameters, mass and density proxies, and demonstrate that most of the local correlations are still observed up to $z \sim 0.65$; at fixed $M_{\rm \star}$, local galaxies have larger DM fraction, on average, than their counterparts at larger redshift. We also interpret these trends with a non universal Initial Mass Function (IMF), finding a strong evolution with redshift, which contrast independent observations and is at odds with the effect of galaxy mergers. For a fixed IMF, the galaxy assembly can be explained, realistically, by mass and size accretion, which can be physically achieved by a series of minor mergers. We reproduce both the $R_{\rm e}-M_{\rm \star}$ and $M_{\rm dyn}/M_{\rm \star}-M_{\rm \star}$ evolution with stellar and dark mass changing at a different rate. This result suggests that the main progenitor galaxy is merging with less massive systems, characterized by a smaller $M_{\rm dyn}/M_{\rm \star}$, consistently with results from halo abundance matching.

The volume of data that will be produced by new-generation surveys requires automatic classification methods to select and analyze sources. Indeed, this is the case for the search for strong gravitational lenses, where the population of the detectable lensed sources is only a very small fraction of the full source population. We apply for the first time a morphological classification method based on a Convolutional Neural Network (CNN) for recognizing strong gravitational lenses in $255$ square degrees of the Kilo Degree Survey (KiDS), one of the current-generation optical wide surveys. The CNN is currently optimized to recognize lenses with Einstein radii $\gtrsim 1.4$ arcsec, about twice the $r$-band seeing in KiDS. In a sample of $21789$ colour-magnitude selected Luminous Red Galaxies (LRG), of which three are known lenses, the CNN retrieves 761 strong-lens candidates and correctly classifies two out of three of the known lenses. The misclassified lens has an Einstein radius below the range on which the algorithm is trained. We down-select the most reliable 56 candidates by a joint visual inspection. This final sample is presented and discussed. A conservative estimate based on our results shows that with our proposed method it should be possible to find $\sim100$ massive LRG-galaxy lenses at $z\lsim 0.4$ in KiDS when completed. In the most optimistic scenario this number can grow considerably (to maximally $\sim$2400 lenses), when widening the colour-magnitude selection and training the CNN to recognize smaller image-separation lens systems.

In recent years, multiple studies have reported substantial populations of large, low-surface-brightness galaxies in local galaxy clusters. Various theories that aim to explain the presence of such ultra-diffuse galaxies (UDGs) have since been proposed. A key question that will help to differentiate between models is whether UDGs have counterparts in lower-mass host haloes, and what their abundance as a function of halo mass is. In this study we extend our previous study of UDGs in galaxy clusters to galaxy groups. We measure the abundance of UDGs in 325 spectroscopically-selected groups from the Galaxy And Mass Assembly (GAMA) survey. We make use of the overlapping imaging from the ESO Kilo-Degree Survey (KiDS), from which we can identify galaxies with mean surface brightnesses within their effective radii down to ~25.5 mag arcsec$^{-2}$ in the r-band. We are able to measure a significant overdensity of UDGs (with sizes r_eff > 1.5 kpc) in galaxy groups down to M200=10^12 Msun, a regime where approximately only 1 in 10 groups contains a UDG that we can detect. We combine measurements of the abundance of UDGs in haloes that cover three orders of magnitude in halo mass, finding that their numbers scale quite steeply with halo mass; N_UDG (R<R200) $\propto$ M200^(1.11+/-0.07). To better interpret this, we also measure the mass-richness relation for brighter galaxies down to $M^*_r$+2.5 in the same GAMA groups, and find a much shallower relation of N_Bright (R<R200) $\propto$ M200^(0.78+/-0.05). This shows that UDGs are relatively more abundant, compared to bright galaxies, in massive clusters than in groups. We discuss implications, but whether this difference is related to a higher destruction rate of UDGs in groups, or whether massive haloes have a positive effect on their formation, is not yet clear.

The Large Synoptic Survey Telescope (LSST) will enable revolutionary studies of galaxies, dark matter, and black holes over cosmic time. The LSST Galaxies Science Collaboration has identified a host of preparatory research tasks required to leverage fully the LSST dataset for extragalactic science beyond the study of dark energy. This Galaxies Science Roadmap provides a brief introduction to critical extragalactic science to be conducted ahead of LSST operations, and a detailed list of preparatory science tasks including the motivation, activities, and deliverables associated with each. The Galaxies Science Roadmap will serve as a guiding document for researchers interested in conducting extragalactic science in anticipation of the forthcoming LSST era.

Verlinde derived gravity as an emergent force from the information flow, through two-dimensional surfaces and recently, by a priori postulating the entanglement of information in 3D space, he derived the effect of the gravitational potential from dark matter (DM) as the entropy displacement of dark energy by baryonic matter. In Emergent Gravity (EG) this apparent DM depends only on the baryonic mass distribution and the present-day value of the Hubble parameter. In this paper we test the EG proposition, formalized by Verlinde for a spherical and isolated mass distribution, using the central velocity dispersion, $\sigma$ and the light distribution in a sample of 4260 massive and local early-type galaxies (ETGs) from the SPIDER sample. Our results remain unaltered if we consider the sample of 807 roundest field galaxies. We derive the predictions by EG for the stellar mass-to-light ratio (M/L) and the Initial Mass Function (IMF), and compare them with the same inferences derived from a) DM-based models, b) MOND and c) stellar population models. We demonstrate that, consistently with a classical Newtonian framework with a DM halo component, or alternative theories of gravity as MOND, the central dynamics can be fitted if the IMF is assumed non-universal. The results can be interpreted with a IMF lighter than a standard Chabrier at low-$\sigma$, and bottom-heavier IMFs at larger $\sigma$. We find lower, but still acceptable, stellar M/L in EG theory, if compared with the DM-based NFW model and with MOND. The results from EG are comparable to what is found if the DM haloes are adiabatically contracted and with expectations from spectral gravity-sensitive features. If the strain caused by the entropy displacement would be not maximal, as adopted in the current formulation, then the dynamics of ETGs could be reproduced with larger M/L. (abridged)

The Kilo-Degree Survey (KiDS) is an ongoing optical wide-field imaging survey with the OmegaCAM camera at the VLT Survey Telescope. It aims to image 1500 square degrees in four filters (ugri). The core science driver is mapping the large-scale matter distribution in the Universe, using weak lensing shear and photometric redshift measurements. Further science cases include galaxy evolution, Milky Way structure, detection of high-redshift clusters, and finding rare sources such as strong lenses and quasars. Here we present the third public data release (DR3) and several associated data products, adding further area, homogenized photometric calibration, photometric redshifts and weak lensing shear measurements to the first two releases. A dedicated pipeline embedded in the Astro-WISE information system is used for the production of the main release. Modifications with respect to earlier releases are described in detail. Photometric redshifts have been derived using both Bayesian template fitting, and machine-learning techniques. For the weak lensing measurements, optimized procedures based on the THELI data reduction and lensfit shear measurement packages are used. In DR3 stacked ugri images, weight maps, masks, and source lists for 292 new survey tiles (~300 sq.deg) are made available. The multi-band catalogue, including homogenized photometry and photometric redshifts, covers the combined DR1, DR2 and DR3 footprint of 440 survey tiles (447 sq.deg). Limiting magnitudes are typically 24.3, 25.1, 24.9, 23.8 (5 sigma in a 2 arcsec aperture) in ugri, respectively, and the typical r-band PSF size is less than 0.7 arcsec. The photometric homogenization scheme ensures accurate colors and an absolute calibration stable to ~2% for gri and ~3% in u. Separately released are a weak lensing shear catalogue and photometric redshifts based on two different machine-learning techniques.

We constrain the average halo ellipticity of ~2 600 galaxy groups from the Galaxy And Mass Assembly (GAMA) survey, using the weak gravitational lensing signal measured from the overlapping Kilo Degree Survey (KiDS). To do so, we quantify the azimuthal dependence of the stacked lensing signal around seven different proxies for the orientation of the dark matter distribution, as it is a priori unknown which one traces the orientation best. On small scales, the major axis of the brightest group/cluster member (BCG) provides the best proxy, leading to a clear detection of an anisotropic signal. In order to relate that to a halo ellipticity, we have to adopt a model density profile. We derive new expressions for the quadrupole moments of the shear field given an elliptical model surface mass density profile. Modeling the signal with an elliptical Navarro-Frenk-White (NFW) profile on scales < 250 kpc, which roughly corresponds to half the virial radius, and assuming that the BCG is perfectly aligned with the dark matter, we find an average halo ellipticity of e_h=0.38 +/- 0.12. This agrees well with results from cold-dark-matter-only simulations, which typically report values of e_h ~ 0.3. On larger scales, the lensing signal around the BCGs does not trace the dark matter distribution well, and the distribution of group satellites provides a better proxy for the halo's orientation instead, leading to a 3--4 sigma detection of a non-zero halo ellipticity at scales between 250 kpc and 750 kpc. Our results suggest that the distribution of stars enclosed within a certain radius forms a good proxy for the orientation of the dark matter within that radius, which has also been observed in hydrodynamical simulations.

We present the VST Optical Imaging of the CDFS and ES1 Fields (VOICE) Survey, a VST INAF Guaranteed Time program designed to provide optical coverage of two 4 deg$^2$ cosmic windows in the Southern hemisphere. VOICE provides the first, multi-band deep optical imaging of these sky regions, thus complementing and enhancing the rich legacy of longer-wavelength surveys with VISTA, Spitzer, Herschel and ATCA available in these areas and paving the way for upcoming observations with facilities such as the LSST, MeerKAT and the SKA. VOICE exploits VST's OmegaCAM optical imaging capabilities and completes the reduction of WFI data available within the ES1 fields as part of the ESO-Spitzer Imaging Extragalactic Survey (ESIS) program providing $ugri$ and $uBVR$ coverage of 4 and 4 deg$^2$ areas within the CDFS and ES1 field respectively. We present the survey's science rationale and observing strategy, the data reduction and multi-wavelength data fusion pipeline. Survey data products and their future updates will be released at http://www.mattiavaccari.net/voice/ and on CDS/VizieR.

Observations of diffuse starlight in the outskirts of galaxies are thought to be a fundamental source of constraints on the cosmological context of galaxy assembly in the $\Lambda$CDM model. Such observations are not trivial because of the extreme faintness of such regions. In this work, we investigate the photometric properties of six massive early type galaxies (ETGs) in the VEGAS sample (NGC 1399, NGC 3923, NGC 4365, NGC 4472, NGC 5044, and NGC 5846) out to extremely low surface brightness levels, with the goal of characterizing the global structure of their light profiles for comparison to state-of-the-art galaxy formation models. We carry out deep and detailed photometric mapping of our ETG sample taking advantage of deep imaging with VST/OmegaCAM in the g and i bands. By fitting the light profiles, and comparing the results to simulations of elliptical galaxy assembly, we identify signatures of a transition between "relaxed" and "unrelaxed" accreted components and can constrain the balance between in situ and accreted stars. The very good agreement of our results with predictions from theoretical simulations demonstrates that the full VEGAS sample of $\sim 100$ ETGs will allow us to use the distribution of diffuse light as a robust statistical probe of the hierarchical assembly of massive galaxies.

As part of the Fornax Deep Survey with the ESO VLT Survey Telescope, we present new $g$ and $r$ bands mosaics of the SW group of the Fornax cluster. It covers an area of $3 \times 2$ square degrees around the central galaxy NGC1316. The deep photometry, the high spatial resolution of OmegaCam and the large covered area allow us to study the galaxy structure, to trace stellar halo formation and look at the galaxy environment. We map the surface brightness profile out to 33arcmin ($\sim 200$kpc $\sim15R_e$) from the galaxy centre, down to $\mu_g \sim 31$ mag arcsec$^{-2}$ and $\mu_r \sim 29$ mag arcsec$^{-2}$. This allow us to estimate the scales of the main components dominating the light distribution, which are the central spheroid, inside 5.5 arcmin ($\sim33$ kpc), and the outer stellar envelope. Data analysis suggests that we are catching in act the second phase of the mass assembly in this galaxy, since the accretion of smaller satellites is going on in both components. The outer envelope of NGC1316 still hosts the remnants of the accreted satellite galaxies that are forming the stellar halo. We discuss the possible formation scenarios for NGC1316, by comparing the observed properties (morphology, colors, gas content, kinematics and dynamics) with predictions from cosmological simulations of galaxy formation. We find that {\it i)} the central spheroid could result from at least one merging event, it could be a pre-existing early-type disk galaxy with a lower mass companion, and {\it ii)} the stellar envelope comes from the gradual accretion of small satellites.

In the modern galaxy surveys photometric redshifts play a central role in a broad range of studies, from gravitational lensing and dark matter distribution to galaxy evolution. Using a dataset of about 25,000 galaxies from the second data release of the Kilo Degree Survey (KiDS) we obtain photometric redshifts with five different methods: (i) Random forest, (ii) Multi Layer Perceptron with Quasi Newton Algorithm, (iii) Multi Layer Perceptron with an optimization network based on the Levenberg-Marquardt learning rule, (iv) the Bayesian Photometric Redshift model (or BPZ) and (v) a classical SED template fitting procedure (Le Phare). We show how SED fitting techniques could provide useful information on the galaxy spectral type which can be used to improve the capability of machine learning methods constraining systematic errors and reduce the occurrence of catastrophic outliers. We use such classification to train specialized regression estimators, by demonstrating that such hybrid approach, involving SED fitting and machine learning in a single collaborative framework, is capable to improve the overall prediction accuracy of photometric redshifts.

In this paper, we present the tools used to search for galaxy clusters in the Kilo Degree Survey (KiDS), and our first results. The cluster detection is based on an implementation of the optimal filtering technique that enables us to identify clusters as over-densities in the distribution of galaxies using their positions on the sky, magnitudes, and photometric redshifts. The contamination and completeness of the cluster catalog are derived using mock catalogs based on the data themselves. The optimal signal to noise threshold for the cluster detection is obtained by randomizing the galaxy positions and selecting the value that produces a contamination of less than 20%. Starting from a subset of clusters detected with high significance at low redshifts, we shift them to higher redshifts to estimate the completeness as a function of redshift: the average completeness is ~ 85%. An estimate of the mass of the clusters is derived using the richness as a proxy. We obtained 1858 candidate clusters with redshift 0 < z_c < 0.7 and mass 13.5 < log(M500/Msun) < 15 in an area of 114 sq. degrees (KiDS ESO-DR2). A comparison with publicly available Sloan Digital Sky Survey (SDSS)-based cluster catalogs shows that we match more than 50% of the clusters (77% in the case of the redMaPPer catalog). We also cross-matched our cluster catalog with the Abell clusters, and clusters found by XMM and in the Planck-SZ survey; however, only a small number of them lie inside the KiDS area currently available.

Photometric redshifts (photo-z's) are fundamental in galaxy surveys to address different topics, from gravitational lensing and dark matter distribution to galaxy evolution. The Kilo Degree Survey (KiDS), i.e. the ESO public survey on the VLT Survey Telescope (VST), provides the unprecedented opportunity to exploit a large galaxy dataset with an exceptional image quality and depth in the optical wavebands. Using a KiDS subset of about 25,000 galaxies with measured spectroscopic redshifts, we have derived photo-z's using i) three different empirical methods based on supervised machine learning, ii) the Bayesian Photometric Redshift model (or BPZ), and iii) a classical SED template fitting procedure (Le Phare). We confirm that, in the regions of the photometric parameter space properly sampled by the spectroscopic templates, machine learning methods provide better redshift estimates, with a lower scatter and a smaller fraction of outliers. SED fitting techniques, however, provide useful information on the galaxy spectral type which can be effectively used to constrain systematic errors and to better characterize potential catastrophic outliers. Such classification is then used to specialize the training of regression machine learning models, by demonstrating that a hybrid approach, involving SED fitting and machine learning in a single collaborative framework, can be effectively used to improve the accuracy of photo-z estimates.

We present a sample of 80 candidate strongly lensed galaxies with flux density above 100mJy at 500{\mu}m extracted from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), over an area of 600 square degrees. Available imaging and spectroscopic data allow us to confirm the strong lensing in 20 cases and to reject it in one case. For other 8 objects the lensing scenario is strongly supported by the presence of two sources along the same line of sight with distinct photometric redshifts. The remaining objects await more follow-up observations to confirm their nature. The lenses and the background sources have median redshifts z_L = 0.6 and z_S = 2.5, respectively, and are observed out to z_L = 1.2 and z_S = 4.2. We measure the number counts of candidate lensed galaxies at 500{\mu}m and compare them with theoretical predictions, finding a good agreement for a maximum magnification of the background sources in the range 10-20. These values are consistent with the magnification factors derived from the lens modelling of individual systems. The catalogue presented here provides sub- mm bright targets for follow-up observations aimed at exploiting gravitational lensing to study with un-precedented details the morphological and dynamical properties of dusty star forming regions in z >~ 1.5 galaxies.

We present cosmological parameter constraints from a tomographic weak gravitational lensing analysis of ~450deg$^2$ of imaging data from the Kilo Degree Survey (KiDS). For a flat $\Lambda$CDM cosmology with a prior on $H_0$ that encompasses the most recent direct measurements, we find $S_8\equiv\sigma_8\sqrt{\Omega_{\rm m}/0.3}=0.745\pm0.039$. This result is in good agreement with other low redshift probes of large scale structure, including recent cosmic shear results, along with pre-Planck cosmic microwave background constraints. A $2.3$-$\sigma$ tension in $S_8$ and `substantial discordance' in the full parameter space is found with respect to the Planck 2015 results. We use shear measurements for nearly 15 million galaxies, determined with a new improved `self-calibrating' version of $lens$fit validated using an extensive suite of image simulations. Four-band $ugri$ photometric redshifts are calibrated directly with deep spectroscopic surveys. The redshift calibration is confirmed using two independent techniques based on angular cross-correlations and the properties of the photometric redshift probability distributions. Our covariance matrix is determined using an analytical approach, verified numerically with large mock galaxy catalogues. We account for uncertainties in the modelling of intrinsic galaxy alignments and the impact of baryon feedback on the shape of the non-linear matter power spectrum, in addition to the small residual uncertainties in the shear and redshift calibration. The cosmology analysis was performed blind. Our high-level data products, including shear correlation functions, covariance matrices, redshift distributions, and Monte Carlo Markov Chains are available at http://kids.strw.leidenuniv.nl.

This is the second paper of a series in which we present measurements of the Supernova (SN) rates from the SUDARE survey. In this paper, we study the trend of the SN rates with the intrinsic colours, the star formation activity and the mass of the parent galaxies. We have considered a sample of about 130000 galaxies and a SN sample of about 50 events. We found that the SN Ia rate per unit mass is higher by a factor of six in the star-forming galaxies with respect to the passive galaxies. The SN Ia rate per unit mass is also higher in the less massive galaxies that are also younger. These results suggest a distribution of the delay times (DTD) less populated at long delay times than at short delays. The CC SN rate per unit mass is proportional to both the sSFR and the galaxy mass. The trends of the Type Ia and CC SN rates as a function of the sSFR and the galaxy mass that we observed from SUDARE data are in agreement with literature results at different redshifts. The expected number of SNe Ia is in agreement with the observed one for all four DTD models considered both in passive and star-forming galaxies so we can not discriminate between different progenitor scenarios. The expected number of CC SNe is higher than the observed one, suggesting a higher limit for the minimum progenitor mass. We also compare the expected and observed trends of the SN Ia rate with the intrinsic U - J colour of the parent galaxy, assumed as a tracer of the age distribution. While the slope of the relation between the SN Ia rate and the U - J color in star-forming galaxies can be reproduced well by all four DTD models considered, only the steepest of them is able to account for the rates and colour in star-forming and passive galaxies with the same value of the SN Ia production efficiency.

We construct a suite of discrete chemo-dynamical models of the giant elliptical galaxy NGC 5846. These models are a powerful tool to constrain both the mass distribution and internal dynamics of multiple tracer populations. We use Jeans models to simultaneously fit stellar kinematics within the effective radius $R_{\rm e}$, planetary nebula (PN) radial velocities out to $3\, R_{\rm e}$, and globular cluster (GC) radial velocities and colours out to $6\,R_{\rm e}$. The best-fitting model is a cored DM halo which contributes $\sim 10\%$ of the total mass within $1\,R_{\rm e}$, and $67\% \pm 10\%$ within $6\,R_{\rm e}$, although a cusped DM halo is also acceptable. The red GCs exhibit mild rotation with $v_{\rm max}/\sigma_0 \sim 0.3$ in the region $R > \,R_{\rm e}$, aligned with but counter-rotating to the stars in the inner parts, while the blue GCs and PNe kinematics are consistent with no rotation. The red GCs are tangentially anisotropic, the blue GCs are mildly radially anisotropic, and the PNe vary from radially to tangentially anisotropic from the inner to the outer region. This is confirmed by general made-to-measure models. The tangential anisotropy of the red GCs in the inner regions could stem from the preferential destruction of red GCs on more radial orbits, while their outer tangential anisotropy -- similar to the PNe in this region -- has no good explanation. The mild radial anisotropy of the blue GCs is consistent with an accretion scenario.

Deep observations of the dwarf elliptical (dE) galaxy NGC 1396 (M$_V = -16.60$, Mass $\sim 4\times10^8$ M$_\odot$), located in the Fornax cluster, have been performed with the VLT/ MUSE spectrograph in the wavelength region from $4750-9350$ \AA{}. In this paper we present a stellar population analysis studying chemical abundances, the star formation history (SFH) and the stellar initial mass function (IMF) as a function of galacto-centric distance. Different, independent ways to analyse the stellar populations result in a luminosity-weighted age of $\sim$ 6 Gyr and a metallicity [Fe/H]$\sim$ $-0.4$, similar to other dEs of similar mass. We find unusually overabundant values of [Ca/Fe] $\sim +0.1$, and under-abundant Sodium, with [Na/Fe] values around $-0.1$, while [Mg/Fe] is overabundant at all radii, increasing from $\sim+0.1$ in the centre to $\sim +0.2$ dex. We notice a significant metallicity and age gradient within this dwarf galaxy. To constrain the stellar IMF of NGC 1396, we find that the IMF of NGC 1396 is consistent with either a Kroupa-like or a top-heavy distribution, while a bottom-heavy IMF is firmly ruled out. An analysis of the abundance ratios, and a comparison with galaxies in the Local Group, shows that the chemical enrichment history of NGC 1396 is similar to the Galactic disc, with an extended star formation history. This would be the case if the galaxy originated from a LMC-sized dwarf galaxy progenitor, which would lose its gas while falling into the Fornax cluster.

Galaxies and their dark matter haloes are part of a complex network of mass structures, collectively called the cosmic web. Using the tidal tensor prescription these structures can be classified into four cosmic environments: voids, sheets, filaments and knots. As the cosmic web may influence the formation and evolution of dark matter haloes and the galaxies they host, we aim to study the effect of these cosmic environments on the average mass of galactic haloes. To this end we measure the galaxy-galaxy lensing profile of 91,195 galaxies, within 0.039 < z < 0.263, from the spectroscopic Galaxy And Mass Assembly (GAMA) survey, using ~100 square degrees of overlapping data from the Kilo-Degree Survey (KiDS). In each of the four cosmic environments we model the contributions from group centrals, satellites and neighbouring groups to the stacked galaxy-galaxy lensing profiles. After correcting the lens samples for differences in the stellar mass distribution, we find no dependence of the average halo mass of central galaxies on their cosmic environment. We do find a significant increase in the average contribution of neighbouring groups to the lensing profile in increasingly dense cosmic environments. We show, however, that the observed effect can be entirely attributed to the galaxy density at much smaller scales (within 4 Mpc/h), which is correlated with the density of the cosmic environments. Within our current uncertainties we find no direct dependence of galaxy halo mass on their cosmic environment.

We report the discovery of a complex extended density enhancement in the Globular Clusters (GCs) in the central $\sim 0.5(^{\circ})^2$ ($\sim 0.06$ Mpc$^2$) of the Fornax cluster, corresponding to $\sim 50\%$ of the area within 1 core radius. This overdensity connects the GC system of NGC1399 to most of those of neighboring galaxies within $\sim 0.6^{\circ}$ ($\sim 210$ kpc) along the W-E direction. The asymmetric density structure suggests that the galaxies in the core of the Fornax cluster experienced a lively history of interactions that have left a clear imprint on the spatial distribution of GCs. The extended central dominant structure is more prominent in the distribution of blue GCs, while red GCs show density enhancements that are more centrally concentrated on the host galaxies. We propose that the relatively small-scale density structures in the red GCs are caused by galaxy-galaxy interactions, while the extensive spatial distribution of blue GCs is due to stripping of GCs from the halos of core massive galaxies by the Fornax gravitational potential. Our investigations is based on density maps of candidate GCs extracted from the multi-band VLT Survey Telescope (VST) survey of Fornax (FDS), identified in a three-dimensional color space and further selected based on their $g$-band magnitude and morphology.

[Abrigded] We have started a new deep, multi-imaging survey of the Fornax cluster, dubbed Fornax Deep Survey (FDS), at the VLT Survey Telescope. In this paper we present the deep photometry inside two square degrees around the bright galaxy NGC1399 in the core of the cluster. We found a very extended and diffuse envelope surrounding the luminous galaxy NGC1399: we map the surface brightness out to 33 arcmin (~ 192 kpc) from the galaxy center and down to about 31 mag/arcsec^2 in the g band. The deep photometry allows us to detect a faint stellar bridge in the intracluster region between NGC1399 and NGC1387. By analyzing the integrated colors of this feature, we argue that it could be due to the ongoing interaction between the two galaxies, where the outer envelope of NGC1387 on its east side is stripped away. By fitting the light profile, we found that it exists a physical break radius in the total light distribution at R=10 arcmin (~58 kpc) that sets the transition region between the bright central galaxy and the outer exponential stellar halo. We discuss the main implications of this work on the build-up of the stellar halo at the center of the Fornax cluster. By comparing with the numerical simulations of the stellar halo formation for the most massive BCGs, we find that the observed stellar halo mass fraction is consistent with a halo formed through the multiple accretion of progenitors with a stellar mass in the range 10^8 - 10^11 M_sun. This might suggest that the halo of NGC1399 has also gone through a major merging event. The absence of a significant number of luminous stellar streams and tidal tails out to 192 kpc suggests that the epoch of this strong interaction goes back to an early formation epoch. Therefore, differently from the Virgo cluster, the extended stellar halo around NGC1399 is characterised by a more diffuse and well-mixed component, including the ICL.

The abundance of compact, massive, early-type galaxies (ETGs) provides important constraints to galaxy formation scenarios. Thanks to the area covered, depth, excellent spatial resolution and seeing, the ESO Public optical Kilo Degree Survey (KiDS), carried out with the VLT Survey Telescope (VST), offers a unique opportunity to conduct a complete census of the most compact galaxies in the Universe. This paper presents a first census of such systems from the first 156 square degrees of KiDS. Our analysis relies on g-, r-, and i-band effective radii ($R_{\rm e}$), derived by fitting galaxy images with PSF-convolved S\'ersic models, high-quality photometric redshifts, $z_{\rm phot}$, estimated from machine learning techniques, and stellar masses, $M_{\rm \star}$, calculated from KiDS aperture photometry. After massiveness ($M_{\rm \star} > 8 \times 10^{10}\, \rm M_{\odot}$) and compactness ($R_{\rm e} < 1.5 \, \rm kpc$ in g-, r- and i-bands) criteria are applied, a visual inspection of the candidates plus near-infrared photometry from VIKING-DR1 are used to refine our sample. The final catalog, to be spectroscopically confirmed, consists of 92 systems in the redshift range $z \sim 0.2-0.7$. This sample, which we expect to increase by a factor of ten over the total survey area, represents the first attempt to select massive super-compact ETGs (MSCGs) in KiDS. We investigate the impact of redshift systematics in the selection, finding that this seems to be a major source of contamination in our sample. A preliminary analysis shows that MSCGs exhibit negative internal colour gradients, consistent with a passive evolution of these systems. We find that the number density of MSCGs is only mildly consistent with predictions from simulations at $z>0.2$, while no such system is found at $z < 0.2$.

In this work we analyze the dark matter (DM) fraction, $f_{DM}$, and mass-to-light ratio mismatch parameter, $\delta_{IMF}$ (computed with respect to a Milky-Way-like IMF), for a sample of 39 dwarf early-type galaxies (dEs) in the Virgo cluster. Both $f_{DM}$ and $\delta_{IMF}$ are estimated within the central (one effective radius) galaxy regions, with a Jeans dynamical analysis that relies on galaxy velocity dispersions, structural parameters, and stellar M/L ratios from the SMAKCED survey. In this first attempt to constrain, simultaneously, the IMF normalization and the DM content, we explore the impact of different assumptions on the DM model profile. On average, for a NFW profile, the $\delta_{IMF}$ is consistent with a Chabrier-like normalization ($\delta_{IMF} \sim 1$), with $f_{DM} \sim 0.35$. One of the main results of the present work is that for at least a few systems the $\delta_{IMF}$ is heavier than the MW-like value (i.e. either top- or bottom-heavy). When introducing tangential anisotropy, larger $\delta_{IMF}$ and smaller $f_{DM}$ are derived. Adopting a steeper concentration-mass relation than that from simulations, we find lower $\delta_{IMF}$ ($< 1$) and larger $f_{DM}$. A constant M/L profile with null $f_{DM}$ gives the heaviest $\delta_{IMF}$ ($\sim 2$). In the MONDian framework, we find consistent results to those for our reference NFW model. If confirmed, the large scatter of $\delta_{IMF}$ for dEs would provide (further) evidence for a non-universal IMF in early-type systems. On average, our reference $f_{DM}$ estimates are consistent with those found for low-$\sigma_{e}$ ($\rm \sim 100 \, \rm km s^{-1}$) early-type galaxies (ETGs). Furthermore, we find $f_{DM}$ consistent with values from the SMAKCED survey, and find a double-value behavior of $f_{DM}$ with stellar mass, which mirrors the trend of dynamical M/L and global star formation efficiency with mass.

The Kilo-Degree Survey (KiDS) is a multi-band imaging survey designed for cosmological studies from weak lensing and photometric redshifts. It uses the ESO VLT Survey Telescope with its wide-field camera OmegaCAM. KiDS images are taken in four filters similar to the SDSS ugri bands. The best-seeing time is reserved for deep r-band observations that reach a median 5-sigma limiting AB magnitude of 24.9 with a median seeing that is better than 0.7arcsec. Initial KiDS observations have concentrated on the GAMA regions near the celestial equator, where extensive, highly complete redshift catalogues are available. A total of 109 survey tiles, one square degree each, form the basis of the first set of lensing analyses, which focus on measurements of halo properties of GAMA galaxies. 9 galaxies per square arcminute enter the lensing analysis, for an effective inverse shear variance of 69 per square arcminute. Accounting for the shape measurement weight, the median redshift of the sources is 0.53. KiDS data processing follows two parallel tracks, one optimized for galaxy shape measurement (for weak lensing), and one for accurate matched-aperture photometry in four bands (for photometric redshifts). This technical paper describes how the lensing and photometric redshift catalogues have been produced (including an extensive description of the Gaussian Aperture and Photometry pipeline), summarizes the data quality, and presents extensive tests for systematic errors that might affect the lensing analyses. We also provide first demonstrations of the suitability of the data for cosmological measurements, and explain how the shear catalogues were blinded to prevent confirmation bias in the scientific analyses. The KiDS shear and photometric redshift catalogues, presented in this paper, are released to the community through http://kids.strw.leidenuniv.nl .

We present the results of our first year of quasar search in the on-going ESO public Kilo Degree Survey (KiDS) and VISTA Kilo-Degree Infrared Galaxy (VIKING) surveys. These surveys are among the deeper wide-field surveys that can be used to uncovered large numbers of z~6 quasars. This allows us to probe a more common population of z~6 quasars that is fainter than the well-studied quasars from the main Sloan Digital Sky Survey. From this first set of combined survey catalogues covering ~250 deg^2 we selected point sources down to Z_AB=22 that had a very red i-Z (i-Z>2.2) colour. After follow-up imaging and spectroscopy, we discovered four new quasars in the redshift range 5.8<z<6.0. The absolute magnitudes at a rest-frame wavelength of 1450 A are between -26.6 < M_1450 < -24.4, confirming that we can find quasars fainter than M^*, which at z=6 has been estimated to be between M^*=-25.1 and M^*=-27.6. The discovery of 4 quasars in 250 deg^2 of survey data is consistent with predictions based on the z~6 quasar luminosity function. We discuss various ways to push the candidate selection to fainter magnitudes and we expect to find about 30 new quasars down to an absolute magnitude of M_1450=-24. Studying this homogeneously selected faint quasar population will be important to gain insight into the onset of the co-evolution of the black holes and their stellar hosts.

We use the first 100 sq. deg. of overlap between the Kilo-Degree Survey (KiDS) and the Galaxy And Mass Assembly (GAMA) survey to determine the galaxy halo mass of ~10,000 spectroscopically-confirmed satellite galaxies in massive ($M > 10^{13}h^{-1}{\rm M}_\odot$) galaxy groups. Separating the sample as a function of projected distance to the group centre, we jointly model the satellites and their host groups with Navarro-Frenk-White (NFW) density profiles, fully accounting for the data covariance. The probed satellite galaxies in these groups have total masses $\log M_{\rm sub} /(h^{-1}{\rm M}_\odot) \approx 11.7 - 12.2$ consistent across group-centric distance within the errorbars. Given their typical stellar masses, $\log M_{\rm \star,sat}/(h^{-2}{\rm M}_\odot) \sim 10.5$, such total masses imply stellar mass fractions of $M_{\rm \star,sat} /M_{\rm sub} \approx 0.04 h^{-1}$ . The average subhalo hosting these satellite galaxies has a mass $M_{\rm sub} \sim 0.015M_{\rm host}$ independent of host halo mass, in broad agreement with the expectations of structure formation in a $\Lambda$CDM universe.

The Shapley Supercluster Survey is a multi-wavelength survey covering an area of ~23 deg^2 (~260 Mpc^2 at z=0.048) around the supercluster core, including nine Abell and two poor clusters, having redshifts in the range 0.045-0.050. The survey aims to investigate the role of the cluster-scale mass assembly on the evolution of galaxies, mapping the effects of the environment from the cores of the clusters to their outskirts and along the filaments. The optical (ugri) imaging acquired with OmegaCAM on the VLT Survey Telescope is essential to achieve the project goals providing accurate multi-band photometry for the galaxy population down to m*+6. We describe the methodology adopted to construct the optical catalogues and to separate extended and point-like sources. The catalogues reach average 5sigma limiting magnitudes within a 3\arcsec diameter aperture of ugri=[24.4,24.6,24.1,23.3] and are 93% complete down to ugri=[23.8,23.8,23.5,22.0] mag, corresponding to ~m*_r+8.5. The data are highly uniform in terms of observing conditions and all acquired with seeing less than 1.1 arcsec full width at half-maximum. The median seeing in r-band is 0.6 arcsec, corresponding to 0.56 kpc h^{-1}_{70} at z=0.048. While the observations in the u, g and r bands are still ongoing, the i-band observations have been completed, and we present the i-band catalogue over the whole survey area. The latter is released and it will be regularly updated, through the use of the Virtual Observatory tools. This includes 734,319 sources down to i=22.0 mag and it is the first optical homogeneous catalogue at such a depth, covering the central region of the Shapley supercluster.

The Kilo-Degree Survey (KiDS) is an optical wide-field imaging survey carried out with the VLT Survey Telescope and the OmegaCAM camera. KiDS will image 1500 square degrees in four filters (ugri), and together with its near-infrared counterpart VIKING will produce deep photometry in nine bands. Designed for weak lensing shape and photometric redshift measurements, the core science driver of the survey is mapping the large-scale matter distribution in the Universe back to a redshift of ~0.5. Secondary science cases are manifold, covering topics such as galaxy evolution, Milky Way structure, and the detection of high-redshift clusters and quasars. KiDS is an ESO Public Survey and dedicated to serving the astronomical community with high-quality data products derived from the survey data, as well as with calibration data. Public data releases will be made on a yearly basis, the first two of which are presented here. For a total of 148 survey tiles (~160 sq.deg.) astrometrically and photometrically calibrated, coadded ugri images have been released, accompanied by weight maps, masks, source lists, and a multi-band source catalog. A dedicated pipeline and data management system based on the Astro-WISE software system, combined with newly developed masking and source classification software, is used for the data production of the data products described here. The achieved data quality and early science projects based on the data products in the first two data releases are reviewed in order to validate the survey data. Early scientific results include the detection of nine high-z QSOs, fifteen candidate strong gravitational lenses, high-quality photometric redshifts and galaxy structural parameters for hundreds of thousands of galaxies. (Abridged)

We estimated photometric redshifts (zphot) for more than 1.1 million galaxies of the ESO Public Kilo-Degree Survey (KiDS) Data Release 2. KiDS is an optical wide-field imaging survey carried out with the VLT Survey Telescope (VST) and the OmegaCAM camera, which aims at tackling open questions in cosmology and galaxy evolution, such as the origin of dark energy and the channel of galaxy mass growth. We present a catalogue of photometric redshifts obtained using the Multi Layer Perceptron with Quasi Newton Algorithm (MLPQNA) model, provided within the framework of the DAta Mining and Exploration Web Application REsource (DAMEWARE). These photometric redshifts are based on a spectroscopic knowledge base which was obtained by merging spectroscopic datasets from GAMA (Galaxy And Mass Assembly) data release 2 and SDSS-III data release 9. The overall 1 sigma uncertainty on Delta z = (zspec - zphot) / (1+ zspec) is ~ 0.03, with a very small average bias of ~ 0.001, a NMAD of ~ 0.02 and a fraction of catastrophic outliers (| Delta z | > 0.15) of ~0.4%.

This work makes use of the VST observations to select variable sources. We use also the IR photometry, SED fitting and X-ray information where available to confirm the nature of the AGN candidates. The IR data, available over the full survey area, allow to confirm the consistency of the variability selection with the IR color selection method, while the detection of variability may prove useful to detect the presence of an AGN in IR selected starburst galaxies.

We present the VST Early-type GAlaxy Survey (VEGAS), which is designed to obtain deep multiband photometry in g, r, i, of about one hundred nearby galaxies down to 27.3, 26.8, and 26 mag/arcsec^2 respectively, using the ESO facility VST/OmegaCAM. The goals of the survey are 1) to map the light distribution up to ten effective radii, r_e, 2) to trace color gradients and surface brightness fluctuation gradients out to a few r_e for stellar population characterization, and 3) to obtain a full census of the satellite systems (globular clusters and dwarf galaxies) out to 20% of the galaxy virial radius. The external regions of galaxies retain signatures of the formation and evolution mechanisms that shaped them, and the study of nearby objects enables a detailed analysis of their morphology and interaction features. To clarify the complex variety of formation mechanisms of early-type galaxies (ETGs), wide and deep photometry is the primary observational step, which at the moment has been pursued with only a few dedicated programs. The VEGAS survey has been designated to provide these data for a volume-limited sample with exceptional image quality. In this commissioning photometric paper we illustrate the capabilities of the survey using g- and i-band VST/OmegaCAM images of the nearby galaxy NGC 4472 and of smaller ETGs in the surrounding field. Our surface brightness profiles reach rather faint levels and agree excellently well with previous literature. Genuine new results concern the detection of an intracluster light tail in NGC 4472 and of various substructures at increasing scales. We have also produced extended (g-i) color profiles. The VST/OmegaCAM data that we acquire in the context of the VEGAS survey provide a detailed view of substructures in the optical emission from extended galaxies, which can be as faint as a hundred times below the sky level.

One of the most peculiar characteristics of Active Galactic Nuclei (AGN) is their variability over all wavelengths. This property has been used in the past to select AGN samples and is foreseen to be one of the detection techniques applied in future multi-epoch surveys, complementing photometric and spectroscopic methods. In this paper, we aim to construct and characterise an AGN sample using a multi-epoch dataset in the r band from the SUDARE-VOICE survey. Our work makes use of the VST monitoring program of an area surrounding the Chandra Deep Field South to select variable sources. We use data spanning a six month period over an area of 2 square degrees, to identify AGN based on their photometric variability. The selected sample includes 175 AGN candidates with magnitude r < 23 mag. We distinguish different classes of variable sources through their lightcurves, as well as X-ray, spectroscopic, SED, optical and IR information overlapping with our survey. We find that 12% of the sample (21/175) is represented by SN. Of the remaining sources, 4% (6/154) are stars, while 66% (102/154) are likely AGNs based on the available diagnostics. We estimate an upper limit to the contamination of the variability selected AGN sample of about 34%, but we point out that restricting the analysis to the sources with available multi-wavelength ancillary information, the purity of our sample is close to 80% (102 AGN out of 128 non-SN sources with multi-wavelength diagnostics). Our work thus confirms the efficiency of the variability selection method in agreement with our previous work on the COSMOS field; in addition we show that the variability approach is roughly consistent with the infrared selection.

We study the total density distribution in the central regions ($<\, 1$ effective radius, $R_{\rm e}$) of early-type galaxies (ETGs), using data from the SPIDER survey. We model each galaxy with two components (dark matter halo + stars), exploring different assumptions for the dark matter (DM) halo profile, and leaving stellar mass-to-light ($M_{\rm \star}/L$) ratios as free fitting parameters to the data. For a Navarro et al. (1996) profile, the slope of the total mass profile is non-universal. For the most massive and largest ETGs, the profile is isothermal in the central regions ($\sim R_{\rm e}/2$), while for the low-mass and smallest systems, the profile is steeper than isothermal, with slopes similar to those for a constant-M/L profile. For a concentration-mass relation steeper than that expected from simulations, the correlation of density slope with mass tends to flatten. Our results clearly point to a "non-homology" in the total mass distribution of ETGs, which simulations of galaxy formation suggest may be related to a varying role of dissipation with galaxy mass.

The Kilo Degree survey (KiDS) is a large-scale optical imaging survey carried out with the VLT Survey Telescope (VST), which is the ideal tool for galaxy evolution studies. We expect to observe millions of galaxies for which we extract the structural parameters in four wavebands (u, g, r and i). This sample will represent the largest dataset with measured structural parameters up to a redshift $z=0.5$. In this paper we will introduce the sample, and describe the 2D fitting procedure using the 2DPHOT environment and the validation of the parameters with an external catalog.

We have started a systematic search of strong lens candidates in the ESO public survey KiDS based on the visual inspection of massive galaxies in the redshift range $0.1<z<0.5$. As a pilot program we have inspected 100 sq. deg., which overlap with SDSS and where there are known lenses to use as a control sample. Taking advantage of the superb image quality of VST/OmegaCAM, the colour information and accurate model subtracted images, we have found 18 new lens candidates, for which spectroscopic confirmation will be needed to confirm their lensing nature and study the mass profile of the lensing galaxies.

The ESO Public Kilo-Degree Survey (KiDS) is an optical wide-field imaging survey carried out with the VLT Survey Telescope and the OmegaCAM camera. KiDS will scan 1500 square degrees in four optical filters (u, g, r, i). Designed to be a weak lensing survey, it is ideal for galaxy evolution studies, thanks to the high spatial resolution of VST, the good seeing and the photometric depth. The surface photometry have provided with structural parameters (e.g. size and S\'ersic index), aperture and total magnitudes have been used to derive photometric redshifts from Machine learning methods and stellar masses/luminositites from stellar population synthesis. Our project aimed at investigating the evolution of the colour and structural properties of galaxies with mass and environment up to redshift $z \sim 0.5$ and more, to put constraints on galaxy evolution processes, as galaxy mergers.

By means of our own cosmological-hydrodynamical simulation and semi-analytical model we studied galaxy population properties in clusters and groups, spanning over 10 different bands from UV to NIR, and their evolution since redshift z=2. We compare our results in terms of galaxy red/blue fractions and luminous-to-faint ratio (LFR) on the Red Sequence (RS) with recent observational data reaching beyond z=1.5. Different selection criteria were tested in order to retrieve galaxies belonging to the RS: either by their quiescence degree measured from their specific SFR ("Dead Sequence"), or by their position in a colour-colour plane which is also a function of sSFR. In both cases, the colour cut and the limiting magnitude threshold were let evolving with redshift, in order to follow the natural shift of the characteristic luminosity in the LF. We find that the Butcher-Oemler effect is wavelength-dependent, with the fraction of blue galaxies increasing steeper in optical colours than in NIR. Besides, only when applying a lower limit in terms of fixed absolute magnitude, a steep BO effect can be reproduced, while the blue fraction results less evolving when selecting samples by stellar mass or an evolving magnitude limit. We then find that also the RS-LFR behaviour, highly debated in the literature, is strongly dependent on the galaxy selection function: in particular its very mild evolution recovered when measured in terms of stellar mass, is in agreement with values reported for some of the highest redshift confirmed (proto)clusters. As to differences through environments, we find that normal groups and (to a lesser extent) cluster outskirts present the highest values of both star forming fraction and LFR at low z, while fossil groups and cluster cores the lowest: this separation among groups begins after z~0.5, while earlier all group star forming properties are undistinguishable.

We present an integral field study of the internal structure, kinematics and stellar population of the almost edge-on, intermediate luminosity ($L_ {*}$) elliptical galaxy NGC 4697. We build extended 2-dimensional (2D) maps of the stellar kinematics and line-strengths of the galaxy up to $\sim 0.7 $ effective radii (R$_{eff}$) using a mosaic of 8 VIMOS (VIsible Multi-Objects Spectrograph on the VLT) integral-field unit pointings. We find clear evidence for a rotation-supported structure along the major axis from the 2D kinematical maps, confirming the previous classification of this system as a `fast-rotator'. We study the correlations between the third and fourth Gauss-Hermite moments of the line-of-sight velocity distribution (LOSVD) $h_3$ and $h_4$ with the rotation parameter ($V/\sigma$), and compare our findings to hydrodynamical simulations. We find remarkable similarities to predictions from gas-rich mergers. Based on photometry, we perform a bulge/disk decomposition and study the stellar population properties of the two components. The bulge and the disk show different stellar populations, with the stars in the bulge being older (age$_{\rm bulge}=13.5^{+1.4}_{-1.4}$ Gyr, age$_{\rm disk}=10.5^{+1.6}_{-2.0}$Gyr) and more metal-poor ($\mathrm{[M/H]_{bulge}} = -0.17^{+0.12}_{-0.1}$, $\mathrm{[M/H]_{disk}}=-0.03^{+0.02}_{-0.1}$). The evidence of a later-formed, more metal-rich disk embedded in an older, more metal-poor bulge, together with the LOSVD structure, supports a mass assembly scenario dominated by gas-rich minor mergers and possibly with a late gas-rich major merger that left a previously rapidly rotating system unchanged. The bulge and the disk do not show signs of different stellar Initial Mass Function slopes, and both match well with a Milky Way-like IMF.

We perform in-depth dynamical modelling of the luminous and dark matter (DM) content of the elliptical galaxy NGC 1407. Our strategy consists of solving the spherical Jeans equations for three independent dynamical tracers: stars, blue GCs and red GCs in a self-consistent manner. We adopt a maximum-likelihood Markov-Chain Monte Carlo fitting technique in the attempt to constrain the inner slope of the DM density profile (the cusp/core problem), and the stellar initial mass function (IMF) of the galaxy. We find the inner logarithmic slope of the DM density profiles to be $\gamma = 0.6\pm0.4$, which is consistent with either a DM cusp ($\gamma = 1$) or with a DM core $(\gamma = 0)$. Our findings are consistent with a Salpeter IMF, and marginally consistent with a Kroupa IMF. We infer tangential orbits for the blue GCs, and radial anisotropy for red GCs and stars. The modelling results are consistent with the virial mass--concentration relation predicted by $\Lambda$CDM simulations. The virial mass of NGC 1407 is $\log$ $M_{\rm vir} = 13.3 \pm 0.2 M_{\odot}$, whereas the stellar mass is $\log M_* = 11.8 \pm 0.1 M_{\odot}$. The overall uncertainties on the mass of NGC 1407 are only 5 per cent at the projected stellar effective radius. We attribute the disagreement between our results and previous X-ray results to the gas not being in hydrostatic equilibrium in the central regions of the galaxy. The halo of NGC 1407 is found be DM dominated, with a dynamical mass-to-light ratio of $M/L=260_{-100} ^{+174} M_{\odot}/L_{\odot, B}$. However, this value can be larger up to a factor of 3 depending on the assumed prior on the DM scale radius.

Optical variability has proven to be an effective way of detecting AGNs in imaging surveys, lasting from weeks to years. In the present work we test its use as a tool to identify AGNs in the VST multi-epoch survey of the COSMOS field, originally tailored to detect supernova events. We make use of the multi-wavelength data provided by other COSMOS surveys to discuss the reliability of the method and the nature of our AGN candidates. Our selection returns a sample of 83 AGN candidates; based on a number of diagnostics, we conclude that 67 of them are confirmed AGNs (81% purity), 12 are classified as supernovae, while the nature of the remaining 4 is unknown. For the subsample of AGNs with some spectroscopic classification, we find that Type 1 are prevalent (89%) compared to Type 2 AGNs (11%). Overall, our approach is able to retrieve on average 15% of all AGNs in the field identified by means of spectroscopic or X-ray classification, with a strong dependence on the source apparent magnitude. In particular, the completeness for Type 1 AGNs is 25%, while it drops to 6% for Type 2 AGNs. The rest of the X-ray selected AGN population presents on average a larger r.m.s. variability than the bulk of non variable sources, indicating that variability detection for at least some of these objects is prevented only by the photometric accuracy of the data. We show how a longer observing baseline would return a larger sample of AGN candidates. Our results allow us to assess the usefulness of this AGN selection technique in view of future wide-field surveys.

We present a study of globular clusters (GCs) and other small stellar systems (SSSs) in the field of NGC3115, observed as part of the VEGAS imaging survey, carried out with the VST telescope. We use deep g and i data of NGC3115, a well-studied lenticular galaxy with excellent scientific literature. This is fundamental to test the methodologies, verify the results, and probe the capabilities of the VEGAS-SSS. Leveraging the large field of view of the VST allow us to accurately study of the distribution and properties of SSSs as a function of galactocentric distance Rgc, well beyond ~20 galaxy effective radii, in a way not often possible. Our analysis of colors, magnitudes and sizes of SSS candidates confirms the results from existing studies, some of which carried out with 8-10m class telescopes, and further extends them to unreached Rgc distances, with similar accuracy. We find a color bimodality for the GC population and a r1/4 profile for the surface density of GCs as for the galaxy light profile. The radial color gradient of blue and red GCs previously found, for instance by the SLUGGS survey, is further extended out to the largest Rgc inspected, ~65 kpc. The surface density profiles of blue and red GCs taken separately are well approximated by a r1/4 density profile, with the fraction of blue GCs being slightly larger at larger Rgc. We do not find hints of a trend for the red GC and for the GC turnover magnitude to vary with radius, but we observe a ~0.2 mag difference in the turnover magnitude of the blue and red GCs subpopulations. Inspecting SSS sizes and colors we obtained a list of UCDs and GC candidates suitable for future spectroscopic follow-up. In conclusion, our study shows the reliability of the methodologies developed to study SSSs in the field of bright early-type galaxies, and the great potential of the VEGAS survey to produce original results on SSSs science.

We present the first deep photometry of a good candidate for a forming polar ring galaxy at redshift z~0.05. This object, named FCSS J033710.0-354727, is a background galaxy in the VST deep field of the Fornax cluster. The deep exposures combined with the high angular resolution of the OmegaCAM at VST allow us to carry out the first detailed photometric analysis for this system in the g and i bands to derive the galaxy structure and colors. Results show that the central object resembles a disk galaxy, surrounded by a ring-like structure 2 times more extended than the central disk. The warped geometry and the presence of bright knots observed along the polar direction, as well as the several debris detected on the NW side with colors comparable to those of the galaxy, suggest that the polar structure is still forming. We argue that the wide polar ring/disk is the result of the ongoing disruption of a companion galaxy in the potential of the central object, which is 2-3 times more massive than the accreting galaxy.

We present an overview of a multi-wavelength survey of the Shapley supercluster (SSC; z~0.05) covering a contiguous area of 260 h^-2_70 Mpc^2 including the supercluster core. The project main aim is to quantify the influence of cluster-scale mass assembly on galaxy evolution in one of the most massive structures in the local Universe. The Shapley supercluster survey (ShaSS) includes nine Abell clusters (A3552, A3554, A3556, A3558, A3559, A3560, A3562, AS0724, AS0726) and two poor clusters (SC1327- 312, SC1329-313) showing evidence of cluster-cluster interactions. Optical (ugri) and near-infrared (K) imaging acquired with VST and VISTA allow us to study the galaxy population down to m*+6 at the supercluster redshift. A dedicated spectroscopic survey with AAOmega on the Anglo-Australian Telescope provides a magnitude-limited sample of supercluster members with 80% completeness at ~m*+3. We derive the galaxy density across the whole area, demonstrating that all structures within this area are embedded in a single network of clusters, groups and filaments. The stellar mass density in the core of the SSC is always higher than 9E09 M_sun Mpc^-3, which is ~40x the cosmic stellar mass density for galaxies in the local Universe. We find a new filamentary structure (~7 Mpc long in projection) connecting the SSC core to the cluster A3559, as well as previously unidentified density peaks. We perform a weak-lensing analysis of the central 1 sqdeg field of the survey obtaining for the central cluster A3558 a mass of M_500=7.63E14 M_sun, in agreement with X-ray based estimates.

We present results from SPH-cosmological simulations, including self-consistent modelling of SN feedback and chemical evolution, of galaxies belonging to two clusters and twelve groups. We reproduce the mass-metallicity (ZM) relation of galaxies classified in two samples according to their star-forming activity, as parametrized by their sSFR, across a redshift range up to z=2. Its slope shows irrelevant evolution in the passive sample, being steeper in groups than in clusters. However, the sub-sample of high-mass passive galaxies only is characterized by a steep increase of the slope with redshift, from which it can be inferred that the bulk of the slope evolution of the ZM relation is driven by the more massive passive objects. (...ABRIDGED...) The ZM relation for the star-forming sample reveals an increasing scatter with redshift, indicating that it is still being built at early epochs. The star-forming galaxies make up a tight sequence in the SFR-M_* plane at high redshift, whose scatter increases with time alongside with the consolidation of the passive sequence. We also confirm the anti-correlation between sSFR and stellar mass, pointing at a key role of the former in determining the galaxy downsizing, as the most significant means of diagnostics of the star formation efficiency. Likewise, an anti-correlation between sSFR and metallicity can be established for the star-forming galaxies, while on the contrary more active galaxies in terms of simple SFR are also metal-richer. We discuss these results in terms of the mechanisms driving the evolution within the high- and low-mass regimes at different epochs: mergers, feedback-driven outflows and the intrinsic variation of the star formation efficiency.

We investigate the evolution of dark and luminous matter in the central regions of early-type galaxies (ETGs) up to z ~ 0.8. We use a spectroscopically selected sample of 154 cluster and field galaxies from the EDisCS survey, covering a wide range in redshifts (z ~ 0.4-0.8), stellar masses ($\log M_{\star}/ M_{\odot}$ ~ 10.5-11.5 dex) and velocity dispersions ($\sigma_{\star}$ ~ 100-300 \, km/s). We obtain central dark matter (DM) fractions by determining the dynamical masses from Jeans modelling of galaxy aperture velocity dispersions and the $M_{\star}$ from galaxy colours, and compare the results with local samples. We discuss how the correlations of central DM with galaxy size (i.e. the effective radius, $R_{\rm e}$), $M_{\star}$ and $\sigma_{\star}$ evolve as a function of redshift, finding clear indications that local galaxies are, on average, more DM dominated than their counterparts at larger redshift. This DM fraction evolution with $z$ can be only partially interpreted as a consequence of the size-redshift evolution. We discuss our results within galaxy formation scenarios, and conclude that the growth in size and DM content which we measure within the last 7 Gyr is incompatible with passive evolution, while it is well reproduced in the multiple minor merger scenario. We also discuss the impact of the IMF on our DM inferences and argue that this can be non-universal with the lookback time. In particular, we find the Salpeter IMF can be better accommodated by low redshift systems, while producing stellar masses at high-$z$ which are unphysically larger than the estimated dynamical masses (particularly for lower-$\sigma_{\star}$ systems).

We study the total density distribution in the central regions (~ 1 effective radius, $R_e$) of early-type galaxies (ETGs), using data from SPIDER and $\rm ATLAS^{3D}$. Our analysis extends the range of galaxy stellar mass ($M_{\star}$) probed by gravitational lensing, down to ~ $10^{10}\, \rm M_{\odot}$. We model each galaxy with two components (dark matter halo + stars), exploring different assumptions for the dark matter (DM) halo profile (i.e. NFW, NFW-contracted, and Burkert profiles), and leaving stellar mass-to-light ($M_{\star}/L$) ratios as free fitting parameters to the data. For all plausible halo models, the best-fitting $M_{\star}/L$, normalized to that for a Chabrier IMF, increases systematically with galaxy size and mass. For an NFW profile, the slope of the total mass profile is non-universal, independently of several ingredients in the modeling (e.g., halo contraction, anisotropy, and rotation velocity in ETGs). For the most massive ($M_{\star}$ ~ $10^{11.5} \, M_{\odot}$) or largest ($R_{\rm e}$ ~ $15 \, \rm kpc$) ETGs, the profile is isothermal in the central regions (~$R_{\rm e}/2$), while for the low-mass ($M_{\star}$ ~ $10^{10.2} \, M_\odot$) or smallest ($R_{\rm e}$ ~ $0.5 \, \rm kpc$) systems, the profile is steeper than isothermal, with slopes similar to those for a constant-$M/L$ profile. For a steeper concentration-mass relation than that expected from simulations, the correlation of density slope with galaxy mass tends to flatten, while correlations with $R_{\rm e}$ and velocity dispersions are more robust. Our results clearly point to a "non-homology" in the total mass distribution of ETGs, which simulations of galaxy formation suggest may be related to a varying role of dissipation with galaxy mass.

STREGA (STRucture and Evolution of the GAlaxy) is a Guaranteed Time survey being performed at the VST (the ESO VLT Survey Telescope) to map about 150 square degrees in the Galactic halo, in order to constrain the mechanisms of galactic formation and evolution. The survey is built as a five-year project, organized in two parts: a core program to explore the surrounding regions of selected stellar systems and a second complementary part to map the southern portion of the Fornax orbit and extend the observations of the core program. The adopted stellar tracers are mainly variable stars (RR~Lyraes and Long Period Variables) and Main Sequence Turn-off stars for which observations in the g,r,i bands are obtained. We present an overview of the survey and some preliminary results for three observing runs that have been completed. For the region centered on $\omega$~Cen (37 deg^2), covering about three tidal radii, we also discuss the detected stellar density radial profile and angular distribution, leading to the identification of extratidal cluster stars. We also conclude that the cluster tidal radius is about 1.2 deg, in agreement with values in the literature based on the Wilson model.

STEP (the SMC in Time: Evolution of a Prototype interacting late-type dwarf galaxy) is a Guaranteed Time Observation survey being performed at the VST (the ESO VLT Survey Telescope). STEP will image an area of 74 deg$^2$ covering the main body of the Small Magellanic Cloud (32 deg$^2$), the Bridge that connects it to the Large Magellanic Cloud (30 deg$^2$) and a small part of the Magellanic Stream (2 deg$^2$). Our $g,r,i,H_{\alpha}$ photometry is able to resolve individual stars down to magnitudes well below the main-sequence turnoff of the oldest populations. In this first paper we describe the observing strategy, the photometric techniques, and the upcoming data products of the survey. We also present preliminary results for the first two fields for which data acquisition is completed, including some detailed analysis of the two stellar clusters IC\,1624 and NGC\,419.

(Abridged) Fossil systems are group- or cluster-sized objects whose luminosity is dominated by a very massive central galaxy. In the current cold dark matter scenario, these objects formed hierarchically at an early epoch of the Universe and then slowly evolved until present day. That is the reason why they are called {\it fossils}. We started an extensive observational program to characterize a sample of 34 fossil group candidates spanning a broad range of physical properties. Deep $r-$band images were taken for each candidate and optical spectroscopic observations were obtained for $\sim$ 1200 galaxies. This new dataset was completed with SDSS DR7 archival data to obtain robust cluster membership and global properties of each fossil group candidate. For each system, we recomputed the magnitude gaps between the two brightest galaxies ($\Delta m_{12}$) and the first and fourth ranked galaxies ($\Delta m_{14}$) within 0.5 $R_{{\rm 200}}$. We consider fossil systems those with $\Delta m_{12} \ge 2$ mag or $\Delta m_{14} \ge 2.5$ mag within the errors. We find that 15 candidates turned out to be fossil systems. Their observational properties agree with those of non-fossil systems. Both follow the same correlations, but fossils are always extreme cases. In particular, they host the brightest central galaxies and the fraction of total galaxy light enclosed in the central galaxy is larger in fossil than in non-fossil systems. Finally, we confirm the existence of genuine fossil clusters. Combining our results with others in the literature, we favor the merging scenario in which fossil systems formed due to mergers of $L^\ast$ galaxies. The large magnitude gap is a consequence of the extreme merger ratio within fossil systems and therefore it is an evolutionary effect. Moreover, we suggest that at least one candidate in our sample could represent a transitional fossil stage.

We study the mass and anisotropy distribution of the giant elliptical galaxy NGC 5846 using stars, as well as the red and blue globular cluster (GC) subpopulations. We break degeneracies in the dynamical models by taking advantage of the different phase space distributions of the two GC subpopulations to unambiguously constrain the mass of the galaxy and the anisotropy of the GC system. Red GCs show the same spatial distribution and behaviour as the starlight, whereas blue GCs have a shallower density profile, a larger velocity dispersion and a lower kurtosis, all of which suggest a different orbital distribution. We use a dispersion-kurtosis Jeans analysis and find that the solutions of separate analyses for the two GC subpopulations overlap in the halo parameter space. The solution converges on a massive dark matter halo, consistent with expectations from $\Lambda$CDM and WMAP7 cosmology in terms of virial mass ($\log M_{DM} \sim13.3 M_{sun}$) and concentration ($c_{vir}\sim8$). This is the first such analysis that solves the dynamics of the different GC subpopulations in a self-consistent manner. Our method improves the uncertainties on the halo parameter determination by a factor of two and opens new avenues for the use of elliptical galaxy dynamics as tests of predictions from cosmological simulations. The implied stellar mass-to-light ratio derived from the dynamical modelling is fully consistent with a Salpeter initial mass function (IMF) and rules out a bottom light IMF. The different GC subpopulations show markedly distinct orbital distributions at large radii, with red GCs having an anisotropy parameter $\beta\sim0.4$ outside $\sim3R_e$, and the blue GCs having $\beta\sim0.15$ at the same radii, while centrally ($\sim1R_e$) they are both isotropic. We discuss the implications of our findings within the two-phase formation scenario for early-type galaxies.

MOdified Newtonian dynamics (MOND) represents a phenomenological alternative to dark matter (DM) for the missing mass problem in galaxies and clusters of galaxies. We analyze the central regions of a local sample of $\sim 220$ early-type galaxies from the $\rm ATLAS^{3D}$ survey, to see if the data can be reproduced without recourse to DM. We estimate dynamical masses in the MOND context through Jeans analysis, and compare to $\rm ATLAS^{3D}$ stellar masses from stellar population synthesis. We find that the observed stellar mass--velocity dispersion relation is steeper than expected assuming MOND with a fixed stellar initial mass function (IMF) and a standard value for the acceleration parameter $a_{\rm 0}$. Turning from the space of observables to model space, a) fixing the IMF, a universal value for $a_{\rm 0}$ cannot be fitted, while, b) fixing $a_{\rm 0}$ and leaving the IMF free to vary, we find that it is "lighter" (Chabrier-like) for low-dispersion galaxies, and "heavier" (Salpeter-like) for high dispersions. This MOND-based trend matches inferences from Newtonian dynamics with DM, and from detailed analysis of spectral absorption lines, adding to the converging lines of evidence for a systematically-varying IMF.

Aims. We investigated the gravitational effects of a scalar field within scalar-tensor gravity as an alternative to dark matter. Motivated by chameleon, symmetron and f(R)-gravity models, we studied a phenomenological scenario where the scalar field has both a mass (i.e. interaction length) and a coupling constant to the ordinary matter which scale with the local properties of the considered astrophysical system. Methods. We analysed the feasibility of this scenario using the modified gravitational potential obtained in its context and applied it to the galactic and hot gas/stellar dynamics in galaxy clusters and elliptical/spiral galaxies respectively. This is intended to be a first step in assessing the viability of this new approach in the context of "alternative gravity" models. Results. The main results are: 1. the velocity dispersion of elliptical galaxies can be fitted remarkably well by the suggested scalar field, with model significance similar to a classical Navarro-Frenk-White dark halo profile; 2. the analysis of the stellar dynamics and the gas equilibrium in elliptical galaxies has shown that the scalar field can couple with ordinary matter with different strengths (different coupling constants) producing and/or depending on the different clustering state of matter components; 3. elliptical and spiral galaxies, combined with clusters of galaxies, show evident correlations among theory parameters which suggest the general validity of our results at all scales and a way toward a possible unification of the theory for all types of gravitational systems we considered. All these results demonstrate that the proposed scalar field scenario can work fairly well as an alternative to dark matter.

We analyze the evolution of colour gradients predicted by the hydrodynamical models of early type galaxies (ETGs) in Pipino et al. (2008), which reproduce fairly well the chemical abundance pattern and the metallicity gradients of local ETGs. We convert the star formation (SF) and metal content into colours by means of stellar population synthetic model and investigate the role of different physical ingredients, as the initial gas distribution and content, and eps_SF, i.e. the normalization of SF rate. From the comparison with high redshift data, a full agreement with optical rest-frame observations at z < 1 is found, for models with low eps_SF, whereas some discrepancies emerge at 1 < z < 2, despite our models reproduce quite well the data scatter at these redshifts. To reconcile the prediction of these high eps_SF systems with the shallower colour gradients observed at lower z we suggest intervention of 1-2 dry mergers. We suggest that future studies should explore the impact of wet galaxy mergings, interactions with environment, dust content and a variation of the Initial Mass Function from the galactic centers to the peripheries.

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.

Given a flurry of recent claims for systematic variations in the stellar initial mass function (IMF), we carry out the first inventory of the observational evidence using different approaches. This includes literature results, as well as our own new findings from combined stellar-populations synthesis (SPS) and Jeans dynamical analyses of data on $\sim$~4500 early-type galaxies (ETGs) from the SPIDER project. We focus on the mass-to-light ratio mismatch relative to the Milky Way IMF, \dimf, correlated against the central stellar velocity dispersion, \sigs. We find a strong correlation between \dimf\ and \sigs, for a wide set of dark matter (DM) model profiles. These results are robust if a uniform halo response to baryons is adopted across the sample. The overall normalization of \dimf, and the detailed DM profile, are less certain, but the data are consistent with standard cold-DM halos, and a central DM fraction that is roughly constant with \sigs. For a variety of related studies in the literature, using SPS, dynamics, and gravitational lensing, similar results are found. Studies based solely on spectroscopic line diagnostics agree on a Salpeter-like IMF at high \sigs, but differ at low \sigs. Overall, we find that multiple independent lines of evidence appear to be converging on a systematic variation in the IMF, such that high-\sigs\ ETGs have an excess of low-mass stars relative to spirals and low-\sigs\ ETGs. Robust verification of super-Salpeter IMFs in the highest-\sigs\ galaxies will require additional scrutiny of scatter and systematic uncertainties. The implications for the distribution of DM are still inconclusive.

There are only a few tracers available to probe the kinematics of individual early-type galaxies beyond one effective radius. Here we directly compare a sample of planetary nebulae (PNe), globular clusters (GCs) and galaxy starlight velocities out to ~4 effective radii, in the S0 galaxy NGC 2768. Using a bulge-to-disk decomposition of a K-band image we assign PNe and starlight to either the disk or the bulge. We show that the bulge PNe and bulge starlight follow the same radial density distribution as the red subpopulation of GCs, whereas the disk PNe and disk starlight are distinct components. We find good kinematic agreement between the three tracers to several effective radii (and with stellar data in the inner regions). Further support for the distinct nature of the two galaxy components come from our kinematic analysis. After separating the tracers into bulge and disk components we find the bulge to be a slowly rotating pressure-supported system, whereas the disk reveals a rapidly rising rotation curve with a declining velocity dispersion profile. The resulting V/sigma ratio for the disk resembles that of a spiral galaxy and hints at an origin for NGC 2768 as a transformed late-type galaxy. A two-component kinematic analysis for a sample of S0s will help to elucidate the nature of this class of galaxy.

We analyze the central dark-matter (DM) content of \sim 4,500 massive (M* \gsim 10^{10} Msun), low-redshift (z<0.1), early-type galaxies (ETGs), with high-quality ugrizYJHK photometry and optical spectroscopy from SDSS and UKIDSS. We estimate the "central" fraction of DM within the K-band effective radius, \Re, using spherically symmetric isotropic galaxy models. We discuss the role of systematics. The main results of the present work are the following: (1) DM fractions increase systematically with both structural parameters and mass proxies, as in previous studies, and decrease with central stellar density. 2) All correlations involving DM fractions are caused by two fundamental ones with galaxy effective radius and central velocity dispersion. These correlations are independent of each other, so that ETGs populate a central-DM plane (DMP), i.e. a correlation among fraction of total-to-stellar mass, effective radius, and velocity dispersion, whose scatter along the total-to-stellar mass axis amounts to \sim 0.15 dex. (3) In general, under the assumption of an isothermal or a constant M/L profile for the total mass distribution, a Chabrier IMF is favoured with respect to a bottom-heavier Salpeter IMF, as the latter produces negative (i.e. unphysical) DM fractions for more than 50% of the galaxies in our sample. For a Chabrier IMF, the DM estimates agree with \LambdaCDM toy-galaxy models based on contracted DM-halo density profiles. We also find agreement with predictions from hydrodynamical simulations. (4) The central DM content of ETGs does not depend significantly on the environment where galaxies reside, with group and field ETGs having similar DM trends. (Abridged)

Within the standard cosmological scenario the Universe is found to be filled by obscure components (dark matter and dark energy) for ~95% of its energy budget. In particular, almost all the matter content in the Universe is given by dark matter, which dominates the mass budget and drives the dynamics of galaxies and clusters of galaxies. Unfortunately, dark matter and dark energy have not been detected and no direct or indirected observations have allowed to prove their existence and amount. For this reason, some authors have suggested that a modification of Einstein Relativity or the change of the Newton's dynamics law (within a relativistic and classical framework, respectively) could allow to replace these unobserved components. We will start discussing the role of dark matter in the early-type galaxies, mainly in their central regions, investigating how its content changes as a function of the mass and the size of each galaxy and few considerations about the stellar Initial mass function have been made. In the second part of the paper we have described, as examples, some ways to overcome the dark matter hypothesis, by fitting to the observations the modified dynamics coming out from general relativistic extended theories and the MOdyfied Newtonian dynamics (MOND).

We present the first analysis of extended stellar kinematics of elliptical galaxies where a Yukawa--like correction to the Newtonian gravitational potential derived from f(R)-gravity is considered as an alternative to dark matter. In this framework, we model long-slit data and planetary nebulae data out to 7 Re of three galaxies with either decreasing or flat dispersion profiles. We use the corrected Newtonian potential in a dispersion-kurtosis Jeans analysis to account for the mass-anisotropy degeneracy. We find that these modified potentials are able to fit nicely all three elliptical galaxies and the anisotropy distribution is consistent with that estimated if a dark halo is considered. The parameter which measures the "strength" of the Yukawa-like correction is, on average, smaller than the one found previously in spiral galaxies and correlates both with the scale length of the Yukawa-like term and the orbital anisotropy.

We investigate the role of the environment on the colour and stellar population gradients in a local sample of ~3500 central and ~1150 satellite SDSS early-type galaxies (ETGs). The environment is parameterized in terms of the number of satellite galaxies, N_gal in each group. For central galaxies, we find that both optical colour and mass-to-light (M/L) ratio gradients are shallower in central galaxies residing in denser environments (higher N_gal). This trend is driven by metallicity gradients, while age gradients appear to be less dependent on the environment and to have a larger scatter. On the other hand, satellites do not show any differences in terms of the environment. The same results are found if galaxies are classified by central age, and both central and satellite galaxies have shallower gradients if they are older and steeper gradients if younger, satellites being independent of ages. In central galaxies, we show that the observed trends can be explained with the occurrence of dry mergings, which are more numerous in denser environments and producing shallower colour gradients because of more uniform metallicity distributions due to the mixing of stellar populations, while no final clues about merging occurrence can be obtained for satellites. Finally we discuss all systematics on stellar population fitting and their impact on the final results.

(Abridged) Fossil systems are galaxy associations dominated by a relatively isolated, bright elliptical galaxy, surrounded by a group of smaller galaxies lacking L* objects. We analyzed the near-infrared photometric and structural properties of a sample of 20 BGGs present in FGs in order to better understand their formation mechanisms. Their surface-brightness distribution was fitted to a Sersic profile using the GASP2D algorithm. Then, the standard scaling relations were derived for the first time for these galaxies and compared with those of normal ellipticals and brightest cluster galaxies in non-fossil systems. The BGGs presented in this study represent a subset of the most massive galaxies in the Universe. We found that their ellipticity profiles are continuously increasing with the galactocentric radius. Our fossil BCGs follow closely the fundamental plane described by normal ellipticals. However, they depart from both the log \sigma_0 vs. log L_{K_{s}} and log r_{\rm e} vs. log L_{K_{s}} relations described by intermediate mass ellipticals. This occurs in the sense that our BGGs have larger effective radii and smaller velocity dispersions than those predicted by these relations. We also found that more elliptical galaxies systematically deviate from the previous relations while more rounder object do not. No similar correlation was found with the Sersic index. The derived scaling relations can be interpreted in terms of the formation scenario of the BGGs. Because our BGGs follow the fundamental plane tilt but they have larger effective radii than expected for intermediate mass ellipticals, we suggest that they only went through dissipational mergers in a early stage of their evolution and then assembled the bulk of their mass through subsequent dry mergers, contrary to previous claims that BGGs in FGs were formed mainly by the merging of gas-rich galaxies.

We analyze the stellar mass-to-light ratio (M/L) gradients in a large sample of local galaxies taken from the Sloan Digital Sky Survey, spanning a wide range of stellar masses and morphological types. As suggested by the well known relationship between M/L ratios and colors, we show that M/L gradients are strongly correlated with colour gradients, which we trace to the effects of age variations. Stellar M/L gradients generally follow patterns of variation with stellar mass and galaxy type that were previous found for colour and metallicty gradients. In late-type galaxies M/L gradients are negative, steepening with increasing mass. In early-type galaxies M/L gradients are shallower while presenting a two-fold trend: they decrease with mass up to a characteristic mass of \M* \sim 10^10.3 M_sun and increase at larger masses. We compare our findings with other analyses and discuss some implications for galaxy formation and for dark matter estimates.

Scaling relations among dark matter (DM) and stellar quantities are a valuable tool to constrain formation scenarios and the evolution of galactic structures. However, most of the DM properties are actually not directly measured, but derived through model dependent mass mapping procedures. It is therefore crucial to adopt theoretically and observationally well founded models. We use here an updated version of the secondary infall model (SIM) to predict the halo density profile, taking into account the effects of angular momentum, dissipative friction and baryons collapse. The resulting family of halo profiles depends on one parameter only, the virial mass, and nicely fits the projected mass and aperture velocity dispersion of a sample of intermediate redshift lens galaxies. We derive DM related quantities (namely the column density and the Newtonian acceleration) and investigate their correlations with stellar mass, luminosity, effective radius and virial mass.

To investigate the origins of S0 galaxies, we present a new method of analyzing their stellar kinematics from discrete tracers such as planetary nebulae. This method involves binning the data in the radial direction so as to extract the most general possible non-parametric kinematic profiles, and using a maximum likelihood fit within each bin in order to make full use of the information in the discrete kinematic tracers. Both disk and spheroid kinematic components are fitted, with a two-dimensional decomposition of imaging data used to attribute to each tracer a probability of membership in the separate components. Likelihood clipping also allows us to identify objects whose properties are not consistent with the adopted model, rendering the technique robust against contaminants and able to identify additional kinematic features. The method is first tested on an N-body simulated galaxy to assess possible sources of systematic error associated with the structural and kinematic decomposition, which are found to be small. It is then applied to the S0 system NGC~1023, for which a planetary nebula catalogue has already been released and analyzed by (Noordermeer et al., 2008). The correct inclusion of the spheroidal component allows us to show that, contrary to previous claims, the stellar kinematics of this galaxy are indistinguishable from those of a normal spiral galaxy, indicating that it may have evolved directly from such a system via gas stripping or secular evolution. The method also successfully identifies a population of outliers whose kinematics are different from those of the main galaxy; these objects can be identified with a stellar stream associated with the companion galaxy NGC~1023A.

The most accepted scenario for the origin of fossil groups (FGs) is that they are galaxy associations in which the merging rate was fast and efficient. These systems have assembled half of their mass at early epoch of the Universe, subsequently growing by minor mergers. They could contain a fossil record of the galaxy structure formation. We have started a project in order to characterize a large sample of FGs. In this paper we present the analysis of the fossil system RX J105453.3+552102. Optical deep images were used for studying the properties of the brightest group galaxy and for computing the photometric luminosity function of the group. We have also performed a detail dynamical analysis of the system based on redshift data for 116 galaxies. This galaxy system is located at z=0.47, and shows a quite large line-of-sight velocity dispersion \sigma_{v}~1000 km/s. Assuming the dynamical equilibrium, we estimated a virial mass of M ~ 10^{15} h_{70} M_{\odot}. No evidence of substructure was found within 1.4 Mpc radius. We found a statistically significant departure from Gaussianity of the group members velocities in the most external regions of the group. This could indicate the presence of galaxies in radial orbits in the external region of the group. We also found that the photometrical luminosity function is bimodal, showing a lack of M_{r} ~ -19.5 galaxies. The brightest group galaxy shows low Sersic parameter (n~2) and a small peculiar velocity. Indeed, our accurate photometry shows that the difference between the brightest and the second brightest galaxies is 1.9 mag in the r-band, while the classical definition of FGs is based on a magnitude gap of 2. We conclude that this fossil system does not follow the empirical definition of FGs. Nevertheless, it is a massive, old and undisturbed galaxy system with little infall of L^{*} galaxies since its initial collapse.

As part of our current programme to test LCDM predictions for dark matter (DM) haloes using extended kinematical observations of early-type galaxies, we present a dynamical analysis of the bright elliptical galaxy NGC 4374 (M84) based on ~450 Planetary Nebulae (PNe) velocities from the PN.Spectrograph, along with extended long-slit stellar kinematics. This is the first such analysis of a galaxy from our survey with a radially constant velocity dispersion profile. We find that the spatial and kinematical distributions of the PNe agree with the field stars in the region of overlap. The velocity kurtosis is consistent with zero at almost all radii. We construct a series of Jeans models, fitting both velocity dispersion and kurtosis to help break the mass-anisotropy degeneracy. Our mass models include DM halos either with shallow cores or with central cusps as predicted by cosmological simulations - along with the novel introduction in this context of adiabatic halo contraction from baryon infall. Both classes of models confirm a very massive dark halo around NGC 4374, demonstrating that PN kinematics data are well able to detect such haloes when present. Considering the default cosmological mass model, we confirm earlier suggestions that bright galaxies tend to have halo concentrations higher than LCDM predictions, but this is found to be solved if either a Salpeter IMF or adiabatic contraction with a Kroupa IMF is assumed. Thus for the first time a case is found where the PN dynamics may well be consistent with a standard dark matter halo. A cored halo can also fit the data, and prefers a stellar mass consistent with a Salpeter IMF. The less dramatic dark matter content found in lower-luminosity "ordinary" ellipticals suggests a bimodality in the halo properties which may be produced by divergent baryonic effects during their assembly histories.

We discuss the colour, age and metallicity gradients in a wide sample of local SDSS early- and late-type galaxies. From the fitting of stellar population models we find that metallicity is the main driver of colour gradients and the age in the central regions is a dominant parameter which rules the scatter in both metallicity and age gradients. We find a consistency with independent observations and a set of simulations. From the comparison with simulations and theoretical considerations we are able to depict a general picture of a formation scenario.

The age and metallicity gradients for a sample of group and cluster galaxies from N-body+hydrodynamical simulation are analyzed in terms of galaxy stellar mass. Dwarf galaxies show null age gradient with a tail of high and positive values for systems in groups and cluster outskirts. Massive systems have generally zero age gradients which turn to positive for the most massive ones. Metallicity gradients are distributed around zero in dwarf galaxies and become more negative with mass; massive galaxies have steeper negative metallicity gradients, but the trend flatten with mass. In particular, fossil groups are characterized by a tighter distribution of both age and metallicity gradients. We find a good agreement with both local observations and independent simulations. The results are also discussed in terms of the central age and metallicity, as well as the total colour, specific star formation and velocity dispersion.

We analyze the correlations between central dark matter (DM) content of early-type galaxies and their sizes and ages, using a sample of intermediate-redshift (z ~ 0.2) gravitational lenses from the SLACS survey, and by comparing them to a larger sample of z ~ 0 galaxies. We decompose the deprojected galaxy masses into DM and stellar components using combinations of strong lensing, stellar dynamics, and stellar populations modeling. For a given stellar mass, we find that for galaxies with larger sizes, the DM fraction increases and the mean DM density decreases, consistently with the cuspy halos expected in cosmological formation scenarios. The DM fraction also decreases with stellar age, which can be partially explained by the inverse correlation between size and age. The residual trend may point to systematic dependencies on formation epoch of halo contraction or stellar initial mass functions. These results are in agreement with recent findings based on local galaxies by Napolitano, Romanowsky & Tortora (2010) and suggest negligible evidence of galaxy evolution over the last ~ 2.5 Gyr other than passive stellar aging.

We analyze the color gradients (CGs) of ~50000 nearby SDSS galaxies. From synthetic spectral models based on a simplified star formation recipe, we derive the mean spectral properties, and explain the observed radial trends of the color as gradients of the stellar population age and metallicity (Z). The most massive ETGs (M_* > 10^{11} Msun) have shallow CGs in correspondence of shallow (negative) Z gradients. In the stellar mass range 10^(10.3-10.5) < M_* < 10^(11) Msun, the Z gradients reach their minimum of ~ -0.5 dex^{-1}. At M_* ~ 10^{10.3-10.5} Msun, color and Z gradient slopes suddenly change. They turn out to anti-correlate with the mass, becoming highly positive at the very low masses. We have also found that age gradients anti-correlate with Z gradients, as predicted by hierarchical cosmological simulations for ETGs. On the other side, LTGs have gradients which systematically decrease with mass (and are always more negative than in ETGs), consistently with the expectation from gas infall and SN feedback scenarios. Z is found to be the main driver of the trend of color gradients, especially for LTGs, but age gradients are not negligible and seem to play a significant role too. We have been able to highlight that older galaxies have systematically shallower age and Z gradients than younger ones. Our results for high-mass galaxies are in perfect agreement with predictions based on the merging scenario, while the evolution of LTGs and younger and less massive ETGs seems to be mainly driven by infall and SN feedback. (Abridged)

We examine correlations between the masses, sizes, and star formation histories for a large sample of low-redshift early-type galaxies, using a simple suite of dynamical and stellar populations models. We confirm an anti-correlation between size and stellar age, and survey for trends with the central content of dark matter (DM). An average relation between central DM density and galaxy size of <rho_DM> ~ Reff^-2 provides the first clear indication of cuspy DM haloes in these galaxies -- akin to standard LCDM haloes that have undergone adiabatic contraction. The DM density scales with galaxy mass as expected, deviating from suggestions of a universal halo profile for dwarf and late-type galaxies. We introduce a new fundamental constraint on galaxy formation by finding that the central DM fraction decreases with stellar age. This result is only partially explained by the size-age dependencies, and the residual trend is in the opposite direction to basic DM halo expectations. Therefore we suggest that there may be a connection between age and halo contraction, and that galaxies forming earlier had stronger baryonic feedback which expanded their haloes, or else lumpier baryonic accretion that avoided halo contraction. An alternative explanation is a lighter initial mass function for older stellar populations.

Dynamical studies of local ETGs and the Fundamental Plane point to a strong dependence of M/L ratio on luminosity (and stellar mass) with a relation of the form $M/L \propto L^{\gamma}$. The "tilt" $\gamma$ may be caused by various factors, including stellar population properties, IMF, rotational support, luminosity profile non-homology and dark matter (DM) fraction. We evaluate the impact of all these factors using a large uniform dataset of local ETGs from Prugniel & Simien (1997). We take particular care in estimating the stellar masses, using a general star formation history, and comparing different population synthesis models. We find that the stellar M/L contributes little to the tilt. We estimate the total M/L using simple Jeans dynamical models, and find that adopting accurate luminosity profiles is important but does not remove the need for an additional tilt component, which we ascribe to DM. We survey trends of the DM fraction within one effective radius, finding it to be roughly constant for galaxies fainter than $M_B \sim -20.5$, and increasing with luminosity for the brighter galaxies; we detect no significant differences among S0s and fast- and slow-rotating ellipticals. We construct simplified cosmological mass models and find general consistency, where the DM transition point is caused by a change in the relation between luminosity and effective radius. A more refined model with varying galaxy star formation efficiency suggests a transition from total mass profiles (including DM) of faint galaxies distributed similarly to the light, to near-isothermal profiles for the bright galaxies. These conclusions are sensitive to various systematic uncertainties which we investigate in detail, but are consistent with the results of dynamics studies at larger radii.

We present new planetary nebulae (PNe) positions, radial velocities, and magnitudes for 6 early-type galaxies obtained with the Planetary Nebulae Spectrograph, their two-dimensional velocity and velocity dispersion fields. We extend this study to include an additional 10 early-type galaxies with PNe radial velocity measurements available from the literature, to obtain a broader description of the outer-halo kinematics in early-type galaxies. These data extend the information derived from stellar kinematics to typically up to ~8 Re. The combination of photometry, stellar and PNe kinematics shows: i) good agreement between the PNe number density and the stellar surface brightness in the region where the two data sets overlap; ii) good agreement between PNe and stellar kinematics; iii) that the mean rms velocity profiles fall into two groups: with of the galaxies characterized by slowly decreasing profiles and the remainder having steeply falling profiles; iv) a larger variety of velocity dispersion profiles; v) that twists and misalignments in the velocity fields are more frequent at large radii, including some fast rotators; vi) that outer haloes are characterised by more complex radial profiles of the specific angular momentum-related lambda_R parameter than observed within 1Re; vii) that many objects are more rotationally dominated at large radii than in their central parts; and viii) that the halo kinematics are correlated with other galaxy properties, such as total luminosity, isophotal shape, total stellar mass, V/sigma, and alpha parameter, with a clear separation between fast and slow rotators.

We present first results of a study of the halo kinematics for a sample of early type galaxies using planetary nebulae (PNe) as kinematical tracers. PNe allow to extend up to several effective radii (Re) the information from absorption line kinematics (confined to within 1 or 2 Re), providing valuable information and constraints for merger simulations and galaxy formation models. We find that the specific angular momentum per unit mass has a more complex radial dependence when the halo region is taken into account and that the halo velocity dispersion is related to the total galaxy luminosity, isophotal shape, and number of PNe per unit of luminosity

We present new Planetary Nebula Spectrograph observations of the ordinary elliptical galaxy NGC 4494, resulting in positions and velocities of 255 PNe out to 7 effective radii (25 kpc). We also present new wide-field surface photometry from MMT/Megacam, and long-slit stellar kinematics from VLT/FORS2. The spatial and kinematical distributions of the PNe agree with the field stars in the region of overlap. The mean rotation is relatively low, with a possible kinematic axis twist outside 1 Re. The velocity dispersion profile declines with radius, though not very steeply, down to ~70 km/s at the last data point. We have constructed spherical dynamical models of the system, including Jeans analyses with multi-component LCDM-motivated galaxies as well as logarithmic potentials. These models include special attention to orbital anisotropy, which we constrain using fourth-order velocity moments. Given several different sets of modelling methods and assumptions, we find consistent results for the mass profile within the radial range constrained by the data. Some dark matter (DM) is required by the data; our best-fit solution has a radially anisotropic stellar halo, a plausible stellar mass-to-light ratio, and a DM halo with an unexpectedly low central density. We find that this result does not substantially change with a flattened axisymmetric model. Taken together with other results for galaxy halo masses, we find suggestions for a puzzling pattern wherein most intermediate-luminosity galaxies have very low concentration halos, while some high-mass ellipticals have very high concentrations. We discuss some possible implications of these results for DM and galaxy formation.

We analyze a new gravitational lens, OAC-GL J1223-1239, serendipitously found in a deep I-band image of the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS). The lens is a L_*, edge-on S0 galaxy at z=0.4656. The gravitational arc has a radius of 0.42 arcsec. We have determined the total mass and the dark matter (DM) fraction within the Einstein radius as a function of the lensed source redshift, which is presently unknown. For z ~ 1.3, which is in the middle of the redshift range plausible for the source according to some external constraints, we find the central velocity dispersion to be ~180 km/s. With this value, close to that obtained by means of the Faber-Jackson relation at the lens redshift, we compute a 30% DM fraction within the Einstein radius (given the uncertainty in the source redshift, the allowed range for the DM fraction is 25-35 % in our lensing model). When compared with the galaxies in the local Universe, the lensing galaxy, OAC-GL J1223-1239 seems to fall in the transition regime between massive DM dominated galaxies and lower-mass, DM deficient systems.

N-body + hydrodynamical simulations of the formation and evolution of galaxy groups and clusters in a LambdaCDM cosmology are used in order to follow the building-up of the colour-magnitude relation in two clusters and in 12 groups. We have found that galaxies, starting from the more massive, move to the Red Sequence (RS) as they get aged over times and eventually set upon a ``dead sequence'' (DS) once they have stopped their bulk star formation activity. Fainter galaxies keep having significant star formation out to very recent epochs and lie broader around the RS. Environment plays a role as galaxies in groups and cluster outskirts hold star formation activity longer than the central cluster regions. However galaxies experiencing infall from the outskirts to the central parts keep star formation on until they settle on to the DS of the core galaxies. Merging contributes to mass assembly until z~1, after which major events only involve the brightest cluster galaxies. The emerging scenario is that the evolution of the colour-magnitude properties of galaxies within the hierarchical framework is mainly driven by star formation activity during dark matter halos assembly. Galaxies progressively quenching their star formation settle to a very sharp ``red and dead'' sequence, which turns out to be universal, its slope and scatter being almost independent of the redshift (since at least z~1.5) and environment. Differently from the DS, the operatively defined RS evolves more evidently with z, the epoch when it changes its slope being closely corresponding to that at which the passive galaxies population takes over the star forming one: this goes from z~1 in clusters down to 0.4 in normal groups.

Recent results from the Planetary Nebula Spectrograph (PN.S) survey have revealed a rapidly falling velocity dispersion profile in the nearby elliptical galaxy NGC 3379, casting doubts on whether this intermediate-luminosity galaxy has the kind of dark matter halo expected in LambdaCDM cosmology. We present a detailed dynamical study of this galaxy, combining long-slit spectroscopy, SAURON integral-field data, and PN.S velocities, reaching to more than seven effective radii (R_e). We construct spherical and axisymmetric dynamical models for these data with the flexible made-to-measure NMAGIC code, in a sequence of gravitational potentials with varying dark halo mass. We find that the data are consistent both with near-isotropic spherical systems dominated by the stellar mass, and with models in massive halos with strongly radially anisotropic outer parts (beta >~ 0.8 at 7R_e). Formal likelihood limits would exclude (at 1 sigma) the model with stars only, as well as halo models with v_circ(7R_e) >~ 250 km/s. A sequence of more realistic axisymmetric models of different inclinations and a small number of triaxial tests confirm the spherical results. All valid models fitting all the data are dynamically stable over Gyrs, including the most anisotropic ones. NGC 3379 may well have a dark matter halo as predicted by recent merger models within LambdaCDM cosmology, provided its outer envelope is strongly radially anisotropic. (abridged)

We investigate the manner in which lenticular galaxies are formed by studying their stellar kinematics: an S0 formed from a fading spiral galaxy should display similar cold outer disc kinematics to its progenitor, while an S0 formed in a minor merger should be more dominated by random motions. In a pilot study to attempt to distinguish between these scenarios, we have measured the planetary nebula (PN) kinematics of the nearby S0 system NGC 1023. Using the Planetary Nebula Spectrograph, we have detected and measured the line-of-sight velocities of 204 candidate PNe in the field of this galaxy. Out to intermediate radii, the system displays the kinematics of a normal rotationally-supported disc system. After correction of its rotational velocities for asymmetric drift, the galaxy lies just below the spiral galaxy Tully-Fisher relation, as one would expect for a fading system. However, at larger radii the kinematics undergo a gradual but major transition to random motion with little rotation. This transition does not seem to reflect a change in the viewing geometry or the presence of a distinct halo component, since the number counts of PNe follow the same simple exponential decline as the stellar continuum with the same projected disc ellipticity out to large radii. The galaxy's small companion, NGC 1023A, does not seem to be large enough to have caused the observed modification either. This combination of properties would seem to indicate a complex evolutionary history in either the transition to form an S0 or in the past life of the spiral galaxy from which the S0 formed. More data sets of this type from both spirals and S0s are needed in order to definitively determine the relationship between these types of system.

We examine the dark matter properties of nearby early-type galaxies using planetary nebulae (PNe) as mass probes. We have designed a specialised instrument, the Planetary Nebula Spectrograph (PN.S) operating at the William Herschel telescope, with the purpose of measuring PN velocities with best efficiency. The primary scientific objective of this custom-built instrument is the study of the PN kinematics in 12 ordinary round galaxies. Preliminary results showing a dearth of dark matter in ordinary galaxies (Romanowsky et al. 2003) are now confirmed by the first complete PN.S datasets. On the other hand early-type galaxies with a "regular" dark matter content are starting to be observed among the brighter PN.S target sample, thus confirming a correlation between the global dark-to-luminous mass virial ratio (f_DM=M_DM/M_star) and the galaxy luminosity and mass.

The "1st Workshop of Astronomy and Astrophysics for Students" was organized with the purpose of bringing togheter scientific Institutions involved in Astrophysical researches, of the Campania region and beyond. In this context the expression "Astronomy and Astrophysics" was used in its broadest meaning, i.e. including all topics related to the study of the Universe: the Sun and the Solar system, Stars, Galaxies, Cosmology, Astroparticles and Astronomical Technologies, with the belief that all scientific areas would benefit from interaction with related, if not strictly connected, research environments. The goal of the Workshop was to offer Students (undergraduate, graduate, PhD and Post-Docs) the opportunity to share their expertise with colleagues studying and working in other Universities and/or Research Institutes involved in research programs of astrophysical interests. At the same time the Workshop was intended to be a place were the partecipants could, possibly for the first time, promote the work done at any level of their educational and/or professional growth. The full Conference Proceedings are available in electronic form at http://www.na.astro.it/astrows06

We present results from Planetary Nebula Spectrograph (PN.S) observations of the elliptical galaxy NGC 3379 and a description of the data reduction pipeline. We detected 214 planetary nebulae of which 191 are ascribed to NGC 3379, and 23 to the companion galaxy NGC 3384. Comparison with data from the literature show that the PN.S velocities have an internal error of <20km/s and a possible offset of similar magnitude. We present the results of kinematic modelling and show that the PN kinematics are consistent with absorption-line data in the region where they overlap. The resulting combined kinematic data set, running from the center of NGC 3379 out to more than seven effective radii (Reff), reveals a mean rotation velocity that is small compared to the random velocities, and a dispersion profile that declines rapidly with radius. From a series of Jeans dynamical models we find the B-band mass-to-light ratio inside 5 Reff to be 8 to 12 in solar units, and the dark matter fraction inside this radius to be less than 40%. We compare these and other results of dynamical analysis with those of dark-matter-dominated merger simulations, finding that significant discrepancies remain, reiterating the question of whether NGC 3379 has the kind of dark matter halo that the current LambdaCDM paradigm requires.

N--body + hydrodynamical simulations of formation and evolution of galaxy groups in a LambdaCDM cosmology have been performed. The properties of the galaxy populations in 12 groups are here discussed, with focus on the colour-magnitude relation in both normal and fossil groups.

Compact groups are associations of a few galaxies in which the environment plays an important role in galaxy evolution. The low group velocity dispersion favors tidal interactions and mergers, which may bring stars from galaxies to the diffuse intragroup light. Numerical simulations of galaxy clusters in hierarchical cosmologies show that the amount of the diffuse light increases with the dynamical evolution of the cluster. We search for diffuse light in the galaxy group HCG 44 in order to determine its luminosity and luminosity fraction. Combining with literature data, we aim to constrain the dynamical status of Hickson compact groups. We use Intra Group planetary nebulae (IGPNe) as tracers of diffuse light. These are detected by the so-called on band-off band technique. We found 12 emission line objects in HCG 44, none of them associated with the galaxies of the group. 6/12 emission line objects are consistent with being IGPNe in HCG 44, but are also consistent with being Ly$\alpha$ background galaxies. Thus we derive an upper limit to the diffuse light fraction in HCG 44 of 4.7%. We find a correlation between the fraction of elliptical galaxies and the amount of diffuse light in Hickson compact groups. Those with large fraction of diffuse light are those with large fractions in number and luminosity of E/S0 galaxies. We propose an evolutionary sequence for Hickson compact groups in which the amount of diffuse light increases with the dynamical evolution of the group.

We present a catalogue of positions, magnitudes and velocities for 3300 emission-line objects found by the Planetary Nebula Spectrograph in a survey of the Andromeda Galaxy, M31. Of these objects, 2615 are found likely to be planetary nebulae (PNe) associated with M31. Initial results from this survey include: the likely non-existence of Andromeda VIII; a universal PN luminosity function, with the exception of a small amount of obscuration, and a small offset in normalization between bulge and disk components; very faint kinematically-selected photometry implying no cut-off in the disk to beyond 4 scalelengths and no halo population in excess of the bulge out to 10 effective bulge radii; disk kinematics that show significant dispersion and asymmetric drift out to large radii, consistent with a warm flaring disk; and no sign of any variation in kinematics with PN luminosity, suggesting that PNe arise from a fairly uniform population of old stars.

We use the mass-to-light gradients in early-type galaxies to infer the global dark matter fraction, f_d=M_d/M_*, for these systems. We discuss implications about the total star formation efficiency in dark-matter halos and show that the trend of $f_{\rm d}$ with mass produces virial mass-to-light ratios which are consistent with semi-analitical models. Preliminary kurtosis analysis of the quasi-constant M/L galaxies in Romanowsky et al. seems at odd with Dekel et al. simulations.

Radial velocities of 40 intracluster planetary nebulae (ICPNe) in the Virgo cluster were obtained with the new multi-fiber FLAMES spectrograph on UT2 at VLT. For the first time, the lambda 4959 AA line of the [OIII] doublet is seen in a large fraction (50%) of ICPNe spectra, and a large fraction of the photometric candidates with m(5007) <\= 27.2 is spectroscopically confirmed. ICPNe with the velocity dispersion of the Virgo cluster are found in our CORE field 1 deg from M87. These may have originated from tidal mass loss of smaller galaxies in the M87 subcluster halo. In a field 0.25 deg from M87, we see an extended stellar halo of M87 in approximate dynamical equilibrium, but with few ICPNe. Finally, in a field near M84/M86, the ICPNe velocities are highly correlated with the galaxy velocities, showing that any well-mixed intracluster population is yet to form. Overall, the measured velocity distributions confirm the non-uniform dynamical structure and on-going assembly of the Virgo cluster.(Based on data collected with the FLAMES spectrograph at the UT2 of the VLT at Cerro Paranal, Chile, operated by ESO, during observing run 71.B-0147(A))

We have investigated the properties of the diffuse light in the Virgo cluster core region, based on the detection of intracluster planetary nebulae (PNe) in four fields. We eliminate the bias from misclassified faint continuum objects, using improved Monte Carlo simulations, and the contaminations by high redshift Ly$\alpha$ galaxies, using the Ly$\alpha$ luminosity function in blank fields. Recent spectroscopic observations confirm that our photometric PN samples are well-understood. We find that the diffuse stellar population in the Virgo core region is inhomogeneous on scales of 30'-90': there exist significant field-to-field variations in the number density of PNe and the inferred amount of intracluster light, with some empty fields, some fields dominated by extended Virgo galaxy halos, and some fields dominated by the true intracluster component. There is no clear trend with distance from M87. The mean surface luminosity density, its rms variation, and the mean surface brightness of diffuse light in our 4 fields are $\Sigma_B = 2.7 x 10^{6}$ L$_{B\odot}$ arcmin$^{-2}$, ${rms} = 2.1 \times 10^{6}$ L$_{B\odot}$ arcmin$^{-2}$, and $\bar{\mu}_{B}=29.0$ mag arcsec$^{-2}$ respectively. Our results indicate that the Virgo cluster is a dynamically young environment, and that the intracluster component is associated at least partially with local physical processes like galaxy interactions or harassment. We also argue, based on kinematic evidence, that the so-called 'over-luminous' PNe in the halo of M84 are dynamically associated with this galaxy, and must thus be brighter than and part of a different stellar population from the normal PN population in elliptical galaxies.

Since the near future should see a rapidly expanding set of probes of the halo masses of individual early-type galaxies, we introduce a convenient parameter for characterising the halo masses from both observational and theoretical results: \dML, the logarithmic radial gradient of the mass-to-light ratio. Using halo density profiles from LCDM simulations, we derive predictions for this gradient for various galaxy luminosities and star formation efficiencies $\epsilon_{SF}$. As a pilot study, we assemble the available \dML\ data from kinematics in early-type galaxies - representing the first unbiassed study of halo masses in a wide range of early-type galaxy luminosities - and find a correlation between luminosity and \dML, such that the brightest galaxies appear the most dark-matter dominated. We find that the gradients in most of the brightest galaxies may fit in well with the LCDM predictions, but that there is also a population of fainter galaxies whose gradients are so low as to imply an unreasonably high star formation efficiency $\epsilon_{SF} > 1$. This difficulty is eased if dark haloes are not assumed to have the standard LCDM profiles, but lower central concentrations.

We present the Capodimonte Deep Field (OACDF), a deep field covering an area of 0.5 square degrees in the B, V, R optical bands plus six medium-band filters in the wavelength range 773--913 nm. The field reaches the following limiting AB magnitudes: B~25.3, V~24.8 and R~25.1 and contains about 50000 extended sources in the AB magnitude range 18 < R < 25 . Hence, it is intermediate between deep pencil beam surveys and very wide but shallow surveys. The main scientific goal of the OACDF is the identification and characterization of early-type field galaxies at different look-back times in order to study different scenarios of galaxy formation. Parallel goals include the search for groups and clusters of galaxies and the search for rare and peculiar objects (gravitational lenses, QSOs, halo White Dwarfs). In this paper we describe the OACDF data reduction, the methods adopted for the extraction of the photometric catalogs, the photometric calibration and the quality assessment of the catalogs by means of galaxy number counts, spectroscopic and photometric redshifts and star colors. We also present the first results of the search for galaxy overdensities. The depth of the OACDF and its relatively large spatial coverage with respect to pencil beam surveys make it a good tool for further studies of galaxy formation and evolution in the redshift range 0--1, as well as for stellar studies.

Using the Planetary Nebula Spectrograph, we have observed and measured the velocities for some 2764 PNe in the disk and halo of the Andromeda galaxy. Preliminary analysis using a basic ring model shows a rotation curve in good agreement with that obtained from HI data out to ~20kpc. Some substructure has also been detected within the velocity field, which can be modeled as the continuation of the tidal-remnant known as the Southern Stream, as it passes through Andromeda's disk.

The Planetary Nebula Spectrograph is a dedicated instrument for measuring radial velocity of individual Planetary Nebulae (PNe) in galaxies. This new instrument is providing crucial data with which to probe the structure of dark halos in the outskirts of elliptical galaxies in particular, which are traditionally lacking of easy interpretable kinematical tracers at large distance from the center. Preliminary results on a sample of intermediate luminosity galaxies have shown little dark matter within 5 ~ R_eff implying halos either not as massive or not as centrally concentrated as CDM predicts (Romanowsky et al. 2003). We briefly discuss whether this is consistent with a systematic trend of the dark matter content with the luminosity as observed in an extended sample of early-type galaxies.

In several nearby intermediate luminosity early-type galaxies, recent observations at large radii have shown a indications of a lack of dark matter, substantially at odds with the prediction from the CDM hierarchical merger models. Here we discuss a pilot observational project for the study of the internal kinematical and dynamical properties of this remarkable sample of galaxies. Using the VIMOS-IFU in its high spectral resolution mode, it would be possible to investigate the regions up to about 1.2 effective radii, taking advantage of the much larger field of view and telescope diameter. This will allow to disclose the presence of any kinematical substructures which could affect the conclusion on the mass modeling and definitely clarify the inner structure of this particular class of early-type galaxies.

We present a possible orbit for the Southern Stream of stars in M31, which connects it to the Northern Spur. Support for this model comes from the dynamics of planetary nebulae (PNe) in the disk of M31: analysis of a new sample of 2611 PNe obtained using the Planetary Nebula Spectrograph reveals ~20 objects whose kinematics are inconsistent with the normal components of the galaxy, but which lie at the right positions and velocities to connect the two photometric features via this orbit. The satellite galaxy M32 is coincident with the stream both in position and velocity, adding weight to the hypothesis that the stream comprises its tidal debris.

A new instrument is providing crucial data with which to probe the structure of dark halos in elliptical galaxies.

Planetary nebulae (PNe) may be the most promising tracers in the halos of early-type galaxies. We have used multi-object spectrographs on the WHT and the VLT, and the new Planetary Nebula Spectrograph on the WHT, to obtain hundreds of PN velocities in a small sample of nearby galaxies. These ellipticals show weak halo rotation, which may be consistent with ab initio models of galaxy formation, but not with more detailed major merger simulations. The galaxies near L* show evidence of a universal declining velocity dispersion profile, and dynamical models indicate the presence of little dark matter within 5 Reff -- implying halos either not as massive or not as centrally concentrated as CDM predicts.

We have found a correlation betwen the M/L global gradients and the structural parameters of the luminous components of a sample of 19 early-type galaxies. Such a correlation supports the hypothesis that there is a connection between the dark matter content and the evolution of the baryonic component in such systems.

The kinematics of the outer parts of three intermediate-luminosity elliptical galaxies have been studied using the Planetary Nebula Spectrograph. The galaxies' velocity dispersion profiles are found to decline with radius; dynamical modeling of the data indicates the presence of little if any dark matter in these galaxies' halos. This surprising result conflicts with findings in other galaxy types, and poses a challenge to current galaxy formation theories.

We use a high resolution collisionless simulation of a Virgo--like cluster in a $\Lambda$CDM cosmology to determine the velocity and clustering properties of the diffuse stellar component in the intracluster region at the present epoch. The simulated cluster builds up hierarchically and tidal interactions between member galaxies and the cluster potential produce a diffuse stellar component free-flying in the intracluster medium. Here we adopt an empirical scheme to identify tracers of the stellar component in the simulation and hence study its properties. We find that at $z=0$ the intracluster stellar light is mostly unrelaxed in velocity space and clustered in structures whose typical clustering radii are about 50 kpc at R=400--500 kpc from the cluster center, and predict the radial velocity distribution expected in spectroscopic follow-up surveys. Finally, we compare the spatial clustering in the simulation with the properties of the Virgo intracluster stellar population, as traced by ongoing intracluster planetary nebulae surveys in Virgo. The preliminary results indicate a substantial agreement with the observed clustering properties of the diffuse stellar population in Virgo.

We have observed emission line objects located in a 0.26 deg^2 field in the M96 (Leo) group, coincident with the intergalactic HI cloud. The emission line objects were selected using the same procedure as used for the search for intracluster planetary nebulae in the Virgo cluster, and their m_5007 luminosity function has a bright cut-off ~1.2 magnitude fainter than for the PNLF associated with the elliptical galaxies in the M96 group.Therefore the vast majority of these emission line objects are compatible with not being intragroup planetary nebulae at the Leo group distance of 10 Mpc. Spectroscopic follow-up of two emission line objects in this Leo field showed that indeed these do not have the [OIII] doublet expected for a real PN. The brighter source is identified as a starburst object at redshift z = 3.128, because of a second emission in the near infrared, identified as FeII (\lambda 2220 AA). From these data we derive three main results: (i) from the absence of PN we can determine a more stringent upper limit to the surface brightness in any old stellar population associated with the Leo HI cloud. (ii) This translates to an upper limit of 1.6% for the fraction of luminosity in a diffuse intragroup component in the densest 3 x 2 deg area of the Leo group,relative to the light in galaxies. (iii) Using this Leo field as a blank field, we derive an average fraction of 13.6% background emission line objects that enter in surveys of Virgo intracluster PN. The small fraction confirms the validity of the selection criteria for intracluster PN candidates in Virgo.

We have identified 38 candidates of intracluster planetary nebulae (ICPNe) in a $34'\times27'$ field in the core of the Virgo cluster based on the Suprime-Cam imaging through two narrow-band filters centered at the redshifted wavelengths of the [OIII]$\lambda=5007\AA$ and the Halpha \lambda=6563\AA lines. Broad-band images in V and R bands are used to check for any emissions in the adjacent continuum. We describe the method briefly and present the list of intracluster planetary nebulae candidates, together with their finding charts. The ICPN candidates show a highly inhomogeneous distribution, which may suggest an association with the M86-M84 subcluster. Fraction of diffuse intracluster light with respect to total light in galaxies is estimated to be about 10%, leading to an estimate of about 20% for the baryon fraction. Spectroscopic follow up and a wider survey are critical to reveal the nature of intracluster stellar population.

We have identified intracluster planetary nebulae in a Virgo cluster core field by imaging with the Subaru Suprime-Cam through two narrow band filters centered at the redshifted wavelengths of the [OIII] $\lambda=5007$ \AA and the H$\alpha$ $\lambda = 6563$ \AA lines; broad-band images in V and R were acquired to check for emission in the adjacent continuum. Emission line objects in Virgo are then selected from the two-color diagram ([OIII] -- Halpha) vs. ([OIII] -- (V+R)), which is calibrated using PNe in M84 (Jacoby et al. 1990). Using both [OIII] and Halpha allows us to distinguish bona-fide planetary nebulae from high redshift emission-line galaxies at the bright end of the [OIII] luminosity function. Spectroscopic observations of a subsample of these objects were made at the TNG and at the VLT, in a region around M84 and in an intracluster field respectively. The observations confirm the efficiency of the combined [OIII]+H$\alpha$ imaging to identify true PNe. We also obtained the first spectrum of an intracluster PN which shows the [OIII] doublet with S/N > 10 and its Halpha emission. From the results based on the spectroscopic follow-up, we derive a lower limit to the fraction of the Virgo cluster light contributed by the intracluster stars at the surveyed position in the cluster core: it amounts to 10%.

The distribution of the radial trends of the mass-to-light ratios (M/L) within an assorted sample of early-type galaxies is discussed. Three classes of galaxies are identified according to their M/L gradients. Two such classes are characterized by the presence or by the absence of a radial gradient of the dark-matter (DM) distribution. A third class contains objects which are likely undergoing interaction; they exhibit steep M/L gradients which are possibly the result of a wrong assumption on their equilibrium conditions. Finally, a possible correlation between DM content and morphological types is briefly discussed.

Evidence for a disturbed velocity structure in the outer regions of the galaxy NGC1399 comes from a re-analysis of the planetary nebulae data from Arnaboldi et al. (1994). We find a strong rotation along a P.A.=140 degrees for R<=140'' which is followed by a rapid drop off at larger radii, where the velocity dispersion starts to increase. This kinematical behavior can be interpreted as an indication for an interaction scenario. Interaction is advocated in previous analysis of the halo regions of this system, based on different dynamical tracers such as globular clusters and X-rays, but in all these studies the mass distribution is derived under the equilibrium hypothesis, which is not appropriate when an interaction takes place. Here we attempt a non-equilibrium dynamical analysis of NGC1399: with a simple model and under the impulse approximation, we show that the observed kinematics is consistent with an energy injection caused by a flyby encounter of NGC 1399 with the nearby system NGC 1404. In this approach, we find a mass-to-light ratio, M/L_B=26 within R=400'', which is about half of that requested when equilibrium is assumed, i.e. M/L_B=56.

We have imaged an empty area of 34'x34' one and a half degree north of the Virgo cluster core to survey for intracluster planetary nebula candidates. We have implemented and tested a fully automatic procedure for the selection of emission line objects in wide-field images, based on the on-off technique from Ciardullo and Jacoby. Freeman et al. have spectroscopically confirmed a sample of intracluster planetary nebulae in one Virgo field. We use the photometric and morphological properties of this sample to test our selection procedure. In our newly surveyed Virgo field, 75 objects were identified as best candidates for intracluster PNe. The luminosity function of the spectroscopically confirmed PNe shows a brighter cut-off than the planetary nebula luminosity function for the inner regions of M87. Such a brighter cut-off is also observed in the newly surveyed field and indicates a smaller distance modulus, implying that the front end of the Virgo cluster is closer to us by a significant amount: 14% closer (2.1 Mpc) than M87 for the spectroscopic field, using the PN luminosity function distance of 14.9 Mpc to M87, and 19% closer (2.8 Mpc) than M87 for the newly surveyed field. Independent distance indicators (Tully-Fisher relation for Virgo spirals and surface brightness fluctuations for Virgo ellipticals) agree with these findings. From these two Virgo cluster fields there is no evidence that the surface luminosity density for the diffuse stellar component in the cluster decreases with radius. The luminosity surface density of the diffuse stellar population is comparable to that of the galaxies.

Multi-object and multi-fiber spectrographs on 4 and 8 meter telescopes make it possible to use extragalactic planetary nebulae (PNe) in the outer halos of early-type galaxies as kinematical tracers, where classical techniques based on integrated stellar light fail. Until now, published PNe radial velocity samples are small, with only a few tens of radial velocity measurements (except for a few cases like NGC 5128 or M31), causing uncertainties in the mass and angular momentum estimates based on these data. To quantify these uncertainties, we have made equilibrium models for spherical galaxies, with and without dark matter, and via Montecarlo simulations we produce radial velocity samples with different sizes. We then apply, to these discrete radial velocity fields, the same standard kinematical analysis as it is commonly done with small samples of observed PNe radial velocities. By comparison of the inferred quantities with those computed from the analytical model, we test for systematic biases and establish a robust procedure to infer the angular momentum distribution and radial velocity dispersion profiles from such samples.

Planetary Nebulae (PNe) are very important kinematical tracers of the outer regions of early-type galaxies, where the integrated light techniques fail. Under ad hoc assumptions, they allow measurements of rotation velocity and velocity dispersion profile from descrete radial velocity fields. We present the results on the precision allowed by different set of radial velocity samples, discuss the hypotheses in the analysis of descrete velocity fields and their impact on the inferred kinematics of the stellar population.