Showing 1–50 of 63 results for author: Michel-Dansac, L

The MUSE eXtremely Deep Field: Detections of circumgalactic SiII* emission at z>~2

Authors: Haruka Kusakabe, Valentin Mauerhofer, Anne Verhamme, Thibault Garel, Jeremy Blaizot, Lutz Wisotzki, Johan Richard, Leindert A. Boogaard, Floriane Leclercq, Yucheng Guo, Adelaide Claeyssens, Thierry Contini, Edmund Christian Herenz, Josephine Kerutt, Michael V. Maseda, Leo Michel-Dansac, Themiya Nanayakkara, Masami Ouchi, Ismael Pessa, Joop Schaye

Abstract: 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 de… ▽ More 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. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 13 pages. 8 figures. Submitted to A&A. Comments are welcome!

Stellar angular momentum of intermediate redshift galaxies in MUSE surveys

Authors: Constanza Muñoz López, Davor Krajnović, B. Epinat, Y. Herrero-Alonso, T. Urrutia, W. Mercier, N. F. Bouché, L. A. Boogaard, T. Contini, L. Michel-Dansac, I. Pessa

Abstract: We quantify the stellar rotation of galaxies by computing the $λ_{R}$ parameter, a proxy for the stellar angular momentum in a sample of 106 galaxies with redshift 0.1 $<$ z $<$ 0.8 and stellar masses from $\sim$10$^{7.5}$ to 10$^{11.8}$ M$_{\odot}$. The sample is located in the CANDELS/GOODS-S and COSMOS fields, and it was observed by various MUSE surveys. We create stellar velocity and velocity… ▽ More We quantify the stellar rotation of galaxies by computing the $λ_{R}$ parameter, a proxy for the stellar angular momentum in a sample of 106 galaxies with redshift 0.1 $<$ z $<$ 0.8 and stellar masses from $\sim$10$^{7.5}$ to 10$^{11.8}$ M$_{\odot}$. The sample is located in the CANDELS/GOODS-S and COSMOS fields, and it was observed by various MUSE surveys. We create stellar velocity and velocity dispersion maps using a full-spectrum fitting technique, covering spatially $\sim$2$R_{e}$ for the galaxies. We study the impact of the atmospheric seeing on the spin parameter and apply corrections when pertinent. Through the analysis of the $λ_{R}-ε$ diagram, we notice that the fraction of round and massive galaxies increases with redshift. We lack galaxies with $λ_{R}$ < 0.1 in the sample and we find only one potential, but uncertain, low-mass slow rotator at z $\sim0.3$. Moreover, we do not see an evident evolution or trend in the stellar angular momentum with redshift. We characterize the sample environment using two indicators: a local estimator based on the Voronoi tesselation method, and a global estimator derived by the use of the Friends-of-Friends algorithm. We find no correlation between the environment and $λ_{R}$ given that we are not probing dense regions or massive galaxy structures. We also analyze the kinematic maps of the sample finding that about 40$\%$ of galaxies are consistent with being regular rotators, having rotating stellar discs with flat velocity dispersion maps, while $\sim20\%$ of galaxies have complex velocity maps and can be identified as non-regular rotators in spite of their $λ_{R}$ values. For the remaining galaxies the classification is uncertain. As we lack galaxies with $λ_{R}$< 0.1, we are not able to identify when galaxies become slow rotators within the surveyed environments, area and redshift range. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: 26 pages, 18 figures, full abstract on the paper

Bipolar Outflows out to 10~kpc for Massive Galaxies at Redshift $z\approx 1$

Authors: Yucheng Guo, Roland Bacon, Nicolas F. Bouché, Lutz Wisotzki, Joop Schaye, Jérémy Blaizot, Anne Verhamme, Sebastiano Cantalupo, Leindert A. Boogaard, Jarle Brinchmann, Maxime Cherrey, Haruka Kusakabe, Ivanna Langan, Floriane Leclercq, Jorryt Matthee, Léo Michel-Dansac, Ilane Schroetter, Martin Wendt

Abstract: 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… ▽ More 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}$. △ Less

Submitted 8 December, 2023; originally announced December 2023.

Comments: published in Nature. 25 pages, 10 figures

Journal ref: 2023, Nature, 624, 53

MAGIC: Muse gAlaxy Groups In Cosmos -- A survey to probe the impact of environment on galaxy evolution over the last 8 Gyr

Authors: B. Epinat, T. Contini, W. Mercier, L. Ciesla, B. C. Lemaux, S. D. Johnson, J. Richard, J. Brinchmann, L. A. Boogaard, D. Carton, L. Michel-Dansac, R. Bacon, D. Krajnovic, H. Finley, I. Schroetter, E. Ventou, V. Abril-Melgarejo, A. Boselli, N. F. Bouché, W. Kollatschny, K. Kovac, M. Paalvast, G. Soucail, T. Urrutia, P. M. Weilbacher

Abstract: We introduce the MUSE gAlaxy Groups in COSMOS (MAGIC) survey, which was built to study the impact of environment on galaxy evolution over the last 8 Gyr. It consists of 17 MUSE fields targeting 14 massive structures at intermediate redshift ($0.3<z<0.8$) in the COSMOS area. We securely measured the redshifts for 1419 sources and identified 76 galaxy pairs and 67 groups of at least 3 members using… ▽ More We introduce the MUSE gAlaxy Groups in COSMOS (MAGIC) survey, which was built to study the impact of environment on galaxy evolution over the last 8 Gyr. It consists of 17 MUSE fields targeting 14 massive structures at intermediate redshift ($0.3<z<0.8$) in the COSMOS area. We securely measured the redshifts for 1419 sources and identified 76 galaxy pairs and 67 groups of at least 3 members using a friends-of-friends algorithm. The environment of galaxies is quantified from group properties, as well as from global and local density estimators. The MAGIC survey has increased the number of objects with a secure spectroscopic redshift over its footprint by a factor of about 5. Most of the new redshifts have apparent magnitudes in the $z^{++}$ band $z_{app}^{++}>21.5$. The spectroscopic redshift completeness is high: in the redshift range of [OII] emitters ($0.25 \le z < 1.5$), where most of the groups are found, it globally reaches a maximum of 80% down to $z_{app}^{++}=25.9$, and locally decreases from $\sim 100$% to $\sim50$% in magnitude bins from $z_{app}^{++}=23-24$ to $z_{app}^{++}=25.5$. We find that the fraction of quiescent galaxies increases with local density and with the time spent in groups. A morphological dichotomy is also found between bulge-dominated quiescent and disk-dominated star-forming galaxies. As environment gets denser, the peak of the stellar mass distribution shifts towards $M_*>10^{10}~M_\odot$, and the fraction of galaxies with $M_*<10^9~M_\odot$ decreases significantly, even for star-forming galaxies. We also highlight peculiar features such as close groups, extended nebulae, and a gravitational arc. Our results suggest that galaxies are preprocessed in groups of increasing mass before entering rich groups and clusters. We publicly release two catalogs containing the properties of galaxies and groups, respectively. △ Less

