Search

We measure the broad impact of galaxy structure on galaxy formation by examining the ongoing star formation and integrated star formation history as revealed through the stellar masses of galaxies at $z < 7$ based on JWST CEERS data from the Extended Groth Strip (EGS). Using the morphological catalog of 3965 visually classified JWST galaxies from Ferreira et al. (2023), we investigate the evolution of stars, and when they form, as a function of morphological type as well as galaxies classified as passive and starburst through spectral energy distributions. Although disk galaxies dominate the structures of galaxies at $z < 7$, we find that these disks are in general either `passive', or on the main-sequence of star formation, and do not contain a large population of starburst galaxies. We also find no significant correlation between morphological type and the star formation rate or colours of galaxies at $z < 7$. In fact, we find that the morphologically classified `spheroids' tend to be blue and are not found to be predominately passive systems at $z > 1.5$. We also find that the stellar mass function for disk galaxies does not evolve significantly during this time, whereas other galaxy types, such as the peculiar population, evolve dramatically, declining at lower redshifts. This indicates that massive peculiars are more common at higher redshifts. We further find that up to $z \sim 7$, the specific star formation rate (sSFR) does not vary with visual morphology, but strongly depends on stellar mass and internal galaxy mass density. This demonstrates that at early epochs galaxy assembly is a mass-driven, rather than a morphologically-driven, process. Quenching of star formation is therefore a mass-dominated process throughout the universe's history, likely due to the presence of supermassive black holes., Comment: MNRAS in press, 20 pages

We present an analysis of rest-frame UV continuum slopes, $\beta$, using a sample of 1011 galaxies at $6.5

We develop a photometric search method for identifying smouldering galaxies at $5< z < 8$, which are defined to have weak emission lines and thus generally have low specific star formation rates and may even be in a state of (temporary) quiescence. The deep NIRCam imaging (${\sim}29.5$ AB mag, 5$\sigma$) from the JADES second data release is essential for finding these systems, as they are faint, relatively quiescent dwarf galaxies ($M_* \sim 10^{8}$-$10^9$ $\mathrm{M}_\odot)$ in the Epoch of Reionisation (EoR). Moreover, medium-band imaging is key, enabling a clear identification of the lack of emission lines in these galaxies, thus betraying their dormant flame. Owing to the young age of the Universe, combined with the likely bursty star formation in these first dwarf galaxies, conventional colour-selection methods like the UVJ diagram likely miss a large fraction of the quiescent population in the EoR. Indeed, we find that smouldering galaxies constitute a considerable fraction (0.10-0.35) of the EoR dwarf galaxy population ($M_* \sim 10^{8}$-$10^{9}$ $\mathrm{M}_\odot$). As predicted by simulations, these first dwarf galaxies are fragile, the star formation in their shallow potential wells easily snuffed out by feedback-driven winds triggered by secular or merger-driven starbursts, with the smouldering fraction increasing with decreasing stellar mass. Finally, we provide observational constraints on the smouldering galaxy comoving number density (${\sim}10^{-4}$-$10^{-5}$ dex$^{-1}$ Mpc$^{-3}$), which, although hampered by incompleteness, should aid in our understanding of the primordial baryon cycle, as current simulations greatly disagree on whether these systems are rare (${\sim}1\%$) or common (${\sim}50\%$) in the EoR., Comment: 22 pages, 17 figures, 2 tables

It is well established that stellar discs are destabilized by sharp features in their phase space, driving recurrent spiral modes. We explore the extent to which surface-density breaks in disc galaxies - which represent sharp changes in the gradient of the disc density - drive new spiral modes. We employ linear perturbation theory to investigate how disc breaks alter the eigenmode spectrum of an otherwise pure exponential disc. We find that the presence of a density break gives rise to a set of new, vigorously growing, modes. For a given multiplicity, these edge modes occur in pairs, with closely separated resonances between each pair. The growth rate of edge modes decreases when the break is weakened or moved outward to lower-density regions of the disc. Both down- and up-bending profiles excite edge modes, whose origin can be best understood via the gravitational torques they exert on the underlying disc. When the profile is down-bending (Type II) the faster growing mode is the inner one while in the up-bending (Type III) case the outer mode is faster growing. In both cases the faster growing mode has a corotation almost coincident with the break. We show that the torques of the edge modes tend to smoothen the break.

We utilize deep JWST NIRCam observations for the first direct constraints on the Galaxy Stellar Mass Function (GSMF) at z>10. Our EPOCHS v1 sample includes 1120 galaxy candidates at 6.5 9 than some studies predict, suggesting a rapid buildup of stellar mass in the early Universe. We find no incompatibility between our results and those of standard cosmological models, as suggested previously, although the most massive galaxies may require a high star formation efficiency. We find that the 'Little Red Dot' galaxies dominate the z=7 GSMF at high-masses, necessitating a better understanding of the relative contributions of AGN and stellar emission. We show that assuming a theoretically motivated top-heavy IMF reduces stellar mass by 0.5 dex without affecting fit quality, but our results remain consistent with existing cosmological models with a standard IMF., Comment: 40 pages, 15 figures. Submitted to ApJ. Data is avaible on Github: https://github.com/tHarvey303/EpochsIV. Comments to corresponding author welcome at thomas.harvey-3@manchester.ac.uk

