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.

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

Every field of Science is undergoing unprecedented changes in the discovery process, and Astronomy has been a main player in this transition since the beginning. The ongoing and future large and complex multi-messenger sky surveys impose a wide exploiting of robust and efficient automated methods to classify the observed structures and to detect and characterize peculiar and unexpected sources. We performed a preliminary experiment on KiDS DR4 data, by applying to the problem of anomaly detection two different unsupervised machine learning algorithms, considered as potentially promising methods to detect peculiar sources, a Disentangled Convolutional Autoencoder and an Unsupervised Random Forest. The former method, working directly on images, is considered potentially able to identify peculiar objects like interacting galaxies and gravitational lenses. The latter instead, working on catalogue data, could identify objects with unusual values of magnitudes and colours, which in turn could indicate the presence of singularities.

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

We use nine different galaxy formation scenarios in ten cosmological simulation boxes from the EAGLE suite of {\Lambda}CDM hydrodynamical simulations to assess the impact of feedback mechanisms in galaxy formation and compare these to observed strong gravitational lenses. To compare observations with simulations, we create strong lenses with $M_\star$ > $10^{11}$ $M_\odot$ with the appropriate resolution and noise level, and model them with an elliptical power-law mass model to constrain their total mass density slope. We also obtain the mass-size relation of the simulated lens-galaxy sample. We find significant variation in the total mass density slope at the Einstein radius and in the projected stellar mass-size relation, mainly due to different implementations of stellar and AGN feedback. We find that for lens selected galaxies, models with either too weak or too strong stellar and/or AGN feedback fail to explain the distribution of observed mass-density slopes, with the counter-intuitive trend that increasing the feedback steepens the mass density slope around the Einstein radius ($\approx$ 3-10 kpc). Models in which stellar feedback becomes inefficient at high gas densities, or weaker AGN feedback with a higher duty cycle, produce strong lenses with total mass density slopes close to isothermal (i.e. -d log({\rho})/d log(r) $\approx$ 2.0) and slope distributions statistically agreeing with observed strong lens galaxies in SLACS and BELLS. Agreement is only slightly worse with the more heterogeneous SL2S lens galaxy sample. Observations of strong-lens selected galaxies thus appear to favor models with relatively weak feedback in massive galaxies.

Metallicity and gas content are intimately related in the baryonic exchange cycle of galaxies, and galaxy evolution scenarios can be constrained by quantifying this relation. To this end, we have compiled a sample of ~400 galaxies in the Local Universe, dubbed "MAGMA" (Metallicity And Gas for Mass Assembly), which covers an unprecedented range in parameter space, spanning more than 5 orders of magnitude in stellar mass (Mstar), star-formation rate (SFR), and gas mass (Mgas), and a factor of ~60 in metallicity [Z, 12+log(O/H)]. Stellar masses and SFRs have been recalculated for all the galaxies using IRAC, WISE and GALEX photometry, and 12+log(O/H) has been transformed, where necessary, to a common metallicity calibration. To assess the true dimensionality of the data, we have applied multi-dimensional principal component analyses (PCAs) to our sample. In confirmation of previous work, we find that even with the vast parameter space covered by MAGMA, the relations between Mstar, SFR, Z and Mgas (MHI+MH2) require only two dimensions to describe the hypersurface. To accommodate the curvature in the Mstar-Z relation, we have applied a piecewise 3D PCA that successfully predicts observed 12+log(O/H) to an accuracy of ~0.1dex. MAGMA is a representative sample of isolated star-forming galaxies in the Local Universe, and can be used as a benchmark for cosmological simulations and to calibrate evolutionary trends with redshift.

