Search

The star formation and metal enrichment histories of galaxies - at any epoch - constitute one of the key properties of galaxies, and their measurement is a core aim of observational extragalactic astronomy. The lack of deep rest-frame optical coverage at high-redshift has made robust constraints elusive, but this is now changing thanks to the \emph{James Webb Space Telescope (JWST)}. In preparation for the constraints provided by \emph{JWST} we explore the star formation and metal enrichment histories of galaxies at $z=5-13$ using the First Light And Reionisation Epoch Simulations (FLARES) suite. Built on the EAGLE model, the unique strategy of FLARES allows us to simulate a wide range of stellar masses (and luminosities) and environments. While we predict significant redshift evolution of average ages and specific star formation rates our core result is a mostly flat relationship of age and specific star formation rate with stellar mass. We also find that galaxies in this epoch predominantly have strongly rising star formation histories, albeit with the magnitude dropping with redshift and stellar mass. In terms of chemical enrichment we predict a strong stellar mass - metallicity relation present at $z=10$ and beyond alongside significant $\alpha$-enhancement. Finally, we find no environmental dependence of the relationship between age, specific star formation rate, or metallicity with stellar mass., Comment: To be submitted to MNRAS. Comments welcome

Full list of authors: Welch, Brian; Coe, Dan; Zackrisson, Erik; de Mink, S. E.; Ravindranath, Swara; Anderson, Jay; Brammer, Gabriel; Bradley, Larry; Yoon, Jinmi; Kelly, Patrick; Diego, Jose M.; Windhorst, Rogier; Zitrin, Adi; Dimauro, Paola; Jimenez-Teja, Yolanda; Abdurro'uf; Nonino, Mario; Acebron, Ana; Andrade-Santos, Felipe; Avila, Roberto J.; Bayliss, Matthew B.; Benitez, Alex; Broadhurst, Tom; Bhatawdekar, Rachana; Bradac, Marusa; Caminha, Gabriel B.; Chen, Wenlei; Eldridge, Jan; Farag, Ebraheem; Florian, Michael; Frye, Brenda; Fujimoto, Seiji; Gomez, Sebastian; Henry, Alaina; Hsiao, Tiger Y-Y; Hutchison, Taylor A.; James, Bethan L.; Joyce, Meridith; Jung, Intae; Khullar, Gourav; Larson, Rebecca L.; Mahler, Guillaume; Mandelker, Nir; McCandliss, Stephan; Morishita, Takahiro; Newshore, Rosa; Norman, Colin; O'Connor, Kyle; Oesch, Pascal A.; Oguri, Masamune; Ouchi, Masami; Postman, Marc; Rigby, Jane R.; Ryan, Russell E., Jr.; Sharma, Soniya; Sharon, Keren; Strait, Victoria; Strolger, Louis-Gregory; Timmes, F. X.; Toft, Sune; Trenti, Michele; Vanzella, Eros; Vikaeus, Anton.--This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited., The gravitationally lensed star WHL 0137–LS, nicknamed Earendel, was identified with a photometric redshift zphot = 6.2 ± 0.1 based on images taken with the Hubble Space Telescope. Here we present James Webb Space Telescope (JWST) Near Infrared Camera images of Earendel in eight filters spanning 0.8–5.0 μm. In these higher-resolution images, Earendel remains a single unresolved point source on the lensing critical curve, increasing the lower limit on the lensing magnification to μ > 4000 and restricting the source plane radius further to r < 0.02 pc, or ∼4000 au. These new observations strengthen the conclusion that Earendel is best explained by an individual star or multiple star system and support the previous photometric redshift estimate. Fitting grids of stellar spectra to our photometry yields a stellar temperature of Teff ≃ 13,000–16,000 K, assuming the light is dominated by a single star. The delensed bolometric luminosity in this case ranges from $\mathrm{log}(L)=5.8$ to 6.6 L⊙, which is in the range where one expects luminous blue variable stars. Follow-up observations, including JWST NIRSpec scheduled for late 2022, are needed to further unravel the nature of this object, which presents a unique opportunity to study massive stars in the first billion years of the universe. © 2022. The Author(s). Published by the American Astronomical Society., E.Z. acknowledges support from the Swedish National Space Board. M.B. acknowledges support from the Slovenian national research agency ARRS through grant N1-0238. A.Z. acknowledges support by grant No. 2020750 from the United States–Israel Binational Science Foundation (BSF) and grant No. 2109066 from the United States National Science Foundation (NSF) and by the Ministry of Science & Technology, Israel. M.T. acknowledges support by the Australian Research Council Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D) through project No. CE170100013. The Cosmic Dawn Center (DAWN) is funded by the Danish National Research Foundation under grant No. 140. J.M.D. acknowledges the support of projects PGC2018-101814-B-100 and MDM-2017-0765. Y.J.-T. acknowledges financial support from the European Union's Horizon 2020 research and innovation program under Marie Skłodowska-Curie grant agreement No. 898633, the MSCA IF Extensions Program of the Spanish National Research Council (CSIC), and the State Agency for Research of the Spanish MCIU through the Center of Excellence Severo Ochoa award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709). P.A.O. acknowledges support by the Swiss National Science Foundation through project grant 200020_207349.

Overdense, metal-rich regions, shielded from stellar radiation might remain neutral throughout reionization and produce metal as well as 21 cm absorption lines. Simultaneous absorption from metals and 21 cm can complement each other as probes of underlying gas properties. We use Aurora, a suite of high resolution radiation-hydrodynamical simulations of galaxy formation, to investigate the occurrence of such "aligned" absorbers. We calculate absorption spectra for 21 cm, OI, CII, SiII and FeII. We find velocity windows with absorption from at least one metal and 21 cm, and classify the aligned absorbers into two categories: 'aligned and cospatial absorbers' and 'aligned but not cospatial absorbers'. While 'aligned and cospatial absorbers' originate from overdense structures and can be used to trace gas properties, 'aligned but not cospatial absorbers' are due to peculiar velocity effects. The incidence of absorbers is redshift dependent, as it is dictated by the interplay between reionization and metal enrichment, and shows a peak at $z \approx 8$ for the aligned and cospatial absorbers. While aligned but not cospatial absorbers disappear towards the end of reionization because of the lack of an ambient 21 cm forest, aligned and cospatial absorbers are associated with overdense pockets of neutral gas which can be found at lower redshift. We produce mock observations for a set of sightlines for the next generation of telescopes like the ELT and SKA1-LOW, finding that given a sufficiently bright background quasar, these telescopes will be able to detect both types of absorbers throughout reionization., Accepted in MNRAS. Comments welcome

We present integral field, far-infrared (FIR) spectroscopy of Mrk 54, a local Lyman Continuum Emitter (LCE), obtained with FIFI-LS on the Stratospheric Observatory for Infrared Astronomy. This is only the second time, after Haro 11, that [C II] 158 $\mu$m and [O III] 88 $\mu$m spectroscopy of the known LCEs have been obtained. We find that Mrk 54 has a strong [C II] emission that accounts for $\sim1$% of the total FIR luminosity, whereas it has only moderate [O III] emission, resulting in the low [O III]/[C II] luminosity ratio of $0.22\pm0.06$. In order to investigate whether [O III]/[C II] is a useful tracer of $f_{\rm esc}$ (LyC escape fraction), we examine the correlations of [O III]/[C II] and (i) the optical line ratio of $\rm O_{32} \equiv$ [O III] 5007 \AA/[O II] 3727 \AA, (ii) specific star formation rate, (iii) [O III] 88 $\mu$m/[O I] 63 $\mu$m ratio, (iv) gas phase metallicity, and (v) dust temperature based on a combined sample of Mrk 54 and the literature data from the Herschel Dwarf Galaxy Survey and the LITTLE THINGS Survey. We find that galaxies with high [O III]/[C II] luminosity ratios could be the result of high ionization (traced by $\rm O_{32}$), bursty star formation, high ionized-to-neutral gas volume filling factors (traced by [O III] 88 $\mu$m/[O I] 63 $\mu$m), and low gas-phase metallicities, which is in agreement with theoretical predictions. We present an empirical relation between the [O III]/[C II] ratio and $f_{\rm esc}$ based on the combination of the [O III]/[C II] and $\rm O_{32}$ correlation, and the known relation between $\rm O_{32}$ and $f_{\rm esc}$. The relation implies that high-redshift galaxies with high [O III]/[C II] ratios revealed by ALMA may have $f_{\rm esc}\gtrsim0.1$, significantly contributing to the cosmic reionization., Comment: 14 pages, 5 figures, Accepted for publication in ApJ

