We present here the Calar Alto Schmidt-Lemaitre Telescope (CASTLE) concept, a technology demonstrator for curved detectors, that will be installed at the Calar Alto Observatory (Spain). This telescope has a wide field of view (2.36x1.56 deg^2) and a design, optimised to generate a Point Spread Function with very low level wings and reduced ghost features, which makes it considerably less susceptible to several systematic effects usually affecting similar systems. These characteristics are particularly suited to study the low surface brightness Universe. CASTLE will be able to reach surface brightness orders of magnitude fainter than the sky background level and observe the extremely extended and faint features around galaxies such as tidal features, stellar halos, intra-cluster light, etc. CASTLE will also be used to search and detect astrophysical transients such as gamma ray bursts (GRB), gravitational wave optical counterparts, neutrino counterparts, etc. This will increase the number of precisely localized GRBs from 20% to 60% (in the case of Fermi/GMB GRBs).

We study the mass distribution and kinematics of the giant elliptical galaxy M87 (NGC 4486) using discrete chemo-dynamical, axisymmetric Jeans equation modelling. Our catalogue comprises 894 globular clusters (GCs) extending to a projected radius of $\sim 430$ kpc with line-of-sight velocities and colours, and Multi Unit Spectroscopic Explorer (MUSE) integral field unit data within the central $2.4$ kpc of the main galaxy. The gravitational potential for our models is a combination of a luminous matter potential with a varying mass-to-light ratio for the main galaxy, a supermassive black hole and a dark matter (DM) potential with a cusped or cored DM halo. The best-fitting models with either a cusped or a cored DM halo show no significant differences and both are acceptable. We obtain a total mass of $(2.16 \pm 0.38) \times 10^{13} M_{\odot}$ within $\sim$ 400 kpc. By including the stellar mass-to-light ratio gradient, the DM fraction increases from $\sim$ 26 percent (with no gradient) to $\sim$ 73 percent within $1\,R_e^{\rm maj}$ (major axis of half-light isophote, 14.2 kpc), and from $\sim$ 84 percent to $\sim$ 94 percent within $5\,R_e^{\rm maj}$ (71.2 kpc). Red GCs have moderate rotation with $V_{\rm max}/\sigma \sim$ 0.4, and blue GCs have weak rotation with $V_{\rm max}/\sigma \sim$ 0.1. Red GCs have tangential velocity dispersion anisotropy, while blue GCs are consistent with being nearly isotropic. Our results suggest that red GCs are more likely to be born in-situ, while blue GCs are more likely to be accreted.

The Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) is a blind narrow-band Halpha+[NII] imaging survey of the Virgo cluster carried out with MegaCam at the Canada-French-Hawaii telescope (CFHT). We use a new set of data extracted from VESTIGE to study the impact of the hostile cluster environment on the star formation process down to the scale of HII regions (~ 50 pc). HII regions are identified and their parameters measured using the HIIphot code on a sample of 114 late-type galaxies spanning a wide range in morphological type (Sa-Sd, Im, BCD), stellar mass (10^6.5 <= M_star <= 10^11 Mo), and star formation activity (10^-3 <= SFR <= 10 Mo yr^-1). Owing to the exquisite average resolution of the VESTIGE data (0.65 arcsec), we detect 11302 HII regions with an Halpha luminosity L(Halpha) >= 10^37 erg s^-1. We show that the typical number of HII regions in gas-stripped objects is significantly lower than in healthy late-types of similar stellar mass. We also show that in these gas-stripped galaxies the number of HII regions significantly drops outside the effective radius, suggesting that the quenching process occurs outside-in, in agreement with other multifrequency observations. These new results consistently confirm that the main mechanism responsible for the decrease of the star formation activity observed in cluster galaxies is ram pressure, allowing us to discard other milder processes such as starvation or strangulation unable to reproduce the observed radially truncated profiles.

