Showing 1–50 of 62 results for author: Melinder, J

MIDIS. Near-infrared rest-frame morphology of massive galaxies at $3<z<5.5$ in the Hubble eXtreme Deep Field

Authors: L. Costantin, S. Gillman, L. A. Boogaard, P. G. Pérez-González, E. Iani, P. Rinaldi, J. Melinder, A. Crespo Gómez, L. Colina, T. R. Greve, G. Östlin, G. Wright, A. Alonso-Herrero, J. Álvarez-Márquez, M. Annunziatella, A. Bik., K. I. Caputi, D. Dicken, A. Eckart, J. Hjorth, O. Ilbert, I. Jermann, A. Labiano, D. Langeroodi, F. Peißker , et al. (7 additional authors not shown)

Abstract: Thanks to decades of observations using HST, the structure of galaxies at redshift $z>2$ has been widely studied in the rest-frame ultraviolet regime, which traces recent star formation from young stellar populations. But, we still have little information about the spatial distribution of the older, more evolved, stellar populations, constrained by the rest-frame infrared portion of galaxies' spec… ▽ More Thanks to decades of observations using HST, the structure of galaxies at redshift $z>2$ has been widely studied in the rest-frame ultraviolet regime, which traces recent star formation from young stellar populations. But, we still have little information about the spatial distribution of the older, more evolved, stellar populations, constrained by the rest-frame infrared portion of galaxies' spectral energy distribution. We present the morphological characterization of a sample of 21 massive galaxies ($\log(M_{\star}/M_{\odot})>9.5$) at redshift $3<z<5.5$. These galaxies are observed as part of the GTO program MIDIS with the Mid-Infrared Instrument (MIRI) onboard JWST. The deep MIRI 5.6~$μ$m imaging allows us to characterize for the first time the rest-frame near-infrared structure of galaxies beyond cosmic noon, at higher redshifts than possible with NIRCam, tracing their older stellar populations. We derive the galaxies' non-parametric morphology and model the galaxies' light distribution with a Sérsic component. We find that at $z>3$ massive galaxies show a smooth distribution of their rest-infrared light, strongly supporting the increasing number of regular disk galaxies already in place at early epochs. On the contrary, the ultraviolet structure obtained from HST observations is generally more irregular, catching the most recent episodes of star formation. Importantly, we find a segregation of morphologies across cosmic time, having massive galaxies at redshift $z>4$ later-type morphologies compared to $z\sim3$ galaxies. These findings suggest a transition phase in galaxy assembly and central mass build up already taking place at $z\sim3-4$. MIRI provides unique information about the structure of the mature stellar population of high-redshift galaxies, unveiling that massive galaxies beyond cosmic noon are prevalently compact disk galaxies with smooth mass distribution. △ Less

Submitted 28 June, 2024; originally announced July 2024.

Comments: Submitted to A&A. 7 pages, 4 figures. Comments welcome

MIDIS: MIRI uncovers Virgil, an extended source at $z\simeq 6.6$ with the photometric properties of Little Red Dots

Authors: Edoardo Iani, Pierluigi Rinaldi, Karina I. Caputi, Marianna Annunziatella, Danial Langeroodi, Jens Melinder, Pablo G. Pérez-González, Javier Álvarez-Márquez, Leindert A. Boogaard, Sarah E. I. Bosman, Luca Costantin, Thibaud Moutard, Luis Colina, Göran Östlin, Thomas R. Greve, Gillian Wright, Almudena Alonso-Herrero, Arjan Bik, Steven Gillman, Alejandro Crespo Gómez, Jens Hjorth, Alvaro Labiano, John P. Pye, Tuomo V. Tikkanen, Paul P. van der Werf

Abstract: We present Virgil, a MIRI extremely red object (MERO) detected with the F1000W filter as part of the MIRI Deep Imaging Survey (MIDIS) observations of the Hubble Ultra Deep Field (HUDF). Virgil is a Lyman-$α$ emitter (LAE) at $z_{spec} = 6.6312\pm 0.0019$ (from VLT/MUSE) with a rest-frame UV-to-optical spectral energy distribution (SED) typical of LAEs at similar redshifts. However, MIRI observatio… ▽ More We present Virgil, a MIRI extremely red object (MERO) detected with the F1000W filter as part of the MIRI Deep Imaging Survey (MIDIS) observations of the Hubble Ultra Deep Field (HUDF). Virgil is a Lyman-$α$ emitter (LAE) at $z_{spec} = 6.6312\pm 0.0019$ (from VLT/MUSE) with a rest-frame UV-to-optical spectral energy distribution (SED) typical of LAEs at similar redshifts. However, MIRI observations reveal an unexpected extremely red color at rest-frame near-infrared wavelengths, $\rm F444W - F1000W = 2.33 \pm 0.06$. Such steep rise in the near-infrared, completely missed without MIRI imaging, is poorly reproduced by models including only stellar populations and hints towards the presence of an Active Galactic Nucleus (AGN). Interestingly, the overall SED shape of Virgil resembles that of the recently discovered population of Little Red Dots (LRDs) but does not meet their compactness criterion: at rest-frame UV-optical wavelengths Virgil's morphology follows a 2D-Sérsic profile with average index $n = 0.93^{+0.85}_{-0.31}$ and $r_e = 0.43$~pkpc. Only at MIRI wavelengths Virgil is unresolved due to the coarser PSF. We also estimate a bolometric luminosity $L_{\rm bol} = (8.4-11.1)\times 10^{44}\rm~erg~s^{-1}$ and a supermassive black hole mass $M_{\rm BH} = (4-7)\times 10^7\rm ~ M_\odot$ in agreement with recently reported values for LRDs. This discovery demonstrates the crucial importance of deep MIRI surveys to find AGN amongst high-$z$ galaxies that otherwise would be completely missed and raises the question of how common Virgil-like objects could be in the early Universe. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 17 pages, 10 figures, 3 tables. Submitted to ApJ

The emergence of the Star Formation Main Sequence with redshift unfolded by JWST

Authors: P. Rinaldi, R. Navarro-Carrera, K. I. Caputi, E. Iani, G. Ostlin, L. Colina, S. Alberts, J. Alvarez-Marquez, M. Annunziatella, L. Boogaard, L. Costantin, J. Hjorth, D. Langeroodi, J. Melinder, T. Moutard, F. Walter

Abstract: We investigate the correlation between stellar mass (M*) and star formation rate (SFR) across the stellar mass range log10(M*/Msun)~6-11. We consider almost 50,000 star-forming galaxies at z~3-7, leveraging data from COSMOS/SMUVS, JADES/GOODS-SOUTH, and MIDIS/XDF. This is the first study spanning such a wide stellar mass range without relying on gravitational lensing effects. We locate our galaxie… ▽ More We investigate the correlation between stellar mass (M*) and star formation rate (SFR) across the stellar mass range log10(M*/Msun)~6-11. We consider almost 50,000 star-forming galaxies at z~3-7, leveraging data from COSMOS/SMUVS, JADES/GOODS-SOUTH, and MIDIS/XDF. This is the first study spanning such a wide stellar mass range without relying on gravitational lensing effects. We locate our galaxies on the SFR-M* plane to assess how the location of galaxies in the star-formation main sequence (MS) and starburst (SB) region evolves with stellar mass and redshift. We find that the two star-forming modes tend to converge at log10(M*/Msun) < 7, with all galaxies found in the SB mode. By dissecting our galaxy sample in stellar mass and redshift, we show that the emergence of the star-formation MS is stellar-mass dependent: while in galaxies with log10(M*/Msun) > 9 the MS is already well in place at z = 5-7, for galaxies with log10(M*/Msun)~7-8 it only becomes significant at z<4. Overall, our results are in line with previous findings that the SB mode dominates amongst low stellar-mass galaxies. The earlier emergence of the MS for massive galaxies is consistent with galaxy downsizing. △ Less

Submitted 19 June, 2024; originally announced June 2024.

Comments: 15 pages, 4 figures. Submitted to ApJ

The Lyman Alpha Reference Sample XV. Relating Ionised Gas Kinematics with Lyman-$α$ observables

Authors: E. C. Herenz, A. Schaible, P. Laursen, A. Runnholm, J. Melinder, A. Le Reste, M. J. Hayes, G. Östlin, J. Cannon, G. Micheva, M. Roth, K. Saha

Abstract: Gas kinematics affect the radiative transfer and escape of hydrogen Lyman-$α$ (Ly$α$) emission from galaxies. We investigate this interplay empirically by relating the ionised gas kinematics of 42 galaxies in the extended Ly$α$ Reference Sample with their Ly$α$ escape fractions, $f_\rm{esc}$, Ly$α$ equivalent widths, $\rm{EW}_\rm{Lyα}$, and Ly$α$ luminosities, $L_\rm{Lyα}$. To this aim we use PMAS… ▽ More Gas kinematics affect the radiative transfer and escape of hydrogen Lyman-$α$ (Ly$α$) emission from galaxies. We investigate this interplay empirically by relating the ionised gas kinematics of 42 galaxies in the extended Ly$α$ Reference Sample with their Ly$α$ escape fractions, $f_\rm{esc}$, Ly$α$ equivalent widths, $\rm{EW}_\rm{Lyα}$, and Ly$α$ luminosities, $L_\rm{Lyα}$. To this aim we use PMAS integral-field spectroscopic observations of the Balmer-$α$ line. We calculate shearing velocities, $v_\rm{shear}$, and intrinsic velocity dispersions, $σ_0^\rm{obs}$ (empirically corrected for beam-smearing effects), as global kinematical measures for each galaxy. The sample is characterised by highly turbulent motions and more than half of the sample shows dispersion dominated kinematics. We uncover clear trends between Ly$α$ observables and global kinematical statistics. We discuss statistically the importance of $v_\rm{shear}$, $σ_0^\rm{obs}$, and $v_\rm{shear}/σ_0^\rm{obs}$ for regulating the Ly$α$ observables in comparison to other galaxy parameters. It emerges that $σ_0^\rm{obs}$ is the dominating parameter for $\rm{EW}_\rm{Lyα}$ and that is as important as nebular extinction, gas covering fraction, and ionising photon production efficiency in regulating $f_\rm{esc}$. A simple scenario where the starburst age is simultaneously regulating turbulence, $\rm{EW}_\rm{Lyα}$, and $f_\rm{esc}$ does not find support by our observations. However, we show that the small scale distribution of dust appears to be influenced by turbulence in some galaxies. In support of our observational result we discuss how turbulence is theoretically expected to play a significant role in modulating $f_\rm{esc}$. (abridged) △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 24 pages (+ 13 pages appendix), 14 figures (+ 5 figures in appendix). Submitted to A&A. Comments welcome!

Feedback in Emerging Extragalactic Star Clusters (JWST--FEAST): Calibration of Star Formation Rates in the Mid-Infrared with NGC 628

Authors: Daniela Calzetti, Angela Adamo, Sean T. Linden, Benjamin Gregg, Mark R. Krumholz, Varun Bajaj, Arjan Bik, Michele Cignoni, Matteo Correnti, Bruce Elmegreen, Helena Faustino Vieira, John S. Gallagher, Kathryn Grasha, Robert A. Gutermuth, Kelsey E. Johnson, Matteo Messa, Jens Melinder, Goran Ostlin, Alex Pedrini, Elena Sabbi, Linda J. Smith, Monica Tosi

Abstract: New JWST near-infrared imaging of the nearby galaxy NGC 628 from the Cycle 1 program JWST-FEAST is combined with archival JWST mid-infrared imaging to calibrate the 21 $μ$m emission as a star formation rate indicator (SFR) at $\sim$120 pc scales. The Pa$α$ ($λ$1.8756 $μ$m) hydrogen recombination emission line targeted by FEAST provides a reference SFR indicator that is relatively insensitive to du… ▽ More New JWST near-infrared imaging of the nearby galaxy NGC 628 from the Cycle 1 program JWST-FEAST is combined with archival JWST mid-infrared imaging to calibrate the 21 $μ$m emission as a star formation rate indicator (SFR) at $\sim$120 pc scales. The Pa$α$ ($λ$1.8756 $μ$m) hydrogen recombination emission line targeted by FEAST provides a reference SFR indicator that is relatively insensitive to dust attenuation, as demonstrated by combining this tracer with the HST H$α$ imaging. Our analysis is restricted to regions that appear compact in nebular line emission and are sufficiently bright to mitigate effects of both age and stochastic sampling of the stellar initial mass function. We find that the 21 $μ$m emission closely correlates with the nebular line emission, with a power-law with exponent=1.07$\pm$0.01, in agreement with past results. We calibrate a hybrid SFR indicator using a combination of H$α$ and 24 $μ$m (extrapolated from 21 $μ$m) tracers and derive the proportionality constant between the two tracers $b=0.095\pm0.007$, which is $\sim$ 3-5 times larger than previous derivations using large regions/entire galaxies. We model these discrepancies as an increasing contribution to the dust heating by progressively older stellar populations for increasing spatial scales, in agreement with earlier findings that star formation is hierarchically distributed in galaxies. Thus, use of hybrid SFR indicators requires prior knowledge of the mean age of the stellar populations dominating the dust heating, which makes their application uncertain. Conversely, non-linear calibrations of SFRs from L(24) alone are more robust, with a factor $\lesssim$2.5 variation across the entire range of L(24) luminosities from HII regions to galaxies. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: 28 pages, 9 figures. Accepted for publication on the Astrophysical Journal

Feedback in Emerging extragAlactic Star clusTers, FEAST: JWST spots PAH destruction in NGC 628 during the emerging phase of star formation

Authors: Alex Pedrini, Angela Adamo, Daniela Calzetti, Arjan Bik, Benjamin Gregg, Sean T. Linden, Varun Bajaj, Jenna E. Ryon, Ahmad A. Ali, Giacomo Bortolini, Matteo Correnti, Bruce G. Elmegreen, Debra Meloy Elmegreen, John S. Gallagher, Kathryn Grasha, Robert A. Gutermuth, Kelsey E. Johnson, Jens Melinder, Matteo Messa, Göran Östlin, Elena Sabbi, Linda J. Smith, Monica Tosi, Helena Faustino Vieira

Abstract: We investigate the emergence phase of young star clusters in the nearby spiral galaxy NGC 628. We use JWST NIRCam and MIRI observations to create spatially resolved maps of the Pa$α$-1.87 $μ$m and Br$α$-4.05 $μ$m hydrogen recombination lines, as well as the 3.3 $μ$m and 7.7 $μ$m emission from polycyclic aromatic hydrocarbons (PAHs). We extract 953 compact HII regions and analyze the PAH emission a… ▽ More We investigate the emergence phase of young star clusters in the nearby spiral galaxy NGC 628. We use JWST NIRCam and MIRI observations to create spatially resolved maps of the Pa$α$-1.87 $μ$m and Br$α$-4.05 $μ$m hydrogen recombination lines, as well as the 3.3 $μ$m and 7.7 $μ$m emission from polycyclic aromatic hydrocarbons (PAHs). We extract 953 compact HII regions and analyze the PAH emission and morphology at $\sim$10 pc scales in the associated photo-dissociation regions (PDRs). While HII regions remain compact, radial profiles help us to define three PAH morphological classes: compact ($\sim$ 42%), extended ($\sim$ 34%) and open ($\sim$ 24%). The majority of compact and extended PAH morphologies are associated with very young star clusters ($<$5 Myr), while open PAH morphologies are mainly associated with star clusters older than 3 Myr. We observe a general decrease in the 3.3 $μ$m and 7.7 $μ$m PAH band emission as a function of cluster age, while their ratio remains constant with age out to 10 Myr and morphological class. The recovered PAH$_{3.3 μ{\rm m}}$/PAH$_{7.7 μ{\rm m}}$ ratio is lower than values reported in the literature for reference models that consider neutral and ionized PAH populations and analyses conducted at galactic scales. The 3.3 $μ$m and 7.7 $μ$m bands are typically associated to neutral and ionised PAHs, respectively. While we expected neutral PAHs to be suppressed in proximity of the ionizing source, the constant PAH$_{3.3 μ{\rm m}}$/PAH$_{7.7 μ{\rm m}}$ ratio would indicate that both families of molecules disrupt at similar rates in proximity of the HII regions. △ Less