We present new observations with the Atacama Large Millimeter/submillimeter Array for a gravitationally-lensed galaxy at $z=9.1$, MACS1149-JD1. [O III] 88-$\mu$m emission is detected at 10$\sigma$ with a spatial resolution of $\sim0.3$ kpc in the source plane, enabling the most distant morpho-kinematic study of a galaxy. The [O III] emission is distributed smoothly without any resolved clumps and shows a clear velocity gradient with $\Delta V_{\rm obs}/2\sigma_{\rm tot}=0.84\pm0.23$, where $\Delta V_{\rm obs}$ is the observed maximum velocity difference and $\sigma_{\rm tot}$ is the velocity dispersion measured in the spatially-integrated line profile, suggesting a rotating system. Assuming a geometrically thin self-gravitating rotation disk model, we obtain $V_{\rm rot}/\sigma_V=0.67_{-0.26}^{+0.73}$, where $V_{\rm rot}$ and $\sigma_V$ are the rotation velocity and velocity dispersion, respectively, still consistent with rotation. The resulting disk mass of $0.65_{-0.40}^{+1.37}\times10^{9}$ M$_\odot$ is consistent with being associated with the stellar mass identified with a 300 Myr-old stellar population independently indicated by a Balmer break in the spectral energy distribution. We conclude that the most of the dynamical mass is associated with the previously-identified mature stellar population that formed at $z\sim15$., Comment: ApJL accepted

We report on the detection of the [CII] 157.7 $\mu$m emission from the Lyman break galaxy (LBG) MACS0416_Y1 at z = 8.3113, by using the Atacama Large Millimeter/submillimeter Array (ALMA). The luminosity ratio of [OIII] 88 $\mu$m (from previous campaigns) to [CII] is 9.31 $\pm$ 2.6, indicative of hard interstellar radiation fields and/or a low covering fraction of photo-dissociation regions. The emission of [CII] is cospatial to the 850 $\mu$m dust emission (90 $\mu$m rest-frame, from previous campaigns), however the peak [CII] emission does not agree with the peak [OIII] emission, suggesting that the lines originate from different conditions in the interstellar medium. We fail to detect continuum emission at 1.5 mm (160 $\mu$m rest-frame) down to 18 $\mu$Jy (3$\sigma$). This nondetection places a strong limit on the dust spectrum, considering the 137 $\pm$ 26 $\mu$Jy continuum emission at 850 $\mu$m. This suggests an unusually warm dust component (T $>$ 80 K, 90% confidence limit), and/or a steep dust-emissivity index ($\beta_{\rm dust}$ $>$ 2), compared to galaxy-wide dust emission found at lower redshifts (typically T $\sim$ 30 - 50 K, $\beta_{\rm dust}$ $\sim$ 1 - 2). If such temperatures are common, this would reduce the required dust mass and relax the dust production problem at the highest redshifts. We therefore warn against the use of only single-wavelength information to derive physical properties, recommend a more thorough examination of dust temperatures in the early Universe, and stress the need for instrumentation that probes the peak of warm dust in the Epoch of Reionization., Comment: 15 pages, 8 figures, Accepted by MNRAS

We give an overview and describe the rationale, methods, and first results from NIRCam images of the JWST "Prime Extragalactic Areas for Reionization and Lensing Science" ("PEARLS") project. PEARLS uses up to eight NIRCam filters to survey several prime extragalactic survey areas: two fields at the North Ecliptic Pole (NEP); seven gravitationally lensing clusters; two high redshift proto-clusters; and the iconic backlit VV 191 galaxy system to map its dust attenuation. PEARLS also includes NIRISS spectra for one of the NEP fields and NIRSpec spectra of two high-redshift quasars. The main goal of PEARLS is to study the epoch of galaxy assembly, AGN growth, and First Light. Five fields, the JWST NEP Time-Domain Field (TDF), IRAC Dark Field (IDF), and three lensing clusters, will be observed in up to four epochs over a year. The cadence and sensitivity of the imaging data are ideally suited to find faint variable objects such as weak AGN, high-redshift supernovae, and cluster caustic transits. Both NEP fields have sightlines through our Galaxy, providing significant numbers of very faint brown dwarfs whose proper motions can be studied. Observations from the first spoke in the NEP TDF are public. This paper presents our first PEARLS observations, their NIRCam data reduction and analysis, our first object catalogs, the 0.9-4.5 $\mu$m galaxy counts and Integrated Galaxy Light. We assess the JWST sky brightness in 13 NIRCam filters, yielding our first constraints to diffuse light at 0.9-4.5 {\mu}m. PEARLS is designed to be of lasting benefit to the community., Comment: Accepted to AJ, comments welcome. We ask anyone who uses our public PEARLS (NEP TDF) data to refer to this overview paper

Dyson spheres are hypothetical megastructures built by advanced extraterrestrial civilizations to harvest radiation energy from stars. Here, we combine optical data from Gaia DR2 with mid-infrared data from AllWISE to set the strongest upper limits to date on the prevalence of partial Dyson spheres within the Milky Way, based on their expected waste-heat signatures. Conservative upper limits are presented on the fraction of stars at G $\leq$ 21 that may potentially host non-reflective Dyson spheres that absorb 1 - 90$\%$ of the bolometric luminosity of their host stars and emit thermal waste-heat in the 100 - 1000 K range. Based on a sample of $\approx$ $2.7\mathrm{e}\,5$ stars within 100 pc, we find that a fraction less than $\approx$ $2\mathrm{e}\,-5$ could potentially host $\sim$300 K Dyson spheres at 90$\%$ completion. These limits become progressively weaker for less complete Dyson spheres due to increased confusion with naturally occurring sources of strong mid-infrared radiation, and also at larger distances, due to the detection limits of WISE. For the $\sim2.9\mathrm{e}\,8$ stars within 5 kpc in our Milky Way sample, the corresponding upper limit on the fraction of stars that could potentially be $\sim$300 K Dyson spheres at 90$\%$ completion is $\leq$ $8\mathrm{e}\,-4$., Accepted for publication in MNRAS

Detailed observations of star forming galaxies at high redshift are critical to understand the formation and evolution of the earliest galaxies. Gravitational lensing provides an important boost, allowing observations at physical scales unreachable in unlensed galaxies. We present three lensed galaxies from the RELICS survey at $z_{phot} = 6 - 10$, including the most highly magnified galaxy at $z_{phot} \sim 6$ (WHL0137-zD1, dubbed the Sunrise Arc), the brightest known lensed galaxy at $z_{phot} \sim 6$ (MACS0308-zD1), and the only spatially resolved galaxy currently known at $z_{phot} \sim 10$ (SPT0615-JD). The Sunrise Arc contains seven star-forming clumps with delensed radii as small as 3 pc, the smallest spatial scales yet observed in a $z>6$ galaxy, while SPT0615-JD contains features measuring a few tens of parsecs. MACS0308-zD1 contains a $r\sim 30$ pc clump with a star formation rate (SFR) of $\sim 3 M_{\odot} \textrm{ yr}^{-1}$, giving it a SFR surface density of $\Sigma_{SFR} \sim 10^3 M_{\odot}\textrm{ yr}^{-1}\textrm{ kpc}^{-2}$. These galaxies provide a unique window into small scale star formation during the Epoch of Reionization. They will be excellent targets for future observations with JWST, including one approved program targeting the Sunrise Arc., Comment: 15 pages, 7 figures, submitted to ApJ