We present the first detection of diffuse dust in the intra-cluster medium of the Virgo cluster out to $\sim$0.4 virial radii, and study the radial variation of its properties on a radial scale of the virial radius. Analysing near-UV - $i$ colours for a sample of $\sim12000$ background galaxies with redshifts $0.02 < z < 0.8$, we find significant colour reddening and relate it to variation in $E(B-V)$ values. The $E(B-V)$ mean profile shows a dust component characterised by an average reddening $E(B-V)\sim0.042 \pm 0.004$ mag within 1.5 degrees ($\sim0.3\, r_{vir}$) from the cluster centre. Assuming a Large Magellanic Cloud extinction law, we derive an average visual extinction $A_{V} = 0.14\pm 0.01$ for a total dust mass, $M_{d} = 2.5\pm0.2\times10^{9}M_{\odot}$, hence a dust-to-gas mass ratio $M_{d}/M_{g} = 3.0\pm 0.3 \times 10^{-4}$. Based on the upper limits on the flux density $\mathrm{I_{250\mu m} = 0.1\, MJy sr^{-1}} $ derived from $Herschel$ data, we estimate an upper limit for the dust temperature of $T_{d} \sim 10\, K$. However, similar densities can be obtained with dust at higher temperatures with lower emissivities. The Virgo cluster has diffuse dust in its intra-cluster medium characterised by different physical properties as those characterising the Milky Way dust. The diffuse dust in Virgo is transported into the cluster space through similar phenomena (stripping) as those building up the optical intra-cluster light, and it constitutes an additional cooling agent of the cluster gas.

The occurrence of planetary nebulae (PNe) in globular clusters (GCs) provides an excellent chance to study low-mass stellar evolution in a special (low-metallicity, high stellar density) environment. We report a systematic spectroscopic survey for the [O{\sc iii}] 5007 emission line of PNe in 1469 Virgo GCs and 121 Virgo ultra-compact dwarfs (UCDs), mainly hosted in the giant elliptical galaxies M87, M49, M86, and M84. We detected zero PNe in our UCD sample and discovered one PN ($M_{5007} = -4.1$ mag) associated with an M87 GC. We used the [O{\sc iii}] detection limit for each GC to estimate the luminosity-specific frequency of PNe, $\alpha$, and measured $\alpha$ in the Virgo cluster GCs to be $\alpha \sim 3.9_{-0.7}^{+5.2}\times 10^{-8}\mathrm{PN}/L_\odot$. $\alpha$ in Virgo GCs is among the lowest values reported in any environment, due in part to the large sample size, and is 5--6 times lower than that for the Galactic GCs. We suggest that $\alpha$ decreases towards brighter and more massive clusters, sharing a similar trend as the binary fraction, and the discrepancy between the Virgo and Galactic GCs can be explained by the observational bias in extragalactic surveys toward brighter GCs. This low but non-zero efficiency in forming PNe may highlight the important role played by binary interactions in forming PNe in GCs. We argue that a future survey of less massive Virgo GCs will be able to determine whether PN production in Virgo GCs is governed by internal process (mass, density, binary fraction), or is largely regulated by external environment.

We have observed the giant elliptical galaxy M87 during the Virgo Environmental Survey Tracing Galaxy Evolution (VESTIGE), a blind narrow-band Halpha+[NII] imaging survey of the Virgo cluster carried out with MegaCam at the Canada French Hawaii Telescope (CFHT). The deep narrow-band image confirmed the presence of a filament of ionised gas extending up to ~ 3 kpc in the north-western direction and ~ 8 kpc to the south-east, with a couple of plumes of ionised gas, the weakest of which, at ~ 18 kpc from the nucleus, was previously unknown. The analysis of deep optical images taken from the NGVS survey confirms that this gas filament is associated with dust seen in absorption which is now detected up to ~ 2.4 kpc from the nucleus. We have also analysed the physical and kinematical properties of the ionised gas filament using deep IFU MUSE data covering the central 4.8 x 4.8 kpc^2 of the galaxy. The spectroscopic data confirms a perturbed kinematics of the ionised gas, with differences in velocity of ~ 700-800 km s^-1 on scales of < 1 kpc. The analysis of 2D diagnostic diagrams and the observed relationship between the shock-sensitive [OI]/Halpha line ratio and the velocity dispersion of the gas suggest that the gas is shock-ionised.

