Researcher profile

Johan Richard

Johan Richard contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
29works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

29 published item(s)

preprint2026arXiv

Caught in Swallowtails: Discovery of Two Swallowtail Image Formations in MS 0451.6-0305

We report the discovery of two swallowtail image formations at $z=2.91$ and $z=6.70$ behind the galaxy cluster MS 0451.6-0305 in JWST-NIRCam imaging. We find that in both of the above lensed systems, the complex image morphology cannot be reproduced by simple fold/cusp caustics, and detailed lens modeling reveals higher-order swallowtail caustic configurations. In the $z=2.91$ lens system, a small part of the source galaxy (which itself is part of a galaxy group) containing atleast two compact knots sits inside the swallowtail caustic, producing a quadruply imaged arc. At two of the image positions of these knots, we infer point source magnifications of $\gtrsim 300$, implying lensing-corrected effective radii of $\lesssim 0.8-1.5$ pc. The $z=6.70$ system exhibits even more complex image formation. We therefore only use the most confidently identified counter-images of knots in this system as constraints in our lens modeling. The resulting model predicts magnifications $\sim20-200$ and lensing-corrected effective radii of $\lesssim 0.8-18.5$ pc for various knots. Together, these two systems represent the first example of observations of multiple swallowtail image formations in a single galaxy cluster and demonstrate the ability of swallowtail caustics to magnify individual substructures at sub-parsec scales, from intermediate redshifts to the first billion years of the Universe.

preprint2026arXiv

Spectroscopic characterization of a remarkable temporally varying, triple-lensed quasar at z=2.67

Gravitationally lensed quasars are viable cosmic tools for constraining a diversity of fundamental astrophysical phenomena; They enable identification of faint, low-mass supermassive black holes, provide unique constraints on the intervening intergalactic or interstellar medium in their sightlines, and can be used to determine key cosmological quantities such as the Hubble constant, $H_0$. However, they are rare phenomena, and it has proven difficult to define efficient, unbiased selection methods.} In this study, we report the independent spectroscopic identification of a remarkable triple-lensed quasar at $z=2.67$, identified based on astrometric measurements from the {\em Gaia} mission, previously identified in Pan-STARRS. Furthermore, a larger spectroscopic follow-up survey of {\em Gaia}-detected candidate lensed quasars. We characterize in detail the three mirror images of the quasar and their spatial and temporal spectroscopic coverage, with focus on the emission-line properties which shows variation across sigthlines and temporal evolution over the $\sim 11$months spectroscopic campaign. We construct a lens model of the foreground source from a combination of the multiple spectra and deep optical imaging, providing a robust halo mass of $M_{\rm h} = (2.78 \pm 0.05)\times 10^{10}M_\odot$. Based on the lens model, the time delay between each sightline is translated into an intrinsic quasar time, allowing us to construct a quasar timeseries over $\sim18$months with monthly cadence. Over months timescales the broad emission lines vary in both velocity offset and equivalent width (EW) as well as an overall increase in ionization. This exemplary triple-lensed quasars demonstrates the viability of identifying such rare lens configurations based purely on the astrometric measurements from the {\em Gaia} mission, which we here provide optimized selection criteria for, for future studies.

preprint2026arXiv

VENUS: Two Faint Little Red Dots Separated by $\sim70\,\mathrm{pc}$ Hidden in a Single Lensed Galaxy at $z\sim7$

We report the identification of a pair of faint little red dots (LRDs), dubbed Red Eyes, in a strongly-lensed galaxy at $z\sim7$ behind the PLCKG004.5-10.5 cluster, identified from the JWST Treasury program VENUS. Red Eyes are spatially resolved on the image plane with distinct colors, while the critical curve lies far north of Red Eyes, clearly requiring two different LRDs rather than a single LRD. Red Eyes is an extremely close pair of LRDs separated by $\sim70\,\mathrm{pc}$ in the source plane with a magnification of $μ\sim20$, which consistently explains another counter-image detected to the north-west. Red Eyes is hosted in a typical star-forming galaxy with $M_{\mathrm{UV,int}}\sim -19$, but its own UV emission is very faint ($M_{\mathrm{UV,int}} \gtrsim -16$). Moreover, Red Eyes does not reside at the galaxy center but lies at an offset position of approximately one effective radius $R_{\mathrm{e}}$ away from the galaxy center. If observed without lensing, Red Eyes would appear as a typical star-forming galaxy at $z\sim 7$ with $M_{\mathrm{UV}}\sim -19$, showing no apparent LRD signatures in either morphology or SED. These results suggest that multiple off-center LRDs, similar to Red Eyes, may be commonly hidden in a typical high-$z$ star-forming galaxy. In this case, various plausible scenarios may emerge, one of which is that intermediate-mass black holes (IMBHs) with $M_\mathrm{BH}\sim10^{4\text{--}6}\,M_\odot$ may form in star clusters on a stellar disk and contribute to the growth of the central supermassive black hole via mergers, with some IMBHs detectable as luminous LRDs in a sufficiently active and massive phase.

preprint2023arXiv

Constraining the physical properties of the first lensed $z\sim10-16$ galaxy candidates with JWST