Submitted 26 June, 2024; v1 submitted 3 June, 2024; originally announced June 2024.

Comments: 25 pages, 14 figures, 3 tables. Accepted for publication in ApJ V2: Minor changes to Figures 7, 8, and 9, and to the text

Feedback in emerging extragalactic star clusters, FEAST: The relation between 3.3 $μ$m PAH emission and Star Formation Rate traced by ionized gas in NGC 628

Authors: Benjamin Gregg, Daniela Calzetti, Angela Adamo, Varun Bajaj, Jenna E. Ryon, Sean T. Linden, Matteo Correnti, Michele Cignoni, Matteo Messa, Elena Sabbi, John S. Gallagher, Kathryn Grasha, Alex Pedrini, Robert A. Gutermuth, Jens Melinder, Ralf Kotulla, Gustavo Pérez, Mark R. Krumholz, Arjan Bik, Göran Östlin, Kelsey E. Johnson, Giacomo Bortolini, Linda J. Smith, Monica Tosi, Subhransu Maji , et al. (1 additional authors not shown)

Abstract: We present maps of ionized gas (traced by Pa$α$ and Br$α$) and 3.3 $μ$m Polycyclic Aromatic Hydrocarbon (PAH) emission in the nearby spiral galaxy NGC 628, derived from new JWST/NIRCam data from the FEAST survey. With this data, we investigate and calibrate the relation between 3.3 $μ$m PAH emission and star formation rate (SFR) in and around emerging young star clusters (eYSCs) on a scale of… ▽ More We present maps of ionized gas (traced by Pa$α$ and Br$α$) and 3.3 $μ$m Polycyclic Aromatic Hydrocarbon (PAH) emission in the nearby spiral galaxy NGC 628, derived from new JWST/NIRCam data from the FEAST survey. With this data, we investigate and calibrate the relation between 3.3 $μ$m PAH emission and star formation rate (SFR) in and around emerging young star clusters (eYSCs) on a scale of ${\sim}40$ pc. We find a tight (correlation coefficient $ρ$${\sim}$0.9) sub-linear (power-law exponent $α$${\sim}$0.75) relation between the 3.3 $μ$m PAH luminosity surface density and SFR traced by Br$α$ for compact, cospatial (within 0.16$''$ or ${\sim}$7 pc) peaks in Pa$α$, Br$α$, and 3.3 $μ$m (eYSC-I). The scatter in the relationship does not correlate well with variations in local interstellar medium (ISM) metallicity due to a radial metallicity gradient, but rather is likely due to stochastic sampling of the stellar initial mass function (IMF) and variations in the PAH heating and age of our sources. The deviation from a linear relation may be explained by PAH destruction in more intense ionizing environments, variations in age, and IMF stochasticity at intermediate to low luminosities. We test our results with various continuum subtraction techniques using combinations of NIRCam bands and find that they remain robust with only minor differences in the derived slope and intercept. An unexpected discrepancy is identified between the relations of hydrogen recombination lines (Pa$α$ versus Br$α$; H$α$ versus Br$α$). △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: 28 pages, 10 figures, 3 tables

Haro 11: The Spatially Resolved Lyman Continuum Sources

Authors: Lena Komarova, M. S. Oey, Svea Hernandez, Angela Adamo, Mattia Sirressi, Claus Leitherer, J. M. Mas-Hesse, Goran Ostlin, Edmund Hodges-Kluck, Arjan Bik, Matthew J. Hayes, Anne E. Jaskot, Daniel Kunth, Peter Laursen, Jens Melinder, T. Emil Rivera-Thorsen

Abstract: As the nearest confirmed Lyman continuum (LyC) emitter, Haro 11 is an exceptional laboratory for studying LyC escape processes crucial to cosmic reionization. Our new HST/COS G130M/1055 observations of its three star-forming knots now reveal that the observed LyC originates in Knots B and C, with $903 - 912~Å$ luminosities of $1.9\pm1.5 \times 10^{40}~\rm erg~s^{-1}$ and… ▽ More As the nearest confirmed Lyman continuum (LyC) emitter, Haro 11 is an exceptional laboratory for studying LyC escape processes crucial to cosmic reionization. Our new HST/COS G130M/1055 observations of its three star-forming knots now reveal that the observed LyC originates in Knots B and C, with $903 - 912~Å$ luminosities of $1.9\pm1.5 \times 10^{40}~\rm erg~s^{-1}$ and $0.9\pm0.7 \times 10^{40}~\rm erg~s^{-1}$, respectively. We derive local escape fractions $f_{\rm{esc, 912}} = 3.4\pm2.9\%$ and $5.1\pm4.3\%$ for Knots B and C, respectively. Our Starburst99 modeling shows dominant populations on the order of $\sim1-4$ Myr and $1-2\times10^7 \rm~M_\odot$ in each knot, with the youngest population in Knot B. Thus, the knot with the strongest LyC detection has the highest LyC production. However, LyC escape is likely less efficient in Knot B than in Knot C due to higher neutral gas covering. Our results therefore stress the importance of the intrinsic ionizing luminosity, and not just the escape fraction, for LyC detection. Similarly, the Ly$α$ escape fraction does not consistently correlate with LyC flux, nor do narrow Ly$α$ red peaks. High observed Ly$α$ luminosity and low Ly$α$ peak velocity separation, however, do correlate with higher LyC escape. Another insight comes from the undetected Knot A, which drives the Green Pea properties of Haro 11. Its density-bounded conditions suggest highly anisotropic LyC escape. Finally, both of the LyC-leaking Knots, B and C, host ultra-luminous X-ray sources (ULXs). While stars strongly dominate over the ULXs in LyC emission, this intriguing coincidence underscores the importance of unveiling the role of accretors in LyC escape and reionization. △ Less

Submitted 1 April, 2024; originally announced April 2024.

Comments: 22 pages, 5 figures. Accepted for publication in The Astrophysical Journal

Glimmers in the Cosmic Dawn: A Census of the Youngest Supermassive Black Holes by Photometric Variability

Authors: Matthew J. Hayes, Jonathan C. Tan, Richard S. Ellis, Alice R. Young, Vieri Cammelli, Jasbir Singh, Axel Runnholm, Aayush Saxena, Ragnhild Lunnan, Benjamin W. Keller, Pierluigi Monaco, Nicolas Laporte, Jens Melinder

Abstract: We report first results from a deep near infrared campaign with the Hubble Space Telescope to obtain late-epoch images of the Hubble Ultra-Deep Field (HUDF), 10-15 years after the first epoch data were obtained. The main objectives are to search for faint active galactic nuclei (AGN) at high redshifts by virtue of their photometric variability, and measure (or constrain) the comoving number densit… ▽ More We report first results from a deep near infrared campaign with the Hubble Space Telescope to obtain late-epoch images of the Hubble Ultra-Deep Field (HUDF), 10-15 years after the first epoch data were obtained. The main objectives are to search for faint active galactic nuclei (AGN) at high redshifts by virtue of their photometric variability, and measure (or constrain) the comoving number density of supermassive black holes (SMBHs), n_SMBH, at early times. In this Letter we present a brief overview of the program and preliminary results regarding eight objects. Three variables are supernovae, two of which are apparently hostless with indeterminable redshifts, although one has previously been recorded at a z\approx 6 galaxy. Two further objects are clear AGN candidates at z=2.0 and 3.2, based on morphology and/or spectroscopy, in particular infrared spectroscopy from JWST. Three variable targets are identified at z=6-7, which are also likely AGN candidates. These sources provide a first measure of n_SMBH in the reionization epoch by photometric variability, which places a firm lower limit of 3x10^{-4} cMpc^{-3}. After accounting for variability and luminosity incompleteness, we estimate n_SMBH \gtrsim 8x10^{-3} cMpc^{-3}, which is the largest value so far reported at these redshifts. This SMBH abundance is also strikingly similar to estimates of n_SMBH in the local Universe. We discuss how these results test various theories for SMBH formation. △ Less

Submitted 24 March, 2024; originally announced March 2024.

Comments: Submitted to ApJ

A NIRCam-dark galaxy detected with the MIRI/F1000W filter in the MIDIS/JADES Hubble Ultra Deep Field

Authors: Pablo G. Pérez-González, Pierluigi Rinaldi, Karina I. Caputi, Javier Álvarez-Márquez, Marianna Annunziatella, Danial Langeroodi, Thibaud Moutard, Leindert Boogaard, Edoardo Iani, Jens Melinder, Luca Costantin, Goran Östlin, Luis Colina, Thomas R. Greve, Gillian Wright, Almudena Alonso-Herrero, Arjan Bik, Sarah E. I. Bosman, Alejandro Crespo Gómez, Daniel Dicken, Andreas Eckart, Macarena García-Marín, Steven Gillman, Manuel Güdel, Thomas Henning , et al. (10 additional authors not shown)

Abstract: We report the discovery of Cerberus, an extremely red object detected with the MIRI Deep Imaging Survey (MIDIS) observations in the F1000W filter of the Hubble Ultra Deep Field. The object is detected at $S/N\sim6$, with $\mathrm{F1000W}\sim27$ mag, and undetected in the NIRCam data gathered by the JWST Advanced Deep Extragalactic Survey, JADES, fainter than the 30.0-30.5 mag $5σ$ detection limits… ▽ More We report the discovery of Cerberus, an extremely red object detected with the MIRI Deep Imaging Survey (MIDIS) observations in the F1000W filter of the Hubble Ultra Deep Field. The object is detected at $S/N\sim6$, with $\mathrm{F1000W}\sim27$ mag, and undetected in the NIRCam data gathered by the JWST Advanced Deep Extragalactic Survey, JADES, fainter than the 30.0-30.5 mag $5σ$ detection limits in individual bands, as well as in the MIDIS F560W ultra-deep data ($\sim$29 mag, $5σ$). Analyzing the spectral energy distribution built with low-$S/N$ ($<5$) measurements in individual optical-to-mid-infrared filters and higher $S/N$ ($\gtrsim5$) in stacked NIRCam data, we discuss the possible nature of this red NIRCam-dark source using a battery of codes. We discard the possibility of Cerberus being a Solar System body based on the $<0.016$" proper motion in the 1-year apart JADES and MIDIS observations. A sub-stellar Galactic nature is deemed unlikely, given that the Cerberus' relatively flat NIRCam-to-NIRCam and very red NIRCam-to-MIRI flux ratios are not consistent with any brown dwarf model. The extragalactic nature of Cerberus offers 3 possibilities: (1) A $z\sim0.4$ galaxy with strong emission from polycyclic aromatic hydrocarbons; the very low inferred stellar mass, $\mathrm{M}_\star=10^{5-6}$ M$_\odot$, makes this possibility highly improbable. (2) A dusty galaxy at $z\sim4$ with an inferred stellar mass $\mathrm{M}_\star\sim10^{8}$ M$_\odot$. (3) A galaxy with observational properties similar to those of the reddest little red dots discovered around $z\sim7$, but Cerberus lying at $z\sim15$, with the rest-frame optical dominated by emission from a dusty torus or a dusty starburst. △ Less

Submitted 26 May, 2024; v1 submitted 26 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJL, results unchanged, manuscript shortened to publish as a letter

SN 2021adxl: A luminous nearby interacting supernova in an extremely low metallicity environment

Authors: S. J. Brennan, S. Schulze, R. Lunnan, J. Sollerman, L. Yan, C. Fransson, I. Irani, J. Melinder, T. -W. Chen, K. De, C. Fremling, Y. -L. Kim, D. Perley, P. J. Pessi, A. J. Drake, M. J. Graham, R. R. Laher, F. J. Masci, J. Purdum, H. Rodriguez

Abstract: SN 2021adxl is a slowly evolving, luminous, Type IIn supernova with asymmetric emission line profiles, similar to the well-studied SN 2010jl. We present extensive optical, near-ultraviolet, and near-infrared photometry and spectroscopy covering ~1.5 years post discovery. SN 2021adxl occurred in an unusual environment, atop a vigorously star-forming region that is offset from its host galaxy core.… ▽ More SN 2021adxl is a slowly evolving, luminous, Type IIn supernova with asymmetric emission line profiles, similar to the well-studied SN 2010jl. We present extensive optical, near-ultraviolet, and near-infrared photometry and spectroscopy covering ~1.5 years post discovery. SN 2021adxl occurred in an unusual environment, atop a vigorously star-forming region that is offset from its host galaxy core. The appearance of Lyα, O II, as well as the compact core, would classify the host of SN 2021adxl as a Blueberry galaxy, analogous to the higher redshift Green Pea galaxies. Using several abundance indicators, we find a metallicity of the explosion environment of only 10% solar, the lowest reported metallicity for a Type IIn SN environment. SN 2021adxl reaches a peak magnitude of r ~ -20.2 mag and since discovery, SN 2021adxl has faded by only ~4 magnitudes in the r band with a cumulative radiated energy of ~1.5e50 erg over 18 months. SN 2021adxl shows strong signs of interaction with a complex circumstellar medium, seen by the detection of X-rays, revealed by the detection of coronal emission lines, and through multi-component hydrogen and helium profiles. In order to further understand this interaction, we model the Hα profile using a Monte-Carlo electron scattering code. The blueshifted high-velocity component is consistent with emission from a radially thin, spherical shell resulting in the broad emission components due to electron scattering. Using the velocity evolution of this emitting shell, we find that the SN ejecta collide with circumstellar material of at least 5 Msun, assuming a steady-state mass-loss rate of 4-6e-3 Msun per year for the first ~200 days of evolution. Continuing the observations of SN 2021adxl may reveal signatures of dust formation or an infrared excess, similar to that seen for SN 2010jl. △ Less

Submitted 20 December, 2023; originally announced December 2023.

Comments: 21 pages with 27 Figures. Submitted to A&A. MC line scattering code can be found at https://github.com/Astro-Sean/escatter

Clumpy star formation and an obscured nuclear starburst in the luminous dusty z=4 galaxy GN20 seen by MIRI/JWST

Authors: A. Bik, J. Álvarez-Márquez, L. Colina, A. Crespo Gómez, F. Peissker, F. Walter, L. A. Boogaard, G. Östlin, T. R. Greve, G. Wright, A. Alonso-Herrero, K. I. Caputi, L. Costantin, A. Eckart, S. Gillman, J. Hjorth, E. Iani, I. Jermann, A. Labiano, D. Langeroodi, J. Melinder, P. G. Pérez-González, J. P. Pye, P. Rinaldi, T. Tikkanen , et al. (6 additional authors not shown)