Submitted 28 March, 2024; v1 submitted 1 December, 2023; originally announced December 2023.

Comments: 27 pages, 22 figures, accepted for publication in A&A

Journal ref: A&A, 683, A205 (2024)

Spatially-resolved Spectroscopic Analysis of Ly$α$ Haloes: Radial Evolution of the Ly$α$ Line Profile out to 60 kpc

Authors: Yucheng Guo, Roland Bacon, Lutz Wisotzki, Thibault Garel, Jérémy Blaizot, Joop Schaye, Jorryt Matthee, Floriane Leclercq, Leindert Boogaard, Johan Richard, Anne Verhamme, Jarle Brinchmann, Léo Michel-Dansac, Haruka Kusakabe

Abstract: Deep MUSE observations have unveiled pervasive Ly$α$ haloes (LAHs) surrounding high-redshift star-forming galaxies. However, the origin of the extended Ly$α$ emission is still a subject of debate. We analyse the average spatial extent and spectral variation of the circumgalactic LAHs by stacking a sample of 155 Ly$α$ emitters (LAEs) at redshift $3<z<4$ in the MUSE Extremely Deep Field (MXDF). With… ▽ More Deep MUSE observations have unveiled pervasive Ly$α$ haloes (LAHs) surrounding high-redshift star-forming galaxies. However, the origin of the extended Ly$α$ emission is still a subject of debate. We analyse the average spatial extent and spectral variation of the circumgalactic LAHs by stacking a sample of 155 Ly$α$ emitters (LAEs) at redshift $3<z<4$ in the MUSE Extremely Deep Field (MXDF). With respect to the Ly$α$ red peak of the target LAE, the Ly$α$ line peak becomes increasingly more blueshifted out to a projected distance of at least 60 kpc, where the velocity offset is $\approx$ 250 km/s. This signal is observed in both the mean and median stacks, and is thus a generic property of the LAE sample with typical Ly$α$ luminosity $\mathrm{\approx 10^{41.1} erg\,s^{-1}}$. We discuss multiple scenarios to explain the blueshift of the circumgalactic Ly$α$ line. The most plausible one is a combination of outflows and inflows. In the inner region of the LAH, the Ly$α$ photons are produced by the central star formation and then scattered within outflows. At larger radii, the infalling cool gas shapes the observed Ly$α$ blueshift. △ Less

Submitted 12 September, 2023; originally announced September 2023.

Comments: Submitted to A&A Letters

The Sphinx Public Data Release: Forward Modelling High-Redshift JWST Observations with Cosmological Radiation Hydrodynamics Simulations

Authors: Harley Katz, Joki Rosdahl, Taysun Kimm, Jeremy Blaizot, Nicholas Choustikov, Marion Farcy, Thibault Garel, Martin G. Haehnelt, Leo Michel-Dansac, Pierre Ocvirk

Abstract: 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 medi… ▽ More 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$α$, 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 △ Less

Submitted 3 December, 2023; v1 submitted 6 September, 2023; originally announced September 2023.

Comments: 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

Simulating the diversity of shapes of the Lyman-$α$ line

Authors: Jeremy Blaizot, Thibault Garel, Anne Verhamme, Harley Katz, Taysun Kimm, Léo Michel-Dansac, Peter D. Mitchell, Joakim Rosdahl, Maxime Trebitsch

Abstract: The Ly$α$ line is a powerful probe of distant galaxies, which contains information about inflowing/outflowing gas through which Ly$α$ 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… ▽ More The Ly$α$ line is a powerful probe of distant galaxies, which contains information about inflowing/outflowing gas through which Ly$α$ 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$α$ 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$α$ profiles and link their properties with galaxy physical parameters offers new perspectives to use Ly$α$ to constrain the mechanisms that regulate galaxy formation and evolution. Notably, our study implies that deeper Ly$α$ surveys may unveil a new population of blue-dominated lines tracing inflowing gas. △ Less

Submitted 17 May, 2023; originally announced May 2023.

Comments: Accepted for publication in MNRAS

The Physics of Indirect Estimators of Lyman Continuum Escape and their Application to High-Redshift JWST Galaxies

Authors: Nicholas Choustikov, Harley Katz, Aayush Saxena, Alex J. Cameron, Julien Devriendt, Adrianne Slyz, Joki Rosdahl, Jeremy Blaizot, Leo Michel-Dansac

Abstract: 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 h… ▽ More 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}$, $Σ_{\rm SFR}$, M$_{\rm UV}$, and H$β$ equivalent width select for one or fewer of our criteria, rendering them either necessary but insufficient or generally poor diagnostics. In contrast, UV slope ($β$), 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$α$ emitters observed with JWST, we find that the majority of them have $f_{\rm esc} \lesssim 10\%$. △ Less

Submitted 13 March, 2024; v1 submitted 17 April, 2023; originally announced April 2023.

Comments: 16 pages, 15 figures, 1 table, amended to reflect version accepted to MNRAS

The MUSE Hubble Ultra Deep Field surveys: Data release II

