Deep imaging of the diffuse light emitted by the stellar fine structures and outer halos around galaxies is now often used to probe their past mass assembly. Because the extended halos survive longer than the relatively fragile tidal features, they trace more ancient mergers. We use images reaching surface brightness limits as low as 28.5-29 mag.arcsec-2 (g-band) to obtain light and color profiles up to 5-10 effective radii of a sample of nearby early-type galaxies. They were acquired with MegaCam as part of the CFHT MATLAS large programme. These profiles may be compared to those produced by simulations of galaxy formation and evolution, once corrected for instrumental effects. Indeed they can be heavily contaminated by the scattered light caused by internal reflections within the instrument. In particular, the nucleus of galaxies generates artificial flux in the outer halo, which has to be precisely subtracted. We present a deconvolution technique to remove the artificial halos that makes use of very large kernels. The technique based on PyOperators is more time efficient than the model-convolution methods also used for that purpose. This is especially the case for galaxies with complex structures that are hard to model. Having a good knowledge of the Point Spread Function (PSF), including its outer wings, is critical for the method. A database of MegaCam PSF models corresponding to different seeing conditions and bands was generated directly from the deep images. It is shown that the difference in the PSFs in different bands causes artificial changes in the color profiles, in particular a reddening of the outskirts of galaxies having a bright nucleus. The method is validated with a set of simulated images and applied to three representative test cases: NGC 3599, NGC 3489, and NGC 4274, and exhibiting for two of them a prominent ghost halo. The method successfully removes it.

Diffuse Galactic light has been observed in the optical since the 1930s. We propose that, when observed in the optical with deep imaging surveys, it can be used as a tracer of the turbulent cascade in the diffuse interstellar medium (ISM), down to scales of about 1 arcsec. Here we present a power spectrum analysis of the dust column density of a diffuse cirrus at high Galactic latitude (l ~ 198 deg, b ~ 32 deg) as derived from the combination of a MegaCam g-band image, obtained as part of the MATLAS large programme at the CFHT, with Planck radiance and Wise 12 micron data. The combination of these three datasets have allowed us to compute the density power spectrum of the HI over scales of more than three orders of magnitude. We found that the density field is well described by a single power law over scales ranging from 0.01 to 50 pc. The exponent of the power spectrum, gamma=-2.9 +- 0.1, is compatible with what is expected for thermally bi-stable and turbulent HI. We did not find any steepening of the power spectrum at small scales indicating that the typical scale at which turbulent energy is dissipated in this medium is smaller than 0.01pc. The ambipolar diffusion scenario that is usually proposed as the main dissipative agent, is consistent with our data only if the density of the cloud observed is higher than the typical values assumed for the cold neutral medium gas. We discuss the new avenue offered by deep optical imaging surveys for the study of the low density ISM structure and turbulence.

NGC 4203 is a nearby early-type galaxy surrounded by a very large, low-column-density HI disc. In this paper we study the star formation efficiency in the gas disc of NGC 4203 by using the UV, deep optical imaging and infrared data. We confirm that the HI disc consists of two distinct components: an inner star forming ring with radius from $\sim$ 1 to $\sim$ 3 R$_{eff}$, and an outer disc. The outer HI disc is 9 times more massive than the inner HI ring. At the location of the inner HI ring we detect spiral-like structure both in the deep $g'-r'$ image and in the 8 $\mu$m $Spitzer$-IRAC image, extending in radius up to $\sim$ 3 R$_{eff}$. These two gas components have a different star formation efficiency likely due to the different metallicity and dust content. The inner component has a star formation efficiency very similar to the inner regions of late-type galaxies. Although the outer component has a very low star formation efficiency, it is similar to that of the outer regions of spiral galaxies and dwarfs. We suggest that these differences can be explained with different gas origins for the two components such as stellar mass loss for the inner HI ring and accretion from the inter galactic medium (IGM) for the outer HI disc. The low level star formation efficiency in the outer HI disc is not enough to change the morphology of NGC 4203, making the depletion time of the HI gas much too long.

We investigate the evolution of stellar population gradients from $z=2$ to $z=0$ in massive galaxies at large radii ($r > 2R_{\mathrm{eff}}$) using ten cosmological zoom simulations of halos with $6 \times 10^{12} M_{\odot} < M_{\mathrm{halo}} < 2 \times 10^{13}M_{\odot}$. The simulations follow metal cooling and enrichment from SNII, SNIa and AGB winds. We explore the differential impact of an empirical model for galactic winds that reproduces the mass-metallicity relation and its evolution with redshift. At larger radii the galaxies, for both models, become more dominated by stars accreted from satellite galaxies in major and minor mergers. In the wind model, fewer stars are accreted, but they are significantly more metal poor resulting in steep global metallicity ($\langle \nabla Z_{\mathrm{stars}} \rangle= -0.35$ dex/dex) and color (e.g. $\langle \nabla g-r \rangle = -0.13$ dex/dex) gradients in agreement with observations. In contrast, colour and metallicity gradients of the models without winds are inconsistent with observations. Age gradients are in general mildly positive at $z=0$ ($\langle \nabla Age_{\mathrm{stars}} \rangle= 0.04$ dex/dex) with significant differences between the models at higher redshift. We demonstrate that for the wind model, stellar accretion is steepening existing in-situ metallicity gradients by about 0.2 dex by the present day and helps to match observed gradients of massive early-type galaxies at large radii. Colour and metallicity gradients are significantly steeper for systems which have accreted stars in minor mergers, while galaxies with major mergers have relatively flat gradients, confirming previous results. This study highlights the importance of stellar accretion for stellar population properties of massive galaxies at large radii, which can provide important constraints for formation models.

Galactic archeology based on star counts is instrumental to reconstruct the past mass assembly of Local Group galaxies. The development of new observing techniques and data-reduction, coupled with the use of sensitive large field of view cameras, now allows us to pursue this technique in more distant galaxies exploiting their diffuse low surface brightness (LSB) light. As part of the Atlas3D project, we have obtained with the MegaCam camera at the Canada-France Hawaii Telescope extremely deep, multi--band, images of nearby early-type galaxies. We present here a catalog of 92 galaxies from the Atlas3D sample, that are located in low to medium density environments. The observing strategy and data reduction pipeline, that achieve a gain of several magnitudes in the limiting surface brightness with respect to classical imaging surveys, are presented. The size and depth of the survey is compared to other recent deep imaging projects. The paper highlights the capability of LSB--optimized surveys at detecting new prominent structures that change the apparent morphology of galaxies. The intrinsic limitations of deep imaging observations are also discussed, among those, the contamination of the stellar halos of galaxies by extended ghost reflections, and the cirrus emission from Galactic dust. The detection and systematic census of fine structures that trace the present and past mass assembly of ETGs is one of the prime goals of the project. We provide specific examples of each type of observed structures -- tidal tails, stellar streams and shells --, and explain how they were identified and classified. We give an overview of the initial results. The detailed statistical analysis will be presented in future papers.