The first deep-field observations of the JWST have immediately yielded a surprisingly large number of very high redshift candidates, pushing the frontier of observability well beyond $z\gtrsim10$. We here present a detailed SED-fitting analysis of the 15 gravitationally lensed $z\sim10-16$ galaxy candidates detected behind the galaxy cluster SMACS J0723.3-7327 in Atek et al. (2022) using the BEAGLE tool. Our analysis makes use of dynamical considerations to place limits on the ages of these galaxies and of all three published SL models of the cluster to account for lensing systematics. We find these galaxies to have relatively low stellar masses $M_{\star}\sim10^7-10^8\,\mathrm{M}_{\odot}$ and young ages $t_{\mathrm{age}}\sim10-100$\,Myr. Due to their very blue UV-slopes, down to $β\sim-3$, all of the galaxies in our sample have extremely low dust attenuations $A_V\lesssim0.02$. Placing the measured parameters into relation, we find a very shallow $M_{\star}-M_{\mathrm{UV}}$-slope and high sSFRs above the main sequence of star-formation with no significant redshift-evolution in either relation. This is in agreement with the bright UV luminosities measured for these objects and indicates that we are naturally selecting galaxies that are currently undergoing a star-bursting episode at the time they are observed. Finally, we discuss the robustness of our high-redshift galaxy sample regarding low-redshift interlopers and conclude that low-redshift solutions can safely be ruled out for roughly half of the sample, including the highest-redshift galaxies at $z\sim12-16$. These objects represent compelling targets for spectroscopic follow-up observations with JWST and ALMA.

preprint2022arXiv

ALMA reveals extended cool gas and hot ionized outflows in a typical star-forming galaxy at $z=7.13$

We present spatially-resolved morphological properties of [CII] 158 $μ$m, [OIII] 88 $μ$m, dust, and rest-frame ultraviolet (UV) continuum emission for A1689-zD1, a strongly lensed, sub-L* galaxy at $z=7.13$, by utilizing deep Atacama Large Millimeter/submillimeter Array (ALMA) and Hubble Space Telescope (HST) observations. While the [OIII] line and UV continuum are compact, the [CII] line is extended up to a radius of $r \sim 12$ kpc. Using multi-band rest-frame far-infrared (FIR) continuum data ranging from 52-400 $μ$m, we find an average dust temperature and emissivity index of $T_{\rm dust} = 41^{+17}_{-14}$ K and $β= 1.7^{+1.1}_{-0.7}$, respectively, across the galaxy. We find slight differences in the dust continuum profiles at different wavelengths, which may indicate that the dust temperature decreases with distance. We map the star-formation rate (SFR) via IR and UV luminosities and determine a total SFR of $37\pm 1~M_\odot~{\rm yr}^{-1}$ with an obscured fraction of $87\%$. While the [OIII] line is a good tracer of the SFR, the [CII] line shows deviation from the local $L_{\rm [CII]}$-SFR relations in the outskirts of the galaxy. Finally, we observe a clear difference in the line profile between [CII] and [OIII], with significant residuals ($\sim 5σ$) in the [OIII] line spectrum after subtracting a single Gaussian model. This suggests a possible origin of the extended [CII] structure from the cooling of hot ionized outflows. The extended [CII] and high-velocity [OIII] emission may both contribute in part to the high $L_{\rm [OIII]}$/$L_{\rm [CII]}$ ratios recently reported in $z>6$ galaxies.

preprint2022arXiv

Multiply lensed star forming clumps in the A521-sys1 galaxy at redshift 1

We study the population of star-forming clumps in A521-sys1, a $\rm z=1.04$ system gravitationally lensed by the foreground ($\rm z=0.25$) cluster Abell 0521. The galaxy presents one complete counter--image with a mean magnification of $\rm μ\sim4$ and a wide arc containing two partial images of A521-sys1 with magnifications reaching $\rm μ>20$, allowing the investigations of clumps down to scales of $\rm R_{eff}<50$ pc. We identify 18 unique clumps with a total of 45 multiple images. Intrinsic sizes and UV magnitudes reveal clumps with elevated surface brightnesses, comparable to similar systems at redshifts $\rm z\gtrsim1.0$. Such clumps account for $\sim40\%$ of the galaxy UV luminosity, implying that a significant fraction of the recent star-formation activity is taking place there. Clump masses range from $\rm 10^6\ M_\odot$ to $\rm 10^9\ M_\odot$ and sizes from tens to hundreds of parsec, resulting in mass surface densities from $10$ to $\rm 10^3\ M_\odot\ pc^{-2}$, with a median of $\rm \sim10^2\ M_\odot\ pc^{-2}$. These properties suggest that we detect star formation taking place across a wide range of scale, from cluster aggregates to giant star-forming complexes. We find ages of less than $100$ Myr, consistent with clumps being observed close to their natal region. The lack of galactocentric trends with mass, mass density, or age and the lack of old migrated clumps can be explained either by dissolution of clumps after few $\sim100$ Myr or by stellar evolution making them fall below the detectability limits of our data.

preprint2022arXiv

Pilot-WINGS: An extended MUSE view of the structure of Abell 370

We investigate the strong-lensing cluster Abell 370 (A370) using a wide Integral Field Unit (IFU) spectroscopic mosaic from the Multi-Unit Spectroscopic Explorer (MUSE). IFU spectroscopy provides significant insight into the structure and mass content of galaxy clusters, yet IFU-based cluster studies focus almost exclusively on the central Einstein-radius region. Covering over 14 arcmin$^2$, the new MUSE mosaic extends significantly beyond the A370 Einstein radius, providing, for the first time, a detailed look at the cluster outskirts. Combining these data with wide-field, multi-band Hubble Space Telescope (HST) imaging from the BUFFALO project, we analyse the distribution of objects within the cluster and along the line of sight. Identifying 416 cluster galaxies, we use kinematics to trace the radial mass profile of the halo, providing a mass estimate independent from the lens model. We also measure radially-averaged properties of the cluster members, tracking their evolution as a function of infall. Thanks to the high spatial resolution of our data, we identify six cluster members acting as galaxy-galaxy lenses, which constrain localized mass distributions beyond the Einstein radius. Finally, taking advantage of MUSE&#39;s 3D capabilities, we detect and analyse multiple spatially extended overdensities outside of the cluster that influence lensing-derived halo mass estimates. We stress that much of this work is only possible thanks to the robust, extended IFU coverage, highlighting its importance even in less optically dense cluster regions. Overall, this work showcases the power of combining HST+MUSE, and serves as the initial step towards a larger and wider program targeting several clusters.