We present an analysis of the quenching of star formation in massive galaxies ($M_* > 10^{9.5} M_\odot$) within the first 0.5 - 3 Gyr of the Universe's history utilizing JWST-CEERS data. We utilize a combination of advanced statistical methods to accurately constrain the intrinsic dependence of quenching in a multi-dimensional and inter-correlated parameter space. Specifically, we apply Random Forest (RF) classification, area statistics, and a partial correlation analysis to the JWST-CEERS data. First, we identify the key testable predictions from two state-of-the-art cosmological simulations (IllustrisTNG & EAGLE). Both simulations predict that quenching should be regulated by supermassive black hole mass in the early Universe. Furthermore, both simulations identify the stellar potential ($\phi_*$) as the optimal proxy for black hole mass in photometric data. In photometric observations, where we have no direct constraints on black hole masses, we find that the stellar potential is the most predictive parameter of massive galaxy quenching at all epochs from $z = 0 - 8$, exactly as predicted by simulations for this sample. The stellar potential outperforms stellar mass, galaxy size, galaxy density, and S\'ersic index as a predictor of quiescence at all epochs probed in JWST-CEERS. Collectively, these results strongly imply a stable quenching mechanism operating throughout cosmic history, which is closely connected to the central gravitational potential in galaxies. This connection is explained in cosmological models via massive black holes forming and growing in deep potential wells, and subsequently quenching galaxies through a mix of ejective and preventative active galactic nucleus (AGN) feedback., Comment: Accepted to ApJ. 32 pages; 10 figures

We present the analysis of the stellar population and star formation history of 181 MIRI selected galaxies at redshift 0-3.5 in the massive galaxy cluster field SMACS J0723.3-7327, commonly referred to as SMACS0723, using the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI). We combine the data with the JWST Near Infrared Camera (NIRCam) catalogue, in conjunction with the Hubble Space Telescope (HST) WFC3/IR and ACS imaging. We find that the MIRI bands capture PAH features and dust emission, significantly enhancing the accuracy of photometric redshift and measurements of the physical properties of these galaxies. The median photo-z's of galaxies with MIRI data are found to have a small 0.1% difference from spectroscopic redshifts and reducing the error by 20 percent. With MIRI data included in SED fits, we find that the measured stellar masses are unchanged, while the star formation rate is systematically lower by 0.1 dex. We also fit the median SED of active galactic nuclei (AGN) and star forming galaxies (SFG) separately. MIRI data provides tighter constraints on the AGN contribution, reducing the typical AGN contributions by ~14 percent. In addition, we also compare the median SED obtained with and without MIRI, and we find that including MIRI data yields steeper optical and UV slopes, indicating bluer colours, lower dust attenuation, and younger stellar populations. In the future, MIRI/MRS will enhance our understanding by providing more detailed spectral information and allowing for the study of specific emission features and diagnostics associated with AGN., Comment: 15 pages, 14 figures and 1 table, submitted to MNRAS

With the successful launch and commissioning of JWST we are now able to routinely spectroscopically probe the rest-frame optical emission of galaxies at $z>6$ for the first time. Amongst the most useful spectral diagnostics used in the optical is the Balmer/4000~\AA\ break; this is, in principle, a diagnostic of the mean ages of composite stellar populations. However, the Balmer break is also sensitive to the shape of the star formation history, the stellar (and gas) metallicity, the presence of nebular continuum emission, and dust attenuation. In this work we explore the origin of the Balmer/4000~\AA\ break using the SYNTHESIZER synthetic observations package. We then make predictions of the Balmer/4000~\AA\ break using the First Light and Reionisation Epoch Simulations (FLARES) at $5

From the TNG50 cosmological simulation we build a sample of 191 well-resolved barred galaxies with a stellar mass $\logMstar > 10$ at $z=0$. We search for box/peanut bulges (BPs) in this sample, finding them in 55\% of cases. We compute $\fbp$, the BP probability for barred galaxies as a function of $\Mstar$, and find that this rises to a plateau, as found in observations of nearby galaxies. The transition mass where $\fbp$ reaches half the plateau value is $\logMstar = 10.14$, consistent with the observational value within measurement errors. We show that this transition in $\fbp$ can be attributed to the youth of the bars at low $\Mstar$, which is a consequence of downsizing of galaxies. Young bars, being generally shorter and weaker, have not yet had time to form BPs. At high mass, while we find a plateau, the value is at $\sim 60\%$ whereas observations saturate at $100\%$. We attribute this difference to excessive heating in TNG50, due to merger activity and to numerical resolution effects. BPs in TNG50 tend to occur in galaxies with more quiescent merger histories. As a result, the main driver of whether a bar hosts a BP in TNG50 is not the galaxy mass, but how long and strong the bar is. Separating the BP sample into those that have visibly buckled and those that have not, we find that fully half of BP galaxies show clear signs of buckling, despite the excessive heating and limited vertical resolution of TNG50., Comment: 20 pages, 24 figures, 5 tables. Submitted to MNRAS