Authors: Roland Bacon, Jarle Brinchmann, Simon Conseil, Michael Maseda, Themiya Nanayakkara, Martin Wendt, Raphael Bacher, David Mary, Peter M. Weilbacher, Davor Krajnovic, Leindert Boogaard, Nicolas Bouche, Thierry Contini, Benoit Epinat, Anna Feltre, Yucheng Guo, Christian Herenz, Wolfram Kollatschny, Haruka Kusakabe, Floriane Leclercq, Leo Michel-Dansac, Roser Pello, Johan Richard, Martin Roth, Gregory Salvignol , et al. (8 additional authors not shown)

Abstract: 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) i… ▽ More 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. △ Less

Submitted 5 December, 2022; v1 submitted 15 November, 2022; originally announced November 2022.

Comments: 46 pages, 48 figures

Journal ref: A&A 670, A4 (2023)

Two Modes of LyC Escape From Bursty Star Formation: Implications for [C II] Deficits and the Sources of Reionization

Authors: Harley Katz, Aayush Saxena, Joki Rosdahl, Taysun Kimm, Jeremy Blaizot, Thibault Garel, Leo Michel-Dansac, Martin Haehnelt, Richard S. Ellis, Laura Penterrici, Julien Devriendt, Adrianne Slyz

Abstract: 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\%$… ▽ More 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}<SFR_{100}}$. The maximum SFRs in these bursts were typically $\sim100$ times greater than the SFR of the galaxy prior to the burst, a rare $2σ$ outlier among the general high-redshift galaxy population. Bursty Leakers are qualitatively similar to ionization-bounded nebulae with holes, exhibiting high ionization parameters and typical HII region gas densities. Remnant Leakers show properties of density-bounded nebulae, having normal ionization parameters but much lower HII region densities. Both types of leaker exhibit [CII]$_{\rm 158μm}$ deficits on the [CII]-SFR$_{100}$ relation, while only Bursty Leakers show deficits when SFR$_{10}$ is used. We predict that [CII] luminosity and SFR indicators such as H$α$ and M$_{\rm 1500Å}$ can be combined to identify both types of LyC leaker and the mode by which photons are escaping. These predictions can be tested with [CII] observations of known $z=3-4$ LyC leakers. Finally, we show that leakers with $f_{\rm esc}>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. △ Less

Submitted 17 October, 2022; originally announced October 2022.

Comments: 16 pages, 14 figures, accepted by MNRAS

Testing SALT Approximations with Numerical Radiation Transfer Code Part 1: Validity and Applicability

Authors: Cody Carr, Leo Michel-Dansac, Jeremy Blaizot, Claudia Scarlata, Alaina Henry, Anne Verhamme

Abstract: 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… ▽ More 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. △ Less

Submitted 28 September, 2022; originally announced September 2022.

LyC escape from SPHINX galaxies in the Epoch of Reionization

Authors: Joakim Rosdahl, Jérémy Blaizot, Harley Katz, Taysun Kimm, Thibault Garel, Martin Haehnelt, Laura C. Keating, Sergion Martin-Alvarez, Léo Michel-Dansac, Pierre Ocvirk

Abstract: 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 ha… ▽ More 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. △ Less

Submitted 7 July, 2022; originally announced July 2022.

Comments: 20 figures, 21 pages plus appendices. Accepted for publication in MNRAS

MgII in the JWST Era: a Probe of Lyman Continuum Escape?

Authors: Harley Katz, Thibault Garel, Joakim Rosdahl, Valentin Mauerhofer, Taysun Kimm, Jérémy Blaizot, Léo Michel-Dansac, Julien Devriendt, Adrianne Slyz, Martin Haehnelt

Abstract: 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 be… ▽ More 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. △ Less

Submitted 23 May, 2022; originally announced May 2022.

Comments: 25 Pages, 20 Figures, MNRAS Accepted

Scaling relations of z~0.25-1.5 galaxies in various environments from the morpho-kinematic analysis of the MAGIC sample

Authors: W. Mercier, B. Epinat, T. Contini, V. Abril-Melgarejo, L. Boogaard, J. Brinchmann, H. Finley, D. Krajnović, L. Michel-Dansac, E. Ventou, N. Bouché, J. Dumoulin, Juan C. B. Pineda

Abstract: The evolution of galaxies is influenced by many physical processes which may vary depending on their environment. We combine Hubble Space Telescope (HST) and Multi-Unit Spectroscopic Explorer (MUSE) data of galaxies at 0.25<z<1.5 to probe the impact of environment on the size-mass relation, the Main Sequence (MS) and the Tully-Fisher relation (TFR). We perform a morpho-kinematic modelling of 593… ▽ More The evolution of galaxies is influenced by many physical processes which may vary depending on their environment. We combine Hubble Space Telescope (HST) and Multi-Unit Spectroscopic Explorer (MUSE) data of galaxies at 0.25<z<1.5 to probe the impact of environment on the size-mass relation, the Main Sequence (MS) and the Tully-Fisher relation (TFR). We perform a morpho-kinematic modelling of 593 [Oii] emitters in various environments in the COSMOS area from the MUSE-gAlaxy Groups In Cosmos (MAGIC) survey. The HST F814W images are modelled with a bulge-disk decomposition to estimate their bulge-disk ratio, effective radius and disk inclination. We use the [Oii]λλ3727, 3729 doublet to extract the ionised gas kinematic maps from the MUSE cubes, and we model them for a sample of 146 [Oii] emitters, with bulge and disk components constrained from morphology and a dark matter halo. We find an offset of 0.03 dex on the size-mass relation zero point between the field and the large structure subsamples, with a richness threshold of N=10 to separate between small and large structures, and of 0.06 dex with N=20. Similarly, we find a 0.1 dex difference on the MS with N=10 and 0.15 dex with N=20. These results suggest that galaxies in massive structures are smaller by 14% and have star formation rates reduced by a factor of 1.3-1.5 with respect to field galaxies at z=0.7. Finally, we do not find any impact of the environment on the TFR, except when using N=20 with an offset of 0.04 dex. We discard the effect of quenching for the largest structures that would lead to an offset in the opposite direction. We find that, at z=0.7, if quenching impacts the mass budget of galaxies in structures, these galaxies would have been affected quite recently, for roughly 0.7-1.5 Gyr. This result holds when including the gas mass, but vanishes once we include the asymmetric drift correction. △ Less