Dark matter halos that reach the HI-cooling mass without prior star formation or external metal pollution represent potential sites for the formation of small - extremely faint - Population III galaxies at high redshifts. Gravitational lensing may in rare cases boost their fluxes to detectable levels, but to find even a small number of such objects in randomly selected regions of the sky requires very large areas to be surveyed. Because of this, a small, wide-field telescope can in principle offer better detection prospects than a large telescope with a smaller field of view. Here, we derive the minimum comoving number density required to allow gravitational lensing to lift such objects at redshift $z=5-16$ above the detection thresholds of blind surveys carried out with the James Webb space telescope (JWST), the Roman space telescope (RST) and Euclid. We find that the prospects for photometric detections of Pop III galaxies is promising, and that they are better for RST than for JWST and Euclid. However, the Pop III galaxies favoured by current simulations have number densities too low to allow spectroscopic detections based on the strength of the HeII1640 emission line in any of the considered surveys unless very high star formation efficiencies ($\epsilon > 0.1$) are envoked. We argue that targeting individual cluster lenses instead of the wide field surveys considered in this paper results in better spectroscopic detection prospects, while for photometric detection, the wide field surveys perform considerably better., Comment: 15 pages, 4 figures

The upcoming James Webb Space Telescope (JWST) will allow observations of high-redshift galaxies at fainter detection levels than ever before, and JWST surveys targeting gravitationally lensed fields are expected to bring z ≳ 6 objects with very low star formation rate (SFR) within reach of spectroscopic studies. As galaxies at lower and lower star formation activity are brought into view, many of the standard methods used in the analysis of integrated galaxy spectra are at some point bound to break down, due to violation of the assumptions of a well-sampled stellar initial mass function (IMF) and a slowly varying SFR. We argue that galaxies with SFR ∼ 0.1 M⊙ yr−1 are likely to turn up at the spectroscopic detection limit of JWST in lensed fields, and investigate to what extent star formation sampling may affect the spectral analysis of such objects. We use the slug spectral synthesis code to demonstrate that such effects are likely to have significant impacts on spectral diagnostics of, for example, the Balmer emission lines. These effects are found to stem primarily from SFRs varying rapidly on short (∼Myr) time-scales due to star formation in finite units (star clusters), whereas the effects of an undersampled IMF is deemed insignificant in comparison. In contrast, the ratio between the He ii- and H i-ionizing flux is found to be sensitive to IMF-sampling as well as ICMF-sampling (sampling of the initial cluster mass function), which may affect interpretations of galaxies containing Population III stars or other sources of hard ionizing radiation.

Recently, spectroscopic detections of O[III] 88 μm and Ly-α emission lines from the z ≍ 9.1 galaxy MACS1149-JD1 have been presented, and with these, some interesting properties of this galaxy were uncovered. One such property is that MACS1149-JD1 exhibits a significant Balmer break at around rest-frame 4000 Å, which may indicate that the galaxy has experienced large variations in star formation rate prior to z ∼ 9, with a rather long period of low star formation activity. While some simulations predict large variations in star formation activity in high-redshift galaxies, it is unclear whether the simulations can reproduce the kind of variations seen in MACS1149-JD1. Here, we utilize synthetic spectra of simulated galaxies from two simulation suites in order to study to what extent these can accurately reproduce the spectral features (specifically the Balmer break) observed in MACS1149-JD1. We show that while the simulations used in this study produce galaxies with varying star formation histories, galaxies such as MACS1149-JD1 would be very rare in the simulations. In principle, future observations with the James Webb Space Telescope may tell us if MACS1149-JD1 represents something rare, or if such galaxies are more common than predicted by current simulations.

A fundamental quest of modern astronomy is to locate the earliest galaxies and study how they influenced the intergalactic medium a few hundred million years after the Big Bang. The abundance of star-forming galaxies is known to decline from redshifts of about 6 to 10, but a key question is the extent of star formation at even earlier times, corresponding to the period when the first galaxies might have emerged. Here we present spectroscopic observations of MACS1149-JD1, a gravitationally lensed galaxy observed when the Universe was less than four per cent of its present age. We detect an emission line of doubly ionized oxygen at a redshift of $9.1096\pm0.0006$, with an uncertainty of one standard deviation. This precisely determined redshift indicates that the red rest-frame optical colour arises from a dominant stellar component that formed about 250 million years after the Big Bang, corresponding to a redshift of about 15. Our results indicate the it may be possible to detect such early episodes of star formation in similar galaxies with future telescopes., Comment: To appear in the 17 May issue of Nature

We present new Atacama Large Millimeter/submillimeter Array Band 7 observational results of a Lyman break galaxy at $ z=7.15 $, B14-65666 ("Big Three Dragons"), which is an object detected in [OIII] 88 $\rm{\mu m}$, [CII] 158 $\rm{\mu m}$, and dust-continuum emission during the epoch of reionization. Our targets are the [NII] 122 $\rm{\mu m}$ fine-structure emission line and underlying 120 $\rm{\mu m}$ dust continuum. The dust continuum is detected with a $ \sim $19$ \sigma $ significance. From far-infrared spectral energy distribution sampled at 90, 120, and 160 $\rm{\mu m}$, we obtaine a best-fit dust temperature of $ 40 $ K ($ 79 $ K) and an infrared luminosity of $ \log_{10}(L_{\rm IR}/{\rm L}_\odot)=11.6$ ($12.1$) at the emissivity index $ \beta = 2.0 $ (1.0). The [NII] 122 $\rm{\mu m}$ line is not detected. The 3$ \sigma $ upper limit of the [NII] luminosity is $ 8.1 \times 10^7\ {\rm L}_\odot$. From the [NII], [OIII], and [CII] line luminosities, we use the Cloudy photoionization code to estimate nebular parameters as functions of metallicity. If the metallicity of the galaxy is high ($ Z > 0.4\ {\rm Z}_\odot$), the ionization parameter and hydrogen density are $ \log_{10} U \simeq -2.7\pm0.1$ and $ n_\text{H} \simeq 50$-$250\ {\rm cm}^{-3}$, respectively, which are comparable to those measured in low-redshift galaxies. The nitrogen-to-oxygen abundance ratio, $\rm{N/O}$, is constrained to be sub-solar. At $ Z < 0.4\ {\rm Z}_\odot$, the allowed $ U $ drastically increases as the assumed metallicity decreases. For high ionization parameters, the $\rm{N/O}$ constraint becomes weak. Finally, our Cloudy models predict the location of B14-65666 on the BPT diagram, thereby allowing a comparison with low-redshift galaxies., Comment: 17 pages, 13 figures, 2 tables. accepted for publication in ApJ

The James Webb Space Telescope (JWST) will observe the rest-frame ultraviolet/optical spectra of galaxies from the epoch of reionization (EoR) in unprecedented detail. While escaping into the intergalactic medium, hydrogen-ionizing (Lyman continuum; LyC) photons from the galaxies will contribute to the bluer end of the UV slope and make nebular emission lines less prominent. We present a method to constrain leakage of the LyC photons using the spectra of high redshift (z ≳ 6) galaxies. We simulate JWST/NIRSpec observations of galaxies at z =6–9 by matching the fluxes of galaxies observed in the Frontier Fields observations of galaxy cluster MACS-J0416. Our method predicts the escape fraction fesc with a mean absolute error Δfesc ≈ 0.14. The method also predicts the redshifts of the galaxies with an error $\left\langle \frac{\Delta z}{(1+z)}\right\rangle \approx 0.0003$.