preprint2022arXiv

Post-starburst galaxies in the centers of intermediate redshift clusters

We present results from MUSE spatially-resolved spectroscopy of 21 post-starburst galaxies in the centers of 8 clusters from $z\sim0.3$ to $z\sim0.4$. We measure spatially resolved star-formation histories (SFHs), the time since quenching ($t_Q$) and the fraction of stellar mass assembled in the past 1.5 Gyr ($μ_{1.5}$). The SFHs display a clear enhancement of star-formation prior to quenching for 16 out of 21 objects, with at least 10% (and up to $>50$%) of the stellar mass being assembled in the past 1.5 Gyr and $t_Q$ ranging from less than 100 Myrs to $\sim800$ Myrs. By mapping $t_Q$ and $μ_{1.5}$, we analyze the quenching patterns of the galaxies. Most galaxies in our sample have quenched their star-formation from the outside-in or show a side-to-side/irregular pattern, both consistent with quenching by ram-pressure stripping. Only three objects show an inside-out quenching pattern, all of which are at the high-mass end of our sample. At least two of them currently host an active galactic nucleus. In two post-starbursts, we identify tails of ionized gas indicating that these objects had their gas stripped by ram pressure very recently. Post-starburst features are also found in the stripped regions of galaxies undergoing ram-pressure stripping in the same clusters, confirming the link between these classes of objects. Our results point to ram-pressure stripping as the main driver of fast quenching in these environments, with active galactic nuclei playing a role at high stellar masses.

preprint2022arXiv

Star formation at the smallest scales; A JWST study of the clump populations in SMACS0723

We present the clump populations detected in 18 lensed galaxies at redshifts 1 to 8.5 within the lensing cluster field SMACS0723. The recent JWST Early Release Observations of this poorly known region of the sky have revealed numerous point-like sources within and surrounding their host galaxies, undetected in the shallower HST images. We use JWST multiband photometry and the lensing model of this galaxy cluster to estimate the intrinsic sizes and magnitudes of the stellar clumps. We derive optical restframe effective radii from $<$10 to 100s pc and masses ranging from $\sim10^5$ to $10^9$ M$_{\odot}$, overlapping with massive star clusters in the local universe. The ages range from 1 Myr to 1 Gyr. We compare the crossing time to the age of the clumps and determine that between 45 and 60 % of the detected clumps are consistent with being gravitationally bound. The lack of Gyr old clumps suggest that the dissolution time scales are shorter than 1 Gyr. We see a significant increase in the luminosity (mass) surface density of the clumps with redshift. Clumps in galaxies at the reionisation era have stellar densities higher than massive clusters in the local universe. We zoom-in into single galaxies at redshift $<6$ and find for two galaxies, the Sparkler and the Firework, that their star clusters/clumps show distinctive colour distributions and location surrounding their host galaxy that are compatible with being accredited or formed during merger events. The ages of some of the compact clusters are between 1 and 4 Gyr, e.g., globular cluster precursors formed around 9-12 Gyr ago. Our study, conducted on a small sample of galaxies, shows the potential of JWST observations for understanding the conditions under which star clusters form in rapidly evolving galaxies.

preprint2022arXiv

The BlueMUSE data reduction pipeline: lessons learned from MUSE and first design choices

BlueMUSE is an integral field spectrograph in an early development stage for the ESO VLT. For our design of the data reduction software for this instrument, we are first reviewing capabilities and issues of the pipeline of the existing MUSE instrument. MUSE has been in operation at the VLT since 2014 and led to discoveries published in more than 600 refereed scientific papers. While BlueMUSE and MUSE have many common properties we briefly point out a few key differences between both instruments. We outline a first version of the flowchart for the science reduction, and discuss the necessary changes due to the blue wavelength range covered by BlueMUSE. We also detail specific new features, for example, how the pipeline and subsequent analysis will benefit from improved handling of the data covariance, and a more integrated approach to the line-spread function, as well as improvements regarding the wavelength calibration which is of extra importance in the blue optical range. We finally discuss how simulations of BlueMUSE datacubes are being implemented and how they will be used to prepare the science of the instrument.

preprint2022arXiv

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

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.

preprint2022arXiv

The MUSE eXtremely Deep Field: Individual detections of Lyα haloes around rest-frame UV-selected galaxies at z~2.9-4.4

Hydrogen Ly$α$ haloes (LAHs) are commonly used as a tracer of the circumgalactic medium (CGM) at high redshifts. In this work, we aim to explore the existence of Ly$α$ haloes around individual UV-selected galaxies, rather than around Ly$α$ emitters (LAEs), at high redshifts. Our sample was continuum-selected with F775W<=27.5, and spectroscopic redshifts were assigned or constrained for all the sources thanks to the deepest (100- to 140-hour) existing Very Large Telescope (VLT)/Multi-Unit Spectroscopic Explorer (MUSE) data with adaptive optics. The final sample includes 21 galaxies that are purely F775W-magnitude selected within the redshift range z=2.9-4.4 and within a UV magnitude range -20<=M1500<= -18, thus avoiding any bias toward LAEs. We tested whether galaxy&#39;s Ly$α$ emission is significantly more extended than the MUSE PSF-convolved continuum component. We find 17 LAHs and four non-LAHs. We report the first individual detections of extended Ly$α$ emission around non-LAEs. The Ly$α$ halo fraction is thus as high as $81.0^{+10.3}_{-11.2}$%, which is close to that for LAEs at z=3-6 in the literature. This implies that UV-selected galaxies generally have a large amount of hydrogen in their CGM. We derived the mean surface brightness (SB) profile for our LAHs with cosmic dimming corrections and find that Ly$α$ emission extends to 5.4 arcsec (~40 physical kpc at the midpoint redshift z=3.6) above the typical 1$σ$ SB limit. The incidence rate of surrounding gas detected in Ly$α$ per one-dimensional line of sight per unit redshift, dn/dz, is estimated to be $0.76^{+0.09}_{-0.09}$ for galaxies with M1500<= -18 mag at z~3.7. Assuming that Ly$α$ emission and absorption arise in the same gas, this suggests, based on abundance matching, that LAHs trace the same gas as damped Ly$α$ systems (DLAs) and sub-DLAs.