JWST has now made it possible to probe the rest-frame optical line emission of high-redshift galaxies extending to z~9, and potentially beyond. To aid in the interpretation of these emerging constraints, in this work we explore predictions for [OIII] emission in high-redshift galaxies using the First Light and Reionisation Epoch Simulations (FLARES). We produce predictions for the [OIII] luminosity function, its correlation with the UV luminosity, and the distribution of equivalent widths (EWs). We also explore how the [OIII] EW correlates with physical properties including specific star formation rate, metallicity, and dust attenuation. Our predictions are largely consistent with recent observational constraints on the luminosity function, average equivalent widths, and line ratios. However, they fail to reproduce the observed tail of high-EW sources and the number density of extreme line emitters. Possibilities to explain these discrepancies include an additional source of ionising photons and/or greater stochasticity in star formation in the model or photometric scatter and/or bias in the observations. With JWST now rapidly building larger samples and a wider range of emission lines the answer to this remaining discrepancy should be available imminently., Comment: 15 pages, accepted for publication in MNRAS, minor changes from original version

We present results on the morphological and structural evolution of a total of 4265 galaxies observed with JWST at $1.5 < z < 8$ in the JWST CEERS observations that overlap with the CANDELS EGS field. This is the biggest visually classified sample observed with JWST yet, $\sim20$ times larger than previous studies, and allows us to examine in detail how galaxy structure has changed over this critical epoch. All sources were classified by six individual classifiers using a simple classification scheme aimed to produce disk/spheroid/peculiar classifications, whereby we determine how the relative number of these morphologies evolves since the Universe's first billion years. Additionally, we explore structural and quantitative morphology measurements using \textsc{Morfometryka}, and show that galaxies at $z > 3$ are not dominated by irregular and peculiar structures, either visually or quantitatively, as previously thought. We find a strong dominance of morphologically selected disk galaxies up to $z = 8$, a far higher redshift than previously thought possible. We also find that the stellar mass and star formation rate densities are dominated by disk galaxies up to $z \sim 6$, demonstrating that most stars in the universe were likely formed in a disk galaxy. We compare our results to theory to show that the fraction of types we find is predicted by cosmological simulations, and that the Hubble Sequence was already in place as early as one billion years after the Big Bang. Additionally, we make our visual classifications public for the community., Comment: 18 pages, 13 figures. Submitted to ApJ; Comments welcome, Classification data will be publicly available at https://github.com/astroferreira/CEERS_JUNE_VISUAL_CLASSIFICATIONS when paper is in press

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

We analyse the photometric and spectroscopic properties of four galaxies in the epoch of reionisation (EoR) within the SMACS 0723 JWST Early Release Observations field. Given the known spectroscopic redshifts of these sources, we investigated the accuracy with which photometric redshifts can be derived using NIRCam photometry alone, finding that F115W imaging is essential to distinguish between z~8 galaxies with high equivalent width (EW) [O III] {\lambda}5007 emission and z~10 Balmer break galaxies. We find that all four sources exhibit strong (\geq 0.6 mag) F356W-F444W colours, which sit at the extreme end of theoretical predictions from numerical simulations. We find that these galaxies deviate (by roughly 0.5 dex) from the local correlation between [O III] {\lambda}5007/H\beta and [Ne III] {\lambda}3869/[O II], which is consistent with the predictions from simulations of high-redshift galaxies having elevated line excitation ratios. We measure the [O III] {\lambda}5007 rest-frame equivalent widths both directly from the spectroscopy, and indirectly as inferred from the strong F356W-F444W colours, finding large [O III] {\lambda}5007 EWs of 225-1740 {\AA}. The [O III] {\lambda}5007 and H\beta EWs are consistent with those seen in extreme, intensely star-forming dwarf galaxies in the local Universe. Our structural analysis indicates that these galaxies are resolved, exhibiting irregular shapes with bright clumps. In line with the predictions from the FLARES hydrodynamic simulations, such intense star formation and extreme nebular conditions are likely the norm, rather than the exception, in the EoR., Comment: 17 pages, 9 figures. Now accepted and published in MNRAS

