Source author record

Brant E. Robertson

Brant E. Robertson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

29works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

29 published item(s)

preprint2026arXiv

On the origins of oxygen: ALMA and JWST characterise the multi-phase, metal-enriched, star-bursting medium within a 'normal' $z > 11$ galaxy

The unexpectedly high abundance of galaxies at $z > 11$ revealed by JWST has sparked a debate on the nature of early galaxies and the physical mechanisms regulating their formation. The Atacama Large Millimeter/submillimeter Array (ALMA) has begun to provide vital insights on their gas and dust content, but so far only for extreme 'blue monsters'. Here we present new, deep ALMA observations of JADES-GS-z11-0, a more typical (sub-$L^*$) $z > 11$ galaxy that bridges the discovery space of JWST and the Hubble Space Telescope. These data confirm the presence of the [O III] 88 $μ$m line at $4.5σ$ significance, precisely at the redshift of several faint emission lines previously seen with JWST/NIRSpec, while the underlying dust continuum remains undetected ($F_ν< 9.0 \, \mathrm{μJy}$), implying an obscured star formation rate (SFR) of $\text{SFR}_\text{IR} \lesssim 6 \, \mathrm{M_\odot \, yr^{-1}}$ and dust mass of $M_\text{dust} \lesssim 1.0 \times 10^{6} \, \mathrm{M_\odot}$ (all $3σ$). The accurate ALMA redshift of $z_\text{[O III]} = 11.1221 \pm 0.0006$ ($\gtrsim \! 5\times$ refined over NIRSpec) helps confirm that redshifts measured purely from the Lyman-$α$ break, even spectroscopically, should properly take into account the effects of potential damped Lyman-$α$ absorption (DLA) systems to avoid systematic overestimates of up to $Δz \approx 0.5$. The [O III] 88 $μ$m luminosity of $L_\text{[O III]} = (1.1 \pm 0.3) \times 10^{8} \, \mathrm{L_\odot}$, meanwhile, agrees well with the scaling relation for local metal-poor dwarfs given the SFR measured by NIRCam, NIRSpec, and MIRI. The spatially resolved MIRI and ALMA emission also underscores that JADES-GS-z11-0 is likely to consist of two low-mass components that are undergoing strong bursts of star formation yet are already pre-enriched in oxygen (~20-30% solar), only 400 Myr after the Big Bang.

preprint2025arXiv

COSMOS2025: The COSMOS-Web galaxy catalog of photometry, morphology, redshifts, and physical parameters from JWST, HST, and ground-based imaging

We present COSMOS2025, the COSMOS-Web catalog of photometry, morphology, photometric redshifts and physical parameters for more than 700,000 galaxies in the Cosmic Evolution Survey (COSMOS) field. This catalog is based on our \textit{James Webb Space Telescope} 255\,h COSMOS-Web program, which provides deep near-infrared imaging in four NIRCam (F115W, F150W, F277W, F444W) and one MIRI (F770W) filter over the central $\sim 0.54 {\, \rm deg}^2$ ($\sim 0.2 {\, \rm deg}^2$ for MIRI) in COSMOS. These data are combined with ground- and space-based data to derive photometric measurements of NIRCam-detected sources using both fixed-aperture photometry (on the space-based bands) and a profile-fitting technique on all 37 bands spanning 0.3-8 micron. We provide morphology for all sources from complementary techniques including profile fitting and machine-learning classification. We derive photometric redshifts, physical parameters and non-parametric star formation histories from spectral energy distribution (SED) fitting. The catalog has been extensively validated against previous COSMOS catalogs and other surveys. Photometric redshift accuracy measured using spectroscopically confirmed galaxies out to $z\sim9$ reaches $σ_{\rm MAD} = 0.012$ at $m_{\rm F444W}<28$ and remains at $σ_{\rm MAD} \lesssim 0.03$ as a function of magnitude, color, and galaxy type. This represents a factor of $\sim 2$ improvement at 26 AB mag compared to COSMOS2020. The catalog is approximately 80\% complete at $\log(M_{\star}/{\rm M}_{\odot}) \sim 9$ at $z \sim 10$ and at $\log(M_{\star}/{\rm M}_{\odot}) \sim 7$ at $z \sim 0.2$, representing a gain of 1\,dex compared to COSMOS2020. COSMOS2025 represents the definitive COSMOS-Web catalog. It is provided with complete documentation, together with redshift probability distributions, and it is ready for scientific exploitation today.

preprint2024arXiv

Inside the bubble: exploring the environments of reionisation-era Lyman-$α$ emitting galaxies with JADES and FRESCO

