We use the SPHINX suite of high-resolution cosmological radiation hydrodynamics simulations to study how spatially and temporally inhomogeneous reionization impacts the baryonic content of dwarf galaxies and cosmic filaments. The SPHINX simulations simultaneously model an inhomogeneous reionization, follow the escape of ionising radiation from thousands of galaxies, and resolve haloes well below the atomic cooling threshold. This makes them an ideal tool for examining how reionization impacts star formation and the gas content of dwarf galaxies. We compare simulations with and without stellar radiation to isolate the effects of radiation feedback from that of supernova, cosmic expansion, and numerical resolution. We find that the gas content of cosmic filaments can be reduced by more than 80% following reionization. The gas inflow rates into haloes with $M_{vir}<10^8M_{\odot}$ are strongly affected and are reduced by more than an order of magnitude compared to the simulation without reionization. A significant increase in gas outflow rates is found for halo masses $M_{vir}<7\times10^7M_{\odot}$. Our simulations show that inflow suppression, rather than photoevaporation, is the dominant mechanism by which the baryonic content of high-redshift dwarf galaxies is regulated. At fixed redshift and halo mass, there is a large scatter in the halo baryon fractions that is entirely dictated by the timing of reionization in the local region surrounding a halo which can change by $\Delta z>3$ at fixed mass. Finally, although the gas content of high-redshift dwarf galaxies is significantly impacted by reionization, we find that most haloes with $M_{vir}<10^8M_{\odot}$ can remain self-shielded and form stars long after reionization, until their local gas reservoir is depleted, suggesting that local group dwarf galaxies do not necessarily exhibit star formation histories that peak prior to $z=6$...

Resonant lines are powerful probes of the interstellar and circumgalactic medium of galaxies. Their transfer in gas being a complex process, the interpretation of their observational signatures, either in absorption or in emission, is often not straightforward. Numerical radiative transfer simulations are needed to accurately describe the travel of resonant line photons in real and in frequency space, and to produce realistic mock observations. This paper introduces RASCAS, a new public 3D radiative transfer code developed to perform the propagation of any resonant line in numerical simulations of astrophysical objects. RASCAS was designed to be easily customisable and to process simulations of arbitrarily large sizes on large supercomputers. RASCAS performs radiative transfer on an adaptive mesh with an octree structure using the Monte Carlo technique. RASCAS features full MPI parallelisation, domain decomposition, adaptive load-balancing, and a standard peeling algorithm to construct mock observations. The radiative transport of resonant line photons through different mixes of species (e.g. \ion{H}{i}, \ion{Si}{ii}, \ion{Mg}{ii}, \ion{Fe}{ii}), including their interaction with dust, is implemented in a modular fashion to allow new transitions to be easily added to the code. RASCAS is very accurate and efficient. It shows perfect scaling up to a minimum of a thousand cores. It has been fully tested against radiative transfer problems with analytic solutions and against various test cases proposed in the literature. Although it was designed to describe accurately the many scatterings of line photons, RASCAS may also be used to propagate photons at any wavelength (e.g. stellar continuum or fluorescent lines), or to cast millions of rays to integrate the optical depths of ionising photons, making it highly versatile.

Several deep observations such as those carried out at the Canada-France-Hawaii Telescope (CFHT) have revealed prominent Low Surface Brightness (LSB) fine structures that change the apparent morphology of galaxies. Previous photometry surveys have developed observational techniques which exploit the diffuse light detected in the external regions of galaxies. In these studies the outer perturbations have been identified and classified like tidal tails, stellar streams, and shells. These structures are tracers of interacting and merging events and they keep a memory of the mass assembly of galaxies. Cosmological numerical simulations are needed to estimate their visibility time-scale (among other properties) in order to reconstruct the past merger history of galaxies. In the present work, we analyze a hydrodynamical cosmological simulation to build a comprehensive interpretation of the properties of fine structures. We make a census of several types of LSB fine structures by visual inspection of individual snapshots at various time. We reconstruct the evolution of the number of fine structures detected around an early-type galaxy and we compare with the merger history of the galaxy. We find that most of fine structures are associated with major and intermediate mass merger events. Their survival time scale ranges between 0.7 and 4 Gyr. Shells and streams remain visible for a longer time than tidal tails. These estimates of survival times provide clues to interpret the shape and frequency of fine structures observed in deep images in terms of mass assembly. We find that the detectability of stellar streams is the most sensitive to the surface brightness limit. We see 2-3 times more streams with a surface brightness cut of 33 mag arcsec$^{-2}$ than with 29 mag arcsec$^{-2}$. The detection of shells display a strong dependence on the projection angle.

Understanding the escape of Lyman continuum (LyC) and Lyman alpha (Lya) photons from molecular clouds is one of the keys to constraining the reionization history of the Universe and the evolution of galaxies at high redshift. Using a set of radiation-hydrodynamic simulations with adaptive mesh refinement, we investigate how photons propagate and escape from turbulent clouds with different masses, star formation efficiencies (SFEs), and metallicities, as well as with different models of stellar spectra and supernova feedback. We find that the escape fractions in both LyC and Lya are generally increasing with time if the cloud is efficiently dispersed by radiation and supernova feedback. When the total SFE is low (1% of the cloud mass), 0.1-5% of LyC photons leave the metal-poor cloud, whereas the fractions increase to 20-70% in clouds with a 10% SFE. LyC photons escape more efficiently if gas metallicity is lower, if the upper mass limit in the stellar initial mass function is higher, if binary interactions are allowed in the evolution of stars, or if additional strong radiation pressure, such as Lya pressure, is present. As a result, the number of escaping LyC photons can easily vary by a factor of $\sim4$ on cloud scales. The escape fractions of Lya photons are systemically higher (60-80%) than those of LyC photons despite large optical depths at line centre ($\tau_0\sim10^6-10^9$). Scattering of Lya photons is already significant on cloud scales, leading to double-peaked profiles with peak separations of $v_{\rm sep}\sim400\,{\rm km\,s^{-1}}$ during the initial stage of the cloud evolution, while it becomes narrower than $v_{\rm sep} \le 150 \, {\rm km\,s^{-1}}$ in the LyC bright phase. Comparisons with observations of low-redshift galaxies suggest that Lya photons require further interactions with neutral hydrogen to reproduce their velocity offset for a given LyC escape fraction.