preprint2022arXiv

The MUSE Hubble Ultra Deep Field Survey XVI. The angular momentum of low-mass star-forming galaxies. A cautionary tale and insights from TNG50

We investigate the specific angular momentum (sAM) $ j(<r)$ profiles of intermediate redshift ($0.4<z<1.4$) star-forming galaxies (SFGs) in the relatively unexplored regime of low masses (down to $M_\star\sim 10^8$M$_{\odot}$), and small sizes (down to $R_{\rm e}\sim 1.5$ kpc) and characterize the sAM scaling relation and its redshift evolution. We have developed a 3D methodology to constrain sAM profiles of the star-forming gas using a forward modeling approach with \galpak{} that incorporates the effects of beam smearing, yielding the intrinsic morpho-kinematic properties even with limited spatial resolution data. Using mock observations from the TNG50 simulation, we find that our 3D methodology robustly recovers the star formation rate (SFR)-weighted $j(<r)$ profiles down to low effective signal-to-noise ratio (SNR) of $\gtrapprox3$. We applied our methodology blindly to a sample of 494 \OII{}-selected SFGs in the MUSE Ultra Deep Field (UDF) 9~arcmin$^2$ mosaic data, covering the unexplored $8<\log M_*/$M$_{\odot}<9$ mass range. We find that the (SFR-weighted) sAM relation follows $j\propto M_\star^α$ with an index $α$ varying from $α=0.3$ to $α=0.5$, from $\log M_\star/$M$_{\odot}=8$ to $\log M_*/$M$_{\odot}=10.5$. The UDF sample supports a redshift evolution consistent with the $(1+z)^{-0.5}$ expectation from a Universe in expansion. The scatter of the sAM sequence is a strong function of the dynamical state with $\log j|_{M_*}\propto 0.65 \times \log(V_{\rm max}/σ)$ where $σ$ is the velocity dispersion at $2 R_{\rm e}$. In TNG50, SFGs also form a $j-M_{\star}-(V/σ)$ plane but it correlates more with galaxy size than with morphological parameters. Our results suggest that SFGs might experience a dynamical transformation before their morphological transformation to becoming passive via either merging or secular evolution.

preprint2021arXiv

ALMA 1.3 mm Survey of Lensed Submillimeter Galaxies (SMGs) Selected by Herschel: Discovery of Spatially Extended SMGs and Implications

We present an ALMA 1.3 mm (Band 6) continuum survey of lensed submillimeter galaxies (SMGs) at $z=1.0\sim3.2$ with an angular resolution of $\sim0.2$&#34;. These galaxies were uncovered by the Herschel Lensing Survey (HLS), and feature exceptionally bright far-infrared continuum emission ($S_\mathrm{peak} \gtrsim 90$ mJy) owing to their lensing magnification. We detect 29 sources in 20 fields of massive galaxy clusters with ALMA. Using both the Spitzer/IRAC (3.6/4.5 $\mathrm{μm}$) and ALMA data, we have successfully modeled the surface brightness profiles of 26 sources in the rest-frame near- and far-infrared. Similar to previous studies, we find the median dust-to-stellar continuum size ratio to be small ($R_\mathrm{e,dust}/R_\mathrm{e,star} = 0.38\pm0.14$) for the observed SMGs, indicating that star formation is centrally concentrated. This is, however, not the case for two spatially extended main-sequence SMGs with a low surface brightness at 1.3 mm ($\lesssim 0.1$ mJy arcsec$^{-2}$), in which the star formation is distributed over the entire galaxy ($R_\mathrm{e,dust}/R_\mathrm{e,star}>1$). As a whole, our SMG sample shows a tight anti-correlation between ($R_\mathrm{e,dust}/R_\mathrm{e,star}$) and far-infrared surface brightness ($Σ_\mathrm{IR}$) over a factor of $\simeq$ 1000 in $Σ_\mathrm{IR}$. This indicates that SMGs with less vigorous star formation (i.e., lower $Σ_\mathrm{IR}$) lack central starburst and are likely to retain a broader spatial distribution of star formation over the whole galaxies (i.e., larger $R_\mathrm{e,dust}/R_\mathrm{e,star}$). The same trend can be reproduced with cosmological simulations as a result of central starburst and potentially subsequent &#34;inside-out&#34; quenching, which likely accounts for the emergence of compact quiescent galaxies at $z\sim2$.

preprint2021arXiv

An Atlas of MUSE Observations towards Twelve Massive Lensing Clusters

