Search

Feedback is widely recognized as an essential condition for Lyman continuum (LyC) escape in star-forming galaxies. However, the mechanisms by which galactic outflows clear neutral gas and dust remain unclear. In this paper, we model the Mg II 2796\r{A}, 2804\r{A} absorption + emission lines in 29 galaxies taken from the Low-z LyC Survey (LzLCS) to investigate the impact of (radiation + mechanical) feedback on LyC escape. Using constraints on Mg$^+$ and photoionization models, we map the outflows' neutral hydrogen content and predict $f_{esc}^{LyC}$ with a multiphase wind model. We measure mass, momentum, and energy loading factors for the neutral winds, which carry up to 10% of the momentum and 1% of the energy in SFR-based deposition rates. We use SED template fitting to determine the relative ages of stellar populations, allowing us to identify radiation feedback dominant systems. We then examine feedback related properties (stellar age, loading factors, etc.) under conditions that optimize feedback efficiency, specifically high star formation rate surface density and compact UV half-light radii. Our findings indicate that the strongest leakers are radiation feedback dominant, lack Mg II outflows, but have extended broad components in higher ionization lines like [O III] 5007\r{A}, as observed by Amor\'in et al. (2024). In contrast, galaxies experiencing supernovae feedback typically exhibit weaker $f_{esc}^{LyC}$ and show evidence of outflows in both Mg II and higher ionization lines. We attribute these findings to rapid or "catastrophic" cooling in the radiation-dominant systems, which, given the low metallicities in our sample, are likely experiencing delayed supernovae., Comment: 34 pages, 16 figures, 7 tables

The vast majority of star-forming galaxies are surrounded by large reservoirs of gas ejected from the interstellar medium. Ultraviolet absorption and emission lines represent powerful diagnostics to constrain the cool phase of these outflows, through resonant transitions of hydrogen and metal ions. The interpretation of these observations is often remarkably difficult as it requires detailed modelling of the propagation of the continuum and emission lines in the gas. To this aim, we present a large public grid of about 20000 simulated spectra which includes HI Lyman-alpha (Lya) and five metal transitions associated with MgII, CII, SiII, and FeII that is accessible online at https://rascas.univ-lyon1.fr/app/idealised_models_grid/. The spectra have been computed with the RASCAS radiative transfer code for 5760 idealised spherical configurations surrounding a central point source emission, and characterised by their column density, Doppler parameter, dust opacity, wind velocity, as well as various density/velocity gradients. Designed to interpret Lya and metal line profiles, our grid exhibits a wide diversity of resonant absorption and emission features, as well as fluorescent lines. We illustrate how it can help better constrain wind properties by performing a joint modelling of observed Lya, CII, and SiII spectra. Using CLOUDY simulations and virial scaling relations, we show that Lya is expected to be a faithful tracer of the gas at T=10^4-10^5 K, even if the medium is highly-ionised. While CII is found to probe the same range of temperatures as Lya, other metal lines merely trace cooler phases (T=10^4 K). As their gas opacity strongly depends on gas temperature, incident radiation field, metallicity and dust depletion, we caution that optically thin metal lines do not necessarily originate from low HI column densities and may not accurately probe Lyman continuum leakage., Comment: 18 pages, 13 figures

The circumgalactic medium (CGM) serves as a baryon reservoir that connects galaxies to the intergalactic medium and fuels star formation. The spatial distribution of the metal-enriched cool CGM has not yet been directly revealed at cosmic noon (z~2-4), as bright emission lines at these redshifts are not covered by optical integral field units. To remedy this situation, we aim for the first-ever detections and exploration of extended SiII* emission (low-ionization state, LIS), referred to as ``SiII* halos'', at redshifts ranging from z=2 to 4 as a means to trace the metal-enriched cool CGM. We use a sample of 39 galaxies with systemic redshifts of z=2.1-3.9 measured with the [CIII] doublet in the MUSE Hubble Ultra Deep Field catalog, which contains integration times spanning from ~30 to 140 hours. We search for extended SiII*1265, 1309, 1533 emission (fluorescent lines) around individual galaxies. We also stack a subsample of 14 UV-bright galaxies. We report five individual detections of SiII*1533 halos. We also confirm the presence of SiII*1533 halos in stacks for the subsample containing UV-bright sources. The other lines do not show secure detections of extended emission in either individual or stacking analyses. These detections may imply that the presence of metal-enriched CGM is a common characteristic for UV-bright galaxies. To investigate whether the origin of SiII* is continuum pumping as suggested in previous studies, we check the consistency of the equivalent width (EW) of SiII* emission and the EW of SiII absorption for the individual halo object with the most reliable detection. We confirm the equivalence, suggesting that photon conservation works for this object and pointing toward continuum pumping as the source of SiII*. We also investigate SiII* lines in a RAMSES-RT zoom-in simulation including continuum pumping and find ubiquitous presence of extended halos., Comment: 16 pages, 10 figures, 1 table, accepted for publication in A&A