We present a kinematic study of a sample of 298 planetary nebulas (PNs) in the outer halo of the central Virgo galaxy M87 (NGC~4486). The line-of-sight velocities of these PNs are used to identify sub-components, to measure the angular momentum content of the smooth M87 halo, and to constrain the orbital distribution of the stars at these large radii. We find that the surface density profile of the ICPNS at $100\, \mathrm{kpc}$ radii has a shallow logarithmic slope, $-\alpha_{\rm ICL}\simeq -0.8$, dominating the light at the largest radii. The angular momentum-related $\lambda_{\mathrm{R}}$ profile for the smooth halo remains in the slow rotator regime, out to 135 kpc average ellipse radius. The velocity dispersione profile is measured to decrease sharply to $\sigma_{\mathrm{halo}} \sim 100\, \mathrm{kms}^{-1}$ at $R_{\rm avg}=135$ kpc. Such a steeply decline of the velocity dispersion and the surface density profile of the smooth halo can be reconciled with the circular velocity curve inferred from assuming hydrostatic equilibrium for the hot X-ray gas if the orbit distribution of the smooth M87 halo changes strongly from approximately isotropic within $R_{\rm avg}\sim 60$ kpc to very radially anisotropic at the largest distances probed. The extended LOSVD of the PNs in the M87 halo allows the identification of several subcomponents: the ICPNs, the "crown" accretion event, and the smooth M87 halo. In galaxies like M87, the presence of these sub-components needs to be taken into account to avoid systematic biases in estimating the total enclosed mass. The dynamical structure inferred from the velocity dispersion profile indicates that the smooth halo of M87 steepens beyond $R_{\rm avg}=60$ kpc and becomes strongly radially anisotropic, and that the velocity dispersion profile is consistent with the X-ray circular velocity curve at these radii without non-thermal pressure effects.

Intra-cluster (IC) populations are expected to be a natural result of the hierarchical assembly of clusters, yet their low space densities make them difficult to detect and study. We present the first definitive kinematic detection of an IC population of globular clusters (GCs) in the Virgo cluster, around the central galaxy, M87. This study focuses on the Virgo core for which the combination of NGVS photometry and follow-up spectroscopy allows us to reject foreground star contamination and explore GC kinematics over the full Virgo dynamical range. The GC kinematics changes gradually with galactocentric distance, decreasing in mean velocity and increasing in velocity dispersion, eventually becoming indistinguishable from the kinematics of Virgo dwarf galaxies at $\mathrm{R>320\, kpc}$. By kinematically tagging M87 halo and intra-cluster GCs we find that 1) the M87 halo has a smaller fraction ($52\pm3\%$) of blue clusters with respect to the IC counterpart ($77\pm10\%$), 2) the $(g'-r')_{0}$ vs $(i'-z')_{0}$ color-color diagrams reveal a galaxy population that is redder than the IC population that may be due to a different composition in chemical abundance and progenitor mass, and 3) the ICGC distribution is shallower and more extended than the M87 GCs, yet still centrally concentrated. The ICGC specific frequency, $S_{N,\mathrm{ICL}}=10.2\pm4.8$, is consistent with what is observed for the population of quenched, low-mass galaxies within 1~Mpc from the cluster's center. The IC population at Virgo's center is thus consistent with being an accreted component from low-mass galaxies tidally stripped or disrupted through interactions, with a total mass of $\mathrm{M_{ICL,tot}=10.8\pm0.1\times10^{11}M_{\odot}}$.

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]

M49 (NGC 4472) is the dominant galaxy in subcluster B of the Virgo Cluster, and a benchmark for studying the build-up of the extended halos of brightest group galaxies in the outskirts of galaxy clusters. We investigate the kinematics in the outer halo of M49, look for substructures, and describe the transition to the surrounding intra-group light. As kinematic tracers we use planetary nebulae (PNe), combining kinematics from the extended Planetary Nebula Spectrograph (PN.S) early-type galaxy survey with our recent deep photometric sample. We study the position-velocity-plane for bright and faint PN populations out to 95 kpc radius, and employ a multi-Gaussian model for the velocity distribution to identify stellar populations with distinct kinematics and histories. We report the detection of stellar-kinematic substructure associated with the interaction of M49 with the dwarf irregular galaxy VCC 1249. We find two kinematically distinct PN populations associated with the main M49 halo and the extended intra-group light (IGL). These have velocity dispersions $\sigma_{halo}\simeq 170$ km/s and $\sigma_{IGL} \simeq400$ km/s at 10-80 kpc radii. The overall luminosity profile and velocity dispersion at $\sim80$ kpc are consistent with a flat circular velocity curve extrapolated from X-ray observations. The dispersion of the PNe associated with the IGL joins onto that of the satellite galaxies in subcluster B at $\sim100$ kpc radius. This is the first time that the transition from halo to IGL is observed based on the velocities of individual stars. Therefore the halo of M49, consisting of at least three distinct components, has undergone an extended accretion history within its parent group potential. The blue colours of the IGL component are consistent with a population of stars formed in low-mass galaxies at redshift $\sim0.5$ that has since evolved passively, as suggested by other data.

The Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) is a blind narrow-band H$\alpha$+[NII] imaging survey carried out with MegaCam at the CFHT. During pilot observations we have observed NGC 4330, an intermediate mass, edge-on star forming spiral currently falling into the core of the Virgo cluster. New deep observations revealed a low surface brightness 10 kpc tail exhibiting a peculiar filamentary structure. The filaments are remarkably parallel one another and clearly indicate the direction of motion of the galaxy in the Virgo potential. Motivated by the detection of these features, indicating ongoing gas stripping, we collected literature photometry in 15 bands from the far-UV to the far-IR and deep optical long slit spectroscopy using the FORS2 instrument at the ESO Very Large Telescope. Using a newly developed Monte Carlo code that jointly fits spectroscopy and photometry, we reconstructed the star formation histories in apertures along the major axis of the galaxy. Our results have been validated against the output of CIGALE, a fitting code which has been previously used for similar studies. We found a clear outside-in gradient with radius of the time when the quenching event started: the outermost radii have been stripped 500 Myr ago, while the stripping has reached the inner 5 kpc from the center in the last 100 Myr. Regions at even smaller radii are currently still forming stars fueled by the presence of HI and H2 gas. When compared to statistical studies of the quenching timescales in the local Universe we find that ram pressure stripping of the cold gas is an effective mechanism to reduce the transformation times for galaxies falling into massive clusters. Future systematic studies of all the active galaxies observed by VESTIGE in the Virgo cluster will extend these results to a robust statistical framework.

The Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) is a blind narrow-band Halpha+[NII] imaging survey carried out with MegaCam at the Canada-France-Hawaii Telescope. The survey covers the whole Virgo cluster region from its core to one virial radius (104 deg^2). The sensitivity of the survey is of f(Halpha) ~ 4 x 10^-17 erg sec-1 cm^-2 (5 sigma detection limit) for point sources and Sigma (Halpha) ~ 2 x 10^-18 erg sec^-1 cm^-2 arcsec^-2 (1 sigma detection limit at 3 arcsec resolution) for extended sources, making VESTIGE the deepest and largest blind narrow-band survey of a nearby cluster. This paper presents the survey in all its technical aspects, including the survey design, the observing strategy, the achieved sensitivity in both the narrow-band Halpha+[NII] and in the broad-band r filter used for the stellar continuum subtraction, the data reduction, calibration, and products, as well as its status after the first observing semester. We briefly describe the Halpha properties of galaxies located in a 4x1 deg^2 strip in the core of the cluster north of M87, where several extended tails of ionised gas are detected. This paper also lists the main scientific motivations of VESTIGE, which include the study of the effects of the environment on galaxy evolution, the fate of the stripped gas in cluster objects, the star formation process in nearby galaxies of different type and stellar mass, the determination of the Halpha luminosity function and of the Halpha scaling relations down to ~ 10^6 Mo stellar mass objects, and the reconstruction of the dynamical structure of the Virgo cluster. This unique set of data will also be used to study the HII luminosity function in hundreds of galaxies, the diffuse Halpha+[NII] emission of the Milky Way at high Galactic latitude, and the properties of emission line galaxies at high redshift.

It is thought that fossil systems are relics of structure formation in the primitive Universe. They are galaxy aggregations that have assembled their mass at high redshift with few or no subsequent accretion. Observationally these systems are selected by large magnitude gaps between their 1st and 2nd ranked galaxies. Nevertheless, there is still debate over whether or not this observational criterium selects dynamically evolved ancient systems. We have studied the properties of the nearby fossil group RXJ075243.6+455653 in order to understand the mass assembly of this system. Deep spectroscopic observations allow us to construct the galaxy luminosity function (LF) of RXJ075243.6+455653 down to M*+ 6. The analysis of the faint-end of the LF in groups and clusters provides valuable information about the mass assembly of the system. In addition, we have analyzed the nearby large-scale structure around this group. We identified 26 group members within r200=0.9 Mpc. The LF of the group shows a flat faint-end slope ( -1.08 +/- 0.33). This low density of dwarf galaxies is confirmed by the low value of the dwarf-to-giant ratio (DGR = 0.99 +/- 0.49) for this system. Both the lack of dwarf galaxies and the low luminosity of the BGG suggests that RXJ075243.6+455653 still has to accrete mass from its nearby environment. This mass accretion will be achieved because it is the dominant structure of a rich environment formed by several groups of galaxies (15) within 7 Mpc from the group center and with +/- 1000$ km/s. RXJ075243.6+455653 is a group of galaxies that has not yet completed the process of its mass assembly. This new mass accretion will change the fossil state of the group. This group is an example of a galaxy aggregation selected by a large magnitude gap but still in the process of the accretion of its mass (Abridged).