The Euclid mission will observe well over a billion galaxies out to $z\sim6$ and beyond. This will offer an unrivalled opportunity to investigate several key questions for understanding galaxy formation and evolution. The first step for many of these studies will be the selection of a sample of quiescent and star-forming galaxies, as is often done in the literature by using well known colour techniques such as the `UVJ' diagram. However, given the limited number of filters available for the Euclid telescope, the recovery of such rest-frame colours will be challenging. We therefore investigate the use of observed Euclid colours, on their own and together with ground-based u-band observations, for selecting quiescent and star-forming galaxies. The most efficient colour combination, among the ones tested in this work, consists of the (u-VIS) and (VIS-J) colours. We find that this combination allows users to select a sample of quiescent galaxies complete to above $\sim70\%$ and with less than 15$\%$ contamination at redshifts in the range $0.75<z<1$. For galaxies at high-z or without the u-band complementary observations, the (VIS-Y) and (J-H) colours represent a valid alternative, with $>65\%$ completeness level and contamination below 20$\%$ at $1<z<2$ for finding quiescent galaxies. In comparison, the sample of quiescent galaxies selected with the traditional UVJ technique is only $\sim20\%$ complete at $z<3$, when recovering the rest-frame colours using mock Euclid observations. This shows that our new methodology is the most suitable one when only Euclid bands, along with u-band imaging, are available.

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 upcoming new generation of optical spectrographs on four-meter-class telescopes will provide invaluable information for reconstructing the history of star formation in individual galaxies up to redshifts of about 0.7. We aim at defining simple but robust and meaningful physical parameters that can be used to trace the coexistence of widely diverse stellar components: younger stellar populations superimposed on the bulk of older ones. We produce spectra of galaxies closely mimicking data from the forthcoming Stellar Populations at intermediate redshifts Survey (StePS), a survey that uses the WEAVE spectrograph on the William Herschel Telescope. First, we assess our ability to reliably measure both ultraviolet and optical spectral indices in galaxies of different spectral types for typically expected signal-to-noise levels. Then, we analyze such mock spectra with a Bayesian approach, deriving the probability density function of r- and u-band light-weighted ages as well as of their difference. We find that the ultraviolet indices significantly narrow the uncertainties in estimating the r- and u-band light-weighted ages and their difference in individual galaxies. These diagnostics, robustly retrievable for large galaxy samples even when observed at moderate signal-to-noise ratios, allow us to identify secondary episodes of star formation up to an age of ~0.1 Gyr for stellar populations older than ~1.5 Gyr, pushing up to an age of ~1 Gyr for stellar populations older than ~5 Gyr. The difference between r-band and u-band light-weighted ages is shown to be a powerful diagnostic to characterize and constrain extended star-formation histories and the presence of young stellar populations on top of older ones. This parameter can be used to explore the interplay between different galaxy star-formation histories and physical parameters such as galaxy mass, size, morphology, and environment.

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.

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

Large scale imaging surveys will increase the number of galaxy-scale strong lensing candidates by maybe three orders of magnitudes beyond the number known today. Finding these rare objects will require picking them out of at least tens of millions of images and deriving scientific results from them will require quantifying the efficiency and bias of any search method. To achieve these objectives automated methods must be developed. Because gravitational lenses are rare objects reducing false positives will be particularly important. We present a description and results of an open gravitational lens finding challenge. Participants were asked to classify 100,000 candidate objects as to whether they were gravitational lenses or not with the goal of developing better automated methods for finding lenses in large data sets. A variety of methods were used including visual inspection, arc and ring finders, support vector machines (SVM) and convolutional neural networks (CNN). We find that many of the methods will be easily fast enough to analyse the anticipated data flow. In test data, several methods are able to identify upwards of half the lenses after applying some thresholds on the lens characteristics such as lensed image brightness, size or contrast with the lens galaxy without making a single false-positive identification. This is significantly better than direct inspection by humans was able to do. (abridged)

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.

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.

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.

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.

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

In this paper we introduce the SEAGLE (i.e. Simulating EAGLE LEnses) program, that approaches the study of galaxy formation through strong gravitational lensing, using a suite of high-resolution hydrodynamic simulations, Evolution and Assembly of GaLaxies and their Environments (EAGLE) project. We introduce the simulation and analysis pipeline and present the first set of results from our analysis of early-type galaxies. We identify and extract an ensemble of simulated lens galaxies and use the GLAMER ray-tracing lensing code to create mock lenses similar to those observed in the SLACS and SL2S surveys, using a range of source parameters and galaxy orientations, including observational effects such as the Point-Spread-Function (PSF), pixelization and noise levels, representative of single-orbit observations with the Hubble Space Telescope (HST) using the ACS-F814W filter. We subsequently model these mock lenses using the code LENSED, treating them in the same way as observed lenses. We also estimate the mass model parameters directly from the projected surface mass density of the simulated galaxy, using an identical mass model family. We perform a three-way comparison of all the measured quantities with real lenses. We find the average total density slope of EAGLE lenses, $t=2.26\; (0.25\; \rm{rms})$ to be higher than SL2S, $t=2.16$ or SLACS, $t=2.08$. We find a very strong correlation between the external shear ($\gamma$) and the complex ellipticity ($\epsilon$), with $\gamma \sim \epsilon/4$. This correlation indicates a degeneracy in the lens mass modeling. We also see a dispersion between lens modeling and direct fitting results, indicating systematical biases.

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.

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.

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)

Photometric redshifts (photo-z's) provide an alternative way to estimate the distances of large samples of galaxies and are therefore crucial to a large variety of cosmological problems. Among the various methods proposed over the years, supervised machine learning (ML) methods capable to interpolate the knowledge gained by means of spectroscopical data have proven to be very effective. METAPHOR (Machine-learning Estimation Tool for Accurate PHOtometric Redshifts) is a novel method designed to provide a reliable PDF (Probability density Function) of the error distribution of photometric redshifts predicted by ML methods. The method is implemented as a modular workflow, whose internal engine for photo-z estimation makes use of the MLPQNA neural network (Multi Layer Perceptron with Quasi Newton learning rule), with the possibility to easily replace the specific machine learning model chosen to predict photo-z's. After a short description of the software, we present a summary of results on public galaxy data (Sloan Digital Sky Survey - Data Release 9) and a comparison with a completely different method based on Spectral Energy Distribution (SED) template fitting.

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 present METAPHOR (Machine-learning Estimation Tool for Accurate PHOtometric Redshifts), a method able to provide a reliable PDF for photometric galaxy redshifts estimated through empirical techniques. METAPHOR is a modular workflow, mainly based on the MLPQNA neural network as internal engine to derive photometric galaxy redshifts, but giving the possibility to easily replace MLPQNA with any other method to predict photo-z's and their PDF. We present here the results about a validation test of the workflow on the galaxies from SDSS-DR9, showing also the universality of the method by replacing MLPQNA with KNN and Random Forest models. The validation test include also a comparison with the PDF's derived from a traditional SED template fitting method (Le Phare).

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.

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.

A variety of fundamental astrophysical science topics require the determination of very accurate photometric redshifts (photo-z's). A wide plethora of methods have been developed, based either on template models fitting or on empirical explorations of the photometric parameter space. Machine learning based techniques are not explicitly dependent on the physical priors and able to produce accurate photo-z estimations within the photometric ranges derived from the spectroscopic training set. These estimates, however, are not easy to characterize in terms of a photo-z Probability Density Function (PDF), due to the fact that the analytical relation mapping the photometric parameters onto the redshift space is virtually unknown. We present METAPHOR (Machine-learning Estimation Tool for Accurate PHOtometric Redshifts), a method designed to provide a reliable PDF of the error distribution for empirical techniques. The method is implemented as a modular workflow, whose internal engine for photo-z estimation makes use of the MLPQNA neural network (Multi Layer Perceptron with Quasi Newton learning rule), with the possibility to easily replace the specific machine learning model chosen to predict photo-z's. We present a summary of results on SDSS-DR9 galaxy data, used also to perform a direct comparison with PDF's obtained by the Le Phare SED template fitting. We show that METAPHOR is capable to estimate the precision and reliability of photometric redshifts obtained with three different self-adaptive techniques, i.e. MLPQNA, Random Forest and the standard K-Nearest Neighbors models.

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

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.

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

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.

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.

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.

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)

We discuss the thermodynamic and dynamical properties of a variable dark energy model with density scaling as $\rho_x \propto (1+z)^{m}$, z being the redshift. These models lead to the creation/disruption of matter and radiation, which affect the cosmic evolution of both matter and radiation components in the Universe. In particular, we have studied the temperature-redshift relation of radiation, which has been constrained using a recent collection of cosmic microwave background (CMB) temperature measurements up to $z \sim 3$. We find that, within the uncertainties, the model is indistinguishable from a cosmological constant which does not exchange any particles with other components. Future observations, in particular measurements of CMB temperature at large redshift, will allow to give firmer bounds on the effective equation of state parameter $w_{eff}$ for such types of dark energy models.

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.

For a general class of analytic f(R)-gravity theories, we discuss the weak field limit in view of gravitational lensing. Though an additional Yukawa term in the gravitational potential modifies dynamics with respect to the standard Newtonian limit of General Relativity, the motion of massless particles results unaffected thanks to suitable cancellations in the post-Newtonian limit. Thus, all the lensing observables are equal to the ones known from General Relativity. Since f(R)-gravity is claimed, among other things, to be a possible solution to overcome for the need of dark matter in virialized systems, we discuss the impact of our results on the dynamical and gravitational lensing analyses. In this framework, dynamics could, in principle, be able to reproduce the astrophysical observations without recurring to dark matter, but in the case of gravitational lensing we find that dark matter is an unavoidable ingredient. Another important implication is that gravitational lensing, in the post-Newtonian limit, is not able to constrain these extended theories, since their predictions do not differ from General Relativity.

The thermodynamic and dynamical properties of a variable dark energy model with density scaling as rho_x \propto (1+z)^m, z being the redshift, are discussed following the outline of Jetzer et al. This kind of models are proven to lead to the creation/disruption of matter and radiation, which affect the cosmic evolution of both matter and radiation components in the Universe. In particular, we have concentrated on the temperature-redshift relation of radiation, which has been constrained using a very recent collection of cosmic microwave background (CMB) temperature measurements up to z ~ 3. For the first time, we have combined this observational probe with a set of independent measurements (Supernovae Ia distance moduli, CMB anisotropy, large-scale structure and observational data for the Hubble parameter), which are commonly adopted to constrain dark energy models. We find that, within the uncertainties, the model is indistinguishable from a cosmological constant which does not exchange any particles with other components. Anyway, while temperature measurements and Supernovae Ia tend to predict slightly decaying models, the contrary happens if CMB data are included. Future observations, in particular measurements of CMB temperature at large redshift, will allow to give firmer bounds on the effective equation of state parameter w_eff of this kind of dark energy models.

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.

Modified Newtonian Dynamics (MOND) has been shown to be able to fit spiral galaxy rotation curves as well as giving a theoretical foundation for empirically determined scaling relations, such as the Tully - Fisher law, without the need for a dark matter halo. As a complementary analysis, one should investigate whether MOND can also reproduce the dynamics of early - type galaxies (ETGs) without dark matter. As a first step, we here show that MOND can indeed fit the observed central velocity dispersion $\sigma_0$ of a large sample of ETGs assuming a simple MOND interpolating functions and constant anisotropy. We also show that, under some assumptions on the luminosity dependence of the Sersic n parameter and the stellar M/L ratio, MOND predicts a fundamental plane for ETGs : a log - linear relation among the effective radius $R_{eff}$, $\sigma_0$ and the mean effective intensity $\langle I_e \rangle$. However, we predict a tilt between the observed and the MOND fundamental planes.

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 investigate scaling relations between the dark matter (DM) halo model parameters for a sample of intermediate redshift early - type galaxies (ETGs) resorting to a combined analysis of Einstein radii and aperture velocity dispersions. Modeling the dark halo with a Navarro - Frenk - White profile and assuming a Salpeter initial mass function (IMF) to estimate stellar masses, we find that the column density ${\cal{S}}$ and the Newtonian acceleration within the halo characteristic radius $r_s$ and effective radius $R_{eff}$ are not universal quantities, but correlate with the luminosity $L_V$, the stellar mass $M_{\star}$ and the halo mass $M_{200}$, contrary to recent claims in the literature. We finally discuss a tight correlation among the DM mass $M_{DM}(R_{eff})$ within the effective radius $R_{eff}$, the stellar mass $M_{\star}(R_{eff})$ and $R_{eff}$ itself. The slopes of the scaling relations discussed here strongly depend, however, on the DM halo model and the IMF adopted so that these ingredients have to be better constrained in order to draw definitive conclusions on the DM scaling relations for ETGs.

The nature of the dark energy is still a mystery and several models have been proposed to explain it. Here we consider a phenomenological model for dark energy decay into photons and particles as proposed by Lima (J. Lima, Phys. Rev. D 54, 2571 (1996)). He studied the thermodynamic aspects of decaying dark energy models in particular in the case of a continuous photon creation and/or disruption. Following his approach, we derive a temperature redshift relation for the CMB which depends on the effective equation of state $w_{eff}$ and on the "adiabatic index" $\gamma$. Comparing our relation with the data on the CMB temperature as a function of the redshift obtained from Sunyaev-Zel'dovich observations and at higher redshift from quasar absorption line spectra, we find $w_{eff}=-0.97 \pm 0.034$, adopting for the adiabatic index $\gamma=4/3$, in good agreement with current estimates and still compatible with $w_{eff}=-1$, implying that the dark energy content being constant in time.

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.

The dark matter content of early,- type galaxies (ETGs) is a hotly debated topic with contrasting results arguing in favour or against the presence of significant dark mass within the effective radius and the change with luminosity and mass. In order to address this question, we investigate here the global mass - to - light ratio $\Upsilon(r) = M(r)/L(r)$ of a sample of 21 lenses observed within the Sloan Lens ACS (SLACS) survey. We follow the usual approach of modeling the galaxy as a two component systems, but we use a phenomenological ansatz for $\Upsilon(r)$, proposed by some of us in Tortora et al. (2007), able to smoothly interpolate between constant $M/L$ models and a wide class of dark matter haloes. The resulting galaxy model is then fitted to the data on the Einstein radius and velocity dispersion. Our phenomenological model turns out to be in well agreement with the data suggesting the presence of massive dark matter haloes in order to explain the lensing and dynamics properties of the SLACS lenses. According to the values of the dark matter mass fraction, we argue that the halo may play a significant role in the inner regions probed by the data, but such a conclusion strongly depends on the adopted initial mass function of the stellar population. Finally, we find that the dark matter mass fraction within $R_{eff}$ scales with both the total luminosity and stellar mass in such a way that more luminous (and hence more massive) galaxies have a larger dark matter content.

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 study the feedback from an AGN on stellar formation within its host galaxy, mainly using one high resolution numerical simulation of the jet propagation within the interstellar medium of an early-type galaxy. In particular, we show that in a realistic simulation where the jet propagates into a two-phase ISM, star formation can initially be slightly enhanced and then, on timescales of few million years, rapidly quenched, as a consequence both of the high temperatures attained and of the reduction of cloud mass (mainly due to Kelvin-Helmholtz instabilities). We then introduce a model of (prevalently) {\em negative} AGN feedback, where an exponentially declining star formation is quenched, on a very short time scale, at a time t_AGN, due to AGN feedback. Using the Bruzual & Charlot (2003) population synthesis model and our star formation history, we predict galaxy colours from this model and match them to a sample of nearby early-type galaxies showing signs of recent episodes of star formation (Kaviraj et al. 2007). We find that the quantity t_gal - t_AGN, where t_gal is the galaxy age, is an excellent indicator of the presence of feedback processes, and peaks significantly around t_gal - t_AGN \approx 0.85 Gyr for our sample, consistent with feedback from recent energy injection by AGNs in relatively bright (M_{B} \lsim -19) and massive nearby early-type galaxies. Galaxies that have experienced this recent feedback show an enhancement of 3 magnitudes in NUV(GALEX)-g, with respect to the unperturbed, no-feedback evolution. Hence they can be easily identified in large combined near UV-optical surveys.

Gravitational lensing represents a powerful tool to estimate the cosmological pa- rameters and the distribution of dark matter. I will describe the main observable quantities, concentrating on strong lensing, that manifests its effect through the formation of spectacular events, like multiple quasars, Einstein rings and arcs in clusters of galaxies. In events where a quasar is lensed by an intervening galaxy, it is possible to give an estimate of the Hubble constant H_0, by choosing a mass density model for the lens, thus allowing to constrain the dark matter content.

Recent observations of the line of sight velocity profile of elliptical galaxies have furnished controversial results with some works favouring the presence of a large amount of dark matter in the outer regions and others arguing in favour of no dark matter at all. In order to shed new light on this controversy, we propose here a new phenomenological description of the total mass profile of galaxies. Under the hypothesis of spherical symmetry, we assume a double power-law expression for the global M/L ratio Upsilon(r)= Upsilon_0(r/r_0) ^{alpha}(1+r/r_0)^{beta}. In particular, Upsilon propto r^{alpha} for r/r_0<<1 so that alpha=0 mimics a constant M/L ratio in the inner regions, while, for (r/r_0>>1), Upsilon propto r^{alpha+beta} thus showing that models with alpha+beta=0 have an asymptotically constant M/L ratio. A wide range of possibilities is obtained by varying the slope parameters in the range we determine on the basis of physical considerations. Choosing a general expression for the luminosity density profile j(r), we work out an effective galaxy model that accounts for all the phenomenology observed in real elliptical galaxies. We derive the main dynamics and lensing properties of such an effective model. We analyze a general class of models, able to take into account different dynamical trends. We are able to obtain analytical expressions for the main dynamical and lensing quantities. We show that constraining the values of alpha+beta makes it possible to analyze the problem of the dark matter in elliptical galaxies. Indeed, positive values of alpha+beta would be a strong evidence for dark matter. Finally we indicate possible future approaches in order to face the observational data, in particular using velocity dispersion profiles and lensed quasar events.

We consider scalar tensor theories of gravity assuming that the scalar field is non minimally coupled with gravity. We use this theory to study evolution of a flat homogeneous and isotropic universe. In this case the dynamical equations can be derived form a point like Lagrangian. We study the general properties of dynamics of this system and show that for a wide range of initial conditions such models lead in a natural way to an accelerated phase of expansion of the universe. Assuming that the point like Lagrangian admits a Noether symmetry we are able to explicitly solve the dynamical equations. We study one particular model and show that its predictions are compatible with observational data, namely the publicly available data on type Ia supernovae, the parameters of large scale structure determined by the 2-degree Field Galaxy Redshift Survey (2dFGRS), the measurements of cosmological distances with the Sunyaev-Zel'dovich effect and the rate of growth of density perturbations}{It turns out that this model have a very interesting feature of producing in a natural way an epoch of accelerated expansion. With an appropriate choice of parameters our model is fully compatible with several observed characteristics of the universe

Abandoning the perfect fluid hypothesis, we investigate here the possibility that the dark energy equation of state (EoS) $w$ is a nonlinear function of the energy density $\rho$. To this aim, we consider four different EoS describing classical fluids near thermodynamical critical points and discuss the main features of cosmological models made out of dust matter and a dark energy term with the given EoS. Each model is tested against the data on the dimensionless coordinate distance to Type Ia Supernovae and radio galaxies, the shift and the acoustic peak parameters and the positions of the first three peaks in the anisotropy spectrum of the comic microwave background radation. We propose a possible interpretation of each model in the framework of scalar field quintessence determining the shape of the self interaction potential $V(\phi)$ that gives rise to each one of the considered thermodynamical EoS. As a general result, we demonstrate that replacing the perfect fluid EoS with more generar expressions gives both the possibility of successfully solving the problem of cosmic acceleration escaping the resort to phantom models.

We present a new family of spherically symmetric models for the luminous components of elliptical and spiral galaxies and their dark matter haloes. Our starting point is a general expression for the logarithmic slope $\alpha(r) = d\log{\rho}/d\log{r}$ from which most of the cuspy models yet available in literature may be derived. We then dedicate our attention to a particular set of models whose logarithmic slope is a power law function of the radius $r$ investigating in detail their dynamics assuming isotropy in the velocity space. While the basic properties (such as the density profile and the gravitational potential) may be expressed analytically, both the distribution function and the observable quantities (surface brightness and line of sight velocity dispersion) have to be evaluated numerically. We also consider the extension to anisotropic models trying two different parameterization. Since the model recently proposed by Navarro et al. (2004) as the best fit to their sample of numerically simulated haloes belongs to the family we present here, analytical approximations are given for the most useful quantities.

We use some of the recently released observational data to test the viability of two classes of minimally coupled scalar field models of quintessence with exponential potentials for which exact solutions of the Einstein equations are known. These models are very sturdy, depending on only one parameter - the Hubble constant. To compare predictions of our models with observations we concentrate on the following data: the power spectrum of the CMBR anisotropy as measured by WMAP, the publicly available data on type Ia supernovae, and the parameters of large scale structure determined by the 2-degree Field Galaxy Redshift Survey (2dFGRS). We use the WMAP data on the age of the universe and the Hubble constant to fix the free parameters in our models. We then show that the predictions of our models are consistent with the observed positions and relative heights of the first 3 peaks in the CMB power spectrum, with the energy density of dark energy as deduced from observations of distant type Ia supernovae, and with parameters of the large scale structure as determined by 2dFGRS, in particular with the average density of dark matter. Our models are also consistent with the results of the Sloan Digital Sky Survey (SDSS). Moreover, we investigate the evolution of matter density perturbations in our quintessential models, solve exactly the evolution equation for the density perturbations, and obtain an analytical expression for the growth index $f$. We verify that the approximate relation f ~ Omega_M^(alpha) also holds in our models.

The structure of the dark matter and the thermodynamical status of the hot gas in galaxy clusters is an interesting and widely discussed topic in modern astrophysics. Recently, Rasia et al. (2004) have proposed a new dynamical model for the mass density profile of clusters of galaxies as a result of a set of high resolution hydrodynamical simulations of structure formation. We investigate the lensing properties of this model evaluating the deflection angle, the lensing potential and the amplification of the images. We reserve particular attention to the structure and position of the critical curves in order to see whether this model is able to produce radial and tangential arcs. To this aim, we also investigate the effect of taking into account the brightest cluster galaxy in the lensing potential and the deviations from spherical symmetry mimicked by an external shear. We also analyze the implication of the gas density and temperature profiles of the Rasia et al.(2004) model on the properties of the X - ray emission and the comptonization parameter that determines the CMBR temperature decrement due to the Sunyaev - Zel'dovich effect.

We present a numerical method to estimate the lensing parameters and the Hubble constant H_0 from quadruply imaged gravitational lens systems. The lens galaxy is modeled using both separable deflection potentials and constant mass-to-light ratio profiles, while possible external perturbations have been taken into account introducing an external shear. The model parameters are recovered inverting the lens and the time delay ratio equations and imposing a set of physically motivated selection criteria. We investigate correlations among the model parameters and the Hubble constant. Finally, we apply the codes to the real lensed quasars PG 1115+080 and RX J0911+0551, and combine the results from these two systems to get H_0 = 56 +/- 23 km/(s Mpc). In addition, we are able to fit to the single systems a general elliptical potential with a non fixed angular part, and then we model the two lens systems with the same potential and a shared H_0: in this last case we estimate H_0=49_(-11)^(+6) km/(s Mpc).