If the dark matter halos of galaxies contain large numbers of subhalos as predicted by the $\Lambda$CDM model, these subhalos are expected to appear in strong galaxy--galaxy lens systems as small--scale perturbations in individual images. We simulate observations of multiply--lensed sub--mm galaxies at $z\sim2$ as a probe of the dark matter halo of a lens galaxy at $z\sim0.5$. We present detection limits for dark substructures based on a visibility plane analysis of simulated Atacama Large Millimeter/submillimeter Array (ALMA) data in bands 7, 8 and 9. We explore two effects: local surface brightness anomalies on angular scales similar to the Einstein radius and the astrometric shift of macroimages. This improves the sensitivity of our lens modeling to the mass of the lens perturber. We investigate the sensitivity of the detection of low--mass subhalos to the projected position of the subhalo on the image plane as well as the source structure and inner density profile of the lens. We demonstrate that, using the most extended ALMA configuration, pseudo-Jaffe subhalos can be detected with 99% confidence down to $M = 10^7 M_\odot$. We show how the detection threshold for the three ALMA bands depends on the projected position of the subhalo with respect to the lensed images and conclude that, despite the highest nominal angular resolution, band 9 provides the poorest sensitivity due to observational noise. All simulations use the \emph{ALMA Full ops most extended} ALMA configuration setup in CASA., Comment: 12 pages, 12 figures - Accepted for publication in MNRAS

Photometric observations of the spectroscopically confirmed $z\approx 9.1$ galaxy MACS1149-JD1 have indicated the presence of a prominent Balmer break in its spectral energy distribution, which may be interpreted as due to very large fluctuations in its past star formation activity. In this paper, we investigate to what extent contemporary simulations of high-redshift galaxies produce star formation rate variations sufficiently large to reproduce the observed Balmer break of MACS1149-JD1. We find that several independent galaxy simulations are unable to account for Balmer breaks of the inferred size, suggesting that MACS1149-JD1 either must be a very rare type of object or that our simulations are missing some key ingredient. We present predictions of spectroscopic Balmer break strength distributions for $z\approx 7-9$ galaxies that may be tested through observations with the upcoming James Webb Space Telescope and also discuss the impact that various assumptions on dust reddening, Lyman continuum leakage and deviations from a standard stellar initial mass function would have on the results., Comment: 11 pages, 8 figures. Accepted for publication in MNRAS

We present new ALMA observations and physical properties of a Lyman Break Galaxy at z=7.15. Our target, B14-65666, has a bright ultra-violet (UV) absolute magnitude, $M_{\rm UV}\approx-22.4$, and has been spectroscopically identified in Ly$\alpha$ with a small rest-frame equivalent width of $\approx4$ \AA. Previous HST image has shown that the target is comprised of two spatially separated clumps in the rest-frame UV. With ALMA, we have newly detected spatially resolved [OIII] 88 $\mu$m, [CII] 158 $\mu$m, and their underlying dust continuum emission. In the whole system of B14-65666, the [OIII] and [CII] lines have consistent redshifts of $7.1520\pm0.0003$, and the [OIII] luminosity, $(34.4\pm4.1)\times10^{8}L_{\rm \odot}$, is about three times higher than the [CII] luminosity, $(11.0\pm1.4)\times10^{8}L_{\rm \odot}$. With our two continuum flux densities, the dust temperature is constrained to be $T_{\rm d}\approx50-60$ K under the assumption of the dust emissivity index of $\beta_{\rm d}=2.0-1.5$, leading to a large total infrared luminosity of $L_{\rm TIR}\approx1\times10^{12}L_{\rm \odot}$. Owing to our high spatial resolution data, we show that the [OIII] and [CII] emission can be spatially decomposed into two clumps associated with the two rest-frame UV clumps whose spectra are kinematically separated by $\approx200$ km s$^{-1}$. We also find these two clumps have comparable UV, infrared, [OIII], and [CII] luminosities. Based on these results, we argue that B14-65666 is a starburst galaxy induced by a major-merger. The merger interpretation is also supported by the large specific star-formation rate (defined as the star-formation rate per unit stellar mass), sSFR$=260^{+119}_{-57}$ Gyr$^{-1}$, inferred from our SED fitting. Probably, a strong UV radiation field caused by intense star formation contributes to its high dust temperature and the [OIII]-to-[CII] luminosity ratio., Comment: 24 pages, 13 figures, 9 tables. Accepted for publication in PASJ. In addition to the dust continuum detection in ALMA Band 6 reported in the original manuscript, we newly report the dust continuum detection in ALMA Band 8 based on our new ALMA Cycle 5 observations

The James Webb Space Telescope (JWST) NIRSpec instrument will allow rest-frame ultraviolet/optical spectroscopy of galaxies in the epoch of reionization (EoR). Some galaxies may exhibit significant leakage of hydrogen-ionizing photons into the intergalactic medium, resulting in faint nebular emission lines. We present a machine learning framework for identifying cases of very high hydrogen-ionizing photon escape from galaxies based on the data quality expected from potential NIRSpec observations of EoR galaxies in lensed fields. We train our algorithm on mock samples of JWST/NIRSpec data for galaxies at redshifts $z=6$--10. To make the samples more realistic, we combine synthetic galaxy spectra based on cosmological galaxy simulations with observational noise relevant for $z\gtrsim 6$ objects of a brightness similar to EoR galaxy candidates uncovered in Frontier Fields observations of galaxy cluster Abell-2744 and MACS-J0416. We find that ionizing escape fractions ($f_\mathrm{esc}$) of galaxies brighter than $m_\mathrm{AB,1500} \approx 27$ mag may be retrieved with mean absolute error $\Delta f_\mathrm{esc}\approx$0.09(0.12) for 24h (1.5h) JWST/NIRSpec exposures at resolution R=100. For 24h exposure time, even fainter galaxies ($m_\mathrm{AB,1500} < 28.5$ mag) can be processed with $\Delta f_\mathrm{esc}\approx$0.14. This framework simultaneously estimates the redshift of these galaxies with a relative error less than 0.03 for both 24h ($m_\mathrm{AB,1500} < 28.5$ mag) and 1.5h ($m_\mathrm{AB,1500} < 27$ mag) exposure times. We also consider scenarios where just a minor fraction of galaxies attain high $f_\mathrm{esc}$ and present the conditions required for detecting a subpopulation of high $f_\mathrm{esc}$ galaxies within the dataset., Comment: 10 pages, 7 figures. Accepted to be published in MNRAS

New challenges in the modelling of galaxy spectra are bound to emerge as upcoming telescopes like the James Webb Space Telescope will allow us to detect galaxies at fainter flux levels and higher redshifts than ever before. Here, we highlight three modelling problems that may become relevant for upcoming observations in the rest-frame ultraviolet/optical at the high-redshift frontier: stellar initial mass function sampling effects, Population III signatures and the leakage of ionizing radiation into the intergalactic medium., Comment: 9 pages, 4 figures, invited review at IAU symposium 341, "Challenges in Panchromatic Galaxy Modelling with Next Generation Facilities"

Context. Characterizing the galaxy population in the early Universe holds the key to understanding the evolution of these objects and the role they played in cosmic reionization. However, there have been very few observations at the very highest redshifts to date. Aims. In order to shed light on the properties of galaxies in the high-redshift Universe and their interstellar media, we observe the Lyman-α emitting galaxy z7_GSD_3811 at z = 7.664 with bands 6 and 8 at the Atacama Large Millimeter/submillimeter Array (ALMA). Methods. We target the far-infrared [O III] 88 μm and [C II] 158 μm emission lines and dust continuum in the star-forming galaxy z7_GSD_3811 with ALMA. We combine these measurements with earlier observations in the rest-frame ultraviolet (UV) in order to characterize the object and compare the results to those of earlier studies that observed [O III] and [C II] emission in high-redshift galaxies. Results. The [O III] 88 μm and [C II] 158 μm emission lines are undetected at the position of z7_GSD_3811, with 3σ upper limits of 1.6 × 108 L⊙ and 4.0 × 107 L⊙, respectively. We do not detect any dust continuum in band 6 nor band 8. The measured rms in the band 8 and band 6 continua are 26 and 9.9 μJy beam−1, respectively. Similar to several other high-redshift galaxies, z7_GSD_3811 exhibits low [C II] emission for its star formation rate compared to local galaxies. Furthermore, our upper limit on the [O III] line luminosity is lower than the previously observed [O III] lines in high-redshift galaxies with similar UV luminosities. Our ALMA band 6 and 8 dust continuum observations imply that z7_GSD_3811 likely has a low dust content, and our non-detections of the [O III] and [C II] lines could indicate that z7_GSD_3811 has a low metallicity (Z ≲ 0.1 Z⊙).