The galaxy M49 (NGC 4472) is the brightest early-type galaxy in the Virgo Cluster. It is located in Subcluster B and has an unusually blue, metal-poor outer halo. Planetary nebulae (PNe) are excellent tracers of diffuse galaxy and intragroup light. We present a photometric survey of PNe in the galaxy's extended halo to characterise its PN population, as well as the surrounding intragroup light (IGL) of the Subcluster B. PNe were identified based on their bright [OIII]5007 \AA\ emission and absence of a broad-band continuum. We identify 738 PNe out to a radius of 155 kpc from M49's centre from which we define a complete sample of 624 PNe within a limiting magnitude of m_5007=28.8. Comparing the PN number density to the broad-band stellar surface brightness profile, we find a variation of the PN-specific frequency (alpha-parameter) with radius. The outer halo beyond 60 kpc has a 3.2 times higher alpha-parameter compared to the main galaxy halo, which is likely due to contribution from the surrounding blue IGL. We use the Planetary Nebulae Luminosity Function (PNLF) as an indicator of distance and stellar population. Its slope, which correlates empirically with galaxy type, varies within the inner halo. In the eastern quadrant of M49, the PNLF slope is shallower, indicating an additional localised, bright PN population following an accretion event, likely that of the dwarf irregular galaxy VCC1249. We also determined a distance modulus of mu = 31.29+/-0.08 for M49, corresponding to a physical distance of 18.1+/-0.6 Mpc, which agrees with a recent surface-brightness fluctuations distance. The PN populations in the outer halo of M49 are consistent with the presence of a main Sersic galaxy halo with a slight (B-V) colour gradient of 10${}^{-4}$ mag/arcsec surrounded by intragroup light with a very blue colour of (B-V)=0.25 and a constant surface brightness mu_V=28.0 mag/arcsec${}^2$.

The diffuse extended outer regions of galaxies are hard to study because they are faint, with typical surface brightness of 1% of the dark night sky. We can tackle this problem by using resolved star tracers which remain visible at large distances from the galaxy centres. This article describes the use of Planetary Nebulae as tracers and the calibration of their properties as indicators of the star formation history, mean age and metallicity of the parent stars in the Milky Way and Local Group galaxies . We then report on the results from a deep, extended, planetary nebulae survey in a 0.5 sqdeg region centred on the brightest cluster galaxy NGC 4486 (M87) in the Virgo cluster core, carried out with SuprimeCam@Subaru and FLAMES-GIRAFFE@VLT. Two PN populations are identified out to 150 kpc distance from the centre of M87. One population is associated with the M87 halo and the second one with the intracluster light in the Virgo cluster core. They have different line-of-sight velocity and spatial distributions, as well as different planetary nebulae specific frequencies and luminosity functions. The intracluster planetary nebulae in the surveyed region correspond to a luminosity of four times the luminosity of the Large Magellanic Cloud. The M87 halo planetary nebulae trace an older, more metal-rich, parent stellar population. A substructure detected in the projected phase-space of the line-of-sight velocity vs. major axis distance for the M87 halo planetary nebulae provides evidence for the recent accretion event of a satellite galaxy with luminosity twice that of M33. The satellite stars were tidally stripped about 1 Gyr ago, and reached apocenter at a major axis distance of 60-90 kpc from the centre of M87. The M87 halo is still growing significantly at the distances where the substructure is detected.