Spectroscopic surveys of massive galaxy clusters reveal the properties of faint background galaxies, thanks to the magnification provided by strong gravitational lensing. We present a systematic analysis of integral-field-spectroscopy observations of 12 massive clusters, conducted with the Multi Unit Spectroscopic Explorer (MUSE). All data were taken under very good seeing conditions (0.6&#34;) in effective exposure times between two and 15 hrs per pointing, for a total of 125 hrs. Our observations cover a total solid angle of ~23 arcmin$^2$ in the direction of clusters, many of which were previously studied by the MACS, Frontier Fields, GLASS and CLASH programs. The achieved emission line detection limit at 5$σ$ for a point source varies between (0.77--1.5)$\times$10$^{-18}$ erg\,s$^{-1}$\,cm$^{-2}$ at 7000Å. We present our developed strategy to reduce these observational data, detect sources and determine their redshifts. We construct robust mass models for each cluster to further confirm our redshift measurements using strong-lensing constraints, and identify a total of 312 strongly lensed sources producing 939 multiple images. The final redshift catalogs contain more than 3300 robust redshifts, of which 40\% are for cluster members and $\sim$30\% for lensed Lyman-$α$ emitters. 14\% of all sources are line emitters not seen in the available HST images, even at the depth of the FFs ($\sim29$ AB). We find that the magnification distribution of the lensed sources in the high-magnification regime ($μ{=}$ 2--25) follows the theoretical expectation of $N(z)\proptoμ^{-2}$. The quality of this dataset, number of lensed sources, and number of strong-lensing constraints enables detailed studies of the physical properties of both the lensing cluster and the background galaxies. The full data products from this work are made available to the community. [abridged]

preprint2021arXiv

Curved detector-based optical design for the VLT/BlueMUSE instrument

BlueMUSE (Blue Multi Unit Spectroscopic Explorer) is a blue-optimised, medium spectral resolution, panoramic integral field spectrograph proposed for the Very Large Telescope (VLT) and based on the MUSE concept. BlueMUSE will open up a new range of galactic and extragalactic science cases allowed by its specific capabilities in the 350 - 580 nm range: an optimised end-to-end transmission down to 350 nm, a larger FoV (up to $1.4 \times 1.4$ arcmin$^2$) sampled at 0.3 arcsec, and a higher spectral resolution ($λ/Δλ\sim 3500$) compared to MUSE. To our knowledge, achieving such capabilities with a comparable mechanical footprint and an identical detector format ($4\text{k} \times 4\text{k}$, 15 $\mathrm{μm}$ CCD) would not be possible with a conventional spectrograph design. In this paper, we present the optomechanical architecture and design of BlueMUSE at pre-phase A level, with a particular attention to some original aspects such as the use of curved detectors.

preprint2021arXiv

Hierarchical fragmentation in high redshift galaxies revealed by hydrodynamical simulations

High-redshift star-forming galaxies have very different morphologies compared to nearby ones. Indeed, they are often dominated by bright star-forming structures of masses up to $10^{8-9}$ $\mathrm{M}_\odot$ dubbed «giant clumps». However, recent observations questioned this result by showing only low-mass structures or no structure at all. We use Adaptative Mesh Refinement hydrodynamical simulations of galaxies with parsec-scale resolution to study the formation of structures inside clumpy high-redshift galaxies. We show that in very gas-rich galaxies star formation occurs in small gas clusters with masses below $10^{7-8}$ $\mathrm{M}_\odot$ that are themselves located inside giant complexes with masses up to $10^8$ and sometimes $10^9$ $\mathrm{M}_\odot$. Those massive structures are similar in mass and size to the giant clumps observed in imaging surveys, in particular with the Hubble Space Telescope. Using mock observations of simulated galaxies, we show that at very high resolution with instruments like the Atacama Large Millimeter Array or through gravitational lensing, only low-mass structures are likely to be detected, and their gathering into giant complexes might be missed. This leads to the non-detection of the giant clumps and therefore introduces a bias in the detection of these structures. We show that the simulated giant clumps can be gravitationally bound even when undetected in mocks representative for ALMA observations and HST observations of lensed galaxies. We then compare the top-down fragmentation of an initially warm disc and the bottom-up fragmentation of an initially cold disc to show that the process of formation of the clumps does not impact their physical properties.

preprint2021arXiv

Radial profiles of lensed $z \sim 1$ galaxies on sub-kiloparsec scales

We study the spatially resolved physical properties of the Cosmic Snake arc in MACS J1206.2-0847 and the arc in Abell 0521 (A521). These are two strongly lensed galaxies at redshifts $z=1.036$ and $z=1.044$. We used observations of the Hubble Space Telescope (HST) and the Atacama Large Millimeter/submillimeter Array (ALMA). The former gives access to the star formation rate (SFR) and stellar mass ($M_\star$), and the latter to the H$_2$ molecular gas mass ($M_{\mathrm{mol}}$). HST and ALMA observations have similar angular resolutions of $0.15^{\prime \prime}-0.2^{\prime \prime}$, which with the help of strong gravitational lensing enable us to reach spatial resolutions down to $\sim 30\,\mathrm{pc}$ and $\sim 100\,\mathrm{pc}$ in these two galaxies, respectively. These resolutions are close to the resolution of observations of nearby galaxies. We study the radial profiles of SFR, $M_\star$, and $M_{\mathrm{mol}}$ surface densities of these high-redshift galaxies and compare the corresponding exponential scale lengths with those of local galaxies. We find that the scale lengths in the Cosmic Snake are about $0.5\,\mathrm{kpc}-1.5\,\mathrm{kpc}$, and they are 3 to 10 times larger in A521. This is a significant difference knowing that the two galaxies have comparable integrated properties. These high-redshift scale lengths are nevertheless comparable to those of local galaxies, which cover a wide distribution. The particularity of our high-redshift radial profiles is the normalisation of the $M_{\mathrm{mol}}$ surface density profiles ($ΣM_{\mathrm{mol}}$), which are offset by up to a factor of 20 with respect to the profiles of $z=0$ counterparts. The SFR surface density profiles are also offset by the same factor as $ΣM_{\mathrm{mol}}$, as expected from the Kennicutt-Schmidt law.

preprint2021arXiv

The MUSE Extremely Deep Field: Evidence for SFR-induced cores in dark-matter dominated galaxies at z=1

Disc-halo decomposition on rotationally supported star-forming galaxies (SFGs) at $z>1$ are often limited to massive galaxies ($M_\star>10^{10}~M_\odot$) and rely on either deep Integral Field Spectroscopy data or stacking analyses. We present a study of the dark matter (DM) content of nine $z\approx1$ SFGs selected Using the brightest [OII] emitters in the deepest Multi-Unit Spectrograph Explorer (MUSE) field to date, namely the 140hr MUSE Extremely Deep Field, we perform disk-halo decompositions on 9 low-mass SFGs (with $10^{8.5}<M_\star<10^{10.5}~M_\odot$) using a novel 3D modeling approach, which together with the exquisite S/N allows us to measure individual rotation curves to $3\times R_e$. The disk-halo decomposition includes a stellar, DM, gas, and occasionally a bulge component. The DM component primarily uses the generalized $α,β,γ$ profile or a Navarro-Frenk-White (NFW) profile. The disk stellar masses $M_\star$ obtained from the [OII] disk-halo decomposition agree with the values inferred from the spectral energy distributions. While the rotation curves show diverse shapes, ranging from rising to declining at large radii, the DM fractions within the half-light radius $f_{\rm DM}(<R_e)$ are found to be 60\% to 95\%, extending to lower masses (densities) recent results on massive SFGs with $M_\star>10^{10}~M_\odot$. The DM halos show constant surface densities of $\sim100~M_\odot$ pc$^{-2}$. Half of the sample shows a strong preference for cored over cuspy DM profiles. The presence of DM cores appears to be related to galaxies with stellar-to-halo mass $\log M_\star/M_{\rm vir}\approx-2.5$. In addition, the cuspiness of the DM profiles is found to be a strong function of the recent star-formation activity. Both of these results are interpreted as evidence for feedback-induced core formation in the Cold Dark Matter context.

preprint2020arXiv

Elevated ionizing photon production efficiency in faint high-equivalent-width Lyman-alpha emitters

While low-luminosity galaxies dominate number counts at all redshifts, their contribution to cosmic Reionization is poorly understood due to a lack of knowledge of their physical properties. We isolate a sample of 35 z~4-5 continuum-faint Lyman-alpha emitters from deep VLT/MUSE spectroscopy and directly measure their Halpha emission using stacked Spitzer/IRAC Ch. 1 photometry. Based on Hubble Space Telescope imaging, we determine that the average UV continuum magnitude is fainter than -16 (~0.01 L_star), implying a median Lyman-alpha equivalent width of 249 Angstroms. By combining the Halpha measurement with the UV magnitude we determine the ionizing photon production efficiency, xi_ion, a first for such faint galaxies. The measurement of log (xi_ion [Hz/erg]) = 26.28 (+0.28; -0.40) is in excess of literature measurements of both continuum- and emission line-selected samples, implying a more efficient production of ionizing photons in these lower-luminosity, Lyman-alpha-selected systems. We conclude that this elevated efficiency can be explained by stellar populations with metallicities between 4e-4 and 0.008, with light-weighted ages less than 3 Myr.

preprint2020arXiv

hybrid-Lenstool: A self-consistent algorithm to model galaxy clusters with strong- and weak-lensing simultaneously

We present a new galaxy cluster lens modeling approach, hybrid-Lenstool, that is implemented in the publicly available modeling software Lenstool. hybrid-Lenstool combines a parametric approach to model the core of the cluster, and a non-parametric (free-form) approach to model the outskirts. hybrid-Lenstool optimizes both strong- and weak-lensing constraints simultaneously (Joint-Fit), providing a self-consistent reconstruction of the cluster mass distribution on all scales. In order to demonstrate the capabilities of the new algorithm, we tested it on a simulated cluster. hybrid-Lenstool yields more accurate reconstructed mass distributions than the former Sequential-Fit approach where the parametric and the non-parametric models are optimized successively. Indeed, we show with the simulated cluster that the mass density profile reconstructed with a Sequential-Fit deviates form the input by $2-3σ$ at all scales while the Joint-Fit gives a profile that is within $1-1.5σ$ of the true value. This gain in accuracy is consequential for recovering mass distributions exploiting cluster lensing and therefore for all applications of clusters as cosmological probes. Finally we found that the Joint-Fit approach yields shallower slope of the inner density profile than the Sequential-Fit approach, thus revealing possible biases in previous lensing studies.

preprint2020arXiv

The BUFFALO HST Survey

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

preprint2020arXiv

The Data Processing Pipeline for the MUSE Instrument

Processing of raw data from modern astronomical instruments is nowadays often carried out using dedicated software, so-called &#34;pipelines&#34; which are largely run in automated operation. In this paper we describe the data reduction pipeline of the Multi Unit Spectroscopic Explorer (MUSE) integral field spectrograph operated at ESO&#39;s Paranal observatory. This spectrograph is a complex machine: it records data of 1152 separate spatial elements on detectors in its 24 integral field units. Efficiently handling such data requires sophisticated software, a high degree of automation and parallelization. We describe the algorithms of all processing steps that operate on calibrations and science data in detail, and explain how the raw science data gets transformed into calibrated datacubes. We finally check the quality of selected procedures and output data products, and demonstrate that the pipeline provides datacubes ready for scientific analysis.

preprint2020arXiv

The Ionizing Photon Production Efficiency ($ξ_{ion}$) Of Lensed Dwarf Galaxies At $z \sim 2 $

We measure the ionizing photon production efficiency ($ξ_{ion}$) of low-mass galaxies ($10^{7.8}$-$10^{9.8}$ $M_{\odot}$) at $1.4<z<2.7$, allowing us to better understand the contribution of dwarf galaxies to the ionizing background and cosmic reionization. We target galaxies that are magnified by the strong lensing galaxy clusters Abell 1689, MACS J0717, and MACS J1149. We utilize Keck/MOSFIRE spectra to measure optical nebular emission line fluxes and HST imaging to measure the rest-UV and rest-optical photometry. We present two methods of stacking. First, we take the average of the log(L$_{Hα}$ /L$_{UV}$) of galaxies in our sample to determine the typical log($ξ_{ion}$). Second, we take the logarithm of the total L$_{Hα}$ over the total L$_{UV}$. We prefer the latter as it provides the total ionizing UV luminosity density of galaxies when multiplied by the non-ionizing UV luminosity density from the UV luminosity function. log($ξ_{ion}$) calculated from the second method is $\sim$ 0.2 dex higher than the first method. We do not find any strong dependence between log($ξ_{ion}$) and stellar mass, M$_{UV}$ or UV spectral slope ($β$). We report a value of log($ξ_{ion}$) $\sim25.47\pm 0.09$ for our UV-complete sample ($-22<M_{UV}<-17.3$) and $\sim25.37\pm0.11$ for our mass-complete sample ($7.8<\log(M_*)<9.8)$. These values are consistent with measurements of more massive, more luminous galaxies in other high-redshift studies that use the same stacking technique. Our log($ξ_{ion}$) is $0.2-0.3$ dex higher than low-redshift galaxies of similar mass, indicating an evolution in the stellar properties, possibly due to metallicity, age, or the prevalence of binary stars. We also find a correlation between log($ξ_{ion}$) and the equivalent widths of H$α$ and [OIII]$λ$5007 fluxes, confirming that these equivalent widths can be used to estimate $ξ_{ion}$.

preprint2020arXiv

The MUSE Hubble Ultra Deep Field Survey XIII. Spatially resolved spectral properties of Lyman alpha haloes around star-forming galaxies at z > 3

We present spatially resolved maps of six individually-detected Lyman alpha haloes (LAHs) as well as a first statistical analysis of the Lyman alpha (Lya) spectral signature in the circum-galactic medium of high-redshift star-forming galaxies using MUSE. Our resolved spectroscopic analysis of the LAHs reveals significant intrahalo variations of the Lya line profile. Using a three-dimensional two-component model for the Lya emission, we measure the full width at half maximum (FWHM), the peak velocity shift and the asymmetry of the Lya line in the core and in the halo of 19 galaxies. We find that the Lya line shape is statistically different in the halo compared to the core for ~40% of our galaxies. Similarly to object-by-object based studies and a recent resolved study using lensing, we find a correlation between the peak velocity shift and the width of the Lya line both at the interstellar and circum-galactic scales. While there is a lack of correlation between the spectral properties and the spatial scale lengths of our LAHs, we find a correlation between the width of the line in the LAH and the halo flux fraction. Interestingly, UV bright galaxies show broader, more redshifted and less asymmetric Lya lines in their haloes. The most significant correlation found is for the FWHM of the line and the UV continuum slope of the galaxy, suggesting that the redder galaxies have broader Lya lines. The generally broad and red line shapes found in the halo component suggests that the Lya haloes are powered either by scattering processes through an outflowing medium, fluorescent emission from outflowing cold clumps of gas, or a mix of both. Considering the large diversity of the Lya line profiles observed in our sample and the lack of strong correlation, the interpretation of our results is still broadly open and underlines the need for realistic spatially resolved models of the LAHs.

preprint2020arXiv

The MUSE Hubble Ultra Deep Field Survey XIV. The evolution of the Lya emitter fraction from z=3 to z=6

The Lya emitter (LAE) fraction, X_LAE, is a potentially powerful probe of the evolution of the intergalactic neutral hydrogen gas fraction. However, uncertainties in the measurement of X_LAE are still debated. Thanks to deep data obtained with MUSE, we can measure the evolution of X_LAE homogeneously over a wide redshift range of z~3-6 for UV-faint galaxies (down to M_1500~-17.75). This is significantly fainter than in former studies, and allows us to probe the bulk of the population of high-z star-forming galaxies. We construct a UV-complete photo-redshift sample following UV luminosity functions and measure the Lya emission with MUSE using the second data release from the MUSE HUDF Survey. We derive the redshift evolution of X_LAE for M_1500 in [-21.75;-17.75] for the first time with a equivalent width range EW(Lya)>=65 A and find low values of X_ LAE<~30% at z<~6. For M_1500 in [-20.25;-18.75] and EW(Lya)<~25 A, our X_LAE values are consistent with those in the literature within 1sigma at z<~5, but our median values are systematically lower than reported values over the whole redshift range. In addition, we do not find a significant dependence of X_LAE on M_1500 for EW(Lya)>~50 A at z~3-4, in contrast with previous work. The differences in X_LAE mainly arise from selection biases for Lyman Break Galaxies (LBGs) in the literature: UV-faint LBGs are more easily selected if they have strong Lya emission, hence X_LAE is biased towards higher values. Our results suggest either a lower increase of X_LAE towards z~6 than previously suggested, or even a turnover of X_LAE at z~5.5, which may be the signature of a late or patchy reionization process. We compared our results with predictions from a cosmological galaxy evolution model. We find that a model with a bursty star formation (SF) can reproduce our observed X_LAE much better than models where SF is a smooth function of time.

preprint2020arXiv

The MUSE Hubble Ultra Deep Field Survey. XV. The mean rest-UV spectra of Ly-alpha emitters at z>3

We investigate the ultraviolet (UV) spectral properties of faint Lyman-$α$ emitters (LAEs) in the redshift range 2.9<z<4.6 and provide material to prepare future observations of the faint Universe. We use data from the MUSE Hubble Ultra Deep Survey to construct mean rest-frame spectra of continuum-faint (median M$_{UV}$ of -18 and down to M$_{UV}$ of -16), low stellar mass (median value of $10^{8.4}$ and down to $10^{7}M_{\odot}$) LAEs at redshift z>3. We compute various averaged spectra of LAEs sub-sampled on the basis of their observational (e.g., Ly$α$ strength, UV magnitude and spectral slope) and physical (e.g., stellar mass and star-formation rate) properties. We search for UV spectral features other than Ly$α$, such as higher-ionization nebular emission lines and absorption features. We successfully observe the OIII]1666 and CIII]909 collisionally excited emission lines and the HeII1640 recombination feature, as well as the resonant CIV1550 doublet either in emission or P-Cygni. We compare the observed spectral properties of the different mean spectra and find the emission lines to vary with the observational and physical properties of the LAEs. In particular, the mean spectra of LAEs with larger Ly$α$ equivalent widths, fainter UV magnitudes, bluer UV spectral slopes and lower stellar masses show the strongest nebular emission. The line ratios of these lines are similar to those measured in the spectra of local metal-poor galaxies, while their equivalent widths are weaker compared to the handful of extreme values detected in individual spectra of z>2 galaxies. This suggests that weak UV features are likely ubiquitous in high z, low-mass and faint LAEs. We publicly release the stacked spectra as they can serve as empirical templates for the design of future observations, such as those with the James Webb Space Telescope and the Extremely Large Telescope.