Recent analysis of strongly-lensed sources in the Hubble Frontier Fields indicates that the rest-frame UV luminosity function of galaxies at $z=$6--8 rises as a power law down to $M_\mathrm{UV}=-15$, and possibly as faint as -12.5. We use predictions from a cosmological radiation hydrodynamic simulation to map these luminosities onto physical space, constraining the minimum dark matter halo mass and stellar mass that the Frontier Fields probe. While previously-published theoretical studies have suggested or assumed that early star formation was suppressed in halos less massive than $10^9$--$10^{11} M_\odot$, we find that recent observations demand vigorous star formation in halos at least as massive as (3.1, 5.6, 10.5)$\times10^9 M_\odot$ at $z=(6,7,8)$. Likewise, we find that Frontier Fields observations probe down to stellar masses of (8.1, 18, 32)$\times10^6 M_\odot$; that is, they are observing the likely progenitors of analogues to Local Group dwarfs such as Pegasus and M32. Our simulations yield somewhat different constraints than two complementary models that have been invoked in similar analyses, emphasizing the need for further observational constraints on the galaxy-halo connection., Comment: 8 pages, 5 figures, re-submitted to MNRAS after incorporating referee's comments

The collapse of supermassive primordial stars in hot, atomically-cooled halos may have given birth to the first quasars at $z \sim$ 15 - 20. Recent numerical simulations of these rapidly accreting stars reveal that they are cool, red hypergiants shrouded by dense envelopes of pristine atomically-cooled gas at 6,000 - 8,000 K, with luminosities $L$ $\gtrsim$ 10$^{10}$ L$_{\odot}$. Could such luminous but cool objects be detected as the first stage of quasar formation in future near infrared (NIR) surveys? We have now calculated the spectra of supermassive primordial stars in their birth envelopes with the Cloudy code. We find that some of these stars will be visible to JWST at $z \lesssim$ 20 and that with modest gravitational lensing Euclid and WFIRST could detect them out to $z \sim$ 10 - 12. Rather than obscuring the star, its accretion envelope enhances its visibility in the NIR today by reprocessing its short-wavelength flux into photons that are just redward of the Lyman limit in the rest frame of the star., Comment: 6 pages, 4 figures, accepted by ApJL

Previous Very Long Baseline Interferometry (VLBI) observations of the quasar B1152+199 at 5GHz has revealed two images of a strongly lensed jet with seemingly discordant morphologies. Whereas the jet appears straight in one of the images, the other exhibits slight curvature on milliarcsecond scales. This is unexpected from the lensing solution and has been interpreted as possible evidence for secondary, small-scale lensing (millilensing) by a compact object with a mass of $~10^5$-$10^7\ M_\odot$ located close to the curved image. The probability for such a superposition is extremely low unless the millilens population has very high surface number density. Here, we revisit the case for millilensing in B1152+199 by combining new global-VLBI data at 8.4GHz with two datasets from the European VLBI Network (EVN) at 5GHz (archival) and at 22GHz (new dataset), and the previously published 5GHz Very Long Baseline Array (VLBA) data. We find that the new data with a more circular synthesized beam, exhibits no apparent milliarcsecond-scale curvature in image B. Various observations of the object spanning $\sim$15 years apart enable us to improve the constraints on lens system (thanks also to the improved astrometry resulting from 22GHz observations) to the point that the only plausible explanation left for the apparent curvature is the artifact due to the shape of the synthesized beam., 9 pages, 4 figures, submitted to MNRAS

We present the Technicolor Dawn simulations, a suite of cosmological radiation-hydrodynamic simulations of the first 1.2 billion years. By modeling a spatially-inhomogeneous UVB on-the-fly with 24 frequencies and resolving dark matter halos down to $10^8 M_\odot$ within 12 $h^{-1}$ Mpc volumes, our simulations unify observations of the intergalactic and circumgalactic media, galaxies, and reionization into a common framework. The only empirically-tuned parameter, the fraction $f_{\mathrm{esc,gal}}(z)$ of ionizing photons that escape the interstellar medium, is adjusted to match observations of the Lyman-$\alpha$ forest and the cosmic microwave background. With this single calibration, our simulations reproduce the history of reionization; the stellar mass-star formation rate relation of galaxies; the number density and metallicity of damped Lyman-$\alpha$ absorbers (DLAs) at $z\sim5$; the abundance of weak metal absorbers; the ultraviolet background (UVB) amplitude; and the Lyman-$\alpha$ flux power spectrum at $z=5.4$. The galaxy stellar mass and UV luminosity functions are underproduced by $\leq2\times$, suggesting an overly vigorous feedback model. The mean transmission in the Lyman-$\alpha$ forest is underproduced at $z4$ Ryd UVB. Solving this problem by increasing metal yields would overproduce both weak absorbers and DLA metallicities. Instead, the observed strength of high-ionization emission from high-redshift galaxies and absorption from their environments suggest that the ionizing flux from conventional stellar population models is too soft., Comment: 24 pages, 17 figures, accepted to MNRAS

We use observations of spatially-aligned C II, C IV, Si II, Si IV, and O I absorbers to probe the slope and intensity of the ultraviolet background (UVB) at $z \sim 6$. We accomplish this by comparing observations with predictions from a cosmological hydrodynamic simulation using three trial UVBs applied in post-processing: a spectrally soft, fluctuating UVB calculated using multi-frequency radiative transfer; a soft, spatially-uniform UVB; and a hard, spatially-uniform "quasars-only" model. When considering our paired high-ionization absorbers (C IV/Si IV), the observed statistics strongly prefer the hard, spatially-uniform UVB. This echoes recent findings that cosmological simulations generically underproduce strong C IV absorbers at $z>5$. A single low/high ionization pair (Si II/Si IV), by contrast, shows a preference for the HM12 UVB, while two more (C II/C IV and O I/C IV) show no preference for any of the three UVBs. Despite this, future observations of specific absorbers, particularly Si IV/C IV, with next-generation telescopes probing to lower column densities should yield tighter constraints on the UVB., 13 pages, 4 figures, accepted by MNRAS

In this paper, we study the effects of a recent drop in star formation rate (SFR) on the spectra of epoch of reionization (EoR) galaxies, and the resulting degeneracy with the spectral features produced by extreme Lyman continuum leakage. In order to study these effects in the wavelength range relevant for the upcoming James Webb Space Telescope (JWST), we utilize synthetic spectra of simulated EoR galaxies from cosmological simulations together with synthetic spectra of partially quenched mock galaxies. We find that rapid declines in the SFR of EoR galaxies could seriously affect the applicability of methods that utilize the equivalent width of Balmer lines and the ultraviolet spectral slope to assess the escape fraction of EoR galaxies. In order to determine if the aforementioned degeneracy can be avoided by using the overall shape of the spectrum, we generate mock NIRCam observations and utilize a classification algorithm to identify galaxies that have undergone quenching. We find that while there are problematic cases, JWST/NIRCam or NIRSpec should be able to reliably identify galaxies with redshifts $z\sim 7$ that have experienced a significant decrease in the SFR (by a factor 10-100) in the past 50-100 Myr with a success rate $\gtrsim 85\%$. We also find that uncertainties in the dust-reddening effects on EoR galaxies significantly affect the performance of the results of the classification algorithm. We argue that studies that aim to characterize the dust extinction law most representative in the EoR would be extremely useful., 10 pages, 7 figures. Accepted for publication in MNRAS