We present a study of the [OIII]5007/[OII]3727 (O32) ratios of star-forming galaxies drawn from MUSE data spanning a redshift range 0.28<z<0.85. Recently discovered Lyman continuum (LyC) emitters have extremely high oxygen line ratios: O32>4. Here we aim to understand the properties and the occurrences of galaxies with such high line ratios. Combining data from several MUSE GTO programmes, we select a population of star-forming galaxies with bright emission lines, from which we draw 406 galaxies for our analysis based on their position in the z-dependent star formation rate (SFR) - stellar mass (M*) plane. Out of this sample 15 are identified as extreme oxygen emitters based on their O32 ratios (3.7%) and 104 galaxies have O32>1 (26%). Our analysis shows no significant correlation between M*, SFR, and the distance from the SFR-M* relation with O32. We find a decrease in the fraction of galaxies with O32>1 with increasing M*, however, this is most likely a result of the relationship between O32 and metallicity, rather than between O32 and M*. We draw a comparison sample of local analogues with <z>~0.03 from SDSS, and find similar incidence rates for this sample. In order to investigate the evolution in the fraction of high O32 emitters with redshift, we bin the sample into three redshift subsamples of equal number, but find no evidence for a dependence on redshift. Furthermore, we compare the observed line ratios with those predicted by nebular models with no LyC escape and find that most of the extreme oxygen emitters can be reproduced by low metallicity models. The remaining galaxies are likely LyC emitter candidates. Finally, based on a comparison between electron temperature estimates from the [OIII4363]/[OIII]5007 ratio of the extreme oxygen emitters and nebular models, we argue that the galaxies with the most extreme O32 ratios have young light-weighted ages.

We present the SPHINX suite of cosmological adaptive mesh refinement simulations, the first radiation-hydrodynamical simulations to simultaneously capture large-scale reionization and the escape of ionizing radiation from thousands of resolved galaxies. Our $5$ and $10$ co-moving Mpc volumes resolve haloes down to the atomic cooling limit and model the inter-stellar medium with better than $\approx10$ pc resolution. The project has numerous goals in improving our understanding of reionization and making predictions for future observations. In this first paper we study how the inclusion of binary stars in computing stellar luminosities impacts reionization, compared to a model that includes only single stars. Owing to the suppression of galaxy growth via strong feedback, our galaxies are in good agreement with observational estimates of the galaxy luminosity function. We find that binaries have a significant impact on the timing of reionization: with binaries, our boxes are $99.9$ percent ionized by volume at $z\approx 7$, while without them our volumes fail to reionize by $z=6$. These results are robust to changes in volume size, resolution, and feedback efficiency. The escape of ionizing radiation from individual galaxies varies strongly and frequently. On average, binaries lead to escape fractions of $\approx 7-10$ percent, about $3.5$ times higher than with single stars only. The higher escape fraction is a result of a shallower decline in ionizing luminosity with age, and is the primary reason for earlier reionization, although the higher integrated luminosity with binaries also plays a sub-dominant role.

The physical origin of the near-ultraviolet MgII emission remains an under-explored domain, contrary to more typical emission lines detected in the spectra of star-forming galaxies. We explore the nebular and physical properties for a sample of 381 galaxies between 0.70 < z < 2.34 drawn from the MUSE Hubble Ultra Deep Survey. The spectra of these galaxies show a wide variety of profiles of the MgII 2796,2803 resonant doublet, from absorption to emission. We present a study on the main drivers for the detection of MgII emission in galaxy spectra. By exploiting photoionization models we verified that the emission-line ratios observed in galaxies with MgII in emission are consistent with nebular emission from HII regions. From a simultaneous analysis of MUSE spectra and ancillary HST information via spectral energy distribution (SED) fitting, we find that galaxies with MgII in emission have lower stellar masses, smaller sizes, bluer spectral slopes and lower optical depth than those with absorption. This leads us to suggest that MgII emission is a potential tracer of physical conditions not merely related to those of the ionized gas. We show that these differences in MgII emission/absorption can be explained in terms of a higher dust and neutral gas content in the interstellar medium (ISM) of galaxies showing MgII in absorption, confirming the extreme sensitivity of MgII to the presence of the neutral ISM. We conclude with an analogy between the MgII doublet and the Ly-alpha line, due to their resonant nature. Further investigations with current and future facilities, including JWST, are promising as the detection of MgII emission and its potential connection with Ly-alpha could provide new insights on the ISM content in the early Universe.

Understanding the origin of strong galactic outflows and the suppression of star formation in dwarf galaxies is a key problem in galaxy formation. Using a set of radiation-hydrodynamic simulations of an isolated dwarf galaxy embedded in a $10^{10}\,M_\odot$ halo, we show that the momentum transferred from resonantly scattered Lyman-$\alpha$ (Lya) photons is an important source of stellar feedback which can shape the evolution of galaxies. We find that Lya feedback suppresses star formation by a factor of two in metal-poor galaxies by regulating the dynamics of star-forming clouds before the onset of supernova explosions (SNe). This is possible because each Lya photon resonantly scatters and imparts 10-300 times greater momentum than in the single scattering limit. Consequently, the number of star clusters predicted in the simulations is reduced by a factor of $\sim 5$, compared to the model without the early feedback. More importantly, we find that galactic outflows become weaker in the presence of strong Lya radiation feedback, as star formation and associated SNe become less bursty. We also examine a model in which radiation field is arbitrarily enhanced by a factor of up to 10, and reach the same conclusion. The typical mass loading factors in our metal-poor dwarf system are estimated to be $\sim5-10$ near the mid plane, while it is reduced to $\sim1$ at larger radii. Finally, we find that the escape of ionizing radiation and hence the reionization history of the Universe is unlikely to be strongly affected by Lya feedback.

We report the discovery of a 10^4 kpc^2 gaseous structure detected in [OII] in an over-dense region of the COSMOS-Gr30 galaxy group at z~0.725 thanks to deep MUSE Guaranteed Time Observations. We estimate the total amount of diffuse ionised gas to be of the order of (~5+-3)x10^10 Msun and explore its physical properties to understand its origin and the source(s) of the ionisation. The MUSE data allow the identification of a dozen of group members embedded in this structure from emission and absorption lines. We extracted spectra from small apertures defined for both the diffuse ionised gas and the galaxies. We investigated the kinematics and ionisation properties of the various galaxies and extended gas regions thanks to line diagnostics (R23, O32 and [OIII]/H\beta) available within the MUSE wavelength range. We compared these diagnostics to photo-ionisation models and shock models. The structure is divided in two kinematically distinct sub-structures. The most extended sub-structure of ionised gas is likely rotating around a massive galaxy and displays filamentary patterns linking some galaxies. The second sub-structure links another massive galaxy hosting an Active Galactic Nucleus to a low mass galaxy but also extends orthogonally to the AGN host disk over ~35 kpc. This extent is likely ionised by the AGN itself. The location of small diffuse regions in the R23 vs. O32 diagram is compatible with photo-ionisation. However, the location of three of these regions in this diagram (low O32, high R23) can also be explained by shocks, which is supported by their large velocity dispersions. One edge-on galaxy shares the same properties and may be a source of shocks. Whatever the hypothesis, the extended gas seems to be non primordial. We favour a scenario where the gas has been extracted from galaxies by tidal forces and AGN triggered by interactions between at least the two sub-structures.