With its exquisite sensitivity, wavelength coverage, and spatial and spectral resolution, the James Webb Space Telescope is poised to revolutionise our view of the distant, high-redshift ($z>5$) Universe. While Webb's spectroscopic observations will be transformative for the field, photometric observations play a key role in identifying distant objects and providing more comprehensive samples than accessible to spectroscopy alone. In addition to identifying objects, photometric observations can also be used to infer physical properties and thus be used to constrain galaxy formation models. However, inferred physical properties from broadband photometric observations, particularly in the absence of spectroscopic redshifts, often have large uncertainties. With the development of new tools for forward modelling simulations it is now routinely possible to predict observational quantities, enabling a direct comparison with observations. With this in mind, in this work, we make predictions for the colour evolution of galaxies at $z=5-15$ using the FLARES: First Light And Reionisation Epoch Simulations cosmological hydrodynamical simulation suite. We predict a complex evolution, driven predominantly by strong nebular line emission passing through individual bands. These predictions are in good agreement with existing constraints from Hubble and Spitzer as well as some of the first results from Webb. We also contrast our predictions with other models in the literature: while the general trends are similar we find key differences, particularly in the strength of features associated with strong nebular line emission. This suggests photometric observations alone should provide useful discriminating power between different models., Comment: 10 pages, 9 figures, accepted for publication in MNRAS