Submitted 22 April, 2022; v1 submitted 19 April, 2022; originally announced April 2022.

Comments: 32 pages, 29 figures, 5 tables, 7 appendices, language checked and appendix table updated (accepted for publication in A&A)

Journal ref: A&A 665, A54 (2022)

The MUSE Extremely Deep Field: a first panoramic view of an Mg II emitting intragroup medium

Authors: Floriane Leclercq, Anne Verhamme, Benoit Epinat, Charlotte Simmonds, Jorryt Matthee, Nicolas F. Bouché, Thibault Garel, Tanya Urrutia, Lutz Wisotzki, Johannes Zabl, Roland Bacon, Valentina Abril-Melgarejo, Leindert Boogaard, Jarle Brinchmann, Sebastiano Cantalupo, Thierry Contini, Josephine Kerutt, Haruka Kusakabe, Michael Maseda, Léo Michel-Dansac, Sowgat Muzahid, Themiya Nanayakkara, Johan Richard, Joop Schaye

Abstract: Using the exquisite MUSE eXtremely Deep Field data, we report the discovery of an MgII emission nebula with an area above a 2$σ$ 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 kp… ▽ More Using the exquisite MUSE eXtremely Deep Field data, we report the discovery of an MgII emission nebula with an area above a 2$σ$ 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. △ Less

Submitted 10 March, 2022; originally announced March 2022.

Comments: Submitted to A&A. Referee's first comments implemented

Journal ref: A&A 663, A11 (2022)

A systematic study of the escape of LyC and Ly$α$ photons from star-forming, magnetized turbulent clouds

Authors: Taysun Kimm, Rebekka Bieri, Sam Geen, Joakim Rosdahl, Jérémy Blaizot, Léo Michel-Dansac, Thibault Garel

Abstract: Understanding the escape of Lyman continuum (LyC) and Lyman $α$ (Ly$α$) 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… ▽ More Understanding the escape of Lyman continuum (LyC) and Lyman $α$ (Ly$α$) 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$α$ shows a trend similar to that of LyC photons, except that the fraction of Ly$α$ photons escaping from the GMCs is larger ($f_{\rm esc}^{\rm Lyα}\approx f_{900}^{0.27}$) and that a GMC with strong turbulence shows larger $f_{\rm esc}^{\rm Lyα}$. The simulated GMCs show a characteristic velocity separation of $Δv\approx 120 \,{\rm km\,s^{-1}}$ in the time-averaged emergent Ly$α$ spectra, suggesting that Ly$α$ could be useful to infer the kinematics of the interstellar and circumgalactic medium. We show that Ly$α$ 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$α$, ultraviolet photons at 1500A, and the Balmer $α$ line are weak. △ Less

Submitted 2 December, 2021; v1 submitted 6 October, 2021; originally announced October 2021.

Comments: 32 pages, 20 figures, Accepted for publication in the ApJS

EMPRESS. IV. Extremely Metal-Poor Galaxies (EMPGs) Including Very Low-Mass Primordial Systems with M*=10^4--10^5 M_sun and 2--3% (O/H)_sun: High (Fe/O) Suggestive of Metal Enrichment by Hypernovae/Pair-Instability Supernovae

Authors: Yuki Isobe, Masami Ouchi, Akihiro Suzuki, Takashi Moriya, Kimihiko Nakajima, Ken'ichi Nomoto, Michael Rauch, Yuichi Harikane, Takashi Kojima, Yoshiaki Ono, Seiji Fujimoto, Akio K. Inoue, Ji Hoon Kim, Yutaka Komiyama, Haruka Kusakabe, Chien-Hsiu Lee, Michael Maseda, Jorryt Matthee, Leo Michel-Dansac, Tohru Nagao, Themiya Nanayakkara, Moka Nishigaki, Masato Onodera, Yuma Sugahara, Yi Xu

Abstract: 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 con… ▽ More 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. △ Less

Submitted 3 November, 2021; v1 submitted 9 August, 2021; originally announced August 2021.

Comments: ApJ in Press

The Nature of High ${\rm [OIII]}_{\rm 88μm}$/${\rm [CII]}_{\rm 158μm}$ Galaxies in the Epoch of Reionization: Low Carbon Abundance and a Top-Heavy IMF?

Authors: Harley Katz, Joakim Rosdahl, Taysun Kimm, Thibault Garel, Jérémy Blaizot, Martin G. Haehnelt, Léo Michel-Dansac, Sergio Martin-Alvarez, Julien Devriendt, Adrianne Slyz, Romain Teyssier, Pierre Ocvirk, Nicolas Laporte, Richard Ellis

Abstract: ALMA observations of $z>6$ galaxies have revealed abnormally high [OIII]$_{\rm 88μm}$/[CII]$_{\rm 158μm}$ ratios and [CII]$_{\rm 158μ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 phenom… ▽ More ALMA observations of $z>6$ galaxies have revealed abnormally high [OIII]$_{\rm 88μm}$/[CII]$_{\rm 158μm}$ ratios and [CII]$_{\rm 158μ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μ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) △ Less

Submitted 5 January, 2022; v1 submitted 2 August, 2021; originally announced August 2021.

Comments: 20 pages, 18 figures, Accepted by MNRAS

Lyman-$α$ as a tracer of cosmic reionisation in the SPHINX radiation-hydrodynamics cosmological simulation

Authors: Thibault Garel, Jérémy Blaizot, Joakim Rosdahl, Léo Michel-Dansac, Martin G. Haehnelt, Harley Katz, Taysun Kimm, Anne Verhamme

Abstract: The Ly$α$ 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$α$ emitters (LAE) during the epoch of reionisation based on a full radiation-hydrodynamics cosmological s… ▽ More The Ly$α$ 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$α$ 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$α$ emission of galaxies in the $10^3$ cMpc$^3$ SPHINX simulation at $6\leq z\leq9$ by computing the full Ly$α$ radiation transfer from ISM to IGM scales. SPHINX is able to reproduce many observational constraints such as the UV/Ly$α$ 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α}\lesssim10^{42}$ erg/s, $M_{\star}\lesssim10^9$ M$_{\odot}$). As intrinsic Ly$α$ emission and internal Ly$α$ escape fractions barely evolve from $z=6$ to 9, the observed suppression of Ly$α$ luminosities with increasing redshift is fully attributed to IGM absorption. For most observable galaxies ($M_{\rm 1500}\lesssim-16$), the Ly$α$ 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$α$ 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$α$ 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\%$. △ Less