We provide, for the first time, robust observational constraints on the galaxy major merger fraction up to $z\approx 6$ using spectroscopic close pair counts. Deep Multi Unit Spectroscopic Explorer (MUSE) observations in the Hubble Ultra Deep Field (HUDF) and Hubble Deep Field South (HDF-S) are used to identify 113 secure close pairs of galaxies among a parent sample of 1801 galaxies spread over a large redshift range ($0.2<z<6$) and stellar masses ($10^7-10^{11} M_\odot$), thus probing about 12 Gyr of galaxy evolution. Stellar masses are estimated from spectral energy distribution (SED) fitting over the extensive UV-to-NIR HST photometry available in these deep Hubble fields, adding Spitzer IRAC bands to better constrain masses for high-redshift ($z\geqslant 3$) galaxies. These stellar masses are used to isolate a sample of 54 major close pairs with a galaxy mass ratio limit of 1:6. Among this sample, 23 pairs are identified at high redshift ($z\geqslant 3$) through their Ly$\alpha$ emission. The sample of major close pairs is divided into five redshift intervals in order to probe the evolution of the merger fraction with cosmic time. Our estimates are in very good agreement with previous close pair counts with a constant increase of the merger fraction up to $z\approx 3$ where it reaches a maximum of 20%. At higher redshift, we show that the fraction slowly decreases down to about 10% at $z\approx6$. The sample is further divided into two ranges of stellar masses using either a constant separation limit of $10^{9.5} M_\odot$ or the median value of stellar mass computed in each redshift bin. Overall, the major close pair fraction for low-mass and massive galaxies follows the same trend. These new, homogeneous, and robust estimates of the major merger fraction since $z\approx6$ are in good agreement with recent predictions of cosmological numerical simulations.

We present spatially resolved stellar kinematic maps, for the first time, for a sample of 17 intermediate redshift galaxies (0.2 < z < 0.8). We used deep MUSE/VLT integral field spectroscopic observations in the Hubble Deep Field South (HDFS) and Hubble Ultra Deep Field (HUDF), resulting from ~30h integration time per field, each covering 1'x1' field of view, with ~0.65" spatial resolution. We selected all galaxies brighter than 25mag in the I band and for which the stellar continuum is detected over an area that is at least two times larger than the spatial resolution. The resulting sample contains mostly late-type disk, main-sequence star-forming galaxies with 10^8.5 - 10^10.5 Msun. Using a full-spectrum fitting technique, we derive two-dimensional maps of the stellar and gas kinematics, including the radial velocity V and velocity dispersion sigma. We find that most galaxies in the sample are consistent with having rotating stellar disks with roughly constant velocity dispersions and that the Vrms=sqrt{V^2+sigma^2} of the gas and stars, a scaling proxy for the galaxy gravitational potential, compare well to each other. These spatially resolved observations of intermediate redshift galaxies suggest that the regular stellar kinematics of disk galaxies that is observed in the local Universe was already in place 4 - 7 Gyr ago and that their gas kinematics traces the gravitational potential of the galaxy, thus is not dominated by shocks and turbulent motions. Finally, we build dynamical axisymmetric Jeans models constrained by the derived stellar kinematics for two specific galaxies and derive their dynamical masses. These are in good agreement (within 25%) with those derived from simple exponential disk models based on the gas kinematics. The obtained mass-to-light ratios hint towards dark matter dominated systems within a few effective radii.

We have conducted a two-layered spectroscopic survey (1'x1' ultra deep and 3'x3' deep regions) in the Hubble Ultra Deep Field (HUDF) with the Multi Unit Spectroscopic Explorer (MUSE). The combination of a large field of view, high sensitivity, and wide wavelength coverage provides an order of magnitude improvement in spectroscopically confirmed redshifts in the HUDF; i.e., 1206 secure spectroscopic redshifts for HST continuum selected objects, which corresponds to 15% of the total (7904). The redshift distribution extends well beyond z>3 and to HST/F775W magnitudes as faint as ~30 mag (AB, 1-sigma). In addition, 132 secure redshifts were obtained for sources with no HST counterparts that were discovered in the MUSE data cubes by a blind search for emission-line features. In total, we present 1338 high quality redshifts, which is a factor of eight increase compared with the previously known spectroscopic redshifts in the same field. We assessed redshifts mainly with the spectral features [OII] at z<1.5 (473 objects) and Lya at 2.9<z<6.7 (692 objects). With respect to F775W magnitude, a 50% completeness is reached at 26.5 mag for ultra deep and 25.5 mag for deep fields, and the completeness remains >~20% up to 28-29 mag and ~27 mag, respectively. We used the determined redshifts to test continuum color selection (dropout) diagrams of high-z galaxies. The selection condition for F336W dropouts successfully captures ~80% of the targeted z~2.7 galaxies. However, for higher redshift selections (F435W, F606W, and F775W dropouts), the success rates decrease to ~20-40%. We empirically redefine the selection boundaries to make an attempt to improve them to ~60%. The revised boundaries allow bluer colors that capture Lya emitters with high Lya equivalent widths falling in the broadbands used for the color-color selection. Along with this paper, we release the redshift and line flux catalog.

We analyse the properties of circumgalactic gas around simulated galaxies in the redshift range z >= 3, utilising a new sample of cosmological zoom simulations. These simulations are intended to be representative of the observed samples of Lyman-alpha emitters recently obtained with the MUSE instrument (halo masses ~10^10-10^11 solar masses). We show that supernova feedback has a significant impact on both the inflowing and outflowing circumgalactic medium by driving outflows, reducing diffuse inflow rates, and by increasing the neutral fraction of inflowing gas. By temporally stacking simulation outputs we find that significant net mass exchange occurs between inflowing and outflowing phases: none of the phases are mass-conserving. In particular, we find that the mass in neutral outflowing hydrogen declines exponentially with radius as gas flows outwards from the halo centre. This is likely caused by a combination of both fountain-like cycling processes and gradual photo/collisional ionization of outflowing gas. Our simulations do not predict the presence of fast-moving neutral outflows in the CGM. Neutral outflows instead move with modest radial velocities (~ 50 kms^-1), and the majority of the kinetic energy is associated with tangential rather than radial motion.