We present early results regarding the morphological and structural properties of galaxies seen with the James Webb Space Telescope at $z > 3$ in the Early Release Observations of SMACS 0723, a galaxy cluster at $z=0.39$. We investigate, for the first time, the optical morphologies of a significant number of $z > 3$ galaxies with accurate photometric redshifts in this field to determine the form of galaxy structure in the relatively early universe. We use visual morphologies and \textsc{Morfometryka} measures to perform quantitative morphology measurements, both parametric with light profile fitting (S\'ersic indices) and non-parametric (CAS values). Using these, we measure the relative fraction of disk, spheroidal, and peculiar galaxies at $3 < z < 8$. We discover the surprising result that at $z > 1.5$ disk galaxies dominate the overall fraction of morphologies, with a factor of $\sim 10$ relative higher number of disk galaxies than seen by the Hubble Space Telescope at these redshifts. Our visual morphological estimates of galaxies align closely with their locations in CAS parameter space and their S\'ersic indices., Comment: 11 pages, 7 figures. Accepted for publication in ApJL

The Azure Window was a natural arch situated in the west coast of Gozo (Maltese Archipelago) that collapsed in March 2017. We employ a Diver Propulsion Vehicle-mounted camera system to capture data for the 3D-modelling of this collapsed arch via photogrammetry. We demonstrate use of this method to document complex underwater geomorphology spread across a large area, and draw up a geomorphic assessment of the site and collapse event on the basis of this 3D-model. The methodology enables a reconstruction and understanding of the collapse event. On account of the high-resolution attained, we are able to cross-match the principal submerged components with their pre-collapse location; this enables an understanding of the dynamics of the collapse, confirmation of rock break-up along existing joints, and mapping the distribution of the rock and debris from the collapse event. We conclude that the key stages in the collapse of the Azure Window entailed erosion at the base of the pillar, leading to the latter's collapse in the southwest direction, breaking into two main sections that separated along the lithological boundary. We also find clear evidence that separation of some sections of the pillar followed pre-existing joints. The bridge collapsed vertically upon loss of support from the pillar, breaking into two main components and many other fragments. We also document further changes at the site post-collapse. We show that this approach can be utilised to understand and characterise such events even when significant time has elapsed since collapse, and rocks have already undergone erosion and significant marine growth., Comment: 33 pages, 14 figures, accepted for publication in Geomorphology

The James Webb Space Telescope (JWST) is set to transform many areas of astronomy, one of the most exciting is the expansion of the redshift frontier to $z>10$. In its first year alone JWST should discover hundreds of galaxies, dwarfing the handful currently known. To prepare for these powerful observational constraints, we use the First Light And Reionisation Epoch (FLARES) simulations to predict the physical and observational properties of the $z>10$ population of galaxies accessible to JWST. This is the first time such predictions have been made using a hydrodynamical model validated at low redshift. Our predictions at $z=10$ are broadly in agreement with current observational constraints on the far-UV luminosity function and UV continuum slope $\beta$, though the observational uncertainties are large. We note tension with recent constraints $z\sim 13$ from Harikane et al. 2022 - compared to these constraints, FLARES predicts objects with the same space density should have an order of magnitude lower luminosity, though this is mitigated slightly if dust attenuation is negligible in these systems. Our predictions suggest that in JWST's first cycle alone, around $600$ galaxies should be identified at $z>10$, with the first small samples available at $z>13$., Comment: Submitted to MNRAS 20/04/2022

The presence of stars on retrograde orbits in disc galaxies is usually attributed to accretion events, both via direct accretion, as well as through the heating of the disc stars. Recent studies have shown that retrograde orbits can also be produced via scattering by dense clumps, which are often present in the early stages of a galaxy's evolution. However, so far it has been unclear whether other internally-driven mechanisms, such as bars, are also capable of driving retrograde motion. Therefore, in this paper, we investigate the efficiencies with which bars and clumps produce retrograde orbits in disc galaxies. We do this by comparing the retrograde fractions and the spatial distributions of the retrograde populations in four $N$-body$+$smooth particle hydrodynamics (SPH) simulations of isolated disc galaxies spanning a range of evolutionary behaviours. We find that both bars and clumps are capable of generating significant retrograde populations of order $\sim 10\%$ of all stars. We also find that while clump-driven retrograde stars may be found at large galactocentric radii, bar-driven retrograde stars remain in the vicinity of the bar, even if the bar dissolves. Consequently, we find that retrograde stars in the Solar Neighbourhood in the clumpy models are exclusively clump-driven, but this is a trace population, constituting $0.01-0.04\%$ of the total stellar population in this region. Finally, we find that neither bars (including dissolving ones) nor clumps in the models are able to produce rotationally supported counter-rotating discs.

Over the years, the Maltese Islands have seen a marked rise in the prevalence of artificial lighting at night. The most evident type of light pollution arising from this evolution in anthropogenic night-time lighting is artificial skyglow via partial back-scattering in the atmosphere, leading to an increase in the Night Sky Brightness (NSB). The importance of understanding and quantifying the geographical distribution of the NSB is underscored by the adverse impact of light pollution on various spheres, from astronomical observation to ecology and human health. For the first time, we present a detailed map of the NSB over the Maltese archipelago carried out with Unihedron Sky Quality Meters. We show that the vast majority of the area of the Maltese Islands is heavily light polluted, with 87% of the area registering a NSB $<$ 20.39~mag$_{\rm SQM}$/arcsec$^2$ (Bortle Class 5 or higher) and 37.3% $<$ 19.09~mag$_{\rm SQM}$/arcsec$^2$ (Bortle Class 6 or higher), with the Milky Way being visible for only 12.8% of the area (adopting a visibility threshold $>$ 20.4 - 21.29~mag$_{\rm SQM}$/arcsec$^2$; Bortle Class 4). Coastal Dark Sky Heritage Areas on the island of Gozo retain generally darker skies than the rest of the islands, but light pollution originating further inland is encroaching upon and adversely affecting these sites. The methodology presented in this study can be adopted for continued future studies in Malta as well as for other regions., Comment: 29 pages, 6 figures, 11 supplementary figures, accepted for publication in Journal of Environmental Management

Compact ellipticals (cEs) are outliers from the scaling relations of early-type galaxies, particularly the mass-metallicity relation which is an important outcome of feedback. The formation of such low-mass, but metal-rich and compact, objects is a long-standing puzzle. Using a pair of high-resolution N-body+gas simulations, we investigate the evolution of a gas-rich low-mass galaxy on a highly radial orbit around a massive host galaxy. As the infalling low-mass galaxy passes through the host's corona at supersonic speeds, its diffuse gas outskirts are stripped by ram pressure, as expected. However, the compactness increases rapidly because of bursty star formation in the gas tidally driven to the centre. The metal-rich gas produced by supernovae and stellar winds is confined by the ram pressure from the surrounding environment, leading to subsequent generations of stars being more metal-rich. After the gas is depleted, tidal interactions enhance the metallicity further via the stripping of weakly bound, old, and metal-poor stars, while the size of the satellite is changed only modestly. The outcome is a metal-rich cE that is consistent with observations. These results argue that classical cEs are neither the stripped remnants of much more massive galaxies nor the merger remnants of normal dwarfs. We present observable predictions that can be used to test our model., Comment: 14 pages, 11 figures; accepted by ApJ. Minor change

We present rest-frame Ly$\alpha$ equivalent widths (EW) of 417 Ly$\alpha$ emitters (LAEs) detected with Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT) at $2.9 < z < 6.6$ in the Hubble Ultra Deep Field. Based on the deep MUSE spectroscopy and ancillary Hubble Space Telescope (HST) photometry data, we carefully measured EW values taking into account extended Ly$\alpha$ emission and UV continuum slopes ($\beta$). Our LAEs reach unprecedented depths, both in Ly$\alpha$ luminosities and UV absolute magnitudes, from log($L_{\rm Ly\alpha}$/erg s$^{-1}$) $\sim$41.0 to 43.0 and from Muv $\sim$ -16 to -21 (0.01-1.0 $L^{*}_{\rm z=3}$). The EW values span the range of $\sim$ 5 to 240 \AA\ or larger, and their distribution can be well fitted by an exponential law $N = N_{\rm 0}$ exp($-$EW/$w_{\rm 0}$). Owing to the high dynamic range in Muv, we find that the scale factor, $w_{\rm 0}$, depends on Muv in the sense that including fainter Muv objects increases $w_{\rm 0}$, i.e., the Ando effect. The results indicate that selection functions affect the EW scale factor. Taking these effects into account, we find that our $w_{\rm 0}$ values are consistent with those in the literature within $1\sigma$ uncertainties at $2.9 < z < 6.6$ at a given threshold of Muv and $L_{\rm Ly\alpha}$. Interestingly, we find 12 objects with EW $>200$ \AA\ above $1\sigma$ uncertainties. Two of these 12 LAEs show signatures of merger or AGN activity: the weak CIV $\lambda 1549$ emission line. For the remaining 10 very large EW LAEs, we find that the EW values can be reproduced by young stellar ages ($< 100$ Myr) and low metallicities ($\lesssim 0.02$ $Z_{\rm \odot}$). Otherwise, at least part of the Ly$\alpha$ emission in these LAEs needs to arise from anisotropic radiative transfer effects, fluorescence by hidden AGN or quasi-stellar object activity, or gravitational cooling., Comment: 22 pages, 12 figures, 9 tables, accepted for publication in A&A (MUSE UDF Series Paper X)