Observations suggest that CII was more abundant than CIV in the intergalactic medium towards the end of the hydrogen reionization epoch. This transition provides a unique opportunity to study the enrichment history of intergalactic gas and the growth of the ionizing background (UVB) at early times. We study how carbon absorption evolves from z=10-5 using a cosmological hydrodynamic simulation that includes a self-consistent multifrequency UVB as well as a well-constrained model for galactic outflows to disperse metals. Our predicted UVB is within 2-4 times that of Haardt & Madau (2012), which is fair agreement given the uncertainties. Nonetheless, we use a calibration in post-processing to account for Lyman-alpha forest measurements while preserving the predicted spectral slope and inhomogeneity. The UVB fluctuates spatially in such a way that it always exceeds the volume average in regions where metals are found. This implies both that a spatially-uniform UVB is a poor approximation and that metal absorption is not sensitive to the epoch when HII regions overlap globally even at column densites of 10^{12} cm^{-2}. We find, consistent with observations, that the CII mass fraction drops to low redshift while CIV rises owing the combined effects of a growing UVB and continued addition of carbon in low-density regions. This is mimicked in absorption statistics, which broadly agree with observations at z=6-3 while predicting that the absorber column density distributions rise steeply to the lowest observable columns. Our model reproduces the large observed scatter in the number of low-ionization absorbers per sightline, implying that the scatter does not indicate a partially-neutral Universe at z=6., 16 pages, 14 figures, accepted to MNRAS

Direct collapse black holes forming in pristine, atomically-cooling haloes at $z \approx 10-20$ may act as the seeds of supermassive black holes (BH) at high redshifts. In order to create a massive BH seed, the host halo needs to be prevented from forming stars. H$_2$ therefore needs to be irradiated by a large flux of Lyman-Werner (LW) UV photons in order to suppress H$_2$ cooling. A key uncertainty in this scenario is the escape fraction of LW radiation from first galaxies, the dominant source of UV photons at this epoch. To better constrain this escape fraction, we have performed radiation-hydrodynamical simulations of the growth of HII regions and their associated photodissociation regions in the first galaxies using the ZEUS-MP code. We find that the LW escape fraction crucially depends on the propagation of the ionisation front (I-front). For an R-type I-front overrunning the halo, the LW escape fraction is always larger than 95%. If the halo recombines later from the outside--in, due to a softened and weakened spectrum, the LW escape fraction in the rest-frame of the halo (the near-field) drops to zero. A detailed and careful analysis is required to analyse slowly moving, D-type I-fronts, where the escape fraction depends on the microphysics and can be as small as 3% in the near-field and 61% in the far-field or as large as 100% in both the near-field and the far-field., 14 pages, 10 figures, 7 tables, accepted for publication in MNRAS

In order to empirically determine the time-scale and environmental dependence of stellar cluster disruption, we have undertaken an analysis of the unprecedented multipointing (seven), multiwaveleng ...

We present the Atacama Large Millimeter/submillimeter Array (ALMA) detection of the [O III] 88 $\mu$m line and rest-frame 90 $\mu$m dust continuum emission in a Y-dropout Lyman break galaxy (LBG), MACS0416_Y1, lying behind the Frontier Field cluster MACS J0416.1-2403. This [O III] detection confirms the LBG with a spectroscopic redshift of $z = 8.3118 \pm 0.0003$, making this object one of the furthest galaxies ever identified spectroscopically. The observed 850 $\mu$m flux density of $137 \pm 26$ $\mu$Jy corresponds to a de-lensed total infrared (IR) luminosity of $L_{\rm IR} = (1.7 \pm 0.3) \times 10^{11} L_{\odot}$ if assuming a dust temperature of $T_{\rm dust} = 50$ K and an emissivity index of $\beta = 1.5$, yielding a large dust mass of $4 \times 10^6 M_{\odot}$. The ultraviolet-to-far IR spectral energy distribution modeling where the [O III] emissivity model is incorporated suggests the presence of a young ($\tau_{\rm age} \approx 4$ Myr), star-forming (SFR $\approx 60 M_{\odot}$ yr$^{-1}$), moderately metal-polluted ($Z \approx 0.2 Z_{\odot}$) stellar component with a mass of $M_{\rm star} = 3 \times 10^8 M_{\odot}$. An analytic dust mass evolution model with a single episode of star-formation does not reproduce the metallicity and dust mass in $\tau_{\rm age} \approx 4$ Myr, suggesting a pre-existing evolved stellar component with $M_{\rm star} \sim 3 \times 10^9 M_{\odot}$ and $\tau_{\rm age} \sim 0.3$ Gyr as the origin of the dust mass., Comment: Accepted by ApJ. 18 pages, 10 figures, 5 tables

We present deep optical and near-infrared (NIR) UBVRIHKs imaging data for 24 blue compact galaxies (BCGs). The individual exposure times are on average similar to 40 min in the optical (B) and simi ...

Using four different suites of cosmological simulations, we generate synthetic spectra for galaxies with different Lyman continuum escape fractions (fesc) at redshifts z=7-9, in the rest-frame wavelength range relevant for the James Webb Space Telescope (JWST) NIRSpec instrument. By investigating the effects of realistic star formation histories and metallicity distributions on the EW(Hb)-beta diagram (previously proposed as a tool for identifying galaxies with very high fesc), we find that neither of these effects are likely to jeopardize the identification of galaxies with extreme Lyman continuum leakage. Based on our models, we expect essentially all z=7-9 galaxies that exhibit rest-frame EW(Hb)< 30 �� to have fesc>0.5. Incorrect assumptions concerning the ionizing fluxes of stellar populations or the dust properties of z>6 galaxies can in principle bias the selection, but substantial model deficiencies of this type will at the same time reveal themselves as an offset between the observed and simulated distribution of z>6 galaxies in the EW(Hb)-beta diagram. Such offsets would thereby allow JWST/NIRSpec measurements of these observables to serve as input for further model refinement., 12 pages, 8 figures, v.2: ApJ, accepted. Model grids are available from http://www.astro.uu.se/~ez/lycan/lycan.html

We propose that one of the sources in the recently detected system CR7 by Sobral et al. (2015) through spectro-photometric measurements at $z = 6.6$ harbors a direct collapse blackhole (DCBH). We argue that the LW radiation field required for direct collapse in source A is provided by sources B and C. By tracing the LW production history and star formation rate over cosmic time for the halo hosting CR7 in a $\Lambda$CDM universe, we demonstrate that a DCBH could have formed at $z\sim 20$. The spectrum of source A is well fit by nebular emission from primordial gas around a BH with MBH $\sim 4.4 \times 10^6 \ \rm M_{\odot}$ accreting at a 40 % of the Eddington rate, which strongly supports our interpretation of the data. Combining these lines of evidence, we argue that CR7 might well be the first DCBH candidate., Comment: 8 pages, 4 figures. Accepted for publication in MNRAS

Recent observations of galaxies at $z \gtrsim 7$, along with the low value of the electron scattering optical depth measured by the Planck mission, make galaxies plausible as dominant sources of ionizing photons during the epoch of reionization. However, scenarios of galaxy-driven reionization hinge on the assumption that the average escape fraction of ionizing photons is significantly higher for galaxies in the reionization epoch than in the local Universe. The NIRSpec instrument on the James Webb Space Telescope (JWST) will enable spectroscopic observations of large samples of reionization-epoch galaxies. While the leakage of ionizing photons will not be directly measurable from these spectra, the leakage is predicted to have an indirect effect on the spectral slope and the strength of nebular emission lines in the rest-frame ultraviolet and optical. Here, we apply a machine learning technique known as lasso regression on mock JWST/NIRSpec observations of simulated $z=7$ galaxies in order to obtain a model that can predict the escape fraction from JWST/NIRSpec data. Barring systematic biases in the simulated spectra, our method is able to retrieve the escape fraction with a mean absolute error of $\Delta f_{\mathrm{esc}} \approx 0.12$ for spectra with $S/N\approx 5$ at a rest-frame wavelength of 1500 {\AA} for our fiducial simulation. This prediction accuracy represents a significant improvement over previous similar approaches., Comment: 13 pages, 11 figures. Accepted for publication in ApJ