preprint2020arXiv

The nature of CR7 revealed with MUSE: a young starburst powering extended Lyman-$α$ emission at z=6.6

CR7 is among the most luminous Lyman-$α$ emitters (LAEs) known at $z = 6.6$ and consists of at least three UV components that are surrounded by Lyman-$α$ (Ly$α$) emission. Previous studies have suggested that it may host an extreme ionising source. Here, we present deep integral field spectroscopy of CR7 with VLT/MUSE. We measure extended emission with a similar halo scale length as typical LAEs at $z\approx5$. CR7&#39;s Ly$α$ halo is clearly elongated along the direction connecting the multiple components, likely tracing the underlying gas distribution. The Ly$α$ emission originates almost exclusively from the brightest UV component, but we also identify a faint kinematically distinct Ly$α$ emitting region nearby a fainter component. Combined with new near-infrared data, the MUSE data show that the rest-frame Ly$α$ equivalent width (EW) is $\approx100$ Å. This is a factor four higher than the EW measured in low-redshift analogues with carefully matched Ly$α$ profiles (and thus arguably HI column density), but this EW can plausibly be explained by star formation. Alternative scenarios requiring AGN powering are also disfavoured by the narrower and steeper Ly$α$ spectrum and much smaller IR to UV ratio compared to obscured AGN in other Ly$α$ blobs. CR7&#39;s Ly$α$ emission, while extremely luminous, resembles the emission in more common LAEs at lower redshifts very well and is likely powered by a young metal poor starburst.