Recent theoretical models suggest that the early phase of galaxy formation could involve an epoch when galaxies are gas-rich but inefficient at forming stars: a "dark galaxy" phase. Here, we report the results of our MUSE (Multi Unit Spectroscopic Explorer) survey for dark galaxies fluorescently illuminated by quasars at $z>3$. Compared to previous studies which are based on deep narrow-band (NB) imaging, our integral field survey provides a nearly uniform sensitivity coverage over a large volume in redshift space around the quasars as well as full spectral information at each location. Thanks to these unique features, we are able to build control samples at large redshift distances from the quasars using the same data taken under the same conditions. By comparing the rest-frame equivalent width (EW$_{0}$) distributions of the Ly$\alpha$ sources detected in proximity to the quasars and in control samples, we detect a clear correlation between the locations of high EW$_{0}$ objects and the quasars. This correlation is not seen in other properties such as Ly$\alpha$ luminosities or volume overdensities, suggesting the possible fluorescent nature of at least some of these objects. Among these, we find 6 sources without continuum counterparts and EW$_{0}$ limits larger than $240\,\mathrm{\AA}$ that are the best candidates for dark galaxies in our survey at $z>3.5$. The volume densities and properties, including inferred gas masses and star formation efficiencies, of these dark galaxy candidates are similar to previously detected candidates at $z\approx2.4$ in NB surveys. Moreover, if the most distant of these are fluorescently illuminated by the quasar, our results also provide a lower limit of $t=60$ Myr on the quasar lifetime., Comment: 28 pages, 13 figures, 5 tables, submitted to ApJ

We present a measurement of the fraction of Lyman $\alpha$ (Ly$\alpha$) emitters ($X_{\rm{Ly} \alpha}$) amongst HST continuum-selected galaxies at $3

Direct Ly $\alpha$ imaging of intergalactic gas at $z\sim2$ has recently revealed giant cosmological structures around quasars, e.g. the Slug Nebula (Cantalupo et al. 2014). Despite their high luminosity, the detection rate of such systems in narrow-band and spectroscopic surveys is less than 10%, possibly encoding crucial information on the distribution of gas around quasars and the quasar emission properties. In this study, we use the MUSE integral-field instrument to perform a blind survey for giant Ly $\alpha$ nebulae around 17 bright radio-quiet quasars at $3

Measurements of the UV-continuum slopes provide valuable information on the physical properties of galaxies forming in the early universe, probing the dust reddening, age, metal content, and even the escape fraction. While constraints on these slopes generally become more challenging at higher redshifts as the UV continuum shifts out of the Hubble Space Telescope bands (particularly at z>7), such a characterisation actually becomes abruptly easier for galaxies in the redshift window z=9.5-10.5 due to the Spitzer/IRAC 3.6um-band probing the rest-UV continuum and the long wavelength baseline between this Spitzer band and the Hubble H-band. Higher S/N constraints on the UV slope are possible at z~10 than at z=8. Here we take advantage of this opportunity and five recently discovered bright z=9.5-10.5 galaxies to present the first measurements of the mean slope for a multi-object sample of galaxy candidates at z~10. We find the measured observed slopes of these candidates are $-2.1\pm0.3\pm0.2$ (random and systematic), only slightly bluer than the measured slopes at 3.5