Galactic outflows are believed to play a critical role in the evolution of galaxies by regulating their mass build-up and star formation. Theoretical models assumes bipolar shapes for the outflows that extends well into the circumgalctic medium (CGM), up to tens of kpc perpendicular to the galaxies. They have been directly observed in the local Universe in several individual galaxies, e.g., around the Milky Way and M82. At higher redshifts, cosmological simulations of galaxy formation predict an increase in the frequency and efficiency of galactic outflows due to the increasing star formation activity. Outflows are responsible for removing potential fuel for star formation from the galaxy, while at the same enriching the CGM and the intergalactic medium. These feedback processes, although incorporated as key elements of cosmological simulations, are still poorly constrained on CGM scales. Here we present an ultra-deep MUSE image of the mean MgII emission surrounding a sample of galaxies at z~1 that strongly suggests the presence of outflowing gas on physical scales of more than 10kpc. We find a strong dependence of the detected signal on the inclination of the central galaxy, with edge-on galaxies clearly showing enhanced MgII emission along the minor axis, while face-on galaxies display much weaker and more isotropic emission. We interpret these findings as supporting the idea that outflows typically have a bipolar cone geometry perpendicular to the galactic disk. We demonstrate that the signal is not dominated by a few outliers. After dividing the galaxy sample in subsamples by mass, the bipolar emission is only detected in galaxies with stellar mass $\mathrm{M_* \gtrsim 10^{9.5} M_\odot}$., Comment: published in Nature. 25 pages, 10 figures

The extended Ly$\alpha$ haloes (LAHs) have been found to be prevalent around high-redshift star-forming galaxies. However, the origin of the LAHs is still a subject of debate. Spatially resolved analysis of Ly$\alpha$ profiles provides an important diagnostic. We analyse the average spatial extent and spectral variation of the circumgalactic LAHs by stacking a sample of 155 Ly$\alpha$ emitters (LAEs) at redshift $3

The recent launch of JWST has ushered in a new era of high-redshift astronomy by providing detailed insights into the gas and stellar populations of galaxies in the epoch of reionization. Interpreting these observations and translating them into constraints on the physics of early galaxy formation is a complex challenge that requires sophisticated models of star formation and the interstellar medium (ISM) in high-redshift galaxies. To this end, we present Version 1 of the Sphinx$^{20}$ public data release. Sphinx$^{20}$ is a full box cosmological radiation hydrodynamics simulation that simultaneously models the large-scale process of cosmic reionization and the detailed physics of a multiphase ISM, providing a statistical sample of galaxies akin to those currently being observed by JWST. The data set contains $\sim14,000$ mock images and spectra of the stellar continuum, nebular continuum, and 52 nebular emission lines, including Ly$\alpha$, for each galaxy in Sphinx$^{20}$ with a star formation rate $\geq0.3\ {\rm M_{\odot}\ yr^{-1}}$. All galaxy emission has been processed with dust radiative transfer and/or resonant line radiative transfer, and data is provided for ten viewing angles for each galaxy. Additionally, we provide a comprehensive set of intrinsic galaxy properties, including halo masses, stellar masses, star formation histories, and ISM characteristics (e.g., metallicity, ISM gas densities, LyC escape fractions). This paper outlines the data generation methods, presents a comparative analysis with JWST ERS and Cycle 1 observations, and addresses data set limitations. The Sphinx$^{20}$ data release can be downloaded at the following URL: https://github.com/HarleyKatz/SPHINX-20-data, Comment: 33 pages, 34 figures, Accepted to The Open Journal of Astrophysics, data can be downloaded at the following URL: https://github.com/HarleyKatz/SPHINX-20-data

