A Direct Measurement of the IGM Opacity to HI Ionizing Photons

We present a new method to directly measure the opacity from HI Lyman limit (LL) absorption k_LL along quasar sightlines by the intergalactic medium (IGM). The approach analyzes the average (``stacked'') spectrum of an ensemble of quasars at a common redshift to infer the mean free path (MFP) to ionizing radiation. We apply this technique to 1800 quasars at z=3.50-4.34 drawn from the Sloan Digital Sky Survey (SDSS), giving the most precise measurements on k_LL at any redshift. From z=3.6 to 4.3, the opacity increases steadily as expected and is well parameterized by MFP = (48.4 +/- 2.1) - (38.0 +/- 5.3)*(z-3.6) h^-1 Mpc (proper distance). The relatively high MFP values indicate that the incidence of systems which dominate k_LL evolves less strongly at z>3 than that of the Lya forest. We infer a mean free path three times higher than some previous estimates, a result which has important implications for the photo-ionization rate derived from the emissivity of star forming galaxies and quasars. Finally, our analysis reveals a previously unreported, systematic bias in the SDSS quasar sample related to the survey's color targeting criteria. This bias potentially affects all z~3 IGM studies using the SDSS database.Comment: 7 pages, 4 figures; Accepted to ApJ

A Vacuum Spectrograph for the Infra-Red

A description is given of a 1-meter focus vacuum grating spectrometer for the infra-red which is distinguished by its versatility and compactness of design. The vacuum container is tubular in form and all optical parts are mounted on a carriage which is readily retractable from the case for adjustment. A novel feature is the use of a sine screw drive which gives a linear wave-length scale. A representative spectrum tracing is given

The origin of ultra diffuse galaxies: stellar feedback and quenching

We test if the cosmological zoom-in simulations of isolated galaxies from the FIRE project reproduce the properties of ultra diffuse galaxies. We show that stellar feedback-generated outflows that dynamically heat galactic stars, together with a passively aging stellar population after imposed quenching (from e.g. infall into a galaxy cluster), naturally reproduce the observed population of red UDGs, without the need for high spin halos or dynamical influence from their host cluster. We reproduce the range of surface brightness, radius and absolute magnitude of the observed z=0 red UDGs by quenching simulated galaxies at a range of different times. They represent a mostly uniform population of dark matter-dominated galaxies with M_star ~1e8 Msun, low metallicity and a broad range of ages. The most massive simulated UDGs require earliest quenching and are therefore the oldest. Our simulations provide a good match to the central enclosed masses and the velocity dispersions of the observed UDGs (20-50 km/s). The enclosed masses of the simulated UDGs remain largely fixed across a broad range of quenching times because the central regions of their dark matter halos complete their growth early. A typical UDG forms in a dwarf halo mass range of Mh~4e10-1e11 Msun. The most massive red UDG in our sample requires quenching at z~3 when its halo reached Mh ~ 1e11 Msun. If it, instead, continues growing in the field, by z=0 its halo mass reaches > 5e11 Msun, comparable to the halo of an L* galaxy. If our simulated dwarfs are not quenched, they evolve into bluer low-surface brightness galaxies with mass-to-light ratios similar to observed field dwarfs. While our simulation sample covers a limited range of formation histories and halo masses, we predict that UDG is a common, and perhaps even dominant, galaxy type around Ms~1e8 Msun, both in the field and in clusters.Comment: 20 pages, 13 figures; match the MNRAS accepted versio

Evolution of giant molecular clouds across cosmic time

Giant molecular clouds (GMCs) are well studied in the local Universe, however, exactly how their properties vary during galaxy evolution is poorly understood due to challenging resolution requirements, both observational and computational. We present the first time-dependent analysis of GMCs in a Milky Way-like galaxy and an Large Magellanic Cloud (LMC)-like dwarf galaxy of the FIRE-2 (Feedback In Realistic Environments) simulation suite, which have sufficient resolution to predict the bulk properties of GMCs in cosmological galaxy formation self-consistently. We show explicitly that the majority of star formation outside the galactic centre occurs within self-gravitating gas structures that have properties consistent with observed bound GMCs. We find that the typical cloud bulk properties such as mass and surface density do not vary more than a factor of 2 in any systematic way after the first Gyr of cosmic evolution within a given galaxy from its progenitor. While the median properties are constant, the tails of the distributions can briefly undergo drastic changes, which can produce very massive and dense self-gravitating gas clouds. Once the galaxy forms, we identify only two systematic trends in bulk properties over cosmic time: a steady increase in metallicity produced by previous stellar populations and a weak decrease in bulk cloud temperatures. With the exception of metallicity, we find no significant differences in cloud properties between the Milky Way-like and dwarf galaxies. These results have important implications for cosmological star and star cluster formation and put especially strong constraints on theories relating the stellar initial mass function to cloud properties

Live Fast, Die Young: GMC lifetimes in the FIRE cosmological simulations of Milky Way-mass galaxies