We present a new integral-field spectroscopic dataset of the central part of the Orion Nebula (M 42), observed with the MUSE instrument at the ESO VLT. We reduced the data with the public MUSE pipeline. The output products are two FITS cubes with a spatial size of ~5.9'x4.9' (corresponding to ~0.76 pc x 0.63 pc) and a contiguous wavelength coverage of 4595...9366 Angstrom, spatially sampled at 0.2". We provide two versions with a sampling of 1.25 Angstrom and 0.85 Angstrom in dispersion direction. Together with variance cubes these files have a size of 75 and 110 GiB on disk. They represent one of the largest integral field mosaics to date in terms of information content. We make them available for use in the community. To validate this dataset, we compare world coordinates, reconstructed magnitudes, velocities, and absolute and relative emission line fluxes to the literature and find excellent agreement. We derive a two-dimensional map of extinction and present de-reddened flux maps of several individual emission lines and of diagnostic line ratios. We estimate physical properties of the Orion Nebula, using the emission line ratios [N II] and [S III] (for the electron temperature $T_e$) and [S II] and [Cl III] (for the electron density $N_e$), and show two-dimensional images of the velocity measured from several bright emission lines., Comment: Resubmitted to A&A after incorporating referee comments; access to full dataset via http://muse-vlt.eu/science/data-releases

Following our previous spectroscopic observations of $z>7$ galaxies with Gemini/GNIRS and VLT/XSHOOTER, which targeted a total of 8 objects, we present here our results from a deeper and larger VLT/FORS2 spectroscopic sample of Wide Field Camera 3 selected $z>7$ candidate galaxies. With our FORS2 setup we cover the 737-1070nm wavelength range, enabling a search for Lyman-$\alpha$ in the redshift range spanning 5.06 - 7.80. We target 22 $z$-band dropouts and find no evidence of Lyman-$\alpha$ emission, with the exception of a tentative detection ($<5\sigma$, which is our adopted criterion for a secure detection) for one object. The upper limits on Lyman-$\alpha$ flux and the broad-band magnitudes are used to constrain the rest-frame Equivalent Widths for this line emission. We analyse our FORS2 observations in combination with our previous GNIRS and XSHOOTER observations, and suggest that a simple model where the fraction of high rest-frame Equivalent Width emitters follows the trend seen at $z=3-6.5$ is inconsistent with our non-detections at $z\sim7.8$ at the 96\% confidence level. This may indicate that a significant neutral HI fraction in the intergalactic medium suppresses Lyman-$\alpha$, with an estimated neutral fraction $\chi_{HI}\sim0.5$, in agreement with other estimates., Comment: accepted for publication in MNRAS

By combining optical and near-IR observations from the Hubble Space Telescope with NIR photometry from the Spitzer Space Telescope it is possible to measure the rest-frame UV-optical colours of galaxies at z=4-8. The UV-optical spectral energy distribution of star formation dominated galaxies is the result of several different factors. These include the joint distribution of stellar masses, ages, and metallicities, and the subsequent reprocessing by dust and gas in the ISM. Using a large cosmological hydrodynamical simulation we investigate the predicted spectral energy distributions of galaxies at high-redshift with a particular emphasis on assessing the potential contribution of nebular emission. We find that the average pure stellar UV-optical colour correlates with both luminosity and redshift such that galaxies at lower-redshift and higher-luminosity are typically redder. Assuming the escape fraction of ionising photons is close to zero, the effect of nebular emission is to redden the UV-optical 1500-V_w colour by, on average, 0.4 mag at z=8 declining to 0.25 mag at z=4. Young and low-metallicity stellar populations, which typically have bluer pure stellar UV-optical colours, produce larger ionising luminosities and are thus more strongly affected by the reddening effects of nebular emission. This causes the distribution of 1500-V_w colours to narrow and the trends with luminosity and redshift to weaken. The strong effect of nebular emission leaves observed-frame colours critically sensitive to the source redshift. For example, increasing the redshift by 0.1 can result in observed frame colours changing by up to ~0.6. These predictions reinforce the need to include nebular emission when modelling the spectral energy distributions of galaxies at high-redshift and also highlight the difficultly in interpreting the observed colours of individual galaxies without precise redshifts., Comment: 13 pages, 17 figures, accepted in MNRAS

We present spectroscopic observations with VLT/XSHOOTER and Subaru/MOIRCS of a relatively bright Y-band drop-out galaxy in the Hubble Ultra Deep Field, first selected by Bunker et al. (2010), McLure et al. (2010) and Bouwens et al. (2010) to be a likely z~8-9 galaxy on the basis of its colours in the HST ACS and WFC3 images. This galaxy, HUDF.YD3 (also known as UDFy-38135539) has been targetted for VLT/SINFONI integral field spectroscopy by Lehnert et al. (2010), who published a candidate Lyman-alpha emission line at z=8.55 from this source. In our independent spectroscopy using two different infrared spectrographs (5 hours with VLT/XSHOOTER and 11 hours with Subaru/MOIRCS) we are unable to reproduce this line. We do not detect any emission line at the spectral and spatial location reported in Lehnert et al. (2010), despite the expected signal in our combined MOIRCS & XSHOOTER data being 5-sigma. The line emission also seems to be ruled out by the faintness of this object in recently extremely deep F105W (Y-band) HST/WFC3 imaging from HUDF12; the line would fall within this filter and such a galaxy should have been detected at Y(AB)=28.6mag (~20 sigma) rather than the marginal Y(AB)~30mag observed in the Y-band image, >3 times fainter than would be expected if the emission lie was real. Hence it appears highly unlikely that the reported Lyman-alpha line emission at z>8 is real, meaning that the highest-redshift sources for which Lyman-alpha emission has been seen are at z=6.9-7.2. It is conceivable that Lyman-alpha does not escape galaxies at higher redshifts, where the Gunn-Peterson absorption renders the Universe optically thick to this line. However, deeper spectroscopy on a larger sample of candidate z>7 galaxies will be needed to test this., Comment: Accepted for publication in MNRAS