We study the oxygen abundance profiles of the gas-phase components in hydrodynamical simulations of pre-prepared disc galaxies including major mergers, close encounters and isolated configurations. We analyse the evolution of the slope of oxygen abundance profiles and the specific star formation rate (sSFR) along their evolution. We find that galaxy-galaxy interactions could generate either positive and negative gas-phase oxygen profiles depending on the state of evolution. Along the interaction, galaxies are found to have metallicity gradients and sSFR consistent with observations, on average. Strong gas inflows produced during galaxy-galaxy interactions or as a result of strong local instabilities in gas-rich discs are able to produce both a quick dilution of the central gas-phase metallicity and a sudden increase of the sSFR. Our simulations show that, during these events, a correlation between the metallicity gradients and the sSFR can be set up if strong gas inflows are triggered in the central regions in short timescales. Simulated galaxies without experiencing strong disturbances evolve smoothly without modifying the metallicity gradients. Gas-rich systems show large dispersion along the correlation. The dispersion in the observed relation could be interpreted as produced by the combination of galaxies with different gas-richness and/or experiencing different types of interactions. Hence, our findings suggest that the observed relation might be the smoking gun of galaxies forming in a hierarchical clustering scenario.

(Abridged) We make use of the deepest VLT/MUSE observations performed so far on the Hubble Deep Field South (HDFS) to characterize the low-mass (< $10^{10}$M$_\odot$) galaxy population at intermediate redshift. We identify a sample of 28 spatially-resolved emission-line galaxies in the deep (27h integration time) MUSE data cube, spread over a redshift interval of 0.2 < z < 1.4. The public HST images and multi-band photometry over the HDFS are used to constrain the stellar mass and star formation rate (SFR) of the galaxies and to perform a morphological analysis. We derive the resolved ionized gas properties of these galaxies from the MUSE data and model the disk (both in 2D and with GalPaK$^{\rm 3D}$) to retrieve their intrinsic gas kinematics. We build a sample of resolved emission-line galaxies of much lower stellar mass and SFR (by $\sim$1-2 orders of magnitude) than previous 3D spectroscopic surveys. Most of the spatially-resolved MUSE-HDFS galaxies have gas kinematics consistent with disk-like rotation, but about 20% have velocity dispersions larger than the rotation velocities, and 30% are part of a close pair and/or show clear signs of recent gravitational interactions. In the high-mass regime, the MUSE-HDFS galaxies follow the Tully-Fisher relation defined from previous surveys in a similar redshift range. This scaling relation extends also to lower masses/velocities but with a higher dispersion. The MUSE-HDFS galaxies follow the scaling relations defined in the local universe between the specific angular momentum and the stellar mass. However, we find that intermediate-redshift star-forming galaxies fill a continuum transition from the spiral to elliptical local scaling relations, according to the dynamical state of the gas. This indicates that some galaxies may lose their angular momentum and become dispersion-dominated prior to becoming passive.