preprint2019arXiv

The Detection of [O III]$λ$4363 in a Lensed, Dwarf Galaxy at $z$ = 2.59: Testing Metallicity Indicators and Scaling Relations at High Redshift and Low Mass

We present Keck/MOSFIRE and Keck/LRIS spectroscopy of A1689-217, a lensed (magnification $\sim7.9$), star-forming (SFR $\sim$ 16 $\rm M_\odot\ yr^{-1}$), dwarf (log($M_\ast/M_\odot$) = $8.07-8.59$), Ly$α$-emitter ($\rm EW_0\sim138$ A) at $z$ = 2.5918. Dwarf galaxies similar to A1689-217 are common at high redshift and likely responsible for reionization, yet few have been studied with detailed spectroscopy. We report a 4.2$σ$ detection of the electron-temperature-sensitive [O III]$λ$4363 emission line and use this line to directly measure an oxygen abundance of 12+$\log$(O/H) = 8.06 $\pm$ 0.12 ($\sim1/4\ \rm Z_\odot$). A1689-217 is the lowest mass galaxy at $z>$ 2 with an [O III]$λ$4363 detection. Using the rest-optical emission lines, we measure A1689-217&#39;s other nebular conditions including electron temperature ($T_e$([O III]) $\sim$ 14,000 K), electron density ($n_e\sim$ 220 $\rm cm^{-3}$) and reddening ($E(B-V)\sim$ 0.39). We study relations between strong-line ratios and direct metallicities with A1689-217 and other galaxies with [O III]$λ$4363 detections at $z\sim0-3.1$, showing that the locally-calibrated, oxygen-based, strong-line relations are consistent from $z\sim0-3.1$. We also show additional evidence that the $\rm O_{32}$ vs. $\rm R_{23}$ excitation diagram can be utilized as a redshift-invariant, direct-metallicity-based, oxygen abundance diagnostic out to $z\sim$ 3.1. From this excitation diagram and the strong-line ratio $-$ metallicity plots, we observe that the ionization parameter at fixed O/H is consistent with no redshift evolution. Although A1689-217 is metal-rich for its $M_\ast$ and SFR, we find it to be consistent within the large scatter of the low-mass end of the Fundamental Metallicity Relation.