Abstract: Dusty star-forming galaxies emit most of their light at far-IR to mm wavelengths as their star formation is highly obscured. Far-IR and mm observations have revealed their dust, neutral and molecular gas properties. The sensitivity of JWST at rest-frame optical and near-infrared wavelengths now allows the study of the stellar and ionized gas content. We investigate the spatially resolved distribut… ▽ More Dusty star-forming galaxies emit most of their light at far-IR to mm wavelengths as their star formation is highly obscured. Far-IR and mm observations have revealed their dust, neutral and molecular gas properties. The sensitivity of JWST at rest-frame optical and near-infrared wavelengths now allows the study of the stellar and ionized gas content. We investigate the spatially resolved distribution and kinematics of the ionized gas in GN20, a dusty star forming galaxy at $z$=4.0548. We present deep MIRI/MRS integral field spectroscopy of the near-infrared rest-frame emission of GN20. We detect spatially resolved \paa, out to a radius of 6 kpc, distributed in a clumpy morphology. The star formation rate derived from \paa\ (144 $\pm$ 9 \msunperyear) is only 7.7 $\pm 0.5 $\% of the infrared star formation rate (1860 $\pm$ 90 \msunperyear). We attribute this to very high extinction (A$_V$ = 17.2 $\pm$ 0.4 mag, or A$_{V,mixed}$ = 44 $\pm$ 3 mag), especially in the nucleus of GN20, where only faint \paa\ is detected, suggesting a deeply buried starburst. We identify four, spatially unresolved, clumps in the \paa\ emission. Based on the double peaked \paa\ profile we find that each clump consist of at least two sub-clumps. We find mass upper limits consistent with them being formed in a gravitationally unstable gaseous disk. The UV bright region of GN20 does not have any detected \paa\ emission, suggesting an age of more than 10 Myrs for this region of the galaxy. From the rotation profile of \paa\ we conclude that the gas kinematics are rotationally dominated and the $v_{rot}/σ_{m} = 3.8 \pm 1.4$ is similar to low-redshift LIRGs. We speculate that the clumps seen in GN20 could contribute to building up the inner disk and bulge of GN20. △ Less

Submitted 2 March, 2024; v1 submitted 5 December, 2023; originally announced December 2023.

Comments: 13 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 686, A3 (2024)

MIDIS: The Relation between Strong (Hb+[OIII]) Emission, Star Formation and Burstiness Around the Epoch of Reionization

Authors: Karina I. Caputi, Pierluigi Rinaldi, Edoardo Iani, Pablo G. Pérez-González, Göran Ostlin, Luis Colina, Thomas R. Greve, Hans-Ulrik Nørgaard-Nielsen, Gillian S. Wright, Javier Alvarez-Márquez, Andreas Eckart, Jens Hjorth, Alvaro Labiano, Olivier Le Fèvre, Fabian Walter, Paul van der Werf, Leindert Boogaard, Luca Costantin, Alejandro Crespo-Gómez, Steven Gillman, Iris Jermann, Danial Langeroodi, Jens Melinder, Florian Peissker, Manuel Güdel , et al. (3 additional authors not shown)

Abstract: We investigate the properties of strong (Hb+[OIII]) emitters before and after the end of the Epoch of Reionization from z=8 to z=5.5. We make use of ultra-deep JWST/NIRCam imaging in the Parallel Field of the MIRI Deep Imaging Survey (MIDIS) in the Hubble eXtreme Deep Field (P2-XDF), in order to select prominent (Hb+[OIII]) emitters (with rest EW_0 > 100 Angstroms) at z=5.5-7, based on their flux… ▽ More We investigate the properties of strong (Hb+[OIII]) emitters before and after the end of the Epoch of Reionization from z=8 to z=5.5. We make use of ultra-deep JWST/NIRCam imaging in the Parallel Field of the MIRI Deep Imaging Survey (MIDIS) in the Hubble eXtreme Deep Field (P2-XDF), in order to select prominent (Hb+[OIII]) emitters (with rest EW_0 > 100 Angstroms) at z=5.5-7, based on their flux density enhancement in the F356W band with respect to the spectral energy distribution continuum. We complement our selection with other (Hb+[OIII]) emitters from the literature at similar and higher (z=7-8) redshifts. We find (non-independent) anti-correlations between EW_0(Hb+[OIII]) and both galaxy stellar mass and age, in agreement with previous studies, and a positive correlation with specific star formation rate (sSFR). On the SFR-M* plane, the (Hb+[OIII]) emitters populate both the star-formation main sequence and the starburst region, which become indistinguishable at low stellar masses (log10(M*) < 7.5). We find tentative evidence for a non-monotonic relation between EW_0(Hb+[OIII]) and SFR, such that both parameters correlate with each other at SFR > 1 Msun/yr, while the correlation flattens out at lower SFRs. This suggests that low metallicities producing high EW_0(Hb+[OIII]) could be important at low SFR values. Interestingly, the properties of the strong emitters and other galaxies (33% and 67% of our z=5.5-7 sample, respectively) are similar, including, in many cases, high sSFR. Therefore, it is crucial to consider both emitters and non-emitters to obtain a complete picture of the cosmic star formation activity around the Epoch of Reionization. △ Less

Submitted 18 June, 2024; v1 submitted 21 November, 2023; originally announced November 2023.

Comments: 14 pages, 8 figures. Accepted for publication at the ApJ. Version v3 includes minor update to match journal version

Nebular C IV 1550 Imaging of the Metal-Poor Starburst Mrk 71: Direct Evidence of Catastrophic Cooling

Authors: M. S. Oey, Amit N. Sawant, Ashkbiz Danehkar, Sergiy Silich, Linda J. Smith, Jens Melinder, Claus Leitherer, Matthew Hayes, Anne E. Jaskot, Daniela Calzetti, You-Hua Chu, Bethan L. James, Goeran Oestlin

Abstract: We use the Hubble Space Telescope ACS camera to obtain the first spatially resolved, nebular imaging in the light of C IV 1548,1551 by using the F150LP and F165LP filters. These observations of the local starburst Mrk 71 in NGC 2366 show emission apparently originating within the interior cavity around the dominant super star cluster (SSC), Knot A. Together with imaging in He II 4686 and supportin… ▽ More We use the Hubble Space Telescope ACS camera to obtain the first spatially resolved, nebular imaging in the light of C IV 1548,1551 by using the F150LP and F165LP filters. These observations of the local starburst Mrk 71 in NGC 2366 show emission apparently originating within the interior cavity around the dominant super star cluster (SSC), Knot A. Together with imaging in He II 4686 and supporting STIS FUV spectroscopy, the morphology and intensity of the C IV nebular surface brightness and the C IV / He II ratio map provide direct evidence that the mechanical feedback is likely dominated by catastrophic radiative cooling, which strongly disrupts adiabatic superbubble evolution. The implied extreme mass loading and low kinetic efficiency of the cluster wind are reasonably consistent with the wind energy budget, which is probably enhanced by radiation pressure. In contrast, the Knot B SSC lies within a well-defined superbubble with associated soft X-rays and He II 1640 emission, which are signatures of adiabatic, energy-driven feedback from a supernova-driven outflow. This system lacks clear evidence of C IV from the limb-brightened shell, as expected for this model, but the observations may not be deep enough to confirm its presence. We also detect a small C IV-emitting object that is likely an embedded compact H II region. Its C IV emission may indicate the presence of very massive stars (> 100 M_sun) or strongly pressure-confined stellar feedback. △ Less

Submitted 20 October, 2023; originally announced October 2023.

Comments: 13 pages, 8 figures. Accepted to ApJ Letters

Journal ref: Astrophys.J.Lett.958:L10,2023

MIDIS: Unveiling the Role of Strong Ha-emitters during the Epoch of Reionization with JWST

Authors: P. Rinaldi, K. I. Caputi, E. Iani, L. Costantin, S. Gillman, P. G. Perez-Gonzalez, G. Ostlin, L. Colina, T. R. Greve, H. U. Noorgard-Nielsen, G. S. Wright, J. Alvarez-Marquez, A. Eckart, M. Garcia-Marin, J. Hjorth, O. Ilbert, S. Kendrew, A. Labiano, O. Le Fevre, J. Pye, T. Tikkanen, F. Walter, P. van der Werf, M. Ward, M. Annunziatella , et al. (18 additional authors not shown)

Abstract: By using the ultra-deep \textit{JWST}/MIRI image at 5.6 $μm$ in the Hubble eXtreme Deep Field, we constrain the role of strong H$α$-emitters (HAEs) during Cosmic Reionization at $z\simeq7-8$. Our sample of HAEs is comprised of young ($<35\;\rm Myr$) galaxies, except for one single galaxy ($\approx 300\;\rm Myr$), with low stellar masses ($\lesssim 10^{9}\;\rm M_{\odot}$). These HAEs show a wide ra… ▽ More By using the ultra-deep \textit{JWST}/MIRI image at 5.6 $μm$ in the Hubble eXtreme Deep Field, we constrain the role of strong H$α$-emitters (HAEs) during Cosmic Reionization at $z\simeq7-8$. Our sample of HAEs is comprised of young ($<35\;\rm Myr$) galaxies, except for one single galaxy ($\approx 300\;\rm Myr$), with low stellar masses ($\lesssim 10^{9}\;\rm M_{\odot}$). These HAEs show a wide range of UV-$β$ slopes, with a median value of $β= -2.15\pm0.21$ which broadly correlates with stellar mass. We estimate the ionizing photon production efficiency ($ξ_{ion,0}$) of these sources (assuming $f_{esc,LyC} = 0\%$), which yields a median value $\rm log_{10}(ξ_{ion,0}/(Hz\;erg^{-1})) = 25.50^{+0.10}_{-0.12}$. We show that $ξ_{ion,0}$ positively correlates with EW$_{0}$(H$α$) and specific star formation rate (sSFR). Instead $ξ_{ion,0}$ weakly anti-correlates with stellar mass and $β$. Based on the $β$ values, we predict $f_{esc, LyC}=4\%^{+3}_{-2}$, which results in $\rm log_{10}(ξ_{ion}/(Hz\;erg^{-1})) = 25.55^{+0.11}_{-0.13}$. Considering this and related findings from the literature, we find a mild evolution of $ξ_{ion}$with redshift. Additionally, our results suggest that these HAEs require only modest escape fractions ($f_{esc, rel}$) of 6$-$15\% to reionize their surrounding intergalactic medium. By only considering the contribution of these HAEs, we estimated their total ionizing emissivity ($\dot{N}_{ion}$) as $\dot{N}_{ion} = 10^{50.53 \pm 0.45}; \text{s}^{-1}\text{Mpc}^{-3}$. When comparing their $\dot{N}_{ion}$ with "non-H$α$ emitter" galaxies across the same redshift range, we find that that strong, young, and low-mass emitters may have played an important role during Cosmic Reionization. △ Less

Submitted 13 June, 2024; v1 submitted 27 September, 2023; originally announced September 2023.

Comments: 21 pages, 11 Figures, 1 table. Accepted for publication in ApJ. Updated version 13/06/2024 -- found a mistake in the reported eq. 2 (now corrected)

MIDIS: JWST NIRCam and MIRI unveil the stellar population properties of Ly$α$-emitters and Lyman-Break galaxies at z ~ 3-7

Authors: Edoardo Iani, Karina I. Caputi, Pierluigi Rinaldi, Marianna Annunziatella, Leindert A. Boogaard, Göran Östlin, Luca Costantin, Steven Gillman, Pablo G. Pérez-González, Luis Colina, Gillian Wright, Almudena Alonso-Herrero, Javier Álvarez-Márquez, Arjan Bik, Sarah E. I. Bosman, Alejandro Crespo-Gómez, Andreas Eckart, Thomas R. Greve, Thomas K. Henning, Jens Hjorth, Iris Jermann, Alvaro Labiano, Danial Langeroodi, Jens Melinder, Thibaud Moutard , et al. (7 additional authors not shown)

Abstract: We study the stellar population properties of 182 spectroscopically-confirmed (MUSE/VLT) Lyman-$α$ emitters (LAEs) and 450 photometrically-selected Lyman-Break galaxies (LBGs) at z = 2.8 - 6.7 in the Hubble eXtreme Deep Field (XDF). Leveraging the combined power of HST and JWST NIRCam and MIRI observations, we analyse their rest-frame UV-through-near-IR spectral energy distributions (SEDs) with MI… ▽ More We study the stellar population properties of 182 spectroscopically-confirmed (MUSE/VLT) Lyman-$α$ emitters (LAEs) and 450 photometrically-selected Lyman-Break galaxies (LBGs) at z = 2.8 - 6.7 in the Hubble eXtreme Deep Field (XDF). Leveraging the combined power of HST and JWST NIRCam and MIRI observations, we analyse their rest-frame UV-through-near-IR spectral energy distributions (SEDs) with MIRI playing a crucial role in robustly assessing the LAE's stellar mass and ages. Our LAEs are low-mass objects (log$_{10}$(M$_\star$[M$_\odot$]) ~ 7.5), with little or no dust extinction (E(B - V) ~ 0.1) and a blue UV continuum slope ($β$ ~ -2.2). While 75% of our LAEs are young (< 100 Myr), the remaining 25% have significantly older stellar populations (> 100 Myr). These old LAEs are statistically more massive, less extinct and have lower specific star formation rate (sSFR) compared to young LAEs. Besides, they populate the M$_\star$ - SFR plane along the main-sequence (MS) of star-forming galaxies, while young LAEs populate the starburst region. The comparison between the LAEs properties to those of a stellar-mass matched sample of LBGs shows no statistical difference between these objects, except for the LBGs redder UV continuum slope and marginally larger E(B - V) values. Interestingly, 48% of the LBGs have ages < 10 Myr and are classified as starbursts, but lack detectable Ly$α$ emission. This is likely due to HI resonant scattering and/or selective dust extinction. Overall, we find that JWST observations are crucial in determining the properties of LAEs and shedding light on the properties and similarities between LAEs and LBGs. △ Less

Submitted 25 January, 2024; v1 submitted 15 September, 2023; originally announced September 2023.

Comments: 22 pages, 15 figures, 5 tables. Accepted for publication in APJ

Spatially-resolved H$α$ and ionizing photon production efficiency in the lensed galaxy MACS1149-JD1 at a redshift of 9.11

Authors: J. Álvarez-Márquez, L. Colina, A. Crespo Gómez, P. Rinaldi, J. Melinder, G. Östlin, M. Annunziatella, A. Labiano, A. Bik, S. Bosman, T. R. Greve, G. Wright, A. Alonso-Herrero, L. Boogaard, R. Azollini, K. I. Caputi, L. Costantin, A. Eckart, M. GarcÍa-MarÍn, S. Gillman, J. Hjorth, E. Iani, O. Ilbert, I. Jermann, D. Langeroodi , et al. (10 additional authors not shown)

Abstract: We present MIRI/JWST medium-resolution spectroscopy (MRS) and imaging (MIRIM) of the lensed galaxy MACS1149-JD1 at a redshift of $z$=9.1092$\pm$0.0002 (Universe age about 530 Myr). We detect, for the first time, spatially resolved H$α$ emission in a galaxy at a redshift above nine. The structure of the H$α$ emitting gas consists of two clumps, S and N. The total H$α$ luminosity implies an instanta… ▽ More We present MIRI/JWST medium-resolution spectroscopy (MRS) and imaging (MIRIM) of the lensed galaxy MACS1149-JD1 at a redshift of $z$=9.1092$\pm$0.0002 (Universe age about 530 Myr). We detect, for the first time, spatially resolved H$α$ emission in a galaxy at a redshift above nine. The structure of the H$α$ emitting gas consists of two clumps, S and N. The total H$α$ luminosity implies an instantaneous star-formation of 5.3$\pm$0.4 $M_{\odot}$ yr$^{-1}$ for solar metallicities. The ionizing photon production efficiency, $\log(ζ_\mathrm{ion})$, shows a spatially resolved structure with values of 25.55$\pm$0.03, 25.47$\pm$0.03, and 25.91$\pm$0.09 Hz erg$^{-1}$ for the integrated galaxy, and clumps S and N, respectively. The H$α$ rest-frame equivalent width, EW$_{0}$(H$α$), is 726$^{+660}_{-182}$ Ángstrom for the integrated galaxy, but presents extreme values of 531$^{+300}_{-96}$ Ángstrom and $\geq$1951 Ángstrom for clumps S and N, respectively. The spatially resolved ionizing photon production efficiency is within the range of values measured in galaxies at redshift above six, and well above the canonical value (25.2$\pm$0.1 Hz erg$^{-1}$). The extreme difference of EW$_{0}$(H$α$) for Clumps S and N indicates the presence of a recent (<5 Myrs) burst in clump N and a star formation over a larger period of time (e.g., $\sim$50 Myr) in clump S. Finally, clump S and N show very different H$α$ kinematics with velocity dispersions of 56$\pm$4 km s$^{-1}$ and 113$\pm$33 km s$^{-1}$, likely indicating the presence of outflows or increased turbulence in the clump N. The dynamical mass, $M_\mathrm{dyn}$= (2.4$\pm$0.5)$\times$10$^{9}$ $M_{\odot}$, is within the range measured with the spatially resolved [OIII]88$μ$m line. △ Less