We present a method to constrain galaxy parameters directly from three-dimensional data cubes. The algorithm compares directly the data with a parametric model mapped in $x,y,\lambda$ coordinates. It uses the spectral lines-spread function (LSF) and the spatial point-spread function (PSF) to generate a three-dimensional kernel whose characteristics are instrument specific or user generated. The algorithm returns the intrinsic modeled properties along with both an `intrinsic' model data cube and the modeled galaxy convolved with the 3D-kernel. The algorithm uses a Markov Chain Monte Carlo (MCMC) approach with a nontraditional proposal distribution in order to efficiently probe the parameter space. We demonstrate the robustness of the algorithm using 1728 mock galaxies and galaxies generated from hydrodynamical simulations in various seeing conditions from 0.6" to 1.2". We find that the algorithm can recover the morphological parameters (inclination, position angle) to within 10% and the kinematic parameters (maximum rotation velocity) to within 20%, irrespectively of the PSF in seeing (up to 1.2") provided that the maximum signal-to-noise ratio (SNR) is greater than $\sim3$ pixel$^{-1}$ and that the ratio of the galaxy half-light radius to seeing radius is greater than about 1.5. One can use such an algorithm to constrain simultaneously the kinematics and morphological parameters of (nonmerging) galaxies observed in nonoptimal seeing conditions. The algorithm can also be used on adaptive-optics (AO) data or on high-quality, high-SNR data to look for nonaxisymmetric structures in the residuals.

The relation between the stellar mass and the star formation rate characterizes how the instantaneous star formation is determined by the galaxy past star formation history and by the growth of the dark matter structures. We deconstruct the M-SFR plane by measuring the specific SFR functions in several stellar mass bins from z=0.2 out to z=1.4. Our analysis is primary based on a MIPS 24$\mu m$ selected catalogue combining the COSMOS and GOODS surveys. We estimate the SFR by combining mid- and far-infrared data for 20500 galaxies. The sSFR functions are derived in four stellar mass bins within the range 9.5<log(M/Msun)<11.5. First, we demonstrate the importance of taking into account selection effects when studying the M-SFR relation. Secondly, we find a mass-dependent evolution of the median sSFR with redshift varying as $sSFR \propto (1+z)^{b}$, with $b$ increasing from $b=2.88$ to $b=3.78$ between $M=10^{9.75}Msun$ and $M=10^{11.1}Msun$, respectively. At low masses, this evolution is consistent with the cosmological accretion rate and predictions from semi-analytical models (SAM). This agreement breaks down for more massive galaxies showing the need for a more comprehensive description of the star-formation history in massive galaxies. Third, we obtain that the shape of the sSFR function is invariant with time at z<1.4 but depends on the mass. We observe a broadening of the sSFR function ranging from 0.28 dex at $M=10^{9.75}Msun$ to 0.46 dex at $M=10^{11.1}Msun$. Such increase in the scatter of the M-SFR relation suggests an increasing diversity of SFHs as the stellar mass increases. Finally, we find a gradual decline of the sSFR with mass as $log(sSFR) \propto -0.17M$. We discuss the numerous physical processes, as gas exhaustion in hot gas halos or secular evolution, which can gradually reduce the sSFR and increase the SFH diversity.

We explore the properties of high-redshift Lyman-alpha emitters (LAE), and their link with the Lyman-Break galaxy population (LBG), using a semi-analytic model of galaxy formation that takes into account resonant scattering of Lya photons in gas outflows. We can reasonably reproduce the abundances of LAEs and LBGs from redshift 3 to 7, as well as most UV LFs of LAEs. The stronger dust attenuation for (resonant) Lya photons compared to UV continuum photons in bright LBGs provides a natural interpretation to the increase of the LAE fraction in LBG samples, X_LAE, towards fainter magnitudes. The redshift evolution of X_LAE seems however very sensitive to UV magnitudes limits and EW cuts. In spite of the apparent good match between the statistical properties predicted by the model and the observations, we find that the tail of the Lya equivalent width distribution (EW > 100 A) cannot be explained by our model, and we need to invoke additional mechanisms. We find that LAEs and LBGs span a very similar dynamical range, but bright LAEs are about 4 times rarer than LBGs in massive halos. Moreover, massive halos mainly contain weak LAEs in our model, which might introduce a bias towards low-mass halos in surveys which select sources with high EW cuts. Overall, our results are consistent with the idea that LAEs and LBGs make a very similar galaxy population. Their apparent differences seem mainly due to EW selections, UV detection limits, and a decreasing Lya-to-UV escape fraction ratio in high SFR galaxies.

We observed the Hubble Deep Field South with the new panoramic integral field spectrograph MUSE that we built and just commissioned at the VLT. The data cube resulting from 27 hours of integration covers one arcmin^2 field of view at an unprecedented depth with a 1 sigma emission line surface brightness limit of 1x$10^{-19}$ erg/s/cm$^2$/arcsec$^2$ and contains ~90,000 spectra. We present the combined and calibrated data cube, and we perform a first-pass analysis of the sources detected in the HDF-S imaging. We measured the redshifts of 189 sources up to a magnitude F814W = 29.5, increasing by more than an order of magnitude the number of known spectroscopic redshifts in this field. We also discovered 26 Lya emitting galaxies which are not detected in the HST WFPC2 deep broad band images. The intermediate spectral resolution of 2.3{\AA} allows us to separate resolved asymmetric Lya emitters, [O II] emitters, and C III] emitters and the large instantaneous wavelength range of 4500{\AA} helps to identify single emission lines. We also show how the three dimensional information of MUSE helps to resolve sources which are confused at ground-based image quality. Overall, secure identifications are provided for 83% of the 227 emission line sources detected in the MUSE data cube and for 32% of the 586 sources identified in the HST catalog of Casertano et al 2000. The overall redshift distribution is fairly flat to z=6.3, with a reduction between z=1.5 to 2.9, in the well-known redshift desert. The field of view of MUSE also allowed us to detect 17 groups within the field. We checked that the number counts of [O II] and Ly-a emitters are roughly consistent with predictions from the literature. Using two examples we demonstrate that MUSE is able to provide exquisite spatially resolved spectroscopic information on intermediate redshift galaxies present in the field.

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.

We trace the formation and advection of several elements within a cosmological adaptive mesh refinement simulation of an L* galaxy. We use nine realisations of the same initial conditions with different stellar Initial Mass Functions (IMFs), mass limits for type-II and type-Ia supernovae (SNII, SNIa) and stellar lifetimes to constrain these sub-grid phenomena. Our code includes self-gravity, hydrodynamics, star formation, radiative cooling and feedback from multiple sources within a cosmological framework. Under our assumptions of nucleosynthesis we find that SNII with progenitor masses of up to 100 Msun are required to match low metallicity gas oxygen abundances. Tardy SNIa are necessary to reproduce the classical chemical evolution knee in [O/Fe]-[Fe/H]: more prompt SNIa delayed time distributions do not reproduce this feature. Within our framework of hydrodynamical mixing of metals and galaxy mergers we find that chemical evolution is sensitive to the shape of the IMF and that there exists a degeneracy with the mass range of SNII. We look at the abundance plane and present the properties of different regions of the plot, noting the distinct chemical properties of satellites and a series of nested discs that have greater velocity dispersions, are more alpha-rich and metal poor with age.

We present a set of hydrodynamical/Nbody controlled simulations of isolated gas rich galaxies that self-consistently include SN feedback and a detailed chemical evolution model, both tested in cosmological simulations. The initial conditions are motivated by the observed star forming galaxies at z ~ 2-3. We find that the presence of a multiphase interstellar media in our models promotes the growth of disc instability favouring the formation of clumps which in general, are not easily disrupted on timescales compared to the migration time. We show that stellar clumps migrate towards the central region and contribute to form a classical-like bulge with a Sersic index, n > 2. Our physically-motivated Supernova feedback has a mild influence on clump survival and evolution, partially limiting the mass growth of clumps as the energy released per Supernova event is increased, with the consequent flattening of the bulge profile. This regulation does not prevent the building of a classical-like bulge even for the most energetic feedback tested. Our Supernova feedback model is able to establish a self-regulated star formation, producing mass-loaded outflows and stellar age spreads comparable to observations. We find that the bulge formation by clumps may coexit with other channels of bulge assembly such as bar and mergers. Our results suggest that galactic bulges could be interpreted as composite systems with structural components and stellar populations storing archaeological information of the dynamical history of their galaxy.

We present an investigation into the origins of a series of interlaced narrow filamentary stellar structures, loops and plumes in the vicinity of the Virgo Cluster, edge-on spiral galaxy, NGC 4216 that were previously identified by the Blackbird Telescope. Using the deeper, higher-resolution and precisely calibrated optical CFHT/MegaCam images obtained as part of the Next Generation Virgo Cluster Survey (NGVS), we confirm the previously identified features and identify a few additional structures. The NGVS data allowed us to make a physical study of these low-surface brightness features and investigate their origin. The likely progenitors of the structures were identified as either already catalogued VCC dwarfs or newly discovered satellites caught in the act of being destroyed. They have the same g-i color index and likely contain similar stellar populations. The alignment of three dwarfs along an apparently single stream is intriguing, and we cannot totally exclude that these are second-generation dwarf galaxies being born inside the filament from the debris of an original dwarf. The observed complex structures, including in particular a stream apparently emanating from a satellite of a satellite, point to a high rate of ongoing dwarf destruction/accretion in the region of the Virgo Cluster where NGC 4216 is located. We discuss the age of the interactions and whether they occurred in a group that is just falling into the cluster and shows signs of so-called "pre-processing" before it gets affected by the cluster environment, or in a group which already ventured towards the central regions of Virgo Cluster.

We test the validity of comparing simulated field disk galaxies with the empirical properties of systems situated within environments more comparable to loose groups, including the Milky Way's Local Group. Cosmological simulations of Milky Way-mass galaxies have been realised in two different environment samples: in the field and in environments with similar properties to the Local Group. Apart from the environments of the galaxies, the samples are kept as homogeneous as possible with equivalent ranges in last major merger time, halo mass and halo spin. Comparison of these two samples allow for systematic differences in the simulations to be identified. Metallicity gradients, disk scale lengths, colours, magnitudes and age-velocity dispersion relations are studied for each galaxy in the suite and the strength of the link between these and environment of the galaxies is studied. The bulge-to-disk ratio of the galaxies show that these galaxies are less spheroid dominated than many other simulated galaxies in literature with the majority of both samples being disk dominated. We find that secular evolution and mergers dominate the spread of morphologies and metallicity gradients with no visible differences between the two environment samples. In contrast with this consistency in the two samples there is tentative evidence for a systematic difference in the velocity dispersion-age relations of galaxies in the different environments. Loose group galaxies appear to have more discrete steps in their velocity dispersion-age relations. We conclude that at the current resolution of cosmological galaxy simulations field environment galaxies are sufficiently similar to those in loose groups to be acceptable proxies for comparison with the Milky Way provided that a similar assembly history is considered.

We report the discovery of a giant HI tail in the intra-group medium of HCG 44 as part of the Atlas3D survey. The tail is ~300 kpc long in projection and contains ~5x10^8 M_sun of HI. We detect no diffuse stellar light at the location of the tail down to ~28.5 mag/arcsec^2 in g band. We speculate that the tail might have formed as gas was stripped from the outer regions of NGC 3187 (a member of HCG 44) by the group tidal field. In this case, a simple model indicates that about 1/3 of the galaxy's HI was stripped during a time interval of <1 Gyr. Alternatively, the tail may be the remnant of an interaction between HCG 44 and NGC 3162, a spiral galaxy now ~650 kpc away from the group. Regardless of the precise formation mechanism, the detected HI tail shows for the first time direct evidence of gas stripping in HCG 44. It also highlights that deep HI observations over a large field are needed to gather a complete census of this kind of events in the local Universe.

By means of high-resolution cosmological hydrodynamical simulations of Milky Way-like disc galaxies, we conduct an analysis of the associated stellar metallicity distribution functions (MDFs). After undertaking a kinematic decomposition of each simulation into spheroid and disc sub-components, we compare the predicted MDFs to those observed in the solar neighbourhood and the Galactic bulge. The effects of the star formation density threshold are visible in the star formation histories, which show a modulation in their behaviour driven by the threshold. The derived MDFs show median metallicities lower by 0.2-0.3 dex than the MDF observed locally in the disc and in the Galactic bulge. Possible reasons for this apparent discrepancy include the use of low stellar yields and/or centrally-concentrated star formation. The dispersions are larger than the one of the observed MDF; this could be due to simulated discs being kinematically hotter relative to the Milky Way. The fraction of low metallicity stars is largely overestimated, visible from the more negatively skewed MDF with respect to the observational sample. For our fiducial Milky Way analog, we study the metallicity distribution of the stars born "in situ" relative to those formed via accretion (from disrupted satellites), and demonstrate that this low-metallicity tail to the MDF is populated primarily by accreted stars. Enhanced supernova and stellar radiation energy feedback to the surrounding interstellar media of these pre-disrupted satellites is suggested as an important regulator of the MDF skewness.

We investigate the accuracy in the recovery of the stellar dynamics of barred galaxies when using axisymmetric dynamical models. We do this by trying to recover the mass-to-light ratio (M/L) and the anisotropy of realistic galaxy simulations using the Jeans Anisotropic Multi-Gaussian Expansion (JAM) method. However, given that the biases we find are mostly due to an application of an axisymmetric modeling algorithm to a non-axisymmetric system and in particular to inaccuracies in the de-projected mass model, our results are relevant for general axisymmetric modelling methods. We run N-body collisionless simulations to build a library with various luminosity distribution, constructed to mimic real individual galaxies, with realistic anisotropy. The final result of our evolved library of simulations contains both barred and unbarred galaxies. The JAM method assumes an axisymmetric mass distribution, and we adopt a spatially constant M/L and anisotropy beta_z=1-sigma_z^2/sigma_R^2 distributions. The models are fitted to two-dimensional maps of the second velocity moments V_rms=sqrt(V^2+sigma^2) of the simulations for various viewing angles (position angle of the bar and inclination of the galaxy). We find that the inclination is generally well recovered by the JAM models, for both barred and unbarred simulations. For unbarred simulations the M/L is also accurately recovered, with negligible median bias and with a maximum one of just Delta(M/L)<1.5% when the galaxy is not too close to face on. At very low inclinations (i<30 deg) the M/L can be significantly overestimated (9% in our tests, but errors can be larger for very face-on views). For barred simulations the M/L is on average (when PA=45 deg) essentially unbiased, but we measure an over/under estimation of up to Delta(M/L)=15% in our tests. The sign of the M/L bias depends on the position angle of the bar as expected. [Abridged]

We examine the chemical properties of 5 cosmological hydrodynamical simulations of an M33-like disc galaxy which have been shown to be consistent with the morphological characteristics and bulk scaling relations expected of late-type spirals. These simulations are part of the Making Galaxies In a Cosmological Context (MaGICC) Project, in which stellar feedback is tuned to match the stellar mass -- halo mass relationship. Each realisation employed identical initial conditions and assembly histories, but differed from one another in their underlying baryonic physics prescriptions, including (a) the efficiency with which each supernova energy couples to the ISM, (b) the impact of feedback associated with massive star radiation pressure, (c) the role of the minimum shut-off time for radiative cooling of Type II SNe remnants, (d) the treatment of metal diffusion, and (e) varying the IMF. Our analysis focusses on the resulting stellar metallicity distribution functions (MDFs) in each simulated (analogous) `solar neighbourhood' and central `bulge' region. We compare the simulated MDFs' skewness, kurtosis, and dispersion (inter-quartile, inter-decile, inter-centile, and inter-tenth-percentile regions) with that of the empirical solar neighbourhood MDF and Local Group dwarfs. We find that the MDFs of the simulated discs are more negatively skewed, with higher kurtosis, than those observed locally. We can trace this difference to the simulations' tight and correlated age-metallicity relations (compared with that of the Milky Way), suggesting that these relations within `dwarf' discs might be steeper than in L* discs and/or the degree of stellar orbital re-distribution and migration inferred locally has not been captured in their entirety, at the resolution of our simulations. The important role of metal diffusion in ameliorating the over-production of extremely metal-poor stars is highlighted.

We examine radial and vertical metallicity gradients using a suite of disk galaxy simulations, supplemented with two classic chemical evolution approaches. We determine the rate of change of gradient and reconcile differences between extant models and observations within the `inside-out' disk growth paradigm. A sample of 25 disks is used, consisting of 19 from our RaDES (Ramses Disk Environment Study) sample, realised with the adaptive mesh refinement code RAMSES. Four disks are selected from the MUGS (McMaster Unbiased Galaxy Simulations) sample, generated with the smoothed particle hydrodynamics (SPH) code GASOLINE, alongside disks from Rahimi et al. (GCD+) and Kobayashi & Nakasato (GRAPE-SPH). Two chemical evolution models of inside-out disk growth were employed to contrast the temporal evolution of their radial gradients with those of the simulations. We find that systematic differences exist between the predicted evolution of radial abundance gradients in the RaDES and chemical evolution models, compared with the MUGS sample; specifically, the MUGS simulations are systematically steeper at high-redshift, and present much more rapid evolution in their gradients. We find that the majority of the models predict radial gradients today which are consistent with those observed in late-type disks, but they evolve to this self-similarity in different fashions, despite each adhering to classical `inside-out' growth. We find that radial dependence of the efficiency with which stars form as a function of time drives the differences seen in the gradients; systematic differences in the sub-grid physics between the various codes are responsible for setting these gradients. Recent, albeit limited, data at redshift z=1.5 are consistent with the steeper gradients seen in our SPH sample, suggesting a modest revision of the classical chemical evolution models may be required.