We test the predictions of spectral synthesis models based on seven different massive-star prescriptions against Legacy ExtraGalactic UV Survey (LEGUS) observations of eight young massive clusters in two local galaxies, NGC 1566 and NGC 5253, chosen because predictions of all seven models are available at the published galactic metallicities. The high angular resolution, extensive cluster inventory and full near-ultraviolet to near-infrared photometric coverage make the LEGUS dataset excellent for this study. We account for both stellar and nebular emission in the models and try two different prescriptions for attenuation by dust. From Bayesian fits of model libraries to the observations, we find remarkably low dispersion in the median E(B-V) (~0.03 mag), stellar masses (~10^4 M_\odot) and ages (~1 Myr) derived for individual clusters using different models, although maximum discrepancies in these quantities can reach 0.09 mag and factors of 2.8 and 2.5, respectively. This is for ranges in median properties of 0.05-0.54 mag, 1.8-10x10^4 M_\odot and 1.6-40 Myr spanned by the clusters in our sample. In terms of best fit, the observations are slightly better reproduced by models with interacting binaries and least well reproduced by models with single rotating stars. Our study provides a first quantitative estimate of the accuracies and uncertainties of the most recent spectral synthesis models of young stellar populations, demonstrates the good progress of models in fitting high-quality observations, and highlights the needs for a larger cluster sample and more extensive tests of the model parameter space., Accepted for publication in MNRAS (14 Jan. 2016). 30 pages, 16 figures, 9 tables

We present IFU observations with MUSE@VLT and deep imaging with FORS@VLT of a dwarf galaxy recently formed within the giant collisional HI ring surrounding NGC 5291. This TDG-like object has the characteristics of typical z=1-2 gas-rich spiral galaxies: a high gas fraction, a rather turbulent clumpy ISM, the absence of an old stellar population, a moderate metallicity and star formation efficiency. The MUSE spectra allow us to determine the physical conditions within the various complex substructures revealed by the deep optical images, and to scrutinize at unprecedented spatial resolution the ionization processes at play in this specific medium. Starburst age, extinction and metallicity maps of the TDG and surrounding regions were determined using the strong emission lines Hbeta, [OIII], [OI], [NII], Halpha and [SII] combined with empirical diagnostics. Discrimination between different ionization mechanisms was made using BPT--like diagrams and shock plus photoionization models. Globally, the physical conditions within the star--forming regions are homogeneous, with in particular an uniform half-solar oxygen abundance. At small scales, the derived extinction map shows narrow dust lanes. Regions with atypically strong [OI] emission line immediately surround the TDG. The [OI] / Halpha ratio cannot be easily accounted for by photoionization by young stars or shock models. At larger distances from the main star--forming clumps, a faint diffuse blue continuum emission is observed, both with the deep FORS images and MUSE data. It does not have a clear counterpart in the UV regime probed by GALEX. A stacked spectrum towards this region does not exhibit any emission line, excluding faint levels of star formation, nor stellar absorption lines that might have revealed the presence of old stars. Several hypotheses are discussed for the origin of these intriguing features., Comment: 13 pages, 15 figures, accepted for publication in A&A

The study of cosmology, galaxy formation and exoplanets has now advanced to a stage where a cosmic inventory of terrestrial planets may be attempted. By coupling semi-analytic models of galaxy formation to a recipe that relates the occurrence of planets to the mass and metallicity of their host stars, we trace the population of terrestrial planets around both solar-mass (FGK type) and lower-mass (M dwarf) stars throughout all of cosmic history. We find that the mean age of terrestrial planets in the local Universe is $7\pm{}1$ Gyr for FGK hosts and $8\pm{}1$ Gyr for M dwarfs. We estimate that hot Jupiters have depleted the population of terrestrial planets around FGK stars by no more than $\approx 10\%$, and that only $\approx 10\%$ of the terrestrial planets at the current epoch are orbiting stars in a metallicity range for which such planets have yet to be confirmed. The typical terrestrial planet in the local Universe is located in a spheroid-dominated galaxy with a total stellar mass comparable to that of the Milky Way. When looking at the inventory of planets throughout the whole observable Universe, we argue for a total of $\approx 1\times 10^{19}$ and $\approx 5\times 10^{20}$ terrestrial planets around FGK and M stars, respectively. Due to light travel time effects, the terrestrial planets on our past light cone exhibit a mean age of just $1.7\pm 0.2$ Gyr. These results are discussed in the context of cosmic habitability, the Copernican principle and searches for extraterrestrial intelligence at cosmological distances., Comment: 11 pages, 8 figures. v.2: Accepted for publication in ApJ. Some changes in quantitative results compared to v.1, mainly due to differences in IMF assumptions

Small galaxies consisting entirely of population III (pop III) stars may form at high redshifts, and could constitute one of the best probes of such stars. Here, we explore the prospects of detecting gravitationally lensed pop III galaxies behind the galaxy cluster J0717.5+3745 (J0717) with both the Hubble Space Telescope (HST) and the upcoming James Webb Space Telescope (JWST). By projecting simulated catalogs of pop III galaxies at z~7-15 through the J0717 magnification maps, we estimate the lensed number counts as a function of flux detection threshold. We find that the ongoing HST survey CLASH, targeting a total of 25 galaxy clusters including J0717, potentially could detect a small number of pop III galaxies if ~1% of the baryons in these systems have been converted into pop III stars. Using JWST exposures of J0717, this limit can be pushed to ~0.1% of the baryons. Ultra-deep JWST observations of unlensed fields are predicted to do somewhat worse, but will be able to probe pop III galaxies with luminosities intermediate between those detectable in HST/CLASH and in JWST observations of J0717. We also explain how current measurements of the galaxy luminosity function at z=7-10 can be used to constrain pop III galaxy models with very high star formation efficiencies (~10% of the baryons converted into pop III stars).

We report the discovery of a ring-like cluster complex in the starburst galaxy NGC 2146. The Ruby Ring, so named due to its appearance, shows a clear ring-like distribution of star clusters around a central object. It is located in one of the tidal streams which surround the galaxy. NGC 2146 is part of the Snapshot Hubble U-band Cluster Survey (SHUCS). The WFC3/F336W data has added critical information to the available archival Hubble Space Telescope imaging set of NGC 2146, allowing us to determine ages, masses, and extinctions of the clusters in the Ruby Ring. These properties have then been used to investigate the formation of this extraordinary system. We find evidence of a spatial and temporal correlation between the central cluster and the clusters in the ring. The latter are about 4 Myr younger than the central cluster, which has an age of 7 Myr. This result is supported by the H alpha emission which is strongly coincident with the ring, and weaker at the position of the central cluster. From the derived total H alpha luminosity of the system we constrain the star formation rate density to be quite high, e.g. ~ 0.47 Msun/yr/kpc^2. The Ruby Ring is the product of an intense and localised burst of star formation, similar to the extended cluster complexes observed in M51 and the Antennae, but more impressive because is quite isolated. The central cluster contains only 5 % of the total stellar mass in the clusters that are determined within the complex. The ring-like morphology, the age spread, and the mass ratio support a triggering formation scenario for this complex. We discuss the formation of the Ruby Ring in a "collect & collapse" framework. The predictions made by this model agree quite well with the estimated bubble radius and expansion velocity produced by the feedback from the central cluster, making the Ruby Ring an interesting case of triggered star formation.

We study the stellar cluster population in two adjacent fields in the nearby, face-on spiral galaxy, M83, using WFC3/HST imaging. The clusters are selected through visual inspection to be centrally concentrated, symmetric, and resolved on the images, which allows us to differentiate between clusters and likely unbound associations. We compare our sample with previous studies and show that the differences between the catalogues are largely due to the inclusion of large numbers of diffuse associations within previous catalogues. The luminosity function of the clusters is well approximated by a power-law with index, -2, over most of the observed range, however a steepening is seen at M_V = -9.3 and -8.8 in the inner and outer fields, respectively. Additionally, we show that the cluster population is inconsistent with a pure power-law mass distribution, but instead exhibits a truncation at the high mass end. If described as a Schechter function, the characteristic mass is 1.6 and 0.5 * 10^5 Msun, for the inner and outer fields, respectively, in agreement with previous estimates of other cluster populations in spiral galaxies. Comparing the predictions of the mass independent disruption (MID) and mass dependent disruption (MDD) scenarios with the observed distributions, we find that both models can accurately fit the data. However, for the MID case, the fraction of clusters destroyed (or mass lost) per decade in age is dependent on the environment, hence, the age/mass distributions of clusters are not universal. In the MDD case, the disruption timescale scales with galactocentric distance (being longer in the outer regions of the galaxy) in agreement with analytic and numerical predictions. Finally, we discuss the implications of our results on other extragalactic surveys, focussing on the fraction of stars that form in clusters and the need (or lack thereof) for infant mortality.