Submitted 24 March, 2024; v1 submitted 12 September, 2023; originally announced September 2023.

Comments: 11 pages, 6 figures, published version in A&A

Journal ref: A&A 686, A85 (2024)

MIDIS: JWST/MIRI reveals the Stellar Structure of ALMA-selected Galaxies in the Hubble-UDF at Cosmic Noon

Authors: Leindert A. Boogaard, Steven Gillman, Jens Melinder, Fabian Walter, Luis Colina, Göran Östlin, Karina I. Caputi, Edoardo Iani, Pablo Pérez-González, Paul van der Werf, Thomas R. Greve, Gillian Wright, Almudena Alonso-Herrero, Javier Álvarez-Márquez, Marianna Annunziatella, Arjan Bik, Sarah Bosman, Luca Costantin, Alejandro Crespo Gómez, Dan Dicken, Andreas Eckart, Jens Hjorth, Iris Jermann, Alvaro Labiano, Danial Langeroodi , et al. (8 additional authors not shown)

Abstract: We present deep James Webb Space Telescope (JWST)/MIRI F560W observations of a flux-limited, ALMA-selected sample of 28 galaxies at z=0.5-3.6 in the Hubble Ultra Deep Field (HUDF). The data from the MIRI Deep Imaging Survey (MIDIS) reveal the stellar structure of the HUDF galaxies at rest-wavelengths of >1 micron for the first time. We revise the stellar mass estimates using new JWST photometry an… ▽ More We present deep James Webb Space Telescope (JWST)/MIRI F560W observations of a flux-limited, ALMA-selected sample of 28 galaxies at z=0.5-3.6 in the Hubble Ultra Deep Field (HUDF). The data from the MIRI Deep Imaging Survey (MIDIS) reveal the stellar structure of the HUDF galaxies at rest-wavelengths of >1 micron for the first time. We revise the stellar mass estimates using new JWST photometry and find good agreement with pre-JWST analysis; the few discrepancies can be explained by blending issues in the earlier lower-resolution Spitzer data. At z~2.5, the resolved rest-frame near-infrared (1.6 micron) structure of the galaxies is significantly more smooth and centrally concentrated than seen by HST at rest-frame 450 nm (F160W), with effective radii of Re(F560W)=1-5 kpc and Sérsic indices mostly close to an exponential (disk-like) profile (n~1), up to n~5 (excluding AGN). We find an average size ratio of Re(F560W)/Re(F160W)~0.7 that decreases with stellar mass. The stellar structure of the ALMA-selected galaxies is indistinguishable from a HUDF reference sample of galaxies with comparable MIRI flux density. We supplement our analysis with custom-made, position-dependent, empirical PSF models for the F560W observations. The results imply that a smoother stellar structure is in place in massive gas-rich, star-forming galaxies at Cosmic Noon, despite a more clumpy rest-frame optical appearance, placing additional constraints on galaxy formation simulations. As a next step, matched-resolution, resolved ALMA observations will be crucial to further link the mass- and light-weighted galaxy structures to the dusty interstellar medium. △ Less

Submitted 26 April, 2024; v1 submitted 31 August, 2023; originally announced August 2023.

Comments: 19 pages, 12 figures, 1 table, Accepted for publication in ApJ

On the evolution of the size of Lyman alpha halos across cosmic time: no change in the circumgalactic gas distribution when probed by line emission

Authors: Axel Runnholm, Matthew J. Hayes, Yu-Heng Lin, Jens Melinder, Claudia Scarlata, Angela Adamo, Ramona Augustin, Arjan Bik, Jérémy Blaizot, John M. Cannon, Sebastiano Cantalupo, Thibault Garel, Max Gronke, Edmund C. Herenz, Floriane Leclercq, Göran Östlin, Celine Peroux, Armin Rasekh, Michael J. Rutkowski, Anne Verhamme, Lutz Wisotzki

Abstract: Lyman $α$ (Ly$α$) is now routinely used as a tool for studying high-redshift galaxies and its resonant nature means it can trace neutral hydrogen around star-forming galaxies. Integral field spectrograph measurements of high-redshift Ly$α$ emitters indicate that significant extended Ly$α$ halo emission is ubiquitous around such objects. We present a sample of redshift 0.23 to 0.31 galaxies observe… ▽ More Lyman $α$ (Ly$α$) is now routinely used as a tool for studying high-redshift galaxies and its resonant nature means it can trace neutral hydrogen around star-forming galaxies. Integral field spectrograph measurements of high-redshift Ly$α$ emitters indicate that significant extended Ly$α$ halo emission is ubiquitous around such objects. We present a sample of redshift 0.23 to 0.31 galaxies observed with the Hubble Space Telescope selected to match the star formation properties of high-$z$ samples while optimizing the observations for detection of low surface brightness Ly$α$ emission. The Ly$α$ escape fractions range between 0.7\% and 37\%, and we detect extended Ly$α$ emission around six out of seven targets. We find Ly$α$ halo to UV scale length ratios around 6:1 which is marginally lower than high-redshift observations, and halo flux fractions between 60\% and 85\% -- consistent with high-redshift observations -- when using comparable methods. However, our targets show additional extended stellar UV emission: we parametrize this with a new double exponential model. We find that this parametrization does not strongly affect the observed Ly$α$ halo fractions. We find that deeper H$α$ data would be required to firmly determine the origin of Ly$α$ halo emission, however, there are indications that H$α$ is more extended than the central FUV profile, potentially indicating conditions favorable for the escape of ionizing radiation. We discuss our results in the context of high-redshift galaxies, cosmological simulations, evolutionary studies of the circumgalactic medium in emission, and the emission of ionizing radiation. △ Less

Submitted 25 April, 2023; originally announced April 2023.

Comments: 20 page, 14 figures, 6 tables. Accepted for publication in MNRAS

Unveiling the gravitationally unstable disc of a massive star-forming galaxy using NOEMA and MUSE

Authors: Johannes Puschnig, Matthew Hayes, Oscar Agertz, Eric Emsellem, John M. Cannon, Alexandra Le Reste, Jens Melinder, Göran Östlin, Christian Herenz, Veronica Menacho

Abstract: Using new high-resolution data of CO (2-1), H-alpha and H-beta obtained with the Northern Extended Millimeter Array (NOEMA) and the Multi-Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope, we have performed a Toomre-Q disc stability analysis and studied star formation, gas depletion times and other environmental parameters on sub-kpc scales within the z~0 galaxy SDSS J125013.84+073444… ▽ More Using new high-resolution data of CO (2-1), H-alpha and H-beta obtained with the Northern Extended Millimeter Array (NOEMA) and the Multi-Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope, we have performed a Toomre-Q disc stability analysis and studied star formation, gas depletion times and other environmental parameters on sub-kpc scales within the z~0 galaxy SDSS J125013.84+073444.5 (LARS 8). The galaxy hosts a massive, clumpy disc and is a proto-typical analogue of main-sequence galaxies at z~1-2. We show that the massive (molecular) clumps in LARS 8 are the result of an extremely gravitationally unstable gas disc, with large scale instabilities found across the whole extent of the rotating disc, with only the innermost 500 pc being stabilized by its bulgelike structure. The radial profiles further reveal that - contrary to typical disc galaxies - the molecular gas depletion time decreases from more than 1 Gyr in the center to less than ~100 Myr in the outskirts of the disc, supporting the findings of a Toomre-unstable disc. We further identified and analysed 12 individual massive molecular clumps. They are virialized and follow the mass-size relation, indicating that on local (cloud/clump) scales the stars form with efficiencies comparable to those in Milky Way clouds. The observed high star formation rate must thus be the result of triggering of cloud/clump formation over large scales due to disc instability. Our study provides evidence that "in-situ" massive clump formation (as also observed at high redshifts) is very efficiently induced by large-scale instabilities. △ Less

Submitted 24 March, 2023; originally announced March 2023.

Comments: Submitted to MNRAS

The Lyman $α$ Reference Sample XIV: Lyman $α$ imaging of 45 low redshift star-forming galaxies and inferences on global emission

Authors: Jens Melinder, Göran Östlin, Matthew Hayes, Armin Rasekh, J. Miguel Mas-Hesse, John M. Cannon, Daniel Kunth, Peter Laursen, Axel Runnholm, E. Christian Herenz, Matteo Messa, Daniel Schaerer, Anne Verhamme, T. Emil Rivera-Thorsen, Lucia Guaita, Thomas Marquart, Johannes Puschnig, Alexandra Le Reste, Andreas Sandberg, Emily Freeland, Joanna Bridge

Abstract: We present Ly $α$ imaging of 45 low redshift star-forming galaxies observed with the Hubble Space Telescope. The galaxies have been selected to have moderate to high star formation rates using far-ultraviolet (FUV) luminosity and \ha equivalent width criteria, but no constraints on Ly $α$ luminosity. We employ a pixel stellar continuum fitting code to obtain accurate continuum subtracted Ly $α$, H… ▽ More We present Ly $α$ imaging of 45 low redshift star-forming galaxies observed with the Hubble Space Telescope. The galaxies have been selected to have moderate to high star formation rates using far-ultraviolet (FUV) luminosity and \ha equivalent width criteria, but no constraints on Ly $α$ luminosity. We employ a pixel stellar continuum fitting code to obtain accurate continuum subtracted Ly $α$, H $α$ and H $β$ maps. We find that Ly $α$ is less concentrated than FUV and optical line emission in almost all galaxies with significant Ly $α$ emission. We present global measurements of Ly $α$ and other quantities measured in apertures designed to capture all of the Ly $α$ emission. We then show how the escape fraction of Ly $α$ relates to a number of other measured quantities (mass, metallicity, star formation, ionisation parameter, and extinction). We find that the escape fraction is strongly anti-correlated with nebular and stellar extinction, weakly anti-correlated with stellar mass, but no conclusive evidence for correlations with other quantities. We show that Ly $α$ escape fractions are inconsistent with common dust extinction laws, and discuss how a combination of radiative transfer effects and clumpy dust models can help resolve the discrepancies. We present a star formation rate calibration based on Ly $α$ luminosity, where the equivalent width of Ly $α$ is used to correct for non-unity escape fraction, and show that this relation provides a reasonably accurate SFR estimate. We also show stacked growth curves of Ly $α$ for the galaxies that can be used to find aperture loss fractions at a given physical radius. △ Less

Submitted 27 February, 2023; originally announced February 2023.