We analyse a set of galaxy interactions performed by using a self-consistent chemo-hydrodynamical model which includes star formation, Supernova feedback and chemical evolution. In agreement with previous works, we find that tidally-induced low-metallicity gas inflows dilute the central oxygen abundance and contribute to the flattening of the metallicity gradients. The tidally-induced inflows trigger starbursts which increase the impact of SN II feedback injecting new chemical elements and driving galactic winds which modulate the metallicity distribution. Although $\alpha$-enhancement in the central regions is detected as a result of the induced starbursts in agreement with previous works, our simulations suggest that this parameter can only provide a timing of the first pericentre mainly for non-retrograde encounters. In order to reproduce wet major mergers at low and high redshifts, we have run simulations with respectively 20 and 50 percent of the disc in form of gas. We find that the more gas-rich encounters behave similarly to the less rich ones, between the first and second pericentre where low-metallicity gas inflows are triggered. However, the higher strength of the inflows triggered in the more gas-rich interactions produces larger metal dilutions factors which are afterward modulated by the new chemical production by Supernova. We find that the more gas-rich interaction develops violent and clumpy star formation triggered by local instabilities all over the disc before the first pericentre, so that if these galaxies were observed at these early stages where no important tidally-induced inflows have been able to develop yet, they would tend to show an excess of oxygen. We find a global mean correlation of both the central abundances and the gradients with the strength of the star formation activity. [abridged]