We present a spectroscopic study of 287 Planetary Nebulas (PNs) in a total area of ~0.4 deg^2 around the BCG M87 in Virgo A. With these data we can distinguish the stellar halo from the co-spatial intracluster light (ICL). PNs were identified from their narrow and symmetric redshifted lambda 5007\4959 Angstrom [OIII] emission lines, and the absence of significant continuum. We implement a robust technique to measure the halo velocity dispersion from the projected phase-space to identify PNs associated with the M87 halo and ICL. The velocity distribution of the spectroscopically confirmed PNs is bimodal, containing a narrow component centred on the systemic velocity of the BCG and an off-centred broader component, that we identify as halo and ICL, respectively. Halo and ICPN have different spatial distributions: the halo PNs follow the galaxy's light, whereas the ICPNs are characterised by a shallower power-law profile. The composite PN number density profile shows the superposition of different PN populations associated with the M87 halo and the ICL, characterised by different PN alpha-parameters, the ICL contributing ~3 times more PNs per unit light. Down to m_5007=28.8, the M87 halo PN luminosity function (PNLF) has a steeper slope towards faint magnitudes than the IC PNLF, and both are steeper than the standard PNLF for the M31 bulge. Moreover, the IC PNLF has a dip at ~1-1.5 mag fainter than the bright cutoff, reminiscent of the PNLFs of systems with extended star formation history. The M87 halo and the Virgo ICL are dynamically distinct components with different density profiles and velocity distribution. The different alpha values and PNLF shapes of the halo and ICL indicate distinct parent stellar populations, consistent with the existence of a gradient towards bluer colours at large radii. These results reflect the hierarchical build-up of the Virgo cluster.

We present kinematic and photometric evidence for an accretion event in the halo of the cD galaxy M87 in the last Gyr. Using velocities for ~300 planetary nebulas (PNs) in the M87 halo, we identify a chevron-like substructure in the PN phase-space. We implement a probabilistic Gaussian mixture model to identify the PNs that belong to the chevron. From analysis of deep V-band images of M87, we find that the region with the highest density of PNs associated to the chevron, is a crown-shaped substructure in the optical light. We assign a total of N_(PN,sub)=54 to the substructure, which extends over ~50 kpc along the major axis where we also observe radial variations of the ellipticity profile and a colour gradient. The substructure has highest surface brightness in a 20kpc x 60kpc region around 70 kpc in radius. In this region, it causes an increase in surface brightness by >60%. The accretion event is consistent with a progenitor galaxy with a V-band luminosity of L=2.8\pm1.0 x 10^9 L_(sun,V), a colour of (B-V)=0.76\pm0.05, and a stellar mass of M=6.4\pm2.3 x 10^9 M_sun. The accretion of this progenitor galaxy has caused an important modification of the outer halo of M87 in the last Gyr. By itself it is strong evidence that the galaxy's cD halo is growing through the accretion of smaller galaxies as predicted by hierarchical galaxy evolution models.

We investigate the diffuse light in the outer regions of the nearby elliptical galaxy M87 in the Virgo cluster, using planetary nebulas (PNs) as tracers. The surveyed areas (0.43 squared degrees) cover M87 up to a radial distance of 150 kpc, in the ransition region between galaxy halo and intracluster light (ICL). All PNs are identified through the on-off band technique using automatic selection criteria based on the distribution of the detected sources in the colour-magnitude diagram and the properties of their point-spread function. We extract a catalogue of 688 objects down to m_5007=28.4, with an estimated residual contamination from foreground stars and background Lyalpha galaxies, which amounts to ~35% of the sample. This is one of the largest extragalactic PN samples in number of candidates, magnitude depth, and radial extent, which allows us to carry out an unprecedented photometric study of the PN population in the outer regions of M87. We find that the logarithmic density profile of the PN distribution is shallower than the surface brightness profile at large radii. This behaviour is consistent with the superposition of two components associated with the halo of M87 and with the ICL, which have different luminosity specific PN numbers, the ICL contributing three times more PNs per unit light. Because of the depth of this survey we are also able to study the shape of the PN luminosity function (PNLF) in the outer regions of M87. We find a slope for the PNLF that is steeper at fainter magnitudes than the standard analytical PNLF formula and adopt a generalised model that treats the slope as a free parameter. Comparing the PNLF of M87 and the M31 bulge, both normalised by the sampled luminosity, the M87 PNLF contains fewer bright PNs and has a steeper slope towards fainter magnitudes.

The luminosity function of planetary nebulae populations in galaxies within 10-15 Mpc distance has a cut-off at bright magnitudes and a functional form that is observed to be invariant in different galaxy morphological types. Thus it is used as a secondary distance indicator in both early and late-type galaxies. Recent deep surveys of planetary nebulae populations in brightest cluster galaxies (BCGs) seem to indicate that their luminosity functions deviate from those observed in the nearby galaxies. We discuss the evidence for such deviations in Virgo, and indicate which physical mechanisms may alter the evolution of a planetary nebula envelope and its central star in the halo of BCGs. We then discuss preliminary results for distances for the Virgo, Hydra I and Coma clusters based on the observed planetary nebulae luminosity functions.