The Ly$\alpha$ line is a powerful probe of distant galaxies, which contains information about inflowing/outflowing gas through which Ly$\alpha$ photons scatter. To develop our understanding of this probe, we post-process a zoom-in radiation-hydrodynamics simulation of a low-mass ($M_* \sim 10^9 M_\odot$) galaxy to construct 22500 mock spectra in 300 directions from $z = 3$ to 4. Remarkably, we show that one galaxy can reproduce the variety of a large sample of spectroscopically observed Ly$\alpha$ line profiles. While most mock spectra exhibit double-peak profiles with a dominant red peak, their shapes cover a large parameter space in terms of peak velocities, peak separation and flux ratio. This diversity originates from radiative transfer effects at ISM and CGM scales, and depends on galaxy inclination and evolutionary phase. Red-dominated lines preferentially arise in face-on directions during post-starburst outflows and are bright. Conversely, accretion phases usually yield symmetric double peaks in the edge-on direction and are fainter. While resonant scattering effects at $< 0.2\times R_{\rm vir}$ are responsible for the broadening and velocity shift of the red peak, the extended CGM acts as a screen and impacts the observed peak separation. The ability of simulations to reproduce observed Ly$\alpha$ profiles and link their properties with galaxy physical parameters offers new perspectives to use Ly$\alpha$ to constrain the mechanisms that regulate galaxy formation and evolution. Notably, our study implies that deeper Ly$\alpha$ surveys may unveil a new population of blue-dominated lines tracing inflowing gas., Comment: Accepted for publication in MNRAS

Reliable indirect diagnostics of LyC photon escape from galaxies are required to understand which sources were the dominant contributors to reionization. While multiple LyC escape fraction ($f_{\rm esc}$) indicators have been proposed to trace favourable conditions for LyC leakage from the interstellar medium of low-redshift ''analog'' galaxies, it remains unclear whether these are applicable at high redshifts where LyC emission cannot be directly observed. Using a library of 14,120 mock spectra of star-forming galaxies with redshifts $4.64 \leq z \leq 10$ from the SPHINX$^{20}$ cosmological radiation hydrodynamics simulation, we develop a framework for the physics that leads to high $f_{\rm esc}$. We investigate LyC leakage from our galaxies based on the criteria that successful LyC escape diagnostics must i) track a high specific star formation rate, ii) be sensitive to stellar population age in the range $3.5-10~$Myr representing the times when supernova first explode to when LyC production significantly drops, and iii) include a proxy for neutral gas content and gas density in the interstellar medium. ${\rm O}_{32}$, $\Sigma_{\rm SFR}$, M$_{\rm UV}$, and H$\beta$ equivalent width select for one or fewer of our criteria, rendering them either necessary but insufficient or generally poor diagnostics. In contrast, UV slope ($\beta$), and ${\rm E(B-V)}$ match two or more of our criteria, rendering them good $f_{\rm esc}$ diagnostics (albeit with significant scatter). Using our library, we build a quantitative model for predicting $f_{\rm esc}$ based on direct observables. When applied to bright $z > 6$ Ly$\alpha$ emitters observed with JWST, we find that the majority of them have $f_{\rm esc} \lesssim 10\%$., Comment: 16 pages, 15 figures, 1 table, amended to reflect version accepted to MNRAS