The mass assembly of galaxies leaves imprints in their outskirts, such as shells and tidal tails. The frequency and properties of such fine structures depend on the main acting mechanisms - secular evolution, minor or major mergers - and on the age of the last substantial accretion event. We use this to constrain the mass assembly history of two apparently relaxed nearby Early-Type Galaxies (ETGs) selected from the Atlas-3D sample, NGC 680 and NGC 5557. Our ultra deep optical images obtained with MegaCam on the Canada-France-Hawaii Telescope reach 29 mag/arcsec^2 in the g-band. They reveal very low-surface brightness (LSB) filamentary structures around these ellipticals. Among them, a gigantic 160 kpc long tail East of NGC 5557 hosts gas-rich star-forming objects. NGC 680 exhibits two major diffuse plumes apparently connected to extended HI tails, as well as a series of arcs and shells. Comparing the outer stellar and gaseous morphology of the two ellipticals with that predicted from models of colliding galaxies, we argue that the LSB features are tidal debris, the star-forming objects near NGC 5557, long lived Tidal Dwarf Galaxies and that each of the two ETGs was assembled during a relatively recent, major wet merger, which likely occurred at a redshift below z = 0.5. The inner kinematics of NGC 680 is typical for fast rotators which make the bulk of nearby ETGs in the Atlas-3D sample. NGC 5557 belongs to the poorly populated class of massive, round, slow rotators that are predicted by semi-analytic models and cosmological simulations to be the end-product of a complex mass accretion history, involving ancient major mergers and more recent minor mergers. Our observations suggest that under specific circumstances a single binary merger may dominate the formation history of such objects and thus that at least some massive ETGs may form at relatively low redshift (abridged).

The Monoceros ring is a collection of stars in nearly-circular orbits at roughly 18 kpc from the Galactic center. It may have originated (i) as the response of the disc to perturbations excited by satellite companions or (ii) from the tidal debris of a disrupted dwarf galaxy. The metallicity of Monoceros stars differs from that of disc stars at comparable Galactocentric distances, an observation that disfavours the first scenario. On the other hand, circular orbits are difficult to accommodate in the tidal-disruption scenario, since it requires a satellite which at the time of disruption was itself in a nearly circular orbit. Such satellite could not have formed at the location of the ring and, given its low mass, dynamical friction is unlikely to have played a major role in its orbital evolution. We search cosmological simulations for low-mass satellites in nearly-circular orbits and find that they result, almost invariably, from orbital changes induced by collisions with more massive satellites: the radius of the circular orbit thus traces the galactocentric distance of the collision. Interestingly, the Sagittarius dwarf, one of the most luminous satellites of the Milky Way, is in a polar orbit that crosses the Galactic plane at roughly the same Galactocentric distance as Monoceros. We use idealized simulations to demonstrate that an encounter with Sagittarius might well have led to the circularization and subsequent tidal demise of the progenitor of the Monoceros ring.

The mass assembly of galaxies leaves various imprints on their surroundings, such as shells, streams and tidal tails. The frequency and properties of these fine structures depend on the mechanism driving the mass assembly: e.g. a monolithic collapse, rapid cold-gas accretion followed by violent disk instabilities, minor mergers or major dry / wet mergers. Therefore, by studying the outskirts of galaxies, one can learn about their main formation mechanism. I present here our on-going work to characterize the outskirts of Early-Type Galaxies (ETGs), which are powerful probes at low redshift of the hierarchical mass assembly of galaxies. This work relies on ultra-deep optical images obtained at CFHT with the wide-field of view MegaCam camera of field and cluster ETGs obtained as part of the Atlas-3D and NGVS projects. State of the art numerical simulations are used to interpret the data. The images reveal a wealth of unknown faint structures at levels as faint as 29 mag arcsec-2 in the g-band. Initial results for two galaxies are presented here.

Extended HI structures around galaxies are of prime importance to probe galaxy formation scenarios. The giant HI ring in the Leo group is one of the largest and most intriguing HI structures in the nearby Universe. Whether it consists of primordial gas, as suggested by the apparent absence of any optical counterpart and the absence of an obvious physical connection to nearby galaxies, or of gas expelled from a galaxy in a collision is actively debated. We present deep wide field-of-view optical images of the ring region obtained with MegaCam on the CFHT. They reveal optical counterparts to several HI and UV condensations along the ring, in the g', r', and i' bands, which likely correspond to stellar associations formed within the gaseous ring. Analyzing the spectral energy distribution of one of these star-forming regions, we found it to be typical for a star-forming region in pre-enriched tidal debris. We then use simulations to test the hypothesis that the Leo ring results from a head-on collision between Leo group members NGC 3384 and M96. According to our model which is able to explain, at least qualitatively, the main observational properties of the system, the Leo ring is consistent with being a collisional ring. It is thus likely another example of extended intergalactic gas made-up of pre-enriched collisional debris.