We present a study of the environments of 17 Lyman-$α$ (Ly$α$) emitting galaxies (LAEs) in the reionisation era ($5.8 < z < 8$) identified by JWST/NIRSpec as part of the JWST Advanced Deep Extragalactic Survey (JADES). Unless situated in sufficiently (re)ionised regions, Ly$α$ emission from these galaxies would be strongly absorbed by neutral gas in the intergalactic medium (IGM). We conservatively estimate sizes of the ionised regions required to reconcile the relatively low Ly$α$ velocity offsets ($Δv_\text{Ly$α$}<300\,\mathrm{km\,s^{-1}}$) with moderately high Ly$α$ escape fractions ($f_\mathrm{esc,\,Lyα}>5\%$) observed in our sample of LAEs, suggesting the presence of ionised hydrogen along the line of sight towards at least eight out of 17 LAEs. We find minimum physical `bubble' sizes of the order of $R_\text{ion}\sim0.1$-$1\,\mathrm{pMpc}$ are required in a patchy reionisation scenario where ionised bubbles containing the LAEs are embedded in a fully neutral IGM. Around half of the LAEs in our sample are found to coincide with large-scale galaxy overdensities seen in FRESCO at $z \sim 5.8$-$5.9$ and $z\sim7.3$, suggesting Ly$α$ transmission is strongly enhanced in such overdense regions, and underlining the importance of LAEs as tracers of the first large-scale ionised bubbles. Considering only spectroscopically confirmed galaxies, we find our sample of UV-faint LAEs ($M_\text{UV}\gtrsim-20\,\mathrm{mag}$) and their direct neighbours are generally not able to produce the required ionised regions based on the Ly$α$ transmission properties, suggesting lower-luminosity sources likely play an important role in carving out these bubbles. These observations demonstrate the combined power of JWST multi-object and slitless spectroscopy in acquiring a unique view of the early Universe during cosmic reionisation via the most distant LAEs.

preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Observational Facilities to Study Dark Matter

We present an overview of future observational facilities that will significantly enhance our understanding of the fundamental nature of dark matter. These facilities span a range of observational techniques including optical/near-infrared imaging and spectroscopy, measurements of the cosmic microwave background, pulsar timing, 21-cm observations of neutral hydrogen at high redshift, and the measurement of gravitational waves. Such facilities are a critical component of a multi-pronged experimental program to uncover the nature of dark matter, while often providing complementary measurements of dark energy, neutrino physics, and inflation.

preprint2021arXiv

Deep Realistic Extragalactic Model (DREaM) Galaxy Catalogs: Predictions for a Roman Ultra-Deep Field

In the next decade, deep galaxy surveys from telescopes such as the James Webb Space Telescope and Roman Space Telescope will provide transformational data sets that will greatly enhance the understanding of galaxy formation during the epoch of reionization (EoR). In this work, we present the Deep Realistic Extragalactic Model (DREaM) for creating synthetic galaxy catalogs. Our model combines dark matter simulations, subhalo abundance matching and empirical models, and includes galaxy positions, morphologies, and spectral energy distributions (SEDs). The resulting synthetic catalog extends to redshifts $z \sim 12 $, and galaxy masses $\log_{10}(M/M_{\odot}) = 5 $ covering an area of $1\, {\rm deg}^2$ on the sky. We use DREaM to explore the science returns of a $1\, {\rm deg}^2$ \emph{Roman} UDF, and to provide a resource for optimizing ultra-deep survey designs. We find that a \emph{Roman} UDF to $\sim 30\,m_{\rm AB}$ will potentially detect more than $10^6$ $M_{\rm UV}<-17$ galaxies, with more than $10^4$ at redshifts $z>7$, offering an unparalleled dataset for constraining galaxy properties during the EoR. Our synthetic catalogs and simulated images are made publicly available to provide the community with a tool to prepare for upcoming data.

preprint2021arXiv

Discovering Supernovae at Epoch of Reionization with Nancy Grace Roman Space Telescope

Massive stars play critical roles for the reionization of the Universe. Individual massive stars at the reionization epoch (z > 6) are too faint to observe and quantify their contributions to reionization. Some massive stars, however, explode as superluminous supernovae (SLSNe) or pair-instability supernovae (PISNe) that are luminous enough to observe even at z > 6 and allow for the direct characterization of massive star properties at the reionization epoch. In addition, hypothetical long-sought-after PISNe are expected to be present preferentially at high redshifts, and their discovery will have a tremendous impact on our understanding of massive star evolution and the formation of stellar mass black holes. The near-infrared Wide Field Instrument on Nancy Grace Roman Space Telescope will excel at discovering such rare high-redshift supernovae. In this work, we investigate the best survey strategy to discover and identify SLSNe and PISNe at z > 6 with Roman. We show that the combination of the F158 and F213 filters can clearly separate both SLSNe and PISNe at z > 6 from nearby supernovae through their colors and magnitudes. The limiting magnitudes are required to be 27.0 mag and 26.5 mag in the F158 and F213 filters, respectively, to identify supernovae at z > 6. If we conduct a 10 deg2 transient survey with these limiting magnitudes for 5 years with a cadence of one year, we expect to discover 22.5 +- 2.8 PISNe and 3.1 +- 0.3 SLSNe at z > 6, depending on the cosmic star-formation history. The same survey is estimated to discover 76.1 +- 8.2 PISNe and 9.1 +- 0.9 SLSNe at 5 < z < 6. Such a supernova survey requires the total observational time of approximately 525 hours in 5 years. The legacy data acquired with the survey will also be beneficial for many different science cases including the study of high-redshift galaxies.

preprint2020arXiv

The Physical Nature of Starburst-Driven Galactic Outflows

We present the fourth of the Cholla Galactic OutfLow Simulations suite (CGOLS). Using a physically-motivated prescription for clustered supernova feedback, we successfully drive a multiphase outflow from a disk galaxy. The high resolution ($< 5\,\mathrm{pc}$) across a relatively large domain ($20\,\mathrm{kpc}$) allows us to capture the hydrodynamic mixing and dynamical interactions between the hot and cool ($T \sim 10^4\,\mathrm{K}$) phases in the outflow, which in turn leads to direct evidence of a qualitatively new mechanism for cool gas acceleration in galactic winds. We show that mixing of momentum from the hot phase to the cool phase accelerates the cool gas to $800\,\mathrm{km}\,\mathrm{s}^{-1}$ on kpc scales, with properties inconsistent with the physical models of ram pressure acceleration or with bulk cooling from the hot phase. The mixing process also affects the hot phase, modifying its radial profiles of temperature, density, and velocity from the expectations of radial supersonic flow. This mechanism provides a physical explanation for the high velocity, blue shifted, low ionization absorption lines often observed in the spectra of starburst and high redshift galaxies.

preprint2020arXiv

The role of galaxies and AGN in reionising the IGM -- III : IGM-galaxy cross-correlations at z~6 from 8 quasar fields with DEIMOS and MUSE

We present improved results of the measurement of the correlation between galaxies and the intergalactic medium (IGM) transmission at the end of reionisation. We have gathered a sample of $13$ spectroscopically confirmed Lyman-break galaxies (LBGs) and $21$ Lyman-$α$ emitters (LAEs) at angular separations $20'' \lesssim θ\lesssim 10'$ ($\sim 0.1-4$ pMpc at $z\sim 6$) from the sightlines to $8$ background $z\gtrsim 6$ quasars. We report for the first time the detection of an excess of Lyman-$α$ transmission spikes at $\sim 10-60$ cMpc from LAEs ($3.2σ$) and LBGs ($1.9σ$). We interpret the data with an improved model of the galaxy-Lyman-$α$ transmission and two-point cross-correlations which includes the enhanced photoionisation due to clustered faint sources, enhanced gas densities around the central bright objects and spatial variations of the mean free path. The observed LAE(LBG)-Lyman-$α$ transmission spike two-point cross-correlation function (2PCCF) constrains the luminosity-averaged escape fraction of all galaxies contributing to reionisation to $\langle f_{\rm esc} \rangle_{M_{\rm UV}<-12} = 0.14_{-0.05}^{+0.28}\,(0.23_{-0.12}^{+0.46})$. We investigate if the 2PCCF measurement can determine whether bright or faint galaxies are the dominant contributors to reionisation. Our results show that a contribution from faint galaxies ($M_{\rm UV} > -20 \, (2σ)$) is necessary to reproduce the observed 2PCCF and that reionisation might be driven by different sub-populations around LBGs and LAEs at $z\sim 6$.

preprint2019arXiv

Clustering with JWST: Constraining Galaxy Host Halo Masses, Satellite Quenching Efficiencies, and Merger Rates at z=4-10

Galaxy clustering measurements can be used to constrain many aspects of galaxy evolution, including galaxy host halo masses, satellite quenching efficiencies, and merger rates. We simulate JWST galaxy clustering measurements at z$\sim$4$-$10 by utilizing mock galaxy samples produced by an empirical model, the UniverseMachine. We also adopt the survey footprints and typical depths of the planned joint NIRCam and NIRSpec Guaranteed Time Observation program planned for Cycle 1 to generate realistic JWST survey realizations and to model high-redshift galaxy selection completeness. We find that galaxy clustering will be measured with $\gtrsim$5$σ$ significance at z$\sim$4$-$10. Halo mass precisions resulting from Cycle 1 angular clustering measurements will be $\sim$0.2 dex for faint (-18 $\gtrsim$ $\mathit{M}_{\mathrm{UV}}^{ }$ $\gtrsim$ -19) galaxies at z$\sim$4$-$10 as well as $\sim$0.3 dex for bright ($\mathit{M}_{\mathrm{UV}}^{ }$ $\sim$ -20) galaxies at z$\sim$4$-$7. Dedicated spectroscopic follow-up over $\sim$150 arcmin$^2$ would improve these precisions by $\sim$0.1 dex by removing chance projections and low-redshift contaminants. Future JWST observations will therefore provide the first constraints on the stellar-halo mass relation in the epoch of reionization and substantially clarify how this relation evolves at z$>$4. We also find that $\sim$1000 individual satellites will be identifiable at z$\sim$4$-$8 with JWST, enabling strong tests of satellite quenching evolution beyond currently available data (z$\lesssim$2). Finally, we find that JWST observations can measure the evolution of galaxy major merger pair fractions at z$\sim$4$-$8 with $\sim$0.1$-$0.2 dex uncertainties. Such measurements would help determine the relative role of mergers to the build-up of stellar mass into the epoch of reionization.

preprint2019arXiv

The Lyman Continuum Escape Survey -- II: Ionizing Radiation as a Function of the [OIII]/[OII] Line Ratio

We discuss the rest-frame optical emission line spectra of a large (~50) sample of z=3.1 Lyman alpha emitters (LAEs) whose physical properties suggest such sources are promising analogs of galaxies in the reionization era. Reliable Lyman continuum escape fractions have now been determined for a large sample of such LAEs from the Lyman Continuum Escape Survey (LACES) undertaken via deep HST imaging in the SSA22 survey area reported in Fletcher et al. (2019). Using new measures of [OII] emission secured from Keck MOSFIRE spectra we re-examine, for a larger sample, earlier claims that Lyman continuum leakages may correlate with the nebular emission line ratio [OIII]/[OII] as expected for density-bound HII regions. We find that a large [OIII]/[OII] line ratio is indeed a necessary condition for Lyman continuum leakage, strengthening earlier claims made using smaller samples at various redshifts. However, not all LAEs with large [OIII]/[OII] line ratios are leakers and leaking radiation appears not to be associated with differences in other spectral diagnostics. This suggests the detection of leaking radiation is modulated by an additional property, possibly the viewing angle for porous HII regions. We discuss our new results in the context of the striking bimodality of LAE leakers and non-leakers found in the LACES program and the implications for the sources of cosmic reionization.

preprint2015arXiv

Cholla : A New Massively-Parallel Hydrodynamics Code For Astrophysical Simulation

We present Cholla (Computational Hydrodynamics On ParaLLel Architectures), a new three-dimensional hydrodynamics code that harnesses the power of graphics processing units (GPUs) to accelerate astrophysical simulations. Cholla models the Euler equations on a static mesh using state-of-the-art techniques, including the unsplit Corner Transport Upwind (CTU) algorithm, a variety of exact and approximate Riemann solvers, and multiple spatial reconstruction techniques including the piecewise parabolic method (PPM). Using GPUs, Cholla evolves the fluid properties of thousands of cells simultaneously and can update over ten million cells per GPU-second while using an exact Riemann solver and PPM reconstruction. Owing to the massively-parallel architecture of GPUs and the design of the Cholla code, astrophysical simulations with physically interesting grid resolutions (> 256^3) can easily be computed on a single device. We use the Message Passing Interface library to extend calculations onto multiple devices and demonstrate nearly ideal scaling beyond 64 GPUs. A suite of test problems highlights the physical accuracy of our modeling and provides a useful comparison to other codes. We then use Cholla to simulate the interaction of a shock wave with a gas cloud in the interstellar medium, showing that the evolution of the cloud is highly dependent on its density structure. We reconcile the computed mixing time of a turbulent cloud with a realistic density distribution destroyed by a strong shock with the existing analytic theory for spherical cloud destruction by describing the system in terms of its median gas density.

preprint2015arXiv

Cosmic Reionization and Early Star-Forming Galaxies: A Joint Analysis of New Constraints from Planck and Hubble Space Telescope

We discuss new constraints on the epoch of cosmic reionization and test the assumption that most of the ionizing photons responsible arose from high redshift star-forming galaxies. Good progress has been made in charting the end of reionization through spectroscopic studies of z~6-8 QSOs, gamma-ray bursts and galaxies expected to host Lyman-alpha emission. However, the most stringent constraints on its duration have come from the integrated optical depth, tau, of Thomson scattering to the cosmic microwave background. Using the latest data on the abundance and luminosity distribution of distant galaxies from Hubble Space Telescope imaging, we simultaneously match the reduced value tau=0.066 +/- 0.012 recently reported by the Planck collaboration and the evolving neutrality of the intergalactic medium with a reionization history within 6 <~ z <~ 10, thereby reducing the requirement for a significant population of very high redshift (z>>10) galaxies. Our analysis strengthens the conclusion that star-forming galaxies dominated the reionization process and has important implications for upcoming 21cm experiments and searches for early galaxies with James Webb Space Telescope.

preprint2015arXiv

New redshift z ~ 9 galaxies in the Hubble Frontier Fields: Implications for early evolution of the UV luminosity density

We present the results of a new search for galaxies at redshift z ~ 9 in the first two Hubble Frontier Fields with completed HST WFC3/IR and ACS imaging. To ensure robust photometric redshift solutions, and to minimize incompleteness, we confine our search to objects with H_{160} < 28.6 (AB mag), consider only image regions with an rms noise sigma_{160} > 30 mag (within a 0.5-arcsec diameter aperture), and insist on detections in both H_{160} and J_{140}. The result is a survey covering an effective area (after accounting for magnification) of 10.9 sq. arcmin, which yields 12 galaxies at 8.4 < z < 9.5. Within the Abell-2744 cluster and parallel fields we confirm the three brightest objects reported by Ishigaki et al. (2014), but recover only one of the four z > 8.4 sources reported by Zheng et al. (2014). In the MACSJ0416.1-240 cluster field we report five objects, and explain why each of these eluded detection or classification as z ~ 9 galaxies in the published searches of the shallower CLASH data. Finally, we uncover four z ~ 9 galaxies from the previously unsearched MACSJ0416.1-240 parallel field. Based on the published magnification maps we find that only one of these 12 galaxies is likely boosted by more than a factor of two by gravitational lensing. Consequently we are able to perform a fairly straightforward reanalysis of the normalization of the z ~ 9 UV galaxy luminosity function as explored previously in the HUDF12 programme. We conclude that the new data strengthen the evidence for a continued smooth decline in UV luminosity density (and hence star-formation rate density) from z ~ 8 to z ~ 9, contrary to recent reports of a marked drop-off at these redshifts. This provides further support for the scenario in which early galaxy evolution is sufficiently extended to explain cosmic reionization.

preprint2014arXiv

Accounting for Cosmic Variance in Studies of Gravitationally-Lensed High-Redshift Galaxies in the Hubble Frontier Field Clusters

Strong gravitational lensing provides a powerful means for studying faint galaxies in the distant universe. By magnifying the apparent brightness of background sources, massive clusters enable the detection of galaxies fainter than the usual sensitivity limit for blank fields. However, this gain in effective sensitivity comes at the cost of a reduced survey volume and, in this {\it Letter}, we demonstrate there is an associated increase in the cosmic variance uncertainty. As an example, we show that the cosmic variance uncertainty of the high redshift population viewed through the Hubble Space Telescope Frontier Field cluster Abell 2744 increases from ~35% at redshift z~7 to >~65% at z~10. Previous studies of high redshift galaxies identified in the Frontier Fields have underestimated the cosmic variance uncertainty that will affect the ultimate constraints on both the faint end slope of the high-redshift luminosity function and the cosmic star formation rate density, key goals of the Frontier Field program.

preprint2013arXiv

New Constraints on Cosmic Reionization from the 2012 Hubble Ultra Deep Field Campaign

Understanding cosmic reionization requires the identification and characterization of early sources of hydrogen-ionizing photons. The 2012 Hubble Ultra Deep Field (UDF12) campaign has acquired the deepest infrared images with the Wide Field Camera 3 aboard Hubble Space Telescope and, for the first time, systematically explored the galaxy population deep into the era when cosmic microwave background (CMB) data indicates reionization was underway. The UDF12 campaign thus provides the best constraints to date on the abundance, luminosity distribution, and spectral properties of early star-forming galaxies. We synthesize the new UDF12 results with the most recent constraints from CMB observations to infer redshift-dependent ultraviolet (UV) luminosity densities, reionization histories, and electron scattering optical depth evolution consistent with the available data. Under reasonable assumptions about the escape fraction of hydrogen ionizing photons and the intergalactic medium clumping factor, we find that to fully reionize the universe by redshift z~6 the population of star-forming galaxies at redshifts z~7-9 likely must extend in luminosity below the UDF12 limits to absolute UV magnitudes of M_UV\sim -13 or fainter. Moreover, low levels of star formation extending to redshifts z~15-25, as suggested by the normal UV colors of z\simeq7-8 galaxies and the smooth decline in abundance with redshift observed by UDF12 to z\simeq10, are additionally likely required to reproduce the optical depth to electron scattering inferred from CMB observations.

preprint2013arXiv

The CASSOWARY spectroscopy survey: A new sample of gravitationally lensed galaxies in SDSS

Bright gravitationally lensed galaxies provide our most detailed view of galaxies at high redshift. Yet as a result of the small number of ultra-bright z~2 lensed systems with confirmed redshifts, most detailed spectroscopic studies have been limited in their scope. With the goal of increasing the number of bright lensed galaxies available for detailed follow-up, we have undertaken a spectroscopic campaign targeting wide separation (>3 arcsec) galaxy-galaxy lens candidates within the Sloan Digital Sky Survey (SDSS). Building on the earlier efforts of our CASSOWARY survey, we target a large sample of candidate galaxy-galaxy lens systems in SDSS using a well-established search algorithm which identifies blue arc-like structures situated around luminous red galaxies. In this paper, we present a new redshift catalog containing 25 lensed sources in SDSS confirmed through spectroscopic follow-up of candidate galaxy-galaxy lens systems. Included in this new sample are two of the brightest galaxies (r=19.6 and 19.7) galaxies known at z~2, a low metallicity (12 + log (O/H)~8.0) extreme nebular line emitting galaxy at z=1.43, and numerous systems for which detailed follow-up will be possible. The source redshifts span 0.9<z<2.5 (median redshift of 1.9), and their optical magnitudes are in the range 19.6<r<22.3. We present a brief source-by-source discussion of the spectroscopic properties extracted from our confirmatory spectra and discuss some initial science results. With more than 50 gravitationally lensed z>1 galaxies now confirmed within SDSS, it will soon be possible for the first time to develop generalized conclusions from detailed spectroscopic studies of the brightest lensed systems at high redshift.

preprint2012arXiv

Evolution of the Sizes of Galaxies over 7<z<12 Revealed by the 2012 Hubble Ultra Deep Field Campaign

We analyze the redshift- and luminosity-dependent sizes of dropout galaxy candidates in the redshift range z~7-12 using deep images from the UDF12 campaign, data which offers two distinct advantages over that used in earlier work. Firstly, we utilize the increased S/N ratio offered by the UDF12 imaging to provide improved size measurements for known galaxies at z=6.5-8 in the HUDF. Specifically, we stack the new deep F140W image with the existing F125W data in order to provide improved measurements of the half-light radii of z-dropouts. Similarly we stack this image with the new deep UDF12 F160W image to obtain new size measurements for a sample of Y-dropouts. Secondly, because the UDF12 data have allowed the construction of the first robust galaxy sample in the HUDF at z>8, we have been able to extend the measurement of average galaxy size out to significantly higher redshifts. Restricting our size measurements to sources which are now detected at >15sigma, we confirm earlier indications that the average half-light radii of z~7-12 galaxies are extremely small, 0.3-0.4 kpc, comparable to the sizes of giant molecular associations in local star-forming galaxies. We also confirm that there is a clear trend of decreasing half-light radius with increasing redshift, and provide the first evidence that this trend continues beyond z~8. Modeling the evolution of the average half-light radius as a power-law (1+z)^s, we obtain a best-fit index of s=-1.28+/-0.13 over the redshift range z~4-12, mid-way between the physically expected evolution for baryons embedded in dark halos of constant mass (s=-1) and constant velocity (s=-1.5). A clear size-luminosity relation, such as that found at lower redshift, is also evident in both our z- and Y-dropout sample. This relation can be interpreted in terms of a constant surface density of star formation over a range in luminosity of 0.05-1.0L*_z=3.(abridged)

preprint2012arXiv

The UV Luminosity Function of star-forming galaxies via dropout selection at redshifts z ~ 7 and 8 from the 2012 Ultra Deep Field campaign

We present a catalog of high redshift star-forming galaxies selected to lie within the redshift range z ~ 7-8 using the Ultra Deep Field 2012 (UDF12), the deepest near-infrared (near-IR) exposures yet taken with the Hubble Space Telescope. As a result of the increased near-infrared exposure time compared to previous HST imaging in this field, we probe 0.65 (0.25) mag fainter in absolute UV magnitude, at z ~ 7 (8), which increases confidence in a measurement of the faint end slope of the galaxy luminosity function. Through a 0.7 mag deeper limit in the key F105W filter that encompasses or lies just longward of the Lyman break, we also achieve a much-refined color-color selection that balances high redshift completeness and a low expected contamination fraction. We improve the number of drop-out selected UDF sources to 47 at z ~ 7 and 27 at z ~ 8. Incorporating brighter archival and ground-based samples, we measure the z ~ 7 UV luminosity function to an absolute magnitude limit of M_UV = -17 and find a faint end Schechter slope of α= -1.87+/- 0.18. Using a similar color-color selection at z ~ 8 that takes account of our newly-added imaging in the F140W filter, and incorporating archival data from the HIPPIES and BoRG campaigns, we provide a robust estimate of the faint end slope at z ~ 8, α= -1.94 +/- 0.23. We briefly discuss our results in the context of earlier work and that derived using the same UDF12 data but with an independent photometric redshift technique (McLure et al 2012).

preprint2011arXiv

Connecting the Gamma Ray Burst Rate and the Cosmic Star Formation History: Implications for Reionization and Galaxy Evolution

(Abridged) The contemporary discoveries of galaxies and gamma ray bursts (GRBs) at high redshift have supplied the first direct information on star formation when the universe was only a few hundred million years old. The probable origin of long duration GRBs in the deaths of massive stars would link the universal GRB rate to the redshift-dependent star formation rate density, although exactly how is currently unknown. As the most distant GRBs and star-forming galaxies probe the reionization epoch, the potential rewards of understanding the redshift-dependent ratio Psi(z) of the GRB rate to star formation rate are significant and include addressing fundamental questions such as incompleteness in rest-frame UV surveys for determining the star formation rate at high redshift and time variations in the stellar initial mass function. Using an extensive sample of 112 GRBs above a fixed luminosity limit drawn from the Second Swift Burst Alert Telescope catalog, we compare the cumulative redshift distribution N(< z) of GRBs with the star formation density rho_sfr(z) measured from UV-selected galaxies over 0 < z < 4. Strong evolution (e.g., Psi(z) \propto (1+z)^{1.5}) is disfavored, while more modest evolution (e.g., Psi(z) \propto (1+z)^{0.5}) is consistent with the data. If such trends continue beyond z ~ 4, we find the discovery rate of distant GRBs implies a star formation rate density much higher than that inferred from UV-selected galaxies. We show that such a star formation history would over-predict the observed stellar mass density at z > 4 measured from rest-frame optical surveys. The resolution of this important disagreement is currently unclear, and the GRB production rate at early times is likely more complex than a simple function of star formation rate and progenitor metallicity.

preprint2011arXiv

Pixel-z: Studying Substructure and Stellar Populations in Galaxies out to z~3 using Pixel Colors I. Systematics

We perform a pixel-by-pixel analysis of 467 galaxies in the GOODS-VIMOS survey to study systematic effects in extracting properties of stellar populations (age, dust, metallicity and SFR) from pixel colors using the pixel-z method. The systematics studied include the effect of the input stellar population synthesis model, passband limitations and differences between individual SED fits to pixels and global SED-fitting to a galaxy's colors. We find that with optical-only colors, the systematic errors due to differences among the models are well constrained. The largest impact on the age and SFR e-folding time estimates in the pixels arises from differences between the Maraston models and the Bruzual&Charlot models, when optical colors are used. This results in systematic differences larger than the 2σ uncertainties in over 10 percent of all pixels in the galaxy sample. The effect of restricting the available passbands is more severe. In 26 percent of pixels in the full sample, passband limitations result in systematic biases in the age estimates which are larger than the 2σ uncertainties. Systematic effects from model differences are reexamined using Near-IR colors for a subsample of 46 galaxies in the GOODS-NICMOS survey. For z > 1, the observed optical/NIR colors span the rest frame UV-optical SED, and the use of different models does not significantly bias the estimates of the stellar population parameters compared to using optical-only colors. We then illustrate how pixel-z can be applied robustly to make detailed studies of substructure in high redshift galaxies such as (a) radial gradients of age, SFR, sSFR and dust and (b) the distribution of these properties within subcomponents such as spiral arms and clumps. Finally, we show preliminary results from the CANDELS survey illustrating how the new HST/WFC3 data can be exploited to probe substructure in z~1-3 galaxies.

preprint2011arXiv

Probing Very Bright End of Galaxy Luminosity Function at z >~ 7 Using Hubble Space Telescope Pure Parallel Observations

We report the first results from the Hubble Infrared Pure Parallel Imaging Extragalactic Survey, which utilizes the pure parallel orbits of the Hubble Space Telescope to do deep imaging along a large number of random sightlines. To date, our analysis includes 26 widely separated fields observed by the Wide Field Camera 3, which amounts to 122.8 sq.arcmin in total area. We have found three bright Y098-dropouts, which are candidate galaxies at z >~ 7.4. One of these objects shows an indication of peculiar variability and its nature is uncertain. The other two objects are among the brightest candidate galaxies at these redshifts known to date L>2L*. Such very luminous objects could be the progenitors of the high-mass Lyman break galaxis (LBGs) observed at lower redshifts (up to z~5). While our sample is still limited in size, it is much less subject to the uncertainty caused by "cosmic variance" than other samples because it is derived using fields along many random sightlines. We find that the existence of the brightest candidate at z~7.4 is not well explained by the current luminosity function (LF) estimates at z~8. However, its inferred surface density could be explained by the prediction from the LFs at z~7 if it belongs to the high-redshift tail of the galaxy population at z~7.

preprint2010arXiv

A Method for Measuring the Bias of High-Redshift Galaxies from Cosmic Variance

As deeper observations discover increasingly distant galaxies, characterizing the properties of high-redshift galaxy populations will become increasingly challenging and paramount. We present a method for measuring the clustering bias of high-redshift galaxies from the field-to-field scatter in their number densities induced by cosmic variance. Multiple widely-separated fields are observed to provide a large number of statistically-independent samples of the high-redshift galaxy population. The expected Poisson uncertainty is removed from the measured dispersion in the distribution of galaxy number counts observed across these many fields, leaving, on average, only the contribution to the scatter expected from cosmic variance. With knowledge of the Lambda Cold Dark Matter power spectrum, the galaxy bias is then calculated from the measured cosmic variance. The results of cosmological N-body simulations can then be used to estimate the halo mass associated with the measured bias. We use Monte Carlo simulations to demonstrate that Hubble Space Telescope pure parallel programs will be able to determine galaxy bias at z>~6 using this method, complementing future measurements from correlation functions.

preprint2010arXiv

Early star-forming galaxies and the reionization of the Universe

Star forming galaxies represent a valuable tracer of cosmic history. Recent observational progress with Hubble Space Telescope has led to the discovery and study of the earliest-known galaxies corresponding to a period when the Universe was only ~800 million years old. Intense ultraviolet radiation from these early galaxies probably induced a major event in cosmic history: the reionization of intergalactic hydrogen. New techniques are being developed to understand the properties of these most distant galaxies and determine their influence on the evolution of the universe.

preprint2010arXiv

Estimating Luminosity Function Constraints from High-Redshift Galaxy Surveys

The installation of the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST) will revolutionize the study of high-redshift galaxy populations. Initial observations of the HST Ultra Deep Field (UDF) have yielded multiple z>~7 dropout candidates. Supplemented by the Great Observatory Origins Deep Survey (GOODS) Early Release Science (ERS) and further UDF pointings, these data will provide crucial information about the most distant known galaxies. However, achieving tight constraints on the z~7 galaxy luminosity function (LF) will require even more ambitious photometric surveys. Using a Fisher matrix approach to fully account for Poisson and cosmic sample variance, as well as covariances in the data, we estimate the uncertainties on LF parameters achieved by surveys of a given area and depth. Applying this method to WFC3 z~7 dropout galaxy samples, we forecast the LF parameter uncertainties for a variety of model surveys. We demonstrate that performing a wide area (~1 deg^2) survey to H_AB~27 depth or increasing the UDF depth to H_AB~30 provides excellent constraints on the high-z LF when combined with the existing UDF GO and GOODS ERS data. We also show that the shape of the matter power spectrum may limit the possible gain of splitting wide area (>~0.5 deg^2) high-redshift surveys into multiple fields to probe statistically independent regions; the increased root-mean-squared density fluctuations in smaller volumes mostly offset the improved variance gained from independent samples.

preprint2010arXiv

The Large Scale Bias of Dark Matter Halos: Numerical Calibration and Model Tests

We measure the clustering of dark matter halos in a large set of collisionless cosmological simulations of the flat LCDM cosmology. Halos are identified using the spherical overdensity algorithm, which finds the mass around isolated peaks in the density field such that the mean density is Delta times the background. We calibrate fitting functions for the large scale bias that are adaptable to any value of Delta we examine. We find a ~6% scatter about our best fit bias relation. Our fitting functions couple to the halo mass functions of Tinker et. al. (2008) such that bias of all dark matter is normalized to unity. We demonstrate that the bias of massive, rare halos is higher than that predicted in the modified ellipsoidal collapse model of Sheth, Mo, & Tormen (2001), and approaches the predictions of the spherical collapse model for the rarest halos. Halo bias results based on friends-of-friends halos identified with linking length 0.2 are systematically lower than for halos with the canonical Delta=200 overdensity by ~10%. In contrast to our previous results on the mass function, we find that the universal bias function evolves very weakly with redshift, if at all. We use our numerical results, both for the mass function and the bias relation, to test the peak-background split model for halo bias. We find that the peak-background split achieves a reasonable agreement with the numerical results, but ~20% residuals remain, both at high and low masses.

preprint2009arXiv

Computational Eulerian Hydrodynamics and Galilean Invariance

Eulerian hydrodynamical simulations are a powerful and popular tool for modeling fluids in astrophysical systems. In this work, we critically examine recent claims that these methods violate Galilean invariance of the Euler equations. We demonstrate that Eulerian hydrodynamics methods do converge to a Galilean-invariant solution, provided a well-defined convergent solution exists. Specifically, we show that numerical diffusion, resulting from diffusion-like terms in the discretized hydrodynamical equations solved by Eulerian methods, accounts for the effects previously identified as evidence for the Galilean non-invariance of these methods. These velocity-dependent diffusive terms lead to different results for different bulk velocities when the spatial resolution of the simulation is kept fixed, but their effect becomes negligible as the resolution of the simulation is increased to obtain a converged solution. In particular, we find that Kelvin-Helmholtz instabilities develop properly in realistic Eulerian calculations regardless of the bulk velocity provided the problem is simulated with sufficient resolution (a factor of 2-4 increase compared to the case without bulk flows for realistic velocities). Our results reiterate that high-resolution Eulerian methods can perform well and obtain a convergent solution, even in the presence of highly supersonic bulk flows.

preprint2009arXiv

LSST Science Book, Version 2.0

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

preprint2009arXiv

Recovering Stellar Population Properties and Redshifts from Broad-Band Photometry of Simulated Galaxies: Lessons for SED Modeling

We present a detailed analysis of our ability to determine stellar masses, ages, reddening and extinction values, and star formation rates of high-redshift galaxies by modeling broad-band SEDs with stellar population synthesis. In order to do so, we computed synthetic optical-to-NIR SEDs for model galaxies taken from hydrodynamical merger simulations placed at redshifts 1.5 < z < 3. Viewed under different angles and during different evolutionary phases, the simulations represent a wide variety of galaxy types (disks, mergers, spheroids). We show that simulated galaxies span a wide range in SEDs and color, comparable to these of observed galaxies. In all star-forming phases, dust attenuation has a large effect on colors, SEDs, and fluxes. The broad-band SEDs were then fed to a standard SED modeling procedure and resulting stellar population parameters were compared to their true values. Disk galaxies generally show a decent median correspondence between the true and estimated mass and age, but suffer from large uncertainties. During the merger itself, we find larger offsets (e.g., log M_recovered - log M_true = -0.13^{+0.10}_{-0.14}). E(B-V) values are generally recovered well, but the estimated total visual absorption Av is consistently too low, increasingly so for larger optical depths. Since the largest optical depths occur during the phases of most intense star formation, it is for the highest SFRs that we find the largest underestimates. The masses, ages, E(B-V), Av, and SFR of merger remnants (spheroids) are very well reproduced. We discuss possible biases in SED modeling results caused by mismatch between the true and template star formation history, dust distribution, metallicity variations and AGN contribution.

preprint2008arXiv

Simulating High-Redshift Disk Galaxies: Applications to Long Duration Gamma-Ray Burst Hosts

The efficiency of star formation governs many observable properties of the cosmological galaxy population, yet many current models of galaxy formation largely ignore the important physics of star formation and the interstellar medium (ISM). Using hydrodynamical simulations of disk galaxies that include a treatment of the molecular ISM and star formation in molecular clouds (Robertson & Kravtsov 2008), we study the influence of star formation efficiency and molecular hydrogen abundance on the properties of high-redshift galaxy populations. In this work, we focus on a model of low-mass, star forming galaxies at 1<~z<~2 that may host long duration gamma-ray bursts (GRBs). Observations of GRB hosts have revealed a population of faint systems with star formation properties that often differ from Lyman-break galaxies (LBGs) and more luminous high-redshift field galaxies. Observed GRB sightlines are deficient in molecular hydrogen, but it is unclear to what degree this deficiency owes to intrinsic properties of the galaxy or the impact the GRB has on its environment. We find that hydrodynamical simulations of low-stellar mass systems at high-redshifts can reproduce the observed star formation rates and efficiencies of GRB host galaxies at redshifts 1<~z<~2. We show that the compact structure of low-mass high-redshift GRB hosts may lead to a molecular ISM fraction of a few tenths, well above that observed in individual GRB sightlines. However, the star formation rates of observed GRB host galaxies imply molecular gas masses of 10^8 - 10^9 M_sun similar to those produced in the simulations, and may therefore imply fairly large average H_2 fractions in their ISM.