We present the second data release of the MUSE Hubble UDF surveys, which includes the deepest spectroscopic survey ever performed. The MUSE data, with their 3D content, amazing depth, wide spectral range, and excellent spatial and medium spectral resolution, are rich in information. This update of the first release incorporates a new 141-hour adaptive-optics-assisted MXDF field (1' diameter FoV) in addition to the reprocessed 10-hour mosaic (3'x3') and the single 31-hour deep field (1'x1'). We have securely identified and measured the redshift of 2221 sources, an increase of 41% compared to the first release. With the exception of 8 stars, the collected sample consists of 25 nearby galaxies (z < 0.25), 677 OII emitters (z=0.25-1.5), 201 galaxies in the MUSE redshift desert range (z=1.5-2.8), and 1308 LAEs (z=2.8-6.7). This represents an order of magnitude more redshifts than the collection of all spectroscopic redshifts obtained before MUSE in the Hubble UDF area (2221 vs 292). At z > 3, the difference is even more striking, with a factor of 65 increase (1308 vs 20). We compared the measured redshifts against three published photometric redshift catalogs and find the photo-z accuracy to be lower than the constraints provided by photo-z fitting codes. 80% of the galaxies have an HST counterpart. They are on average faint, with a median magnitude of 25.7 and 28.7 for the OII and Ly-alpha emitters, respectively. SED fits show that these galaxies tend to be low-mass star-forming galaxies, with a median stellar mass of 6.2 10**8 M and a median SFR of 0.4 M/yr. 20% of our catalog, or 424 galaxies, have no HST counterpart. The vast majority of these new sources are high EQW z>2.8 LAEs that are detected by MUSE thanks to their bright and asymmetric broad Ly-alpha line. We release advanced data products, specific software, and a web interface to select and download data sets., Comment: 46 pages, 48 figures

We use the SPHINX$^{20}$ cosmological radiation hydrodynamics simulation to study how Lyman Continuum (LyC) photons escape from galaxies and the observational signatures of this escape. We define two classes of LyC leaker: Bursty Leakers and Remnant Leakers, based on their star formation rates (SFRs) that are averaged over 10 Myr (SFR$_{10}$) or 100 Myr (SFR$_{100}$). Both have $f_{\rm esc}>20\%$ and experienced an extreme burst of star formation, but Bursty Leakers have ${\rm SFR_{10}>SFR_{100}}$, while Remnant Leakers have ${\rm SFR_{10}20\%$ dominate the ionizing photon budget at $z\gtrsim7.5$ but the contribution from galaxies with $f_{\rm esc}<5\%$ becomes significant at the tail-end of reionization., Comment: 16 pages, 14 figures, accepted by MNRAS

Absorption line spectroscopy offers one of the best opportunities to constrain the properties of galactic outflows and the environment of the circumgalactic medium. Extracting physical information from line profiles is difficult, however, for the physics governing the underlying radiation transfer is complicated and depends on many different parameters. Idealized analytical models are necessary to constrain the large parameter spaces efficiently, but are typically plagued by model degeneracy and systematic errors. Comparison tests with idealized numerical radiation transfer codes offer an excellent opportunity to confront both of these issues. In this paper, we present a detailed comparison between SALT, an analytical radiation transfer model for predicting UV spectra of galactic outflows, with the numerical radiation transfer software, RASCAS. Our analysis has lead to upgrades to both models including an improved derivation of SALT and a customizable adaptive mesh refinement routine for RASCAS. We explore how well SALT, when paired with a Monte Carlo fitting procedure, can recover flow parameters from non-turbulent and turbulent flows. When the velocity and density gradients are excluded, we find that flow parameters are well recovered from high resolution (20 $\rm{km}$ $\rm{s}^{-1}$) data and moderately well from medium resolution (100 $\rm{km}$ $\rm{s}^{-1}$) data without turbulence at a S/N = 10, while derived quantities (e.g., mass outflow rates, column density, etc.) are well recovered at all resolutions. In the turbulent case, biased errors emerge in the recovery of individual parameters, but derived quantities are still well recovered.

We measure escape fractions, $f_{\rm esc}$, of ionizing radiation from galaxies in the SPHINX suite of cosmological radiation-hydrodynamical simulations of reionization, resolving halos with $M_{\rm vir} \gtrapprox 7.5 \times 10^7 \ M_{\odot}$ with a minimum cell width of $\approx 10$ pc. Our new and largest $20$ co-moving Mpc wide volume contains tens of thousands of star-forming galaxies with halo masses up to a few times $10^{11} \ M_{\odot}$. The simulated galaxies agree well with observational constraints of the UV luminosity function in the Epoch of Reionization. The escape fraction fluctuates strongly in individual galaxies over timescales of a few Myrs, due to its regulation by supernova and radiation feedback, and at any given time a tiny fraction of star-forming galaxies emits a large fraction of the ionizing radiation escaping into the inter-galactic medium. Statistically, $f_{\rm esc}$ peaks in intermediate-mass, intermediate-brightness, and low-metallicity galaxies ($M_{*} \approx 10^7 \ M_{\odot}$, $M_{1500} \approx -17$, $Z\lesssim 5 \times 10^{-3} \ Z_{\odot}$), dropping strongly for lower and higher masses, brighter and dimmer galaxies, and more metal-rich galaxies. The escape fraction correlates positively with both the short-term and long-term specific star formation rate. According to SPHINX, galaxies too dim to be yet observed, with $M_{1500} \gtrapprox -17$, provide about $55$ percent of the photons contributing to reionization. The global averaged $f_{\rm esc}$ naturally decreases with decreasing redshift, as predicted by UV background models and low-redshift observations. This evolution is driven by decreasing specific star formation rates over cosmic time., Comment: 20 figures, 21 pages plus appendices. Accepted for publication in MNRAS

Limited constraints on the evolution of the Lyman Continuum (LyC) escape fraction represent one of the primary uncertainties in the theoretical determination of the reionization history. Due to the intervening intergalactic medium (IGM), the possibility of observing LyC photons directly in the epoch of reionization is highly unlikely. For this reason, multiple indirect probes of LyC escape have been identified, some of which are used to identify low-redshift LyC leakers (e.g. O32), while others are primarily useful at $z>6$ (e.g. [OIII]/[CII] far infrared emission). The flux ratio of the resonant MgII doublet emission at 2796$\dot{\rm A}$ and 2803$\dot{\rm A}$ as well as the MgII optical depth have recently been proposed as ideal diagnostics of LyC leakage that can be employed at $z>6$ with JWST. Using state-of-the-art cosmological radiation hydrodynamics simulations post-processed with CLOUDY and resonant-line radiative transfer, we test whether MgII is indeed a useful probe of LyC leakage. Our simulations indicate that the majority of bright, star-forming galaxies with high LyC escape fractions are expected to be MgII emitters rather than absorbers at $z=6$. However, we find that the MgII doublet flux ratio is a more sensitive indicator of dust rather than neutral hydrogen, limiting its use as a LyC leakage indicator to only galaxies in the optically thin regime. Given its resonant nature, we show that MgII will be an exciting probe of the complex kinematics in high-redshift galaxies in upcoming JWST observations., Comment: 25 Pages, 20 Figures, MNRAS Accepted

Using the exquisite MUSE eXtremely Deep Field data, we report the discovery of an MgII emission nebula with an area above a 2$\sigma$ significance level of 1000 proper kpc$^2$, providing the first panoramic view of the spatial distribution of magnesium in the intragroup medium of a low mass group of five star-forming galaxies at z=1.31. The galaxy group members are separated by less than 50 physical kpc in projection and $\approx$120 km/s in velocity space. The most massive galaxy has a stellar mass of 10$^{9.35}$ M$_\odot$ and shows an MgII P-Cygni line profile indicating the presence of an outflow, which is consistent with the spatially resolved spectral analysis showing $\approx+$120 km/s shift of the MgII emission lines with respect to the systemic redshift. The other galaxies are less massive and only show MgII in emission. The detected MgII nebula has a maximal projected extent of $\approx$70 kpc including a low surface brightness (2 $\times$ 10$^{-19}$ erg/s/cm$^{2}$/arcsec$^{2}$) gaseous bridge between two subgroups of galaxies. The presence of absorption features in the spectrum of a background galaxy located at an impact parameter of 19 kpc from the closest galaxy of the group indicates the presence of gas enriched in magnesium even beyond the detected nebula seen in emission, suggesting that we are observing the tip of a larger intragroup medium. The observed MgII velocity gradient suggests an overall rotation of the structure along the major axis of the most massive galaxy. Our MUSE data also reveal extended Fe II* emission in the vicinity of the most massive galaxy, aligned with its minor axis. Extended [OII] emission is found around the galaxy group members and at the location of the MgII bridge. Our results suggest that both tidal stripping effects from galaxy interactions and outflows are enriching the intragroup medium of this system., Comment: Submitted to A&A. Referee's first comments implemented

Understanding the escape of Lyman continuum (LyC) and Lyman $\alpha$ (Ly$\alpha$) photons from giant molecular clouds (GMCs) is crucial if we are to study the reionization of the Universe and to interpret spectra of observed galaxies at high redshift. To this end, we perform high-resolution, radiation-magneto-hydrodynamic simulations of GMCs with self-consistent star formation and stellar feedback. We find that a significant fraction (15-70%) of ionizing radiation escapes from the simulated GMCs with different masses ($10^5$ and $10^6\,M_\odot$), as the clouds are dispersed within about $2$-$5\,{\rm Myr}$ from the onset of star formation. The fraction of LyC photons leaked is larger when the GMCs are less massive, metal-poor, less turbulent, and less dense. The most efficient leakage of LyC radiation occurs when the total star formation efficiency of a GMC is about 20%. The escape of Ly$\alpha$ shows a trend similar to that of LyC photons, except that the fraction of Ly$\alpha$ photons escaping from the GMCs is larger ($f_{\rm esc}^{\rm Ly\alpha}\approx f_{900}^{0.27}$) and that a GMC with strong turbulence shows larger $f_{\rm esc}^{\rm Ly\alpha}$. The simulated GMCs show a characteristic velocity separation of $\Delta v\approx 120 \,{\rm km\,s^{-1}}$ in the time-averaged emergent Ly$\alpha$ spectra, suggesting that Ly$\alpha$ could be useful to infer the kinematics of the interstellar and circumgalactic medium. We show that Ly$\alpha$ luminosities are a useful indicator of the LyC escape, provided the number of LyC photons can be deduced through stellar population modeling. Finally, we find that the correlations between the escape fractions of Ly$\alpha$, ultraviolet photons at 1500A, and the Balmer $\alpha$ line are weak., Comment: 32 pages, 20 figures, Accepted for publication in the ApJS

We present Keck/LRIS follow-up spectroscopy for 13 photometric candidates of extremely metal poor galaxies (EMPGs) selected by a machine-learning technique applied to the deep (~26 AB mag) optical and wide-area (~500 deg^2) Subaru imaging data in the EMPRESS survey. Nine out of the 13 candidates are EMPGs with an oxygen abundance (O/H) less than ~10% solar value (O/H)_sun, and four sources are contaminants of moderately metal-rich galaxies or no emission-line objects. Notably, two out of the nine EMPGs have extremely-low stellar masses and oxygen abundances of 5*10^4--7*10^5 M_sun and 2--3% (O/H)_sun, respectively. With a sample of five EMPGs with (Fe/O) measurements, two (three) of which are taken from this study (the literature), we confirm that two EMPGs with the lowest (O/H) ratios of ~2% (O/H)_sun show high (Fe/O) ratios of ~0.1, close to the solar abundance ratio. Comparing galaxy chemical enrichment models, we find that the two EMPGs cannot be explained by a scenario of metal-poor gas accretion/episodic star-formation history due to their low (N/O) ratios. We conclude that the two EMPGs can be reproduced by an inclusion of bright hypernovae and/or hypothetical pair-instability supernovae (SNe) preferentially produced in a metal-poor environment. This conclusion implies that primordial galaxies at z~10 could have a high abundance of Fe that is not originated from Type Ia SNe with delays, and that Fe may not serve as a cosmic clock for primordial galaxies., Comment: ApJ in Press

ALMA observations of $z>6$ galaxies have revealed abnormally high [OIII]$_{\rm 88\mu m}$/[CII]$_{\rm 158\mu m}$ ratios and [CII]$_{\rm 158\mu m}$ deficits compared to local galaxies. Numerous solutions have been proposed including differences in C and O abundance ratios, observational bias, and differences in ISM properties. In order to elucidate the underlying physics that drives this high-redshift phenomenon, we employ SPHINX$^{20}$, a state-of-the-art, cosmological radiation-hydrodynamics simulation, that resolves detailed ISM properties of thousands of galaxies in the epoch of reionization. We find that the observed $z>6$ [OIII]-SFR and [CII]-SFR relations can only be reproduced when the C/O abundance ratio is $\sim8\times$ lower than Solar and the total metal production is $\sim4\times$ higher than that of a Kroupa IMF. This implies that high-redshift galaxies are potentially primarily enriched by low-metallicity core-collapse supernovae with a more top-heavy IMF. As AGB stars and type-Ia supernova begin to contribute to the galaxy metallicity, both the [CII]-SFR and [CII] luminosity functions are predicted to converge to observed values at $z\sim4.5$. While we demonstrate that ionisation parameter, LyC escape fraction, ISM gas density, and CMB attenuation all drive galaxies towards higher [OIII]/[CII], observed values at $z>6$ can only be reproduced with substantially lower C/O abundances compared to Solar. The combination of [CII] and [OIII] can be used to predict the values of ionisation parameter, ISM gas density, and LyC escape fraction and we provide estimates of these quantities for nine observed $z>6$ galaxies. Finally, we demonstrate that [OI]$_{\rm 63\mu m}$ can be used as a replacement for [CII] in high-redshift galaxies where [CII] is unobserved and argue that more observation time should be used to target [OI] at $z>6$. (Abridged), Comment: 20 pages, 18 figures, Accepted by MNRAS

The Ly$\alpha$ emission line is one of the most promising probes of cosmic reionisation but isolating the signature of a change in the ionisation state of the IGM is challenging because of intrinsic evolution and internal radiation transfer effects. We present the first study of the evolution of Ly$\alpha$ emitters (LAE) during the epoch of reionisation based on a full radiation-hydrodynamics cosmological simulation that is able to capture both the large-scale process of reionisation and the small-scale properties of galaxies. We predict the Ly$\alpha$ emission of galaxies in the $10^3$ cMpc$^3$ SPHINX simulation at $6\leq z\leq9$ by computing the full Ly$\alpha$ radiation transfer from ISM to IGM scales. SPHINX is able to reproduce many observational constraints such as the UV/Ly$\alpha$ luminosity functions and stellar mass functions at z $\geq$ 6 for the dynamical range probed by our simulation ($M_{\rm 1500}\gtrsim-18$, $L_{\rm Ly\alpha}\lesssim10^{42}$ erg/s, $M_{\star}\lesssim10^9$ M$_{\odot}$). As intrinsic Ly$\alpha$ emission and internal Ly$\alpha$ escape fractions barely evolve from $z=6$ to 9, the observed suppression of Ly$\alpha$ luminosities with increasing redshift is fully attributed to IGM absorption. For most observable galaxies ($M_{\rm 1500}\lesssim-16$), the Ly$\alpha$ line profiles are slightly shifted to the red due to internal radiative transfer effects which mitigates the effect of IGM absorption. Overall, the enhanced Ly$\alpha$ suppression during reionisation traces the IGM neutral fraction $x_{\rm HI}$ well but the predicted amplitude of this reduction is a strong function of the Ly$\alpha$ peak shift, which is set at ISM/CGM scales. We find that a large number of LAEs could be detectable in very deep surveys during reionisation when $x_{\rm HI}$ is still $\approx 50\%$., Comment: MNRAS accepted, 20 pages, 19 figures, plus appendices

We report the discovery of diffuse extended Ly-alpha emission from redshift 3.1 to 4.5, tracing cosmic web filaments on scales of 2.5-4 comoving Mpc. These structures have been observed in overdensities of Ly-alpha emitters in the MUSE Extremely Deep Field, a 140 hour deep MUSE observation located in the Hubble Ultra Deep Field. Among the 22 overdense regions identified, 5 are likely to harbor very extended Ly-alpha emission at high significance with an average surface brightness of $\mathrm{5 \times 10^{-20} erg s^{-1} cm^{-2} arcsec^{-2}}$. Remarkably, 70% of the total Ly-alpha luminosity from these filaments comes from beyond the circumgalactic medium of any identified Ly-alpha emitters. Fluorescent Ly-alpha emission powered by the cosmic UV background can only account for less than 34% of this emission at z$\approx$3 and for not more than 10% at higher redshift. We find that the bulk of this diffuse emission can be reproduced by the unresolved Ly-alpha emission of a large population of ultra low luminosity Ly-alpha emitters ($\mathrm{<10^{40} erg s^{-1}}$), provided that the faint end of the Ly-alpha luminosity function is steep ($\alpha \lessapprox -1.8$), it extends down to luminosities lower than $\mathrm{10^{38} - 10^{37} erg s^{-1}}$ and the clustering of these Ly-alpha emitters is significant (filling factor $< 1/6$). If these Ly-alpha emitters are powered by star formation, then this implies their luminosity function needs to extend down to star formation rates $\mathrm{< 10^{-4} M_\odot yr^{-1}}$. These observations provide the first detection of the cosmic web in Ly-alpha emission in typical filamentary environments and the first observational clue for the existence of a large population of ultra low luminosity Ly-alpha emitters at high redshift., Comment: 28 pages, 19 figures

We present the first results from SPHINX-MHD, a suite of cosmological radiation-magnetohydrodynamics simulations designed to study the impact of primordial magnetic fields (PMFs) on galaxy formation and the evolution of the intergalactic medium during the epoch of reionization. The simulations are among the first to employ multi-frequency, on-the-fly radiation transfer and constrained transport ideal MHD in a cosmological context to simultaneously model the inhomogeneous process of reionization as well as the growth of PMFs. We run a series of $(5\,\text{cMpc})^3$ cosmological volumes, varying both the strength of the seed magnetic field ($B_0$) and its spectral index ($n_B$). We find that PMFs that have $n_B > -0.562\log_{10}\left(\frac{B_0}{1{\rm n}G}\right) - 3.35$ produce electron optical depths ($\tau_e$) that are inconsistent with CMB constraints due to the unrealistically early collapse of low-mass dwarf galaxies. For $n_B\geq-2.9$, our constraints are considerably tighter than the $\sim{\rm n}G$ constraints from Planck. PMFs that do not satisfy our constraints have little impact on the reionization history or the shape of the UV luminosity function. Likewise, detecting changes in the Lya forest due to PMFs will be challenging because photoionisation and photoheating efficiently smooth the density field. However, we find that the first absorption feature in the global 21cm signal is a sensitive indicator of the properties of the PMFs, even for those that satisfy our $\tau_e$ constraint. Furthermore, strong PMFs can marginally increase the escape of LyC photons by up to 25\% and shrink the effective radii of galaxies by $\sim44\%$ which could increase the completeness fraction of galaxy surveys. Finally, our simulations show that surveys with a magnitude limit of ${\rm M_{UV,1500}=-13}$ can probe the sources that provide the majority of photons for reionization out to $z=12$., Comment: 29 pages, 27 figures, accepted MNRAS

Galaxies in dense environments are subject to interactions and mechanisms which directly affect their evolution by lowering their gas fractions and reducing their star-forming capacity earlier than their isolated counterparts. The aim of our project is to get new insights about the role of environment on the stellar and baryonic content of galaxies using a kinematic approach, through the study of the Tully-Fisher relation (TFR). We study a sample of galaxies in 8 groups spanning a redshift range of $0.5

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., Comment: Accepted for publication in A&A. Code available at http://rascas.univ-lyon1.fr/

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$..., Comment: 21 pages, 11 figures, MNRAS Accepted

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., Comment: 23 pages, 20 figures, accepted for publication in MNRAS

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., Comment: 17 pages, 12 figures, 2 tables, accepted for publication in A&A

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., Comment: 18 pages, 15 figures. Accepted for publication in MNRAS

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., Comment: 19 pages, 16 figures, Accepted for publication in MNRAS

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., Comment: 21 pages, 12 figures, 3 tables, accepted for publication in A&A (27 October 2017)

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., Comment: 25 pages, 27 figures, A&A for accepted (26 June 2017)

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., Comment: 24 pages, 17 figures, submitted to MNRAS after responding to a first referee report

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., Comment: 11 pages, 7 figures

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., Comment: MNRAS in press, On-line version of the catalog of galaxies with available deep images and figures (including navigation tools) are available from http://irfu.cea.fr/Projets/matlas/public/Atlas3D/atlas3d-pXXIX.html

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., Comment: 11 Pages, 7 Figures, 3 Tables, Accepted for publication in ApJ

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., Comment: MNRAS, accepted

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], Comment: 29 pages, 29 figures, LaTeX. Accepted for publication in MNRAS

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], Comment: 13 pages, 13 figures, accepted for publication in MNRAS

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)., Comment: 20 pages, 14 figures, MNRAS in press

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., Comment: 6 pages, 4 figures, to match version published in MNRAS Letters (http://onlinelibrary.wiley.com/doi/10.1111/j.1745-3933.2011.01035.x/abstract)

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., Comment: Published in ApJ Letters

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., Comment: 10 pages, 11 Postscript figures, published in A&A

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., Comment: Major revision, A&A in press

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., Comment: 14 pages, 16 figures (13 in eps, 3 in jpg format). Accepted for publication in A&A