We focus on an analysis of galaxies of high stellar mass and low metallicity. We cross-correlated the Millenium Galaxy Catalogue (MGC) and the Sloan Digital Sky Survey (SDSS) galaxy catalogue to provide a sample of MGC objects with high resolution imaging and both spectroscopic and photometric information available in the SDSS database. For each galaxy in our sample, we conducted a systematic morphological analysis by visual inspection of MGC images using their luminosity contours. The galaxies are classified as either disturbed or undisturbed objects. We divide the sample into three metallicity regions, within wich we compare the properties of disturbed and undisturbed objects. We find that the fraction of galaxies that are strongly disturbed, indicative of being merger remnants, is higher when lower metallicity objects are considered. The three bins analysed consist of approximatively 15%, 20%, and 50% disturbed galaxies (for high, medium, and low metallicity, respectively). Moreover, the ratio of the disturbed to undisturbed relative distributions of the population age indicator, Dn(4000), in the low metallicity bin, indicates that the disturbed objects have substantially younger stellar populations than their undisturbed counterparts. In addition, we find that an analysis of colour distributions provides similar results, showing that low metallicity galaxies with a disturbed morphology are bluer than those that are undisturbed. The bluer colours and younger populations of the low metallicity, morphologically disturbed objects suggest that they have experienced a recent merger with an associated enhanced star formation rate. [abridged]

We investigate whether the formation mechanism of boxy and peanut-shaped (B/PS) bulges could depend on the gas content of the galaxy. We have performed N-body simulations with and without a gaseous component. In the second case star formation/feedback recipes have also been implemented to create new stellar populations. As in many previous studies, in our N-body collisionless simulation, the B/PS is due to the classical break in the z mirror symmetry lasting roughly 200 Myr. When a gaseous component and star formation recipes are added to the simulation, the bulge-growing mechanism is quite different. The young stellar population that is born in the thin gaseous disc rapidly populates vertical resonant orbits triggered by the combined effects of the linear horizontal and vertical ILRs. This leads to a B/PS bulge mainly made of stellar material younger than the surrounding population. The non-linear analysis of the orbital structure shows that the main orbit family responsible for the B/PS is not the same in the two cases. The 2:2:1 orbits prevail in the collisionless simulation whereas additional asymmetrical families contribute to the B/PS if a dissipative component is present and can form new stars. We found that 2:3:1 and 2:5:1 orbits trap a significant fraction of the mass. A flat ringed discy stellar component also appears simultaneously with the thickening of the young population. It is due to the star formation in a nuclear gaseous disc located in the central kpc, inside the ILR, and accumulated there by the torques exerted by the large-scale bar. Remarkably, it remains flat throughout the simulation although it develops a nuclear bar, leading to a double-barred galaxy. We predict that two populations of B/PS bulges could exist and even coexist in the same galaxy.

We study the mass-metallicity relation of galaxies in pairs and in isolation taken from the SDSS-DR4 using the stellar masses and oxygen abundances derived by Tremonti et al. (2004). Close galaxy pairs, defined by projected separation r_p < 25kpc/h and radial velocity Delta_V < 350 km/s, are morphologically classified according to the strength of the interaction signs. We find that only for pairs showing signs of strong interactions, the mass-metallicity relation differs significantly from that of galaxies in isolation. In such pairs, the mean gas-phase oxygen abundances of galaxies with low stellar masses (Mstar ~< 10^9 Msun/h) exhibit an excess of 0.2 dex. Conversely, at larger masses (Mstar >~ 10^10 Msun/h) galaxies have a systematically lower metallicity, although with a smaller difference (-0.05 dex). Similar trends are obtained if g-band magnitudes are used instead of stellar masses. In minor interactions, we find that the less massive member is systematically enriched, while a galaxy in interaction with a comparable stellar mass companion shows a metallicity decrement with respect to galaxies in isolation. We argue that metal-rich starbursts triggered by a more massive component, and inflows of low metallicity gas induced by comparable or less massive companion galaxies, provide a natural scenario to explain our findings.

We report about the fact that the stellar population that is born in the gas inflowing towards the central regions can be vertically unstable leading to a B/PS feature remarkably bluer that the surrounding bulge. Using new chemodynamical simulations we show that this young population does not remain as flat as the gaseous nuclear disc and buckles out of the plane to form a new boxy bulge. We show that such a young B/PS bulge can be detected in colour maps.

We have investigated the accuracy and reliability of six methods used to determine the length of stellar bars in galaxies or N-body simulations. All these methods use ellipse fitting and Fourier decomposition of the surface brightness. We have applied them to N-body simulations that include stars, gas, star formation, and feedback. Stellar particles were photometrically calibrated to make B and K-band mock images. Dust absorption is also included. We discuss the advantages and drawbacks of each method, the effects of projection and resolution, as well as the uncertainties introduced by the presence of dust. The use of N-body simulations allows us to compare the location of Ultra Harmonic Resonance (UHR or 4/1) and corotation (CR) with measured bar lengths. We show that the minimum of ellipticity located just outside the bulk of the bar is correlated with the corotation, whereas the location of the UHR can be approximated using the phase of the fitted ellipses or the phase of the m=2 Fourier development of the surface brightness. We give evidence that the classification of slow/fast bars, based on the ratio R = Rcr/Rbar could increase from 1 (fast bar) to 1.4 (slow bar) just by a change of method. We thus conclude that one has to select the right bar-length estimator depending on the application, since these various estimators do not define the same physical area.

We present a detailed analysis of the evolution of a simulated isolated disc galaxy. The simulation includes stars, gas, star formation and simple chemical yields. Stellar particles are split in two populations: the old one is present at the beginning of the simulation and is calibrated according to various ages and metallicities; the new population borns in the course of the simulation and inherits the metallicity of the gas particles. The results have been calibrated in four wavebands with the spectro-photometric evolutionary model GISSEL2000 (Bruzual & Charlot 1993). Dust extinction has also been taken into account. A rest-frame morphological and bidimensional photometric analysis has been performed on simulated images, with the same tools as for observations. The effects of the stellar bar formation and the linked star formation episode on the global properties of the galaxy (mass and luminosity distribution, colours, isophotal radii) have been analysed. In particular, we have disentangled the effects of stellar evolution from dynamic evolution to explain the cause of the isophotal radii variations. We show that the dynamic properties (e.g. mass) of the area enclosed by any isophotal radius depends on the waveband and on the level of star formation activity. It is also shown that the bar isophotes remain thinner than mass isodensities a long time (> 0.7 Gyr) after the maximum of star formation rate. We show that bar ellipticity is very wavelength dependent as suggested by real observations. Effects of dust extinction on photometric and morphological measurements are systematically quantified.