Population III galaxies, made partly or exclusively of metal-free stars, are predicted to exist at high redshifts and may produce very strong Lya emission. A substantial fraction of these Lya photons are likely absorbed in the intergalactic medium at z>6, but recent simulations suggest that significant Lya emission may be detectable up to z~8.5, i.e. well into the reionization epoch. Here, we argue that high-redshift population III galaxies with strong Lya emission can be identified in Hubble Space Telescope imaging data because of their unusual colours. We quantify this effect in some of the filters used in Y-band dropout searches for galaxies at z~8 and find that population III galaxies with high Lya fluxes may exhibit much bluer J-H colours at z=8-10 than any normal type of galaxy at these redshifts. This colour signature can arise even if pop III stars account for as little as ~1e-3 to ~1e-2 of the stellar mass in these galaxies. Some of the anomalously blue objects reported in current Y-band dropout samples do in fact meet the colour criteria for Lya-emitting population III galaxies.

The numerous and massive young star clusters in blue compact galaxies (BCGs) are used to investigate the properties of their hosts. We test whether BCGs follow claimed relations between cluster populations and their hosts, such as the the fraction of the total luminosity contributed by the clusters as function of the mean star formation rate density; the $V$ band luminosity of the brightest youngest cluster as related to the mean host star formation rate; and the cluster formation efficiency (i.e., the fraction of star formation happening in star clusters) versus the density of the SFR. We find that BCGs follow the trends, supporting a scenario where cluster formation and environmental properties of the host are correlated. They occupy, in all the diagrams, the regions of higher SFRs, as expected by the extreme nature of the starbursts operating in these systems. We find that the star clusters contribute almost to the 20 % of the UV luminosity of the hosts. We suggest that the BCG starburst environment has most likely favoured the compression and collapse of the giant molecular clouds, enhancing the local star formation efficiency, so that massive clusters have been formed. The estimated cluster formation efficiency supports this scenario. BCGs have a cluster formation efficiency comparable to luminous IR galaxies and spiral starburst nuclei (the averaged value is about 35 %) which is much higher than the 8 - 10 % reported for quiescent spirals and dwarf star-forming galaxies.

Using multi-wavelength imaging from the Wide Field Camera 3 on the Hubble Space Telescope we study the stellar cluster populations of two adjacent fields in the nearby face-on spiral galaxy, M83. The observations cover the galactic centre and reach out to ~6 kpc, thereby spanning a large range of environmental conditions, ideal for testing empirical laws of cluster disruption. The clusters are selected by visual inspection to be centrally concentrated, symmetric, and resolved on the images. We find that a large fraction of objects detected by automated algorithms (e.g. SExtractor or Daofind) are not clusters, but rather are associations. These are likely to disperse into the field on timescales of tens of Myr due to their lower stellar densities and not due to gas expulsion (i.e. they were never gravitationally bound). We split the sample into two discrete fields (inner and outer regions of the galaxy) and search for evidence of environmentally dependent cluster disruption. Colour-colour diagrams of the clusters, when compared to simple stellar population models, already indicate that a much larger fraction of the clusters in the outer field are older by tens of Myr than in the inner field. This impression is quantified by estimating each cluster's properties (age, mass, and extinction) and comparing the age/mass distributions between the two fields. Our results are inconsistent with "universal" age and mass distributions of clusters, and instead show that the ambient environment strongly affects the observed populations.

We have investigated the starburst properties of the luminous blue compact galaxy ESO 185-IG13. The galaxy has been imaged with the high resolution cameras onboard to the Hubble Space Telescope. From the UV to the IR, the data reveal a system shaped by hundreds of young star clusters, and fine structures, like a tidal stream and a shell. The presence of numerous clusters and the perturbed morphology indicate that the galaxy has been involved in a recent merger event. Using previous simulations of shell formation in galaxy mergers we constrain potential progenitors of ESO 185-IG13. The analysis of the star cluster population is used to investigate the properties of the present starburst and to date the final merger event, which has produced hundreds of clusters younger than 100 Myr. We have found a peak of cluster formation only 3.5 Myr old. A large fraction of these clusters will not survive after 10-20 Myr, due to the "infant mortality" caused by gas expulsion. However, this sample of clusters represents an unique chance to investigate the youngest phases of cluster evolution. As already observed in the analog blue compact galaxy Haro 11, a fraction of young clusters are affected by a flux excess at wavelengths longer than 8000 \AA. Ages, masses, and extinctions of clusters with this NIR excess are estimated from UV and optical data. We discuss similarities and differences of the observed NIR excess in ESO 185-IG13 clusters with other cases in the literature. The cluster ages and masses are used to distinguish among the potential causes of the excess. We observe, as in Haro 11, that the use of the IR and the (commonly used) I band data results in overestimates of age and mass in clusters affected by the NIR excess. This has important implications for a number of related studies of star clusters.

We present a HST multiband analysis of the large population of star clusters in the blue compact galaxy (BCG) Mrk 930. We have reconstructed the spectral energy distributions of the star clusters and estimated age, mass, and extinction for a representative sample. We observe a very young cluster population with 70% of the systems formed less than 10 Myr ago. The peak in the star cluster age distribution at 4 Myr is corroborated by the presence of Wolf-Rayet spectral features, and by the observed optical and IR lines ratios [OIII]/H_beta and [NeIII]/[NeII]. The recovered extinction in these very young clusters shows large variations, with a decrease at older ages. It is likely that our analysis is limited to the optically brightest objects (i.e. systems only partially embedded in their natal cocoons). We map the extinction across the galaxy using low-resolution spectra and the H_alpha/H_beta ratio, as obtained from ground-based narrow band imaging. We find that the mean optical extinction derived in the starburst regions is close to the averaged value observed in the clusters, but locally, do not trace the more extinguished clusters. Previous HST studies of BCGs have revealed a population of young and extremely red super star clusters. We detect a considerable fraction of clusters affected by a red excess also in Mrk 930. The nature of the red excess, which turns up at near-IR wavelengths remains unknown. We compare the cluster and the star formation history, the latter derived from the fit of spectral population synthesis models to the spectra. We find a general agreement between the two independently estimated quantities. Using the cluster properties we perform a study of the host environmental properties. We find that the cluster formation efficiency is significantly higher, suggesting a key role of the environment for the formation of these massive objects.[Abridged]

A star enshrouded in a Dyson sphere with high covering fraction may manifest itself as an optically subluminous object with a spectrophotometric distance estimate significantly in excess of its parallax distance. Using this criterion, the Gaia mission will in coming years allow for Dyson-sphere searches that are complementary to searches based on waste-heat signatures at infrared wavelengths. A limited search of this type is also possible at the current time, by combining Gaia parallax distances with spectrophotometric distances from ground-based surveys. Here, we discuss the merits and shortcomings of this technique and carry out a limited search for Dyson-sphere candidates in the sample of stars common to Gaia Data Release 1 and RAVE Data Release 5. We find that a small fraction of stars indeed display distance discrepancies of the type expected for nearly complete Dyson spheres. To shed light on the properties of objects in this outlier population, we present follow-up high-resolution spectroscopy for one of these stars, the late F-type dwarf TYC 6111-1162-1. The spectrophotometric distance of this object is about twice that derived from its Gaia parallax, and there is no detectable infrared excess. While our analysis largely confirms the stellar parameters and the spectrophotometric distance inferred by RAVE, a plausible explanation for the discrepant distance estimates of this object is that the astrometric solution has been compromised by an unseen binary companion, possibly a rather massive white dwarf ($\approx 1\ M_\odot$). This scenario can be further tested through upcoming Gaia data releases., Comment: 9 pages, 4 figures, v.2: ApJ, in press