We present the first measurement of the lifetimes of giant molecular clouds (GMCs) in cosmological simulations at z = 0, using the Latte suite of FIRE-2 simulations of Milky Way (MW) mass galaxies. We track GMCs with total gas mass ≳10⁵ M⊙ at high spatial (∼1 pc), mass (7100 M⊙), and temporal (1 Myr) resolution. Our simulated GMCs are consistent with the distribution of masses for massive GMCs in the MW and nearby galaxies. We find GMC lifetimes of 5–7 Myr, or 1–2 freefall times, on average, with less than 2 per cent of clouds living longer than 20 Myr. We find decreasing GMC lifetimes with increasing virial parameter, and weakly increasing GMC lifetimes with galactocentric radius, implying that environment affects the evolutionary cycle of GMCs. However, our GMC lifetimes show no systematic dependence on GMC mass or amount of star formation. These results are broadly consistent with inferences from the literature and provide an initial investigation into ultimately understanding the physical processes that govern GMC lifetimes in a cosmological setting

The Great Observatories Origins Deep Survey: Constraints on the Lyman Continuum Escape Fraction Distribution of Lyman--Break Galaxies at 3.4<z<4.5

We use ultra-deep ultraviolet VLT/VIMOS intermediate-band and VLT/FORS1 narrow-band imaging in the GOODS Southern field to derive limits on the distribution of the escape fraction (f_esc) of ionizing radiation for L >~ L*(z=3) Lyman Break Galaxies (LBGs) at redshift 3.4--4.5. Only one LBG, at redshift z=3.795, is detected in its Lyman continuum (LyC; S/N~5.5), the highest redshift galaxy currently known with a direct detection. Its ultraviolet morphology is quite compact (R_eff=0.8, kpc physical). Three out of seven AGN are also detected in their LyC, including one at redshift z=3.951 and z850 = 26.1. From stacked data (LBGs) we set an upper limit to the average f_esc in the range 5%--20%, depending on the how the data are selected (e.g., by magnitude and/or redshift). We undertake extensive Monte Carlo simulations that take into account intergalactic attenuation, stellar population synthesis models, dust extinction and photometric noise in order to explore the moments of the distribution of the escaping radiation. Various distributions (exponential, log-normal and Gaussian) are explored. We find that the median f_esc is lower than ~6% with an 84% percentile limit not larger than 20%. If this result remains valid for fainter LBGs down to current observational limits, then the LBG population might be not sufficient to account for the entire photoionization budget at the redshifts considered here, with the exact details dependent upon the assumed ionizing background and QSO contribution thereto. It is possible that f_esc depends on the UV luminosity of the galaxies, with fainter galaxies having higher f_esc, and estimates of f_esc from a sample of faint LBG from the HUDF (i775<28.5) are in broad quantitative agreement with such a scenario.Comment: 58 pages, 23 figures; submitted to ApJ, revised version in response to referee's comment

A Definitive Survey for Lyman Limit Systems at z~3.5 with the Sloan Digital Sky Survey

We perform a semi-automated survey for tau>=2 Lyman Limit systems (LLSs) in quasar spectra from the Sloan Digital Sky Survey, Data Release 7. From a starting sample of 2473 quasars with zem=3.6-4.4, we analyze 469 spectra meeting strict seletion criteria for a total redshift path Dz=93.8 and identify 192 intervening systems at z>3.3. The incidence of tau>=2 LLSs per unit redshift, l(z), is well described by a single-power law at these redshifts: l(z) = C_LLS [(1+z)/(1+z_*)]^gamma, with z_*=3.7, C_LLS = 1.9+/-0.2, and gamma = 5.2+/-1.5 (68% c.l.). These values are systematically lower than previous estimates (especially at z<4) but are consistent with recent measurements of the mean free path to ionizing radiation. Extrapolations of this power-law to z=0 are inconsistent with previous estimations of l(z) at z<1 and suggest a break at z~2, similar to that observed for the Lya forest. Our results also indicate that the systems giving rise to LLS absorption decrease by ~50% in comoving number density and/or physical size from z=4 to 3.3, perhaps due to an enhanced extragalactic ultraviolet background. The observations place an integral constraint on the HI frequency distribution f(N_HI,X) and indicate that the power-law slope beta= dln[f(N,X)]/dln[N] is likely shallower than beta = -1 at N_HI=10^18 cm^-2. Including other constraints on f(N_HI,X) from the literature, we infer that beta is steeper than beta = -1.7 at N_HI~10^15 cm^-2, implying at least two inflections in f(N_HI,X). We also perform a survey for proximate LLSs (PLLSs) and find that l(z)_PLLS is systematically lower ~25% than intervening systems. Finally, we estimate that systematic effects impose an uncertainty of 10-20% in the l(z) measurements; these effects may limit the precision of all future surveys.Comment: 26 pages, 17 figures (most in color). Submitted to Ap

Reproducing the CO-to-H₂ conversion factor in cosmological simulations of Milky-Way-mass galaxies

We present models of CO(1–0) emission from Milky-Way-mass galaxies at redshift zero in the FIRE-2 cosmological zoom-in simulations. We calculate the molecular abundances by post-processing the simulations with an equilibrium chemistry solver while accounting for the effects of local sources, and determine the emergent CO(1–0) emission using a line radiative transfer code. We find that the results depend strongly on the shielding length assumed, which, in our models, sets the attenuation of the incident UV radiation field. At the resolution of these simulations, commonly used choices for the shielding length, such as the Jeans length, result in CO abundances that are too high at a given H₂ abundance. We find that a model with a distribution of shielding lengths, which has a median shielding length of ∼3 pc in cold gas (T < 300 K) for both CO and H₂, is able to reproduce both the observed CO(1–0) luminosity and inferred CO-to-H₂ conversion factor at a given star formation rate compared with observations. We suggest that this short shielding length can be thought of as a subgrid model, which controls the amount of radiation that penetrates giant molecular clouds