Submitted 7 April, 2021; originally announced April 2021.

Comments: MNRAS accepted, 20 pages, 19 figures, plus appendices

The MUSE Extremely Deep Field: the Cosmic Web in Emission at High Redshift

Authors: Roland Bacon, David Mary, Thibault Garel, Jeremy Blaizot, Michael Maseda, Joop Schaye, Lutz Wisotzki, Simon Conseil, Jarle Brinchmann, Floriane Leclercq, Valentina Abril-Melgarejo, Leindert Boogaard, Nicolas Bouché, Thierry Contini, Anna Feltre, Bruno Guiderdoni, Christian Herenz, Wolfram Kollatschny, Haruka Kusakabe, Jorryt Matthee, Léo Michel-Dansac, Themiya Nanayakkara, Johan Richard, Martin Roth, Kasper B. Schmidt , et al. (7 additional authors not shown)

Abstract: 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 ver… ▽ More 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 ($α\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. △ Less

Submitted 10 February, 2021; originally announced February 2021.

Comments: 28 pages, 19 figures

Journal ref: A&A 647, A107 (2021)

Introducing SPHINX-MHD: The Impact of Primordial Magnetic Fields on the First Galaxies, Reionization, and the Global 21cm Signal

Authors: Harley Katz, Sergio Martin-Alvarez, Joakim Rosdahl, Taysun Kimm, Jérémy Blaizot, Martin G. Haehnelt, Léo Michel-Dansac, Thibault Garel, Jose Oñorbe, Julien Devriendt, Adrianne Slyz, Omar Attia, Romain Teyssier

Abstract: 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 id… ▽ More 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 ($τ_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 $τ_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$. △ Less

Submitted 22 July, 2021; v1 submitted 27 January, 2021; originally announced January 2021.

Comments: 29 pages, 27 figures, accepted MNRAS

The Tully-Fisher relation in dense groups at $z \sim 0.7$ in the MAGIC survey

Authors: Valentina Abril-Melgarejo, Benoît Epinat, Wilfried Mercier, Thierry Contini, Leindert A. Boogaard, Jarle Brinchmann, Hayley Finley, Léo Michel-Dansac, Emmy Ventou, Philipe Amram, Davor Krajnović, Guillaume Mahler, Juan C. B. Pineda, Johan Richard

Abstract: 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… ▽ More 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<z<0.8$ and located in 10 pointings of the MAGIC MUSE Guaranteed Time Observations program. We perform a morpho-kinematics analysis of this sample and set up a selection based on galaxy size, [OII] emission line doublet signal-to-noise ratio, bulge-to-disk ratio and nuclear activity to construct a robust kinematic sample of 67 star-forming galaxies. This selection considerably reduces the number of outliers in the TFR, which are predominantly dispersion-dominated galaxies. Our results suggest a significant offset of the TFR zero-point between galaxies in low- and high-density environments, whatever kinematics estimator is used. This can be interpreted as a decrease of either stellar mass by $\sim 0.05 - 0.3$ dex or an increase of rotation velocity by $\sim 0.02 - 0.06$ dex for galaxies in groups, depending on the samples used for comparison. We also studied the stellar and baryon mass fractions within stellar disks and found they both increase with stellar mass, the trend being more pronounced for the stellar component alone. These fractions do not exceed 50%. We show that this evolution of the TFR is consistent either with a decrease of star formation or with a contraction of the mass distribution due to the environment. These two effects probably act together with their relative contribution depending on the mass regime. △ Less

Submitted 15 May, 2023; v1 submitted 20 January, 2021; originally announced January 2021.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 647, A152 (2021)

UV absorption lines and their potential for tracing the Lyman continuum escape fraction

Authors: V. Mauerhofer, A. Verhamme, J. Blaizot, T. Garel, T. Kimm, L. Michel-Dansac, J. Rosdahl

Abstract: The neutral intergalactic medium above redshift 6 is opaque to ionizing radiation, therefore one needs indirect measurements of the escape fraction of ionizing photons from galaxies of this epoch. Low-ionization state absorption lines are a common feature in the spectrum of galaxies, showing a diversity of strengths and shapes. Since these lines indicate the presence of neutral gas in front of the… ▽ More The neutral intergalactic medium above redshift 6 is opaque to ionizing radiation, therefore one needs indirect measurements of the escape fraction of ionizing photons from galaxies of this epoch. Low-ionization state absorption lines are a common feature in the spectrum of galaxies, showing a diversity of strengths and shapes. Since these lines indicate the presence of neutral gas in front of the stars, they have been proposed to carry information on the escape of ionizing radiation from galaxies. We study which processes are responsible for the shape of the absorption lines, to better understand their origin. We then explore whether the absorption lines can be used to predict the escape fractions. Using a radiation-hydrodynamical zoom-in simulation and the radiative transfer code RASCAS, we generate mock CII 1334 and LyB lines of a virtual galaxy at redshift 3 as seen from many directions of observation. We also compute the escape fraction of ionizing photons in those directions and look for correlations between the lines and the escape fractions. We find that the resulting mock absorption lines are comparable to observations and that the lines and the escape fractions vary strongly depending on the direction of observation. Gas velocity and dust always affect the absorption profile significantly. We find no strong correlations between observable LyB or CII 1334 and the escape fraction. After correcting the continuum for attenuation by dust to recover the intrinsic continuum, the residual flux of CII 1334 correlates well with the escape fraction for directions with a dust corrected residual flux larger than 30%. For other directions, the relations have a strong dispersion, and the residual flux overestimates the escape fraction for most cases. Concerning LyB, the residual flux after dust correction does not correlate with the escape fraction but can be used as a lower limit. (abridged) △ Less

Submitted 15 December, 2020; v1 submitted 7 December, 2020; originally announced December 2020.

Comments: 25 pages, 22 figures, accepted for publication in Astronomy and Astrophysics. New title, shorter abstract, corrected typos and url for Krome website

Journal ref: A&A 646, A80 (2021)

RASCAS: RAdiation SCattering in Astrophysical Simulations

Authors: Leo Michel-Dansac, Jeremy Blaizot, Thibault Garel, Anne Verhamme, Taysun Kimm, Maxime Trebitsch

Abstract: 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 sp… ▽ More 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. △ Less

Submitted 30 January, 2020; originally announced January 2020.

Comments: Accepted for publication in A&A. Code available at http://rascas.univ-lyon1.fr/

Journal ref: A&A 635, A154 (2020)

Probing the merger history of red early-type galaxies with their faint stellar substructures

Authors: B. Mancillas, P. -A. Duc, F. Combes, F. Bournaud, E. Emsellem, M. Martig, L. Michel-Dansac

Abstract: 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 perturbation… ▽ More 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. △ Less

Submitted 25 October, 2019; v1 submitted 16 September, 2019; originally announced September 2019.

Comments: 16 pages, 14 figures, accepted in A&A (v2 after language editing)

Journal ref: A&A 632, A122 (2019)

How to Quench a Dwarf Galaxy: The Impact of Inhomogeneous Reionization on Dwarf Galaxies and Cosmic Filaments

Authors: Harley Katz, Marius Ramsoy, Joakim Rosdahl, Taysun Kimm, Jeremy Blaizot, Martin G. Haehnelt, Leo Michel-Dansac, Thibault Garel, Clotilde Laigle, Julien Devriendt, Adrianne Slyz

Abstract: 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 t… ▽ More 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 $Δ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$... △ Less

Submitted 26 February, 2020; v1 submitted 27 May, 2019; originally announced May 2019.

Comments: 21 pages, 11 figures, MNRAS Accepted

Understanding the escape of LyC and Ly$α$ photons from turbulent clouds

Authors: Taysun Kimm, Jeremy Blaizot, Thibault Garel, Leo Michel-Dansac, Harley Katz, Joakim Rosdahl, Anne Verhamme, Martin Haehnelt

Abstract: 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,… ▽ More 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 ($τ_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. △ Less

Submitted 4 April, 2019; v1 submitted 17 January, 2019; originally announced January 2019.

Comments: 23 pages, 20 figures, accepted for publication in MNRAS

Properties and redshift evolution of star-forming galaxies with high [OIII]/[OII] ratios with MUSE at 0.28<z<0.85

Authors: M. Paalvast, A. Verhamme, L. A. Straka, J. Brinchmann, E. C. Herenz, D. Carton, M. L. P. Gunawardhana, L. A. Boogaard, S. Cantalupo, T. Contini, B. Epinat, H. Inami, R. A. Marino, M. V. Maseda, L. Michel-Dansac, S. Muzahid, T. Nanayakkara, G. Pezzulli, J. Richard, J. Schaye, M. C. Segers, T. Urrutia, M. Wendt, L. Wisotzki

Abstract: 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 sele… ▽ More 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. △ Less

Submitted 14 August, 2018; originally announced August 2018.

Comments: 17 pages, 23 figures, accepted for publication in A&A

Journal ref: A&A 618, A40 (2018)

The MUSE Hubble Ultra Deep Field Survey: XII. MgII emission and absorption in star-forming galaxies

Authors: Anna Feltre, Roland Bacon, Laurence Tresse, Hayley Finley, David Carton, Jérémy Blaizot, Nicolas Bouché, Thibault Garel, Hanae Inami, Leindert A. Boogaard, Jarle Brinchmann, Stéphane Charlot, Jacopo Chevallard, Thierry Contini, Leo Michel-Dansac, Guillaume Mahler, Raffaella A. Marino, Michael V. Maseda, Johan Richard, Kasper B. Schmidt, Anne Verhamme

Abstract: 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 MgI… ▽ More 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. △ Less

Submitted 5 June, 2018; originally announced June 2018.

Comments: 17 pages, 12 figures, 2 tables, accepted for publication in A&A

Journal ref: A&A 617, A62 (2018)

The SPHINX Cosmological Simulations of the First Billion Years: the Impact of Binary Stars on Reionization

Authors: Joakim Rosdahl, Harley Katz, Jeremy Blaizot, Taysun Kimm, Leo Michel-Dansac, Thibault Garel, Martin Haehnelt, Pierre Ocvirk, Romain Teyssier

Abstract: 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… ▽ More 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. △ Less

Submitted 25 June, 2018; v1 submitted 22 January, 2018; originally announced January 2018.

Comments: 18 pages, 15 figures. Accepted for publication in MNRAS

Impact of Lyman alpha pressure on metal-poor dwarf galaxies

Authors: Taysun Kimm, Martin Haehnelt, Jeremy Blaizot, Harley Katz, Leo Michel-Dansac, Thibault Garel, Joakim Rosdahl, Romain Teyssier

Abstract: 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-$α$ (Lya) photons is an important source of stellar feedback which… ▽ More 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-$α$ (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. △ Less

Submitted 15 January, 2018; originally announced January 2018.

Comments: 19 pages, 16 figures, Accepted for publication in MNRAS

The MUSE Hubble Ultra Deep Field Survey. IX. Evolution of galaxy merger fraction since z~6

Authors: E. Ventou, T. Contini, N. Bouché, B. Epinat, J. Brinchmann, R. Bacon, H. Inami, D. Lam, A. Drake, T. Garel, L. Michel-Dansac, R. Pello, M. Steinmetz, P. M. Weilbacher, L. Wisotzki, M. Carollo

Abstract: 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… ▽ More 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$α$ 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. △ Less

Submitted 1 November, 2017; originally announced November 2017.

Journal ref: A&A 608, A9 (2017)

Ionised gas structure of 100 kpc in an over-dense region of the galaxy group COSMOS-Gr30 at z ~ 0.7

Authors: Benoît Epinat, Thierry Contini, Hayley Finley, Leindert Boogaard, Adrien Guérou, Jarle Brinchmann, David Carton, Léo Michel-Dansac, Roland Bacon, Sebastiano Cantalupo, Marcella Carollo, Stephen Hamer, Wolfram Kollatschny, Davor Krajnović, Raffaella Anna Marino, Johan Richard, Geneviève Soucail, Peter M. Weilbacher, Lutz Wisotzki

Abstract: 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 al… ▽ More 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β) 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. △ Less

Submitted 21 November, 2017; v1 submitted 30 October, 2017; originally announced October 2017.

Comments: 21 pages, 12 figures, 3 tables, accepted for publication in A&A (27 October 2017)

Journal ref: A&A 609, A40 (2018)

The MUSE Hubble Ultra Deep Field Survey: V. Spatially resolved stellar kinematics of galaxies at redshift $0.2\lesssim z \lesssim 0.8$

Authors: Adrien Guerou, Davor Krajnovic, Benoit Epinat, Thierry Contini, Eric Emsellem, Nicolas Bouche, Roland Bacon, Leo Michel-Dansac, Johan Richard, Peter M. Weilbacher, Joop Schaye, Raffaella Anna Marino, Mark den Brok, Santiago Erroz-Ferrer

Abstract: 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 sel… ▽ More 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. △ Less

Submitted 20 October, 2017; originally announced October 2017.

Comments: 25 pages, 27 figures, A&A for accepted (26 June 2017)

The MUSE Hubble Ultra Deep Field Survey: II. Spectroscopic redshifts and comparisons to color selections of high-redshift galaxies

Authors: H. Inami, R. Bacon, J. Brinchmann, J. Richard, T. Contini, S. Conseil, S. Hamer, M. Akhlaghi, N. Bouche, B. Clement, G. Desprez, A. B. Drake, T. Hashimoto, F. Leclercq, M. Maseda, L. Michel-Dansac, M. Paalvast, L. Tresse, E. Ventou, W. Kollatschny, L. A. Boogaard, H. Finley, R. A. Marino, J. Schaye, L. Wisotzki

Abstract: 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 spectroscopi… ▽ More 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. △ Less

Submitted 10 October, 2017; originally announced October 2017.

Comments: 26 pages, 26 figures, 4 tables, accepted for publication in A&A (MUSE UDF Series Paper II), the redshift catalogs will be available at the CDS

Gas flows in the circumgalactic medium around simulated high-redshift galaxies

Authors: Peter Mitchell, Jeremy Blaizot, Julien Devriendt, Taysun Kimm, Leo Michel-Dansac, Joakim Rosdahl, Adrianne Slyz

Abstract: 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 impac… ▽ More 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. △ Less

Submitted 10 October, 2017; originally announced October 2017.

Comments: 24 pages, 17 figures, submitted to MNRAS after responding to a first referee report

The evolution of the metallicity gradient and the star formation efficiency in disc galaxies

Authors: Emanuel Sillero, Patricia B. Tissera, Diego G. Lambas, Leo Michel-Dansac

Abstract: 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 positi… ▽ More 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. △ Less

Submitted 1 September, 2017; originally announced September 2017.

Comments: 11 pages, 7 figures

Deep MUSE observations in the HDFS. Morpho-kinematics of distant star-forming galaxies down to $10^8$M$_\odot$

Authors: T. Contini, B. Epinat, N. Bouché, J. Brinchmann, L. A. Boogaard, E. Ventou, R. Bacon, J. Richard, P. M. Weilbacher, L. Wisotzki, D. Krajnović, J-B. Vielfaure, E. Emsellem, H. Finley, H. Inami, J. Schaye, M. Swinbank, A. Guérou, T. Martinsson, L. Michel-Dansac, I. Schroetter, M. Shirazi, G. Soucail

Abstract: (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 imag… ▽ More (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. △ Less

Submitted 29 September, 2016; v1 submitted 1 December, 2015; originally announced December 2015.

Comments: 26 pages. Published in A&A

Journal ref: A&A 591, A49 (2016)

The UV, Lyman-alpha, and dark matter halo properties of high redshift galaxies

Authors: T. Garel, J. Blaizot, B. Guiderdoni, L. Michel-Dansac, M. Hayes, A. Verhamme

Abstract: 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… ▽ More 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. △ Less

Submitted 23 March, 2015; originally announced March 2015.

Comments: accepted in MNRAS

GalPak3D: A Bayesian parametric tool for extracting morpho-kinematics of galaxies from 3D data

Authors: N. Bouché, H. Carfantan, I. Schroetter, L. Michel-Dansac, T. Contini

Abstract: 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,λ$ 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 algorith… ▽ More 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,λ$ 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. △ Less

Submitted 16 August, 2015; v1 submitted 26 January, 2015; originally announced January 2015.

Comments: 16 pages, 10 figures, accepted to publication in AJ, revised version after proofs corrections. Algorithm available at http://galpak.irap.omp.eu

Journal ref: 2015, AJ 150, 92

The MUSE 3D view of the Hubble Deep Field South

Authors: R. Bacon, J. Brinchmann, J. Richard, T. Contini, A. Drake, M. Franx, S. Tacchella, J. Vernet, L. Wisotzki, J. Blaizot, N. Bouché, R. Bouwens, S. Cantalupo, C. M. Carollo, D. Carton, J. Caruana, B. Clément, S. Dreizler, B. Epinat, B. Guiderdoni, C. Herenz, T. -O. Husser, S. Kamann, J. Kerutt, W. Kollatschny , et al. (13 additional authors not shown)

Abstract: 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 comb… ▽ More 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Å allows us to separate resolved asymmetric Lya emitters, [O II] emitters, and C III] emitters and the large instantaneous wavelength range of 4500Å 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. △ Less

Submitted 27 November, 2014; originally announced November 2014.

Comments: 21 pages, 20 figures. Note that the catalogues, data cube and associated spectra will be released upon acceptance of the paper

Journal ref: A&A 575, A75 (2015)

Evolution of the specific Star Formation Rate Function at z<1.4 - Dissecting the mass-SFR plane in COSMOS and GOODS

Authors: O. Ilbert, S. Arnouts, E. Le Floc'h, H. Aussel, H. Bethermin, P. Capak, B. C. Hsieh, M. Kajisawa, A. Karim, O. Le Fevre, N. Lee, S. Lilly, H. J. McCracken, L. Michel-Dansac, T. Moutard, M. A. Renzini, M. Salvato, D. B. Sanders, N. Scoville, K. Sheth, J. D. Silverman, V. Smolcic, Y. Taniguchi, L. Tresse

Abstract: 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$μm$ select… ▽ More 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$μ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. △ Less

Submitted 30 March, 2015; v1 submitted 17 October, 2014; originally announced October 2014.

Comments: 24 pages, 16 figures, 3 tables, published version in A&A

The Atlas3D project -- XXIX. The new look of early-type galaxies and surrounding fields disclosed by extremely deep optical images

Authors: Pierre-Alain Duc, Jean-Charles Cuillandre, Emin Karabal, Michele Cappellari, Katherine Alatalo, Leo Blitz, Frederic Bournaud, Martin Bureau, Alison F. Crocker, Roger L. Davies, Timothy A. Davis, P. T. de Zeeuw, Eric Emsellem, Sadegh Khochfar, Davor Krajnovic, Harald Kuntschner, Richard M. McDermid, Leo Michel-Dansac, Raffaella Morganti, Thorsten Naab, Tom Oosterloo, Sanjaya Paudel, Marc Sarzi, Nicholas Scott, Paolo Serra , et al. (2 additional authors not shown)

Abstract: 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 proje… ▽ More 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. △ Less

Submitted 3 October, 2014; originally announced October 2014.

Comments: 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

Chemodynamics of a Simulated Disc Galaxy: Initial Mass Functions and SNIa Progenitors

Authors: C. G. Few, S. Courty, B. K. Gibson, L. Michel-Dansac, F. Calura

Abstract: 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, hydro… ▽ More 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. △ Less

Submitted 21 August, 2014; originally announced August 2014.

Comments: 20 pages, 9 figures, accepted for publication in MNRAS

Clumpy Disc and Bulge Formation

Authors: J. Perez, O. Valenzuela, P. Tissera, L. Michel-Dansac

Abstract: 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… ▽ More 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. △ Less

Submitted 20 August, 2013; originally announced August 2013.

Comments: Accepted for publication in MNRAS - Aug. 20, 2013

The Next Generation Virgo Cluster Survey. IV. NGC 4216: A Bombarded Spiral in the Virgo Cluster

Authors: Sanjaya Paudel, Pierre-Alain Duc, Patrick Cote, Jean-Charles Cuillandre, Laura Ferrarese, Etienne Ferriere, Stephen D. J. Gwyn, J. Christopher Mihos, Bernd Vollmer, Michael L. Balogh, Ray G. Carlberg, Samuel Boissier, Alessandro Boselli, Patrick R. Durrell, Eric Emsellem, Lauren A. MacArthur, Simona Mei, Leo Michel-Dansac, Wim van Driel

Abstract: 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… ▽ More 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. △ Less

Submitted 26 February, 2013; originally announced February 2013.

Comments: 11 Pages, 7 Figures, 3 Tables, Accepted for publication in ApJ

Properties of simulated Milky Way-mass galaxies in loose group and field environments

Authors: C. G. Few, B. K. Gibson, S. Courty, L. Michel-Dansac, C. B. Brook, G. S. Stinson

Abstract: 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 fr… ▽ More 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. △ Less

Submitted 3 October, 2012; originally announced October 2012.

Comments: 16 pages, 11 figures, abstract abridged for arXiv. Accepted for publication in Astronomy & Astrophysics

Discovery of a giant HI tail in the galaxy group HCG 44

Authors: Paolo Serra, Baerbel Koribalski, Pierre-Alain Duc, Tom Oosterloo, Richard M. McDermid, Leo Michel-Dansac, Eric Emsellem, Jean-Charles Cuillandre, Katherine Alatalo, Leo Blitz, Maxime Bois, Frederic Bournaud, Martin Bureau, Michele Cappellari, Alison F. Crocker, Roger L. Davies, Timothy A. Davis, P. T. de Zeeuw, Sadegh Khochfar, Davor Krajnovic, Harald Kuntschner, Pierre-Yves Lablanche, Raffaella Morganti, Thorsten Naab, Marc Sarzi , et al. (3 additional authors not shown)

Abstract: 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 HC… ▽ More 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. △ Less

Submitted 18 September, 2012; originally announced September 2012.

Comments: MNRAS, accepted

The Atlas3D project - XII. Recovery of the mass-to-light ratio of simulated early-type barred galaxies with axisymmetric dynamical models

Authors: Pierre-Yves Lablanche, Michele Cappellari, Eric Emsellem, Frederic Bournaud, Leo Michel-Dansac, Katherine Alatalo, Leo Blitz, Maxime Bois, Martin Bureau, Roger L. Davies, Timothy A. Davis, P. T. de Zeeuw, Pierre-Alain Duc, Sadegh Khochfar, Davor Krajnovic, Harald Kuntschner, Raffaella Morganti, Richard M. McDermid, Thorsten Naab, Tom Oosterloo, Marc Sarzi, Nicholas Scott, Paolo Serra, Anne-Marie Weijmans, Lisa M. Young

Abstract: 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… ▽ More 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] △ Less

Submitted 1 June, 2012; originally announced June 2012.

Comments: 29 pages, 29 figures, LaTeX. Accepted for publication in MNRAS

The Distribution of Metals in Cosmological Hydrodynamical Simulations of Dwarf Disk Galaxies

Authors: K. Pilkington, B. K. Gibson, C. B. Brook, F. Calura, G. S. Stinson, R. J. Thacker, L. Michel-Dansac, J. Bailin, H. M. P. Couchman, J. Wadsley, T. R. Quinn, A. Maccio

Abstract: 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 --… ▽ More 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. △ Less

Submitted 21 May, 2012; originally announced May 2012.

Comments: Accepted for publication in MNRAS; 12 pages; 7 figures