Comments: 48 pages, 21 figures, accepted for publication in The Astrophysical Journal Supplement Series. Multipanel figures for all 45 galaxies will be available in the online version of the published paper (and here: https://stockholmuniversity.box.com/s/ghvnvjyhxts64hpbqh90uwcj5dlczlhm ). Data products will be available as a MAST HLSP soon (placeholder DOI: 10.17909/pbe1-m743)

Life beyond 30: probing the -20<M_UV<-17 luminosity function at 8<z<13 with the NIRCam parallel field of the MIRI Deep Survey

Authors: Pablo G. Pérez-González, Luca Costantin, Danial Langeroodi, Pierluigi Rinaldi, Marianna Annunziatella, Olivier Ilbert, Luis Colina, Hans Ulrik Noorgaard-Nielsen, Thomas Greve, Göran Ostlin, Gillian Wright, Almudena Alonso-Herrero, Javier Álvarez-Márquez, Karina I. Caputi, Andreas Eckart, Olivier Le Fèvre, Álvaro Labiano, Macarena García-Marín, Jens Hjorth, Sarah Kendrew, John P. Pye, Tuomo Tikkanen, Paul van der Werf, Fabian Walter, Martin Ward , et al. (19 additional authors not shown)

Abstract: We present the ultraviolet luminosity function and an estimate of the cosmic star formation rate density at $8<z<13$ derived from deep NIRCam observations taken in parallel with the MIRI Deep Survey (MDS) of the Hubble Ultra Deep Field (HUDF), NIRCam covering the parallel field 2 (HUDF-P2). Our deep (40 hours) NIRCam observations reach a F277W magnitude of 30.8 ($5σ$), more than 2 magnitudes deepe… ▽ More We present the ultraviolet luminosity function and an estimate of the cosmic star formation rate density at $8<z<13$ derived from deep NIRCam observations taken in parallel with the MIRI Deep Survey (MDS) of the Hubble Ultra Deep Field (HUDF), NIRCam covering the parallel field 2 (HUDF-P2). Our deep (40 hours) NIRCam observations reach a F277W magnitude of 30.8 ($5σ$), more than 2 magnitudes deeper than JWST public datasets already analyzed to find high redshift galaxies. We select a sample of 44 $z>8$ galaxy candidates based on their dropout nature in the F115W and/or F150W filters, a high probability for their photometric redshifts, estimated with three different codes, being at $z>8$, good fits based on $χ^2$ calculations, and predominant solutions compared to $z<8$ alternatives. We find mild evolution in the luminosity function from $z\sim13$ to $z\sim8$, i.e., only a small increase in the average number density of $\sim$0.2 dex, while the faint-end slope and absolute magnitude of the knee remain approximately constant, with values $α=-2.2\pm0.1$ and $M^*=-20.8\pm0.2$ mag. Comparing our results with the predictions of state-of-the-art galaxy evolution models, we find two main results: (1) a slower increase with time in the cosmic star formation rate density compared to a steeper rise predicted by models; (2) nearly a factor of 10 higher star formation activity concentrated in scales around 2 kpc in galaxies with stellar masses $\sim10^8$ M$_\odot$ during the first 350 Myr of the Universe, $z\sim12$, with models matching better the luminosity density observational estimations $\sim$150 Myr later, by $z\sim9$. △ Less

Submitted 22 May, 2023; v1 submitted 5 February, 2023; originally announced February 2023.

Comments: Accepted for publication in ApJL, in memoriam Hans Ulrik Noorgaard-Nielsen and Olivier Le Fèvre

MIDIS: Strong (Hb + [OIII]) and Ha emitters at redshift $z \simeq 7-8$ unveiled with JWST/NIRCam and MIRI imaging in the Hubble eXtreme Deep Field (XDF)

Authors: P. Rinaldi, K. I. Caputi, L. Costantin, S. Gillman, E. Iani, P. G. Perez Gonzalez, G. Oestlin, L. Colina, T. Greve, H. U. Noorgard-Nielsen, G. S. Wright, A. Alonso-Herrero, J. Alvarez-Marquez, A. Eckart, M. Garcia-Marin, J. Hjorth, O. Ilbert, S. Kendrew, A. Labiano, O. Le Fevre, J. Pye, T. Tikkanen, F. Walter, P. van der Werf, M. Ward , et al. (22 additional authors not shown)

Abstract: We make use of \textit{JWST} medium and broad-band NIRCam imaging, along with ultra-deep MIRI $5.6 \rm μm$ imaging, in the Hubble eXtreme Deep Field (XDF) to identify prominent line emitters at $z\simeq 7-8$. Out of a total of 58 galaxies at $z\simeq 7-8$, we find 18 robust candidates ($\simeq$31\%) for (H$β$ + [OIII]) emitters, based on their enhanced fluxes in the F430M and F444W filters, with E… ▽ More We make use of \textit{JWST} medium and broad-band NIRCam imaging, along with ultra-deep MIRI $5.6 \rm μm$ imaging, in the Hubble eXtreme Deep Field (XDF) to identify prominent line emitters at $z\simeq 7-8$. Out of a total of 58 galaxies at $z\simeq 7-8$, we find 18 robust candidates ($\simeq$31\%) for (H$β$ + [OIII]) emitters, based on their enhanced fluxes in the F430M and F444W filters, with EW$_{0}$(H$β$ +[OIII]) $\simeq 87 - 2100$ Å. Among these emitters, 16 lie in the MIRI coverage area and 12 exhibit a clear flux excess at $5.6 \, \rm μm$, indicating the simultaneous presence of a prominent H$α$ emission line with EW$_{0}$(H$α$) $\simeq 200-3000$ Å. This is the first time that H$α$ emission can be detected in individual galaxies at $z>7$. The H$α$ line, when present, allows us to separate the contributions of H$β$ and [OIII] to the (H$β$ +[OIII]) complex, and derive H$α$-based star formation rates (SFRs). We find that in most cases [OIII]/H$β> 1$. Instead, two galaxies have [OIII]/H$β< 1$, indicating that the NIRCam flux excess is mainly driven by H$β$. This could potentially imply extremely low metallicities. Most prominent line emitters are very young starbursts or galaxies on their way to/from the starburst cloud. They make for a cosmic SFR density $\rm log_{10}(ρ_{SFR_{Hα}}) \simeq -2.35$, which is about a quarter of the total value ($\rm log_{10}(ρ_{SFR_{tot}}) \simeq -1.76$) at $z\simeq 7-8$. Therefore, the strong H$α$ emitters likely had a significant role in reionization. △ Less

Submitted 5 June, 2023; v1 submitted 25 January, 2023; originally announced January 2023.

Comments: 18 pages, 11 figures. The paper has been accepted for publication in ApJ

Tidally offset neutral gas in Lyman continuum emitting galaxy Haro 11

Authors: Alexandra Le Reste, John M. Cannon, Matthew J. Hayes, John L. Inoue, Amanda A. Kepley, Jens Melinder, Veronica Menacho, Angela Adamo, Arjan Bik, Timmy Ejdetjärn, Gyula I. G. Józsa, Göran Östlin, Sarah H. Taft

Abstract: Around 400 million years after the Big Bang, the ultraviolet emission from star-forming galaxies reionized the Universe. Ionizing radiation (Lyman Continuum, LyC) is absorbed by cold neutral hydrogen gas (HI) within galaxies, hindering the escape of LyC photons. Since the HI reservoir of LyC emitters has never been mapped, major uncertainties remain on how LyC photons escape galaxies and ionize th… ▽ More Around 400 million years after the Big Bang, the ultraviolet emission from star-forming galaxies reionized the Universe. Ionizing radiation (Lyman Continuum, LyC) is absorbed by cold neutral hydrogen gas (HI) within galaxies, hindering the escape of LyC photons. Since the HI reservoir of LyC emitters has never been mapped, major uncertainties remain on how LyC photons escape galaxies and ionize the intergalactic medium. We have directly imaged the neutral gas in the nearby reionization-era analog galaxy Haro 11 with the 21cm line to identify the mechanism enabling ionizing radiation escape. We find that merger-driven interactions have caused a bulk offset of the neutral gas by about 6 kpc from the center of the galaxy, where LyC emission production sites are located. This could facilitate the escape of ionizing radiation into our line of sight. Galaxy interactions can cause both elevated LyC production and large-scale displacement of HI from the regions where these photons are produced. They could contribute to the anisotropic escape of LyC radiation from galaxies and the reionization of the Universe. We argue for a systematic assessment of the effect of environment on LyC production and escape. △ Less

Submitted 19 December, 2023; v1 submitted 6 January, 2023; originally announced January 2023.

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

Spatially resolved gas and stellar kinematics in compact starburst galaxies

Authors: Arjan Bik, Göran Östlin, Matthew Hayes, Jens Melinder, Veronica Menacho

Abstract: The kinematics of galaxies provide valuable insights in their physics and assembly history. Kinematics are governed not only by the gravitational potential, but also by merger events and stellar feedback processes such as stellar winds and supernova explosions. Aims. We aim at identifying what governs the kinematics in a sample of SDSS selected nearby starburst galaxies, by obtaining spatially res… ▽ More The kinematics of galaxies provide valuable insights in their physics and assembly history. Kinematics are governed not only by the gravitational potential, but also by merger events and stellar feedback processes such as stellar winds and supernova explosions. Aims. We aim at identifying what governs the kinematics in a sample of SDSS selected nearby starburst galaxies, by obtaining spatially resolved measurements of the gas and stellar kinematics. We obtain near-infrared integral field K-band spectroscopy with VLT/SINFONI of 15 compact starburst galaxies. We derive the integrated as well as spatially resolved stellar and gas kinematics. The stellar kinematics are derived from the CO absorption bands, and Pa$α$ and Br$γ$ emission lines are used for the gas kinematics. Based on the integrated spectra we find that the majority of galaxies have gas and stellar velocity dispersion that are comparable. A spatially resolved comparison shows that the six galaxies that deviate show evidence for a bulge or stellar feedback. Two galaxies are identified as mergers based on their double peaked emission lines. In our sample, we find a negative correlation between the ratio of the rotational velocity over the velocity dispersion (vrot/$σ$) and the star formation rate surface density. We propose a scenario where the global kinematics of the galaxies are determined by gravitational instabilities that affect both the stars and gas. This process could be driven by mergers or accretion events. Effects of stellar feedback on the ionised gas are more localised and detected only in the spatially resolved analysis. The mass derived from the velocity dispersion provides a reliable mass even if the galaxy cannot be spatially resolved. The technique used in this paper is applicable to galaxies at low and high redshift with the next generation of infrared focussed telescopes (JWST, ELT). △ Less

Submitted 1 September, 2022; v1 submitted 23 August, 2022; originally announced August 2022.

Comments: 26 pages, 31 figures. Accepted for publication in Astronomy and Astrophysics, corrected by language editor

Journal ref: A&A 666, A161 (2022)

A bottom-up search for Lyman-continuum leakage in the Hubble Ultra Deep Field

Authors: T. Emil Rivera-Thorsen, Matthew Hayes, Jens Melinder

Abstract: Context: Studies of the production and escape of Lyman Continuum from galaxies often rely on array of indirect observational tracers in preselection of candidate leakers. Aims: Here, we investigate how much ionizing radiation might be missed due to these selection criteria by completely removing them and performing a search selected purely from rest-frame LyC emission; and how that affects our e… ▽ More Context: Studies of the production and escape of Lyman Continuum from galaxies often rely on array of indirect observational tracers in preselection of candidate leakers. Aims: Here, we investigate how much ionizing radiation might be missed due to these selection criteria by completely removing them and performing a search selected purely from rest-frame LyC emission; and how that affects our estimates of the ionizing background. Methods: We invert the conventional method and perform a bottom-up search for Lyman-continuum leaking galaxies at redshifts $2 < z < 3.5$. Using archival data from HST and VLT/MUSE, we run source finding software on UV-filter HST images from the HUDF, and subject all detected sources to a series of tests to eliminate those that are inconsistent with being ionizing sources. Results: We find 6 new and one previously identified candidate leakers with absolute escape fractions ranging from 36% to 100%. Our filtering criteria eliminate one object previously reported as a candidate ionizing emitter in the literature, while we report non-detection in the rest frame Lyman continuum of two other previously reported sources. We find that our candidates make a contribution to the metagalactic ionizing field of $\log_{10}(ε_ν) = 25.32(+0.25)(-0.21)$ and $25.29(+0.27)(-0.22)$ erg/s/Hz/cMpc$^3$ for the full set of candidates and for the 4 strongest candidates only; both values are higher than but consistent with other recent figures in the literature. Conclusions: Our findings suggest that galaxies that do not meet the usual selection criteria may make a non-negligible contribution to the cosmic ionizing field. We recommend that similar searches be carried out on a larger scale in well-studied fields with both UV and large ancillary data coverage, for example in the full set of CANDELS fields. △ Less

Submitted 5 August, 2022; v1 submitted 21 June, 2022; originally announced June 2022.

Comments: 22 pages, 8 tables, 12 figures. This is the version accepted for publication in A&A. Contains only minor edits from v1

Journal ref: A&A 666, A145 (2022)

The Lyman Alpha Reference Sample XIII: High-Angular Resolution 21cm HI observations of Ly$α$ emitting galaxies

Authors: A. Le Reste, M. Hayes, J. M. Cannon, E. C. Herenz, J. Melinder, V. Menacho, G. Östlin, J. Puschnig, T. E. Rivera-Thorsen, D. Kunth, N. Velikonja

Abstract: The Ly$α$ emission line is one of the main observables of galaxies at high redshift, but its output depends strongly on the neutral gas distribution and kinematics around the star-forming regions where UV photons are produced. We present observations of Ly$α$ and 21-cm HI emission at comparable scales with the goal to qualitatively investigate how the neutral interstellar medium (ISM) properties i… ▽ More The Ly$α$ emission line is one of the main observables of galaxies at high redshift, but its output depends strongly on the neutral gas distribution and kinematics around the star-forming regions where UV photons are produced. We present observations of Ly$α$ and 21-cm HI emission at comparable scales with the goal to qualitatively investigate how the neutral interstellar medium (ISM) properties impact Ly$α$ transfer in galaxies. We have observed 21-cm HI at the highest angular resolution possible (~ 3" beam) with the VLA in two local galaxies from the Lyman Alpha Reference Sample. We contrast this data with HST Ly$α$ imaging and spectroscopy, and MUSE and PMAS ionized gas observations. In LARS08, high intensity Ly$α$ emission is co-spatial with high column density HI where dust content is the lowest. The Ly$α$ line is strongly redshifted, consistent with velocity redistribution which allows Ly$α$ escape from high column density neutral medium with low dust content. In eLARS01, high intensity Ly$α$ emission is located in regions of low column density HI, below the HI data sensitivity limit ($<2\times10^{20}\,$cm$^{-2}$). The perturbed ISM distribution with low column density gas in front of the Ly$α$ emission region plays an important role in the escape. In both galaxies, the faint Ly$α$ emission ($\sim 1\times10^{-16}$erg.s$^{-1}$cm$^{-2}$arcsec$^{-2}$) traces intermediate H$α$ emission regions where HI is found, regardless of the dust content. Dust seems to modulate, but not prevent, the formation of a faint Ly$α$ halo. This study suggests the existence of scaling relations between dust, H$α$, HI, and Ly$α$ emission in galaxies. △ Less

Submitted 13 June, 2022; originally announced June 2022.

Comments: 28 pages, 14 figures, accepted for publication in ApJ

The Lyman Alpha Reference Sample XII: Morphology of extended Lyman alpha emission in star-forming galaxies

Authors: Armin Rasekh, Jens Melinder, Göran Östlin, Matthew Hayes, Edmund. C. Herenz, Axel Runnholm, Daniel Kunth, J. Miguel Mas Hesse, Anne Verhamme, John M. Cannon

Abstract: We use Hubble space telescope data of 45 nearby star-forming galaxies to investigate properties of Lyman-alpha (Ly$α$) halos, Ly$α$ morphology, and the star-forming characteristics of galaxies. We study how the morphology of Ly$α$ emission is related to other Ly$α$ observables. Furthermore, we study the interdependencies of Ly$α$ morphological quantities. We studied the spatial extent of Ly$α$ usi… ▽ More We use Hubble space telescope data of 45 nearby star-forming galaxies to investigate properties of Lyman-alpha (Ly$α$) halos, Ly$α$ morphology, and the star-forming characteristics of galaxies. We study how the morphology of Ly$α$ emission is related to other Ly$α$ observables. Furthermore, we study the interdependencies of Ly$α$ morphological quantities. We studied the spatial extent of Ly$α$ using surface brightness profiles in i) using circular apertures and ii) within faint Ly$α$ isophotes. We measured the average intensity and the size of the regions with a high star formation rate density. The morphology of the galaxies was quantified by computing centroid position, axis ratio, and position angle in the Ly$α$, ultraviolet continuum, and I band maps. We found that galaxies with more extended star-forming regions possess larger Ly$α$ halos. Furthermore, galaxies with more elongated Ly$α$ morphology are also more extended in Ly$α$. Our data suggest that Ly$α$ bright galaxies appear rounder in their Ly$α$ morphology, and there is less of a contribution from their Ly$α$ halo to their overall luminosity. We compared our results with studies at high redshift and found that whilst the Ly$α$ extent in the inner regions of the galaxies in our sample are similar to the high$-z$ Ly$α$ emitters (LAEs), Ly$α$ halos are more extended in high$-z$ LAEs. Our analysis suggests that the Ly$α$ morphology affects the measurement of other observable quantities concerning Ly$α$ emission, and some of the conclusions drawn from high redshift LAEs might be biased towards galaxies with specific Ly$α$ shapes. In particular, faint Ly$α$ emitters have larger Ly$α$ scale lengths and halo fractions. This implies that faint Ly$α$ emitters are harder to detect at high redshift than previously believed. △ Less

Submitted 4 October, 2021; originally announced October 2021.

Journal ref: A&A 662, A64 (2022)

Ionized gas properties of the extreme starburst galaxy Haro 11. Temperature and metal abundance discrepancies

Authors: V. Menacho, G. Östlin A. Bik, A. Adamo, N. Bergvall, L. Della Bruna, M. Hayes, J. Melinder, E. T. Rivera-Thorsen

Abstract: We use high quality VLT/MUSE data to study the kinematics and the ionized gas properties of Haro 11, a well known starburst merger system and the closest confirmed Lyman continuum leaking galaxy. We present results from integrated line maps, and from maps in three velocity bins comprising the blueshifted, systemic and redshifted emission. The kinematic analysis reveals complex velocities resulting… ▽ More We use high quality VLT/MUSE data to study the kinematics and the ionized gas properties of Haro 11, a well known starburst merger system and the closest confirmed Lyman continuum leaking galaxy. We present results from integrated line maps, and from maps in three velocity bins comprising the blueshifted, systemic and redshifted emission. The kinematic analysis reveals complex velocities resulting from the interplay of virial motions and momentum feedback. Star formation happens intensively in three compact knots (knots A, B and C), but one, knot C, dominates the energy released in supernovae. The halo is characterised by low gas density and extinction, but with large temperature variations, coincident with fast shock regions. Moreover, we find large temperature discrepancies in knot C, when using different temperature-sensitive lines. The relative impact of the knots in the metal enrichment differs. While knot B is strongly enriching its closest surrounding, knot C is likely the main distributor of metals in the halo. In knot A, part of the metal enriched gas seems to escape through low density channels towards the south. We compare the metallicities from two methods and find large discrepancies in knot C, a shocked area, and the highly ionized zones, that we partially attribute to the effect of shocks. This work shows, that traditional relations developed from averaged measurements or simplified methods, fail to probe the diverse conditions of the gas in extreme environments. We need robust relations that include realistic models where several physical processes are simultaneously at work. △ Less

Submitted 23 May, 2021; originally announced May 2021.

Comments: Paper accepted for publication in MNRAS

Spatially resolved CIII]$λ$1909 emission in Haro 11

Authors: Genoveva Micheva, Göran Östlin, Jens Melinder, Matthew Hayes, M. S. Oey, Akio K. Inoue, Ikuru Iwata, Angela Adamo, Lutz Wisotzki, Kimihiko Nakajima