The recent Hubble Space Telescope near-infrared imaging with the Wide-Field Camera #3 (WFC3) of the GOODS-South field in the CANDELS program covering nearly 100arcmin^2, along with already existing Advanced Camera for Surveys optical data, makes possible the search for bright galaxy candidates at redshift z ~ 7 - 9 using the Lyman-break technique. We present the first analysis of z'-drop z ~ 7 candidate galaxies in this area, finding 19 objects. We also analyse Y-drops at z ~ 8, trebling the number of bright (H_AB < 27 mag) Y-drops from our previous work, and compare our results with those of other groups based on the same data. The bright high redshift galaxy candidates we find serve to better constrain the bright end of the luminosity function at those redshift, and may also be more amenable to spectroscopic confirmation than the fainter ones presented in various previous work on the smaller fields (the Hubble Ultra Deep Field and the WFC 3 Early Release Science observations). We also look at the agreement with previous luminosity functions derived from WFC3 drop-out counts, finding a generally good agreement, except for the luminosity function of Yan et al. (2010) at z ~ 8, which is strongly ruled out., Comment: Accepted for publication in MNRAS. Originally submitted on 2012 June 18. arXiv admin note: text overlap with arXiv:1006.3545

We present Gemini/GNIRS spectroscopic observations of 4 z-band (z~7) dropout galaxies and VLT/XSHOOTER observations of one z-band dropout and 3 Y-band (z~8-9) dropout galaxies in the Hubble Ultra Deep Field, which were selected with Wide Field Camera 3 imaging on the Hubble Space Telescope. We find no evidence of Lyman-alpha emission with a typical 5-sigma sensitivity of 5X10^-18erg/cm^2/s, and we use the upper limits on Lyman-alpha flux and the broad-band magnitudes to constrain the rest-frame equivalent widths for this line emission. Accounting for incomplete spectral coverage, we survey 3.0 z-band dropouts and 2.9 Y-band dropouts to a Lyman-alpha rest-frame equivalent width limit > 120Ang (for an unresolved emission line); for an equivalent width limit of 50Ang the effective numbers of drop-outs surveyed fall to 1.2 z-band drop-outs and 1.5 Y-band drop-outs. A simple model where the fraction of high rest-frame equivalent width emitters follows the trend seen at z=3-6.5 is inconsistent with our non-detections at z=7-9 at the ~ 1-sigma level for spectrally unresolved lines, which may indicate that a significant neutral HI fraction in the intergalactic medium suppresses the Lyman-alpha line in z-drop and Y-drop galaxies at z > 7., Comment: Accepted for publication in MNRAS

The acquisition of deep Near-IR imaging with Wide Field Camera 3 on the Hubble Space Telescope has provided the opportunity to study the very-high redshift Universe. For galaxies up to $z\approx 7.7$ sufficient wavelength coverage exists to probe the rest-frame ultraviolet (UV) continuum without contamination from either Lyman-$\alpha$ emission or the Lyman-$\alpha$ break. In this work we use Near-IR imaging to measure the rest-frame UV continuum colours of galaxies at $4.7

We present a search for galaxies at 7.6

The addition of Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST) has led to a dramatic increase in our ability to study the z>6 Universe. The increase in the near-infrared (NIR) sensitivity of WFC3 over previous instruments has enabled us to reach apparent magnitudes approaching 29 (AB). This allows us to probe the rest-frame ultraviolet (UV) continuum, redshifted into the NIR at $z>6$. Taking advantage of the large optical depths at this redshift, resulting in the Lyman-alpha break, we use a combination of WFC3 imaging and pre-existing Advanced Camera for Surveys (ACS) imaging to search for z approx 7 over 4 fields. Our analysis reveals 29 new z approx 7 star forming galaxy candidates in addition to 16 pre-existing candidates already discovered in these fields. The improved statistics from our doubling of the robust sample of z-drop candidates confirms the previously observed evolution of the bright end of the luminosity function., Comment: 15 pages, accepted in MNRAS