Abstract: The CIII]1909 (hereafter, C III]) line is the strongest ultraviolet emission line after Ly$α$ and is therefore of interest to high redshift studies of star-forming (SF) galaxies near the epoch of reionization. It is thought that C III] emission is strongest in galaxies with subsolar metallicity and low mass, however, spectral observations of numerous such galaxies at high and low redshift produce… ▽ More The CIII]1909 (hereafter, C III]) line is the strongest ultraviolet emission line after Ly$α$ and is therefore of interest to high redshift studies of star-forming (SF) galaxies near the epoch of reionization. It is thought that C III] emission is strongest in galaxies with subsolar metallicity and low mass, however, spectral observations of numerous such galaxies at high and low redshift produce inconclusive or even contradictory results. We present the first ever C III ] imaging, obtained with HST/STIS for the low-redshift SF galaxy Haro 11. Cluster parameters like stellar mass, dust fraction and dust attenuation, and ionization parameter, obtained through spectral energy distribution fitting, show no correlation with the CIII] equivalent width (EW), which may be due to a combination of the limitation of the models and the age-homogeneity of the cluster population. Comparing the ratio of C III] emission line flux from individual clusters to that of H$α$, [O III], and [O II] we find that the clusters with the highest EW(C III]) can be reconciled only with Cloudy models with an extremely high C/O ratio of $\ge$ 1.4(C/O)$_\odot$ for an ionizing population of single stars, binary stars, or a mixture of binary stars and active galactic nuclei. Given the point-like nature of strong C III], the integrated total strength of EW(C III]) becomes dependent on the morphology of the galaxy, which would explain the large scatter in EW(C III]) strengths, observed in galaxies with otherwise similar SF properties, and of similarly low metallicity and stellar mass. △ Less

Submitted 9 October, 2020; originally announced October 2020.

Comments: Accepted for publication in ApJ

Star cluster formation in the most extreme environments: Insights from the HiPEEC survey

Authors: A. Adamo, K. Hollyhead, M. Messa, J. E. Ryon, V. Bajaj, A. Runnholm, S. Aalto, D. Calzetti, J. S. Gallagher, M. J. Hayes, J. M. D. Kruijssen, S. König, S. S. Larsen, J. Melinder, E. Sabbi, L. J. Smith, G. Östlin

Abstract: We present the Hubble imaging Probe of Extreme Environments and Clusters (HiPEEC) survey. We fit HST NUV to NIR broadband and H$α$ fluxes, to derive star cluster ages, masses, extinctions and determine the star formation rate (SFR) of 6 merging galaxies. These systems are excellent laboratories to trace cluster formation under extreme gas physical conditions, rare in the local universe, but typica… ▽ More We present the Hubble imaging Probe of Extreme Environments and Clusters (HiPEEC) survey. We fit HST NUV to NIR broadband and H$α$ fluxes, to derive star cluster ages, masses, extinctions and determine the star formation rate (SFR) of 6 merging galaxies. These systems are excellent laboratories to trace cluster formation under extreme gas physical conditions, rare in the local universe, but typical for star-forming galaxies at cosmic noon. We detect clusters with ages of 1-500 Myr and masses that exceed $10^7$ M$_\odot$. The recent cluster formation history and their distribution within the host galaxies suggest that systems like NGC34, NGC1614, NGC4194 are close to their final coalescing phase, while NGC3256, NGC3690, NGC6052 are at an earlier/intermediate stage. A Bayesian analysis of the cluster mass function in the age interval 1-100 Myr provides strong evidence in 4 of the 6 galaxies that an exponentially truncated power law better describes the observed mass distributions. For two galaxies, the fits are inconclusive due to low number statistics. We determine power-law slopes $β\sim-1.5$ to $-2.0$, and truncation masses, M$_c$, between $10^6$ and a few times $10^7$ M$_\odot$, among the highest values reported in the literature. Advanced mergers have higher M$_c$ than early/intermediate merger stage galaxies, suggesting rapid changes in the dense gas conditions during the merger. We compare the total stellar mass in clusters to the SFR of the galaxy, finding that these systems are among the most efficient environments to form star clusters in the local universe. △ Less

Submitted 28 August, 2020; originally announced August 2020.

Comments: Accepted for publication in MNRAS; data and catalogues will be released at this URL http://dx.doi.org/10.17909/t9-cn0b-ht83

Lyman Alpha Reference Sample: X. Predicting Lyman alpha output from starforming galaxies using multivariate regression

Authors: Axel Runnholm, Matthew Hayes, Jens Melinder, Emil Rivera-Thorsen, Göran Östlin, John Cannon, Daniel Kunth

Abstract: Understanding the production and escape of Lyman $α$ (Ly$α$) radiation from star-forming galaxies is a long standing problem in astrophysics. The ability to predict the Ly$α$ luminosity of galaxies would open up new ways of exploring the Epoch of Reionization (EoR), and to estimate Ly$α$ emission from galaxies in cosmological simulations where radiative transfer calculations cannot be done. We app… ▽ More Understanding the production and escape of Lyman $α$ (Ly$α$) radiation from star-forming galaxies is a long standing problem in astrophysics. The ability to predict the Ly$α$ luminosity of galaxies would open up new ways of exploring the Epoch of Reionization (EoR), and to estimate Ly$α$ emission from galaxies in cosmological simulations where radiative transfer calculations cannot be done. We apply multivariate regression methods to the Lyman Alpha Reference Sample dataset to obtain a relation between the galaxy properties and the emitted Ly$α$. The derived relation predicts the Ly$α$ luminosity of our galaxy sample to good accuracy, regardless of whether we consider only direct observables (root-mean-square (RMS) dispersion around the relation of $\sim 0.19$ dex) or derived physical quantities (RMS $\sim 0.27$ dex). We confirm the predictive ability on a separate sample of compact star-forming galaxies and find that the prediction works well, but that aperture effects on measured Ly$α$ luminosity may be important, depending on the redshift of the galaxy. We apply statistical feature selection techniques to determine an order of importance of the variables in our dataset, enabling future observations to be optimized for predictive ability. When using physical variables, we are able to determine that the most important predictive parameters are, in order, star formation rate, dust extinction, compactness and the gas covering fraction. We discuss the application of our results in terms of studying the EoR and intensity mapping experiments. △ Less

Submitted 27 February, 2020; originally announced February 2020.

Comments: Accepted for publication in ApJ. 23 pages, 6 figures

Accretion-Driven Sources in Spatially Resolved Ly$α$ Emitters

Authors: B. Dittenber, M. S. Oey, E. Hodges-Kluck, E. Gallo, M. Hayes, G. Oestlin, J. Melinder

Abstract: Ly$α$ emission is a standard tracer of starburst galaxies at high redshift. However, a number of local Ly$α$ emitters (LAEs) are X-ray sources, suggesting a possible origin of Ly$α$ photons other than young, hot stars, and which may be active at much later ages relative to the parent starburst. Resolved, nearby LAEs offer the opportunity to discriminate between diffuse X-ray emission arising from… ▽ More Ly$α$ emission is a standard tracer of starburst galaxies at high redshift. However, a number of local Ly$α$ emitters (LAEs) are X-ray sources, suggesting a possible origin of Ly$α$ photons other than young, hot stars, and which may be active at much later ages relative to the parent starburst. Resolved, nearby LAEs offer the opportunity to discriminate between diffuse X-ray emission arising from supernova-heated gas, high-mass X-ray binaries (HMXBs), or low-luminosity active galactic nuclei (LLAGN). We examine archival X-ray imaging from Chandra and XMM-Newton for 11 galaxies with spatially resolved Ly$α$ imaging to determine the luminosity, morphology, and spectral hardness of the X-ray sources. The data are consistent with 9 of the 12, bright Ly$α$ sources being driven by luminous, $10^{40}$ erg s$^{-1}$ X-ray sources. Half of the 8 Chandra sources are unresolved. The data suggest that nuclear activity, whether from LLAGN or nuclear starbursts, may play an important role in Ly$α$ emission. Our results also suggest a significant link between Ly$α$ emission and HMXBs, ULXs, and/or LLAGN, which would imply that Ly$α$ may be generated over timescales 1 - 2 orders of magnitude longer than produced by photoionization from OB stars. This highlights a critical need to quantify the relative contributions of different sources across cosmic time, to interpret Ly$α$ observations and the resulting properties of distant galaxies. △ Less

Submitted 3 February, 2020; v1 submitted 21 January, 2020; originally announced January 2020.

Comments: 10 pages, 3 figures, 2 tables. Accepted by ApJ Letters

Star-forming clumps in the Lyman Alpha Reference Sample of galaxies -- I. Photometric analysis and clumpiness

Authors: Matteo Messa, Angela Adamo, Göran Östlin, Jens Melinder, Matthew Hayes, Johanna S. Bridge, John Cannon

Abstract: We study young star-forming clumps on physical scales of 10-500 pc in the Lyman-Alpha Reference Sample (LARS), a collection of low-redshift (z = 0.03-0.2) UV-selected star-forming galaxies. In each of the 14 galaxies of the sample, we detect clumps for which we derive sizes and magnitudes in 5 UV-optical filters. The final sample includes $\sim$1400 clumps, of which $\sim$600 have magnitude uncert… ▽ More We study young star-forming clumps on physical scales of 10-500 pc in the Lyman-Alpha Reference Sample (LARS), a collection of low-redshift (z = 0.03-0.2) UV-selected star-forming galaxies. In each of the 14 galaxies of the sample, we detect clumps for which we derive sizes and magnitudes in 5 UV-optical filters. The final sample includes $\sim$1400 clumps, of which $\sim$600 have magnitude uncertainties below 0.3 in all filters. The UV luminosity function for the total sample of clumps is described by a power-law with slope $α= -2.03^{+0.11}_{-0.13}$. Clumps in the LARS galaxies have on average $Σ_{SFR}$ values higher than what observed in HII regions of local galaxies and comparable to typical SFR densities of clumps in z = 1-3 galaxies. We derive the clumpiness as the relative contribution from clumps to the UV emission of each galaxy, and study it as a function of galactic-scale properties, i.e. $Σ_{SFR}$ and the ratio between rotational and dispersion velocities of the gas ($v_s/σ_0$). We find that in galaxies with higher $Σ_{SFR}$ or lower $v_s/σ_0$, clumps dominate the UV emission of their host systems. All LARS galaxies with Ly$α$ escape fractions larger than 10% have more than 50% of the UV luminosity from clumps. We tested the robustness of these results against the effect of different physical resolutions. At low resolution, the measured clumpiness appears more elevated than if we could resolve clumps down to single clusters. This effect is small in the redshift range covered by LARS, thus our results are not driven by the physical resolution. △ Less

Submitted 11 June, 2019; originally announced June 2019.

Comments: 20 pages, 11 figures (+appendix: 6 pages, 6 figures). Accepted for publication in MNRAS

The impact of Stellar feedback from velocity-dependent ionised gas maps. -- A MUSE view of Haro 11

Authors: V. Menacho, G. Östlin, A. Bik, L. Della Bruna, J. Melinder, A. Adamo, M. Hayes, E. C. Herenz, N. Bergvall

Abstract: We have used the capability of the MUSE instrument to explore the impact of stellar feedback at large scales in Haro 11, a galaxy under extreme starburst condition and one of the first galaxies where Lyman continuum (LyC) has been detected. Using Ha, [OIII] and [OI] emission lines from deep MUSE observations, we have constructed a sequence of velocity-dependent maps of the Ha emission, the state o… ▽ More We have used the capability of the MUSE instrument to explore the impact of stellar feedback at large scales in Haro 11, a galaxy under extreme starburst condition and one of the first galaxies where Lyman continuum (LyC) has been detected. Using Ha, [OIII] and [OI] emission lines from deep MUSE observations, we have constructed a sequence of velocity-dependent maps of the Ha emission, the state of the ionised gas and a tracer of fast shocks. These allowed us to investigate the ionisation structure of the galaxy in 50 kms^2 bins over a velocity range of -400 to 350 kms. The ionised gas in Haro 11 is assembled by a rich arrangement of structures, such as superbubbles, filaments, arcs and galactic ionised channels, whose appearances change drastically with velocity. The central star forming knots and the star forming dusty arm are the main engines that power the strong mechanical feedback in this galaxy, although with different impact on the ionisation structure. Haro 11 appears to leak LyC radiation in many directions. We found evidence of a kpc-scale fragmented superbubble, that may have cleared galactic-scale channels in the ISM. Additionally, the southwestern hemisphere is highly ionised in all velocities, hinting at a density bound scenario. A compact kpc-scale structure of lowly ionised gas coincides with the diffuse Lya emission and the presence of fast shocks. Finally, we find evidence that a significant fraction of the ionised gas mass may escape the gravitational potential of the galaxy. △ Less

Submitted 15 July, 2019; v1 submitted 27 March, 2019; originally announced March 2019.

Comments: Pubisched version. Ionisation values were corrected after discovering a bug in a code used to generate the maps

Super star cluster feedback driving ionization, shocks and outflows in the halo of the nearby starburst ESO 338-IG04

Authors: A. Bik, G. Östlin, V. Menacho, A. Adamo, M. Hayes, E. C. Herenz, J. Melinder

Abstract: Stellar feedback strongly affects the interstellar medium (ISM) of galaxies. Stellar feedback in the first galaxies likely plays a major role in enabling the escape of LyC photons, which contribute to the re-ionization of the Universe. Nearby starburst galaxies serve as local analogues allowing for a spatially resolved assessment of the feedback processes in these galaxies. We characterize the fee… ▽ More Stellar feedback strongly affects the interstellar medium (ISM) of galaxies. Stellar feedback in the first galaxies likely plays a major role in enabling the escape of LyC photons, which contribute to the re-ionization of the Universe. Nearby starburst galaxies serve as local analogues allowing for a spatially resolved assessment of the feedback processes in these galaxies. We characterize the feedback effects from the star clusters in the local high-redshift analogue ESO 338-IG04 on the ISM and compare the results with the properties of the most massive clusters. We use high quality VLT/MUSE optical integral field data to derive the physical properties of the ISM such as ionization, density, shocks, and perform new fitting of the spectral energy distributions of the brightest clusters in ESO 338-IG04 from HST imaging. ESO 338-IG04 has a large ionized halo which we detect to a distance of 9 kpc. We identify 4 Wolf-Rayet (WR) clusters based on the blue and red WR bump. We follow previously identified ionization cones and find that the ionization of the halo increases with distance. Analysis of the galaxy kinematics shows two complex outflows driven by the numerous young clusters in the galaxy. We find a ring of shocked emission traced by an enhanced [OI]/H$α$ ratio surrounding the starburst and at the end of the outflow. Finally we detect nitrogen enriched gas associated with the outflow, likely caused by the WR stars in the massive star clusters. Photo-ionization dominates the central starburst and sets the ionization structure of the entire halo, resulting in a density bounded halo, facilitating the escape of LyC photons. Outside the central starburst, shocks triggered by an expanding super bubble become important. The shocks at the end of the outflow suggest interaction between the hot outflowing material and the more quiescent halo gas. △ Less

Submitted 10 September, 2018; originally announced September 2018.

Comments: Accepted for publication in Astronomy and Astrophysics, 22 pages, 15 figures

Journal ref: A&A 619, A131 (2018)

The Lyman Alpha Reference Sample IX: Revelations from deep surface photometry

Authors: Genoveva Micheva, Göran Östlin, Erik Zackrisson, Matthew Hayes, Jens Melinder, Lucia Guaita, John M. Cannon, Joanna S. Bridge, Daniel Kunth, Andreas Sandberg

Abstract: The Lyman Alpha Reference Sample (LARS) of 14 star-forming galaxies offers a wealth of insight into the workings of these local analogs to high-redshift star-forming galaxies. The sample has been well-studied in terms of LyA and other emission line properties, such as HI mass, gas kinematics, and morphology. We analyze deep surface photometry of the LARS sample in UBIK broadband imaging obtained a… ▽ More The Lyman Alpha Reference Sample (LARS) of 14 star-forming galaxies offers a wealth of insight into the workings of these local analogs to high-redshift star-forming galaxies. The sample has been well-studied in terms of LyA and other emission line properties, such as HI mass, gas kinematics, and morphology. We analyze deep surface photometry of the LARS sample in UBIK broadband imaging obtained at the Nordic Optical Telescope and the Canada-France-Hawaii Telescope, and juxtaposition their derived properties with a sample of local high-redshift galaxy analogs, namely, with blue compact galaxies (BCGs). We construct radial surface brightness and color profiles with both elliptical and isophotal integration, as well as RGB images, deep contours, color maps, a burst fraction estimate, and a radial mass-to-light ratio profile for each LARS galaxy. Standard morphological parameters like asymmetry, clumpiness, the Gini and M20 coefficients are [...] analyzed, as well as isophotal asymmetry profiles for each galaxy. [...] We compare the LARS to the properties of the BCG sample and highlight the differences. Several diagnostics indicate that the LARS galaxies have highly disturbed morphologies even at the level of the faintest isophotes [...]. The ground-based photometry [...] reveals previously unexplored isophotes [...]. The burst fraction estimate suggests a spatially more extended burst region in LARS than in the BCGs. [...] The galaxies in the LARS sample appear to be in earlier stages of a merger event compared to the BCGs. Standard morphological diagnostics like asymmetry, clumpiness, Gini and M20 coefficients cannot separate the two samples, although an isophotal asymmetry profile successfully captures the average difference in morphology. These morphological diagnostics do not show any correlation with the equivalent width or the escape fraction of Lyman Alpha. [abridged] △ Less

Submitted 19 March, 2018; originally announced March 2018.

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

Journal ref: A&A 615, A46 (2018)

The Lyman Alpha Reference Sample. VIII. Characterizing Lyman-Alpha Scattering in Nearby Galaxies

Authors: Joanna S. Bridge, Matthew Hayes, Jens Melinder, Göran Östlin, Caryl Gronwall, Robin Ciardullo, Hakim Atek, John M. Cannon, Max Gronke, Lucia Guaita, Alex Hagen, Edmund Christian Herenz, Daniel Kunth, Peter Laursen, J. Miguel Mas-Hesse, Stephen A. Pardy

Abstract: We examine the dust geometry and Lyα scattering in the galaxies of the Lyman Alpha Reference Sample (LARS), a set of 14 nearby (0.02 < $z$ < 0.2) Lyα emitting and starbursting systems with Hubble Space Telescope Lyα, Hα, and Hβ imaging. We find that the global dust properties determined by line ratios are consistent with other studies, with some of the LARS galaxies exhibiting clumpy dust media wh… ▽ More We examine the dust geometry and Lyα scattering in the galaxies of the Lyman Alpha Reference Sample (LARS), a set of 14 nearby (0.02 < $z$ < 0.2) Lyα emitting and starbursting systems with Hubble Space Telescope Lyα, Hα, and Hβ imaging. We find that the global dust properties determined by line ratios are consistent with other studies, with some of the LARS galaxies exhibiting clumpy dust media while others of them show significantly lower Lyα emission compared to their Balmer decrement. With the LARS imaging, we present Lyα/Hα and Hα/Hβ maps with spatial resolutions as low as $\sim$ 40 pc, and use these data to show that in most galaxies, the dust geometry is best modeled by three distinct regions: a central core where dust acts as a screen, an annulus where dust is distributed in clumps, and an outer envelope where Lyα photons only scatter. We show that the dust that affects the escape of Lyα is more restricted to the galaxies' central regions, while the larger Lyα halos are generated by scattering at large radii. We present an empirical modeling technique to quantify how much Lyα scatters in the halo, and find that this "characteristic" scattering distance correlates with the measured size of the Lyα halo. We note that there exists a slight anti-correlation between the scattering distance of Lyα and global dust properties. △ Less

Submitted 8 November, 2017; originally announced November 2017.

Comments: 32 pages, 51 figures, accepted to ApJ

MUSE Illuminates Channels for Lyman Continuum Escape in the Halo of SBS 0335-52E

Authors: E. C. Herenz, M. Hayes, P. Papaderos, J. M. Cannon, A. Bik, J. Melinder, G. Östlin

Abstract: We report on the discovery of ionised gas filaments in the circum-galactic halo of the extremely metal-poor compact starburst SBS 0335-052E in a 1.5h integration with the MUSE integral-field spectrograph. We detect these features in H$α$ and [OIII] emission down to surface-brightness levels of $5 \times 10^{-19}$erg s$^{-1}$cm$^{-2}$arcsec$^{-2}$. The filaments have projected diameters of 2.1 kpc… ▽ More We report on the discovery of ionised gas filaments in the circum-galactic halo of the extremely metal-poor compact starburst SBS 0335-052E in a 1.5h integration with the MUSE integral-field spectrograph. We detect these features in H$α$ and [OIII] emission down to surface-brightness levels of $5 \times 10^{-19}$erg s$^{-1}$cm$^{-2}$arcsec$^{-2}$. The filaments have projected diameters of 2.1 kpc and extend more than 9 kpc to the north and north-west from the main stellar body. We also detect extended nebular HeII $λ$4686 emission that brightens towards the north-west at the rim of a star-burst driven super-shell, suggestive of a locally enhanced UV radiation field due to shocks. We also present a velocity field of the ionised gas. The filaments appear to connect seamlessly in velocity space to the kinematical disturbances caused by the shell. Similar to high-$z$ star-forming galaxies, the ionised gas in this galaxy is dispersion dominated. We argue that the filaments were created via feedback from the starburst and that these ionised structures in the halo may act as escape channels for Lyman continuum radiation in this gas-rich system. △ Less

Submitted 5 October, 2017; v1 submitted 23 August, 2017; originally announced August 2017.

Comments: Revised version after peer review. Accepted for publication in A&A letters

The Lyman Continuum Escape and ISM properties in Tololo 1247-232 - New Insights from HST and VLA

Authors: Johannes Puschnig, Matthew Hayes, Göran Östlin, Thøger E. Rivera-Thorsen, Jens Melinder, John M. Cannon, Veronica Menacho, Erik Zackrisson, Nils Bergvall, Elisabet Leitet

Abstract: Low- and intermediate mass galaxies are widely discussed as cause of reionization at redshift $z\sim10-6$. However, observational proof of galaxies that are leaking ionizing radiation (Lyman continuum; LyC) is a currently ongoing challenge and the list of LyC emitting candidates is still short. Tololo 1247-232 is among those very few galaxies with recently reported leakage. We performed intermedia… ▽ More Low- and intermediate mass galaxies are widely discussed as cause of reionization at redshift $z\sim10-6$. However, observational proof of galaxies that are leaking ionizing radiation (Lyman continuum; LyC) is a currently ongoing challenge and the list of LyC emitting candidates is still short. Tololo 1247-232 is among those very few galaxies with recently reported leakage. We performed intermediate resolution ultraviolet (UV) spectroscopy with the Cosmic Origins Spectrograph (COS) onboard the Hubble Space Telescope and confirm ionizing radiation emerging from Tololo 1247-232. Adopting an improved data reduction procedure, we find that LyC escapes from the central stellar clusters, with an escape fraction of 1.5$\pm$0.5% only, i.e. the lowest value reported for the galaxy so far. We further make use of FUV absorption lines of Si II and Si IV as a probe of the neutral and ionized interstellar medium. We find that most of the ISM gas is ionized, likely facilitating LyC escape from density bounded regions. Neutral gas covering as a function of line-of-sight velocity is derived using the apparent optical depth method. The ISM is found to be sufficiently clumpy, supporting the direct escape of LyC photons. We further report on broadband UV and optical continuum imaging as well as narrowband imaging of Ly$α$, H$α$ and H$β$. Using stellar population synthesis, a Ly$α$ escape fraction of 8% was derived. We also performed VLA 21cm imaging. The hydrogen hyperfine transition was not detected, but a deep upper limit atomic gas mass of $\lesssim10^9 M_{\odot}$ could be derived. The upper limit gas fraction defined as $\frac{M_{H I}}{M_*}$ is only 20%. Evidence is found that the H I gas halo is relatively small compared to other Lyman Alpha emitters. △ Less

Submitted 19 April, 2017; originally announced April 2017.

A Multiply-Imaged z $\sim$ 6.3 Lyman Alpha Emitter candidate behind Abell 2261

Authors: Claes-Erik Rydberg, Adi Zitrin, Erik Zackrisson, Jens Melinder, Daniel J. Whalen, Ralf S. Klessen, Juan Gonzalez, Göran Östlin, Daniela Carollo

Abstract: While the Lyman-$α$ ($\mathrm{Ly}α$) emission line serves as an important tool in the study of galaxies at $z\lesssim 6$, finding Ly$α$ emitters (LAE) at significantly higher redshifts has been more challenging, probably because of the increasing neutrality of the intergalactic medium above $z\sim6$. Galaxies with extremely high rest-frame Ly$α$ equivalent widths, EW(Ly$α$) $\gtrsim 150$ Å, at… ▽ More While the Lyman-$α$ ($\mathrm{Ly}α$) emission line serves as an important tool in the study of galaxies at $z\lesssim 6$, finding Ly$α$ emitters (LAE) at significantly higher redshifts has been more challenging, probably because of the increasing neutrality of the intergalactic medium above $z\sim6$. Galaxies with extremely high rest-frame Ly$α$ equivalent widths, EW(Ly$α$) $\gtrsim 150$ Å, at $z>6$ are good candidates for Ly$α$ follow-up observations, and can stand out in multiband imaging surveys because of their unusual colors. We have conducted a photometric search for such objects in the Cluster Lensing And Supernova survey with Hubble (CLASH), and report here the identification of three likely gravitationally-lensed images of a single LAE candidate at $z\sim6.3$, behind the galaxy cluster Abell 2261($z = 0.225$). In the process, we also measured with Keck/MOSFIRE the first spectroscopic redshift of a multiply-imaged galaxy behind Abell 2261, at $z = 3.337$. This allows us to calibrate the lensing model, which in turn is used to study the properties of the candidate LAE. Population III galaxy spectral energy distribution (SED) model fits to the CLASH broadband photometry of the possible LAE provide a slightly better fit than Population I/II models. The best fitted model suggests intrinsic EW(Ly$α$) $\approx 160$ Å after absorption in the interstellar and intergalactic medium. Future spectroscopic observations will examine this prediction as well as shed more light on the morphology of this object, which indicates it may be a merger of two smaller galaxies. △ Less

Submitted 20 October, 2016; originally announced October 2016.

Comments: 13 pages, 6 figures

O VI Emission Imaging of a Galaxy with the Hubble Space Telescope: a Warm Gas Halo Surrounding the Intense Starburst SDSS J115630.63+500822.1

Authors: Matthew Hayes, Jens Melinder, Göran Östlin, Claudia Scarlata, Matthew D. Lehnert, Gustav Mannerström-Jansson

Abstract: We report results from a new HST study of the OVI 1032,1038Å doublet in emission around intensely star-forming galaxies. The programme aims to characterize the energy balance in starburst galaxies and gas cooling in the difficult-to-map coronal temperature regime of 2-5 x $10^5$K. We present the first resolved image of gas emission in the OVI line. Our target, SDSS J1156+5008, is very compact in t… ▽ More We report results from a new HST study of the OVI 1032,1038Å doublet in emission around intensely star-forming galaxies. The programme aims to characterize the energy balance in starburst galaxies and gas cooling in the difficult-to-map coronal temperature regime of 2-5 x $10^5$K. We present the first resolved image of gas emission in the OVI line. Our target, SDSS J1156+5008, is very compact in the continuum but displays OVI emission to radii of 23 kpc. The surface brightness profile is well fit by an exponential with a scale of 7.5kpc. This is ten times the size of the photoionized gas, and we estimate that 1/6 the total OVI luminosity comes from resonantly scattered continuum radiation. Spectroscopy - which closely resembles a stacked sample of archival spectra - confirms the OVI emission, and determines the column density and outflow velocity from blueshifted absorption. The combination of measurements enables several new calculations with few assumptions. The OVI regions fill only ~$10^{-3}$ of the volume. By comparing the cooling time with the cloud sound-crossing time, the cooling distance with the size, and the pressure in the OVI and nebular gas, we conclude that the OVI-bearing gas cannot have been lifted to the scale height at this temperature, and must be cooling in situ through this coronal temperature regime. The coronal phase contains ~1% of the ionized mass, and its kinetic energy is currently ~1% of the budget set by supernova feedback. However a much larger amount of the gas must have cooled through this phase during the star formation episode. The outflow exceeds the escape velocity and the gas may become unbound, but it will recombine before it escapes and become visible to Lyman (and OI) spectroscopy. The mapping of this gas represents a crucial step in further constraining galaxy formation scenarios and guiding the development of future satellites. △ Less

Submitted 16 June, 2016; v1 submitted 14 June, 2016; originally announced June 2016.

Comments: Accepted by the Astrophysical Journal. 25 pages, 11 figures. Section 7 presents calculated properties of warm halo gas. Version 2 fixes PDF compatibility issue for some PDF viewers

The Lyman alpha Reference Sample VI: Lyman alpha escape from the edge-on disk galaxy Mrk1486

Authors: Florent Duval, Göran Östlin, Matthew Hayes, Erik Zackrisson, Anne Verhamme, Ivana Orlitova, Angela Adamo, Lucia Guaita, Jens Melinder, John M. Cannon, Peter Laursen, Thoger Rivera-Thorsen, E. Christian Herenz, Pieter Gruyters, J. Miguel Mas-Hesse, Daniel Kunth, Andreas Sandberg, Daniel Schaerer, Tore Mansson

Abstract: While numerical simulations suggest that the strength of the Lyman alpha (Lya) line of star-forming disk galaxies strongly depends on the inclination at which they are observed (i.e. from edge-on to face-on, we expect to see a change from an attenuated Lya line to a strong Lya emission line), recent observations with the Hubble space telescope (HST) have highlighted few low-redshift highly incline… ▽ More While numerical simulations suggest that the strength of the Lyman alpha (Lya) line of star-forming disk galaxies strongly depends on the inclination at which they are observed (i.e. from edge-on to face-on, we expect to see a change from an attenuated Lya line to a strong Lya emission line), recent observations with the Hubble space telescope (HST) have highlighted few low-redshift highly inclined (edge-on) disk galaxies that breaks this trend. We aim to understand how a strong Lya emission line is able to escape from one of those inclined disk galaxies, named Mrk1486 (z=0.0338). For that purpose we used a large set of HST imaging and spectroscopic data to investigate both the ISM structure and the dominant source of Lya radiation inside Mrk1486. Moreover, we used a 3D Monte Carlo Lya radiation transfer code to study the radiative transfer of Lya and UV continuum photons inside a 3D geometry of neutral hydrogen (HI) and dust that models the ISM structure at the galaxy center. The analysis of IFU Halpha spectroscopic data of Mrk1486 indicates the presence of two bipolar galactic winds of HI gas above and bellow the disk plane of Mrk1486. Furthermore, comparing different diagnostic diagrams (such as [OIII]5007/Hbeta versus [OI]6300/Halpha) to photo- and shock-ionization models, we find that the Lya production of Mrk1486 is dominated by photoionization inside the galaxy disk. From this perspective, our numerical simulations succeed in reproducing the strength and spectral shape of the observed Lya line of Mrk1486 by assuming a scenario in which the Lya photons are produced inside the disk, travel along the galactic winds and scatter on cool HI materials toward the observer. As bipolar galactic winds are ubiquitous in star-forming disk galaxies, this mechanism may explain the origin of strong Lya emission lines recently observed from highly inclined galaxies at high-redshift. △ Less

Submitted 2 December, 2015; originally announced December 2015.

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

The Lyman alpha reference sample. VII. Spatially resolved H$α$ kinematics

Authors: Edmund Christian Herenz, Pieter Gruyters, Ivana Orlitova, Matthew Hayes, Göran Östlin, John M. Cannon, Martin M. Roth, Arjan Bik, Stephen Pardy, Héctor Otí-Floranes, J. Miguel Mas-Hesse, Angela Adamo, Hakim Atek, Florent Duval, Lucia Guaita, Daniel Kunth, Peter Laursen, Jens Melinder, Johannes Puschnig, Thøger E. Rivera-Thorsen, Daniel Schaerer, Anne Verhamme

Abstract: We present integral field spectroscopic observations with the Potsdam Multi Aperture Spectrophotometer of all 14 galaxies in the $z\sim 0.1$ Lyman Alpha Reference Sample (LARS). We produce 2D line of sight velocity maps and velocity dispersion maps from the Balmer $α$ (H$α$) emission in our data cubes. These maps trace the spectral and spatial properties of the LARS galaxies' intrinsic Ly$α$ radia… ▽ More We present integral field spectroscopic observations with the Potsdam Multi Aperture Spectrophotometer of all 14 galaxies in the $z\sim 0.1$ Lyman Alpha Reference Sample (LARS). We produce 2D line of sight velocity maps and velocity dispersion maps from the Balmer $α$ (H$α$) emission in our data cubes. These maps trace the spectral and spatial properties of the LARS galaxies' intrinsic Ly$α$ radiation field. We show our kinematic maps spatially registered onto the Hubble Space Telescope H$α$ and Lyman $α$ (Ly$α$) images. Only for individual galaxies a causal connection between spatially resolved H$α$ kinematics and Ly$α$ photometry can be conjectured. However, no general trend can be established for the whole sample. Furthermore, we compute non-parametric global kinematical statistics -- intrinsic velocity dispersion $σ_0$, shearing velocity $v_\mathrm{shear}$, and the $v_\mathrm{shear}/σ_0$ ratio -- from our kinematic maps. In general LARS galaxies are characterised by high intrinsic velocity dispersions (54\,km\,s$^{-1}$ median) and low shearing velocities (65\,km\,s$^{-1}$ median). $v_\mathrm{shear}/σ_0$ values range from 0.5 to 3.2 with an average of 1.5. Noteworthy, five galaxies of the sample are dispersion dominated systems with $v_\mathrm{shear}/σ_0 <1$ and are thus kinematically similar to turbulent star forming galaxies seen at high redshift. When linking our kinematical statistics to the global LARS Ly$α$ properties, we find that dispersion dominated systems show higher Ly$α$ equivalent widths and higher Ly$α$ escape fractions than systems with $v_\mathrm{shear}/σ_0 > 1$. Our result indicates that turbulence in actively star-forming systems is causally connected to interstellar medium conditions that favour an escape of Ly$α$ radiation. △ Less

Submitted 17 November, 2015; originally announced November 2015.

Comments: 26 pages, 15 figures, accepted for publication in A&A

Journal ref: A&A 587, A78 (2016)

Limits on Lyman Continuum escape from z=2.2 H-alpha emitting galaxies

Authors: A. Sandberg, G. Östlin, J. Melinder, A. Bik, L. Guaita

Abstract: The leakage of Lyman continuum photons from star forming galaxies is an elusive parameter. When observed, it provides a wealth of information on star formation in galaxies and the geometry of the interstellar medium, and puts constraints on the role of star forming galaxies in the reionization of the universe. H-alpha-selected galaxies at z~2 trace the highest star formation population at the peak… ▽ More The leakage of Lyman continuum photons from star forming galaxies is an elusive parameter. When observed, it provides a wealth of information on star formation in galaxies and the geometry of the interstellar medium, and puts constraints on the role of star forming galaxies in the reionization of the universe. H-alpha-selected galaxies at z~2 trace the highest star formation population at the peak of cosmic star formation history, providing a base for directly measuring Lyman continuum escape. Here we present this method, and highlight its benefits as well as caveats. We also use the method on 10 H-alpha emitters in the Chandra Deep Field South at z=2.2, also imaged with the Hubble Space Telescope in the ultraviolet. We find no individual Lyman continuum detections, and our stack puts a 5 sigma upper limit on the average absolute escape fraction of <24%, consistent with similar studies. With future planned observations, the sample sizes would rapidly increase and the method presented here would provide very robust constraints on the escape fraction. △ Less

Submitted 2 November, 2015; originally announced November 2015.

Comments: 6 pages, 2 figures. Accepted for publication in Astrophysical Journal Letters

The super star cluster driven feedback in ESO338-IG04 and Haro 11

Authors: Arjan Bik, Göran Östlin, Veronica Menacho, Angela Adamo, Matthew Hayes, Jens Melinder, Philippe Amram

Abstract: The stellar content of young massive star clusters emit large amounts of Lyman continuum photons and inject momentum into the inter stellar medium (ISM) by the strong stellar winds of the most massive stars in the cluster. When the most massive stars explode as supernovae, large amounts of mechanical energy are injected in the ISM. A detailed study of the ISM around these massive cluster provides… ▽ More The stellar content of young massive star clusters emit large amounts of Lyman continuum photons and inject momentum into the inter stellar medium (ISM) by the strong stellar winds of the most massive stars in the cluster. When the most massive stars explode as supernovae, large amounts of mechanical energy are injected in the ISM. A detailed study of the ISM around these massive cluster provides insights on the effect of cluster feedback. We present high quality integral field spectroscopy taken with VLT/MUSE of two starburst galaxies: ESO 338-IG04 and Haro 11. Both galaxies contain a significant number of super star clusters. The MUSE data provide us with an unprecedented view of the state and kinematics of the ionized gas in the galaxy allowing us to study the effect of stellar feedback on small and large spatial scales. We present our recent results on studying the ISM state of these two galaxies. The data of both galaxies show that the mechanical and ionization feedback of the super star clusters in the galaxy modify the state and kinematics of the ISM substancially by creating highly ionized bubbles around the cluster, making the central part of the galaxy highly ionized. This shows that the HII regions around the individual clusters are density bounded, allowing the ionizing photons to escape and ionize the ISM further out. △ Less

Submitted 7 October, 2015; originally announced October 2015.

Comments: 6 pages, 4 figures, to appear in proceedings, IAU Symposium 316, Formation, Evolution, and Survival of Massive Clusters, eds. C. Charbonnel & A. Nota

VLT/MUSE view of the highly ionized outflow cones in the nearby starburst ESO338-IG04

Authors: A. Bik, G. Östlin, M. Hayes, A. Adamo, J. Melinder, P. Amram

Abstract: The Ly$α$ line is an important diagnostic for star formation at high redshift, but interpreting its flux and line profile is difficult because of the resonance nature of Ly$α$. Trends between the escape of Ly$α$ photons and dust and properties of the interstellar medium (ISM) have been found, but detailed comparisons between Ly$α$ emission and the properties of the gas in local high-redshift analo… ▽ More The Ly$α$ line is an important diagnostic for star formation at high redshift, but interpreting its flux and line profile is difficult because of the resonance nature of Ly$α$. Trends between the escape of Ly$α$ photons and dust and properties of the interstellar medium (ISM) have been found, but detailed comparisons between Ly$α$ emission and the properties of the gas in local high-redshift analogs are vital for understanding the relation between Ly$α$ emission and galaxy properties. For the first time, we can directly infer the properties of the ionized gas at the same location and similar spatial scales of the extended Ly$α$ halo around ESO 338-IG04. We obtained VLT/MUSE integral field spectra. We used ionization parameter mapping of the [SII]/[OIII] line ratio and the kinematics of H$α$ to study the ionization state and kinematics of the ISM of ESO338-IG04. The velocity map reveals two outflows. The entire central area of the galaxy is highly ionized by photons leaking from the HII regions around the youngest star clusters. Three highly ionized cones have been identified, of which one is associated with an outflow. We propose a scenario where the outflows are created by mechanical feedback of the older clusters, while the highly ionized gas is caused by the hard ionizing photons emitted by the youngest clusters. A comparison with the Ly$α$ map shows that the (approximately bipolar) asymmetries observed in the Ly$α$ emission are consistent with the base of the outflows detected in H$α$. No clear correlation with the ionization cones is found. The mechanical and ionization feedback of star clusters significantly changes the state of the ISM by creating ionized cones and outflows. The comparison with Ly$α$ suggests that especially the outflows could facilitate the escape of Ly$α$ photons [Abridged]. △ Less

Submitted 31 March, 2015; v1 submitted 23 March, 2015; originally announced March 2015.

Comments: Accepted for publication in A&A Letters, 4 pages, 2 figures

Journal ref: A&A 576, L13 (2015)

The Lyman Alpha Reference Sample: V. The impact of neutral ISM kinematics and geometry on Lyman Alpha escape

Authors: Thøger E. Rivera-Thorsen, Matthew Hayes, Göran Östlin, Florent Duval, Ivana Orlitová, Anne Verhamme, J. Miguel Mas-Hesse, Daniel Schaerer, John M. Cannon, Héctor Otí-Floranes, Andreas Sandberg, Lucia Guaita, Angela Adamo, Hakim Atek, E. Christian Herenz, Daniel Kunth, Peter Laursen, Jens Melinder

Abstract: We present high-resolution far-UV spectroscopy of the 14 galaxies of the Lyman Alpha Reference Sample; a sample of strongly star-forming galaxies at low redshifts ($0.028 < z < 0.18$). We compare the derived properties to global properties derived from multi band imaging and 21 cm HI interferometry and single dish observations, as well as archival optical SDSS spectra. Besides the Lyman $α$ line… ▽ More We present high-resolution far-UV spectroscopy of the 14 galaxies of the Lyman Alpha Reference Sample; a sample of strongly star-forming galaxies at low redshifts ($0.028 < z < 0.18$). We compare the derived properties to global properties derived from multi band imaging and 21 cm HI interferometry and single dish observations, as well as archival optical SDSS spectra. Besides the Lyman $α$ line, the spectra contain a number of metal absorption features allowing us to probe the kinematics of the neutral ISM and evaluate the optical depth and and covering fraction of the neutral medium as a function of line-of-sight velocity. Furthermore, we show how this, in combination with precise determination of systemic velocity and good Ly$α$ spectra, can be used to distinguish a model in which separate clumps together fully cover the background source, from the "picket fence" model named by Heckman et al. (2011). We find that no one single effect dominates in governing Ly$α$ radiative transfer and escape. Ly$α$ escape in our sample coincides with a maximum velocity-binned covering fraction of $\lesssim 0.9$ and bulk outflow velocities of $\gtrsim 50$ km s$^{-1}$, although a number of galaxies show these characteristics and yet little or no Ly$α$ escape. We find that Ly$α$ peak velocities, where available, are not consistent with a strong backscattered component, but rather with a simpler model of an intrinsic emission line overlaid by a blueshifted absorption profile from the outflowing wind. Finally, we find a strong anticorrelation between H$α$ equivalent width and maximum velocity-binned covering factor, and propose a heuristic explanatory model. △ Less

Submitted 5 March, 2015; v1 submitted 3 March, 2015; originally announced March 2015.

Comments: 28 pages, 19 figures, 5 tables

Journal ref: 2015 ApJ 805 14

The Lyman alpha reference sample: IV. Morphology at low and high redshift

Authors: L. Guaita, J. Melinder, M. Hayes, G. Ostlin, J. Gonzalez, LARS team

Abstract: We measured the sizes and morphological parameters of LARS galaxies in the continuum, Lya, and Ha images. We studied morphology by using the Gini coefficient vs M20 and asymmetry vs concentration diagrams. We then simulated LARS galaxies at z~2 and 5.7, performing the same morphological measurements. We also investigated the detectability of LARS galaxies in current deep field observations. The su… ▽ More We measured the sizes and morphological parameters of LARS galaxies in the continuum, Lya, and Ha images. We studied morphology by using the Gini coefficient vs M20 and asymmetry vs concentration diagrams. We then simulated LARS galaxies at z~2 and 5.7, performing the same morphological measurements. We also investigated the detectability of LARS galaxies in current deep field observations. The subsample of LAEs within LARS (LARS-LAEs) was stacked to provide a comparison to stacking studies performed at high redshift. LARS galaxies have continuum size, stellar mass, and rest-frame absolute magnitude typical of Lyman break analogues in the local Universe and also similar to 2<z<3 star-forming galaxies and massive LAEs. LARS optical morphology is consistent with the one of merging systems, and irregular or starburst galaxies. For the first time we quantify the morphology in Lya images: even if a variety of intrinsic conditions of the interstellar medium can favour the escape of Lya photons, LARS-LAEs appear small in the continuum, and their Lya is compact. LARS galaxies tend to be more extended in Lya than in the rest-frame UV. It means that Lya photons escape by forming haloes around HII regions of LARS galaxies. The stack of LARS-LAE Lya images is peaked in the centre, indicating that the conditions, which make a galaxy an LAE, tend to produce a concentrated surface brightness profile. On the other hand, the stack of all LARS galaxies is shallower and more extended. This can be caused by the variety of dust and HI amount and distribution, which produces a more complex, patchy, and extended profile, like the one observed for Lyman break galaxies that can contribute to the stack. We cannot identify a single morphological property that controls whether a galaxy emits a net positive Lya flux. However, the LARS-LAEs have continuum properties consistent with merging systems. △ Less

Submitted 10 January, 2015; originally announced January 2015.

Comments: Accepted date: Dec 19th, 2014. arXiv admin note: text overlap with arXiv:astro-ph/0311352

A search for Population III galaxies in CLASH. I. Singly-imaged candidates at high redshift

Authors: Claes-Erik Rydberg, Erik Zackrisson, Adi Zitrin, Lucia Guaita, Jens Melinder, Saghar Asadi, Juan Gonzalez, Göran Östlin, Tina Ström

Abstract: Population III galaxies are predicted to exist at high redshifts and may be rendered sufficiently bright for detection with current telescopes when gravitationally lensed by a foreground galaxy cluster. Population III galaxies that exhibit strong Lya emission should furthermore be identifiable from broadband photometry because of their unusual colors. Here, we report on a search for such objects a… ▽ More Population III galaxies are predicted to exist at high redshifts and may be rendered sufficiently bright for detection with current telescopes when gravitationally lensed by a foreground galaxy cluster. Population III galaxies that exhibit strong Lya emission should furthermore be identifiable from broadband photometry because of their unusual colors. Here, we report on a search for such objects at z > 6 in the imaging data from the Cluster Lensing And Supernova survey with Hubble (CLASH), covering 25 galaxy clusters in 16 filters. Our selection algorithm returns five singly-imaged candidates with Lya-like color signatures, for which ground-based spectroscopy with current 8-10 m class telescopes should be able to test the predicted strength of the Lya line. None of these five objects have been included in previous CLASH compilations of high-redshift galaxy candidates. However, when large grids of spectral synthesis models are applied to the study of these objects, we find that only two of these candidates are significantly better fitted by Population III models than by more mundane, low-metallicity stellar populations. △ Less

Submitted 19 February, 2015; v1 submitted 20 November, 2014; originally announced November 2014.

Comments: 18 pages, 10 figures, accepted by ApJ