Showing 1–50 of 79 results for author: Finlator, K

A massive, neutral gas reservoir permeating a galaxy proto-cluster after the reionization era

Authors: Kasper E. Heintz, Jake S. Bennett, Pascal A. Oesch, Albert Sneppen, Douglas Rennehan, Joris Witstok, Renske Smit, Simone Vejlgaard, Chamilla Terp, Umran S. Koca, Gabriel B. Brammer, Kristian Finlator, Matthew J. Hayes, Debora Sijacki, Rohan P. Naidu, Jorryt Matthee, Francesco Valentino, Nial R. Tanvir, Páll Jakobsson, Peter Laursen, Darach J. Watson, Romeel Davé, Laura C. Keating, Alba Covelo-Paz

Abstract: Galaxy clusters are the most massive, gravitationally-bound structures in the Universe, emerging through hierarchical structure formation of large-scale dark matter and baryon overdensities. Early galaxy ``proto-clusters'' are believed to be important physical drivers of the overall cosmic star-formation rate density and serve as ``hotspots'' for the reionization of the intergalactic medium. Our u… ▽ More Galaxy clusters are the most massive, gravitationally-bound structures in the Universe, emerging through hierarchical structure formation of large-scale dark matter and baryon overdensities. Early galaxy ``proto-clusters'' are believed to be important physical drivers of the overall cosmic star-formation rate density and serve as ``hotspots'' for the reionization of the intergalactic medium. Our understanding of the formation of these structures at the earliest cosmic epochs is, however, limited to sparse observations of their galaxy members, or based on phenomenological models and cosmological simulations. Here we report the detection of a massive neutral, atomic hydrogen (HI) gas reservoir permeating a galaxy proto-cluster at redshift $z=5.4$, observed one billion years after the Big Bang. The presence of this cold gas is revealed by strong damped Lyman-$α$ absorption features observed in several background galaxy spectra taken with JWST/NIRSpec in close on-sky projection. While overall the sightlines probe a large range in HI column densities, $N_{\rm HI} = 10^{21.7}-10^{23.5}$ cm$^{-2}$, they are similar across nearby sightlines, demonstrating that they probe the same dense, neutral gas. This observation of a massive, large-scale overdensity of cold neutral gas challenges current large-scale cosmological simulations and has strong implications for the reionization topology of the Universe. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: Submitted

Galaxy-Absorber Association in the Epoch of Reionization: Galactic Population Luminosity Distribution for Different Absorbers at $10 \geq z \geq 5.5$

Authors: Samir Kušmić, Kristian Finlator, Ezra Huscher, Maya Steen

Abstract: How do galaxies of different luminosities contribute to the metal absorber populations of varying species and strength? We present our analysis of the predicted metal contributions from galaxies as observed in quasar absorption line spectra during the end of the Epoch of Reionization (EoR; $10 \geq z \geq 5.5$). This was done by implementing on-the-fly particle tracking into the latest \textsc{Tec… ▽ More How do galaxies of different luminosities contribute to the metal absorber populations of varying species and strength? We present our analysis of the predicted metal contributions from galaxies as observed in quasar absorption line spectra during the end of the Epoch of Reionization (EoR; $10 \geq z \geq 5.5$). This was done by implementing on-the-fly particle tracking into the latest \textsc{Technicolor Dawn} simulation and then linking CII, CIV, SiII, SiIV, OI, and MgII absorbers to host galaxies in post-processing. We define the Host Galaxy Luminosity Distribution (HGLD) as the rest-frame ultraviolet luminosity distribution of galaxies contributing ions to an absorber, weighted by the fractional contribution, and compute its dependence on ion and absorber strength. The HGLD shape is predicted to be indistinguishable from the field luminosity function, indicating that there is no relationship between the absorber strength or ion and the luminosity of the dominant contributing galaxy. Switching from galaxy luminosity to stellar mass, the predicted host galaxy mass distributions (HGMD) indicate that more-massive galaxies contribute a higher fraction of metal ions to absorbers of each species, with the HGMD of stronger absorbers extending out to higher masses. We conclude that the fraction of absorbing metal ions contributed by galaxies increases weakly with stellar mass, but the scatter in luminosity at fixed stellar mass obscures this relationship. For the same reason, we predict that observational analyses of the absorber-galaxy relationship will uncover stronger trends with stellar mass than with luminosity. △ Less

Submitted 30 April, 2024; originally announced May 2024.

Comments: 18 pages, 9 figures, 6 tables. Submitted to ApJ

The JWST-PRIMAL Legacy Survey. A JWST/NIRSpec reference sample for the physical properties and Lyman-$α$ absorption and emission of $\sim 500$ galaxies at $z=5.5-13.4$

Authors: K. E. Heintz, G. B. Brammer, D. Watson, P. A. Oesch, L. C. Keating, M. J. Hayes, Abdurro'uf, K. Z. Arellano-Córdova, A. C. Carnall, C. R. Christiansen, F. Cullen, R. Davé, P. Dayal, A. Ferrara, K. Finlator, J. P. U. Fynbo, S. R. Flury, V. Gelli, S. Gillman, R. Gottumukkala, K. Gould, T. R. Greve, S. E. Hardin, T. Y. -Y Hsiao, A. Hutter , et al. (23 additional authors not shown)

Abstract: One of the surprising early findings with JWST has been the discovery of a strong "roll-over" or a softening of the absorption edge of Ly$α$ in a large number of galaxies at ($z\gtrsim 6$), in addition to systematic offsets from photometric redshift estimates and fundamental galaxy scaling relations. This has been interpreted as damped Ly$α$ absorption (DLA) wings from high column densities of neu… ▽ More One of the surprising early findings with JWST has been the discovery of a strong "roll-over" or a softening of the absorption edge of Ly$α$ in a large number of galaxies at ($z\gtrsim 6$), in addition to systematic offsets from photometric redshift estimates and fundamental galaxy scaling relations. This has been interpreted as damped Ly$α$ absorption (DLA) wings from high column densities of neutral atomic hydrogen (HI), signifying major gas accretion events in the formation of these galaxies. To explore this new phenomenon systematically, we assemble the JWST/NIRSpec PRImordial gas Mass AssembLy (PRIMAL) legacy survey of 494 galaxies at $z=5.5-13.4$. We characterize this benchmark sample in full and spectroscopically derive the galaxy redshifts, metallicities, star-formation rates, and ultraviolet slopes. We define a new diagnostic, the Ly$α$ damping parameter $D_{\rm Lyα}$ to measure and quantify the Ly$α$ emission strength, HI fraction in the IGM, or local HI column density for each source. The JWST-PRIMAL survey is based on the spectroscopic DAWN JWST Archive (DJA-Spec). All the software, reduced spectra, and spectroscopically derived quantities and catalogs are made publicly available in dedicated repositories. The fraction of strong galaxy DLAs are found to be in the range $65-95\%$ at $z>5.5$. The fraction of strong Ly$α$ emitters (LAEs) is found to increase with decreasing redshift, in qualitative agreement with previous observational results, and are predominantly associated with low-metallicity and UV faint galaxies. By contrast, strong DLAs are observed in galaxies with a variety of intrinsic physical properties. Our results indicate that strong DLAs likely reflect a particular early assembly phase of reionization-era galaxies, at which point they are largely dominated by pristine HI gas accretion. [abridged] △ Less

Submitted 2 April, 2024; originally announced April 2024.

Comments: 20 pages, 18 figures, 3 tables. Submitted to A&A. Comments welcome! All data and catalogs are available through the DAWN JWST Archive (DJA): https://dawn-cph.github.io/dja/ and https://github.com/keheintz/jwst-primal

Impact of Helium II resonant absorption on the UVB modeled in three dimensions

Authors: Ezra Huscher, Kristian Finlator, Samir Kušmić, Maya Steen

Abstract: We implement a treatment of Helium {\sc ii} absorption and re-emission into the \TD cosmological simulations to study its impact on the metagalactic ultraviolet background (UVB) in three dimensions. By comparing simulations with and without He {\sc ii} reprocessing, we show that it weakens the mean UVB by $\sim$3 dex from $z = 10$ to $z = 5$ between 3.5 and 4 Ryd, where the He {\sc ii} Lyman-serie… ▽ More We implement a treatment of Helium {\sc ii} absorption and re-emission into the \TD cosmological simulations to study its impact on the metagalactic ultraviolet background (UVB) in three dimensions. By comparing simulations with and without He {\sc ii} reprocessing, we show that it weakens the mean UVB by $\sim$3 dex from $z = 10$ to $z = 5$ between 3.5 and 4 Ryd, where the He {\sc ii} Lyman-series resonance occurs. In overdense regions, the overall UVB amplitude is higher and the impact of \HeII reprocessing is weaker, qualitatively indicating an early start to \HeII reionization near galaxies. Comparing our simulations to two popular one-dimensional UV models, we find good agreement up to 3 Ryd at $z = 5$. At higher energies, our simulation shows significantly greater He {\sc ii} absorption because it accounts for He {\sc ii} arising both in diffuse regions and in Lyman limit systems. By contrast, the comparison models account only for He {\sc ii} in Lyman limit systems, which are subdominant prior to the completion of He {\sc ii} reionization. The H {\sc i} and He {\sc ii} reionization histories are nearly unaffected by He {\sc ii} reprocessing although the cosmic star formation rate density is altered by up to $4\%$. The cosmic mass density of C {\sc iv} is reduced by $\sim$2 dex when He {\sc ii} is accounted for while Si {\sc iv}, C {\sc ii}, Mg {\sc ii}, Si {\sc ii}, and O {\sc i} are unaffected. △ Less

Submitted 29 March, 2024; originally announced April 2024.

Comments: 10 pages, 8 figures

Searching for [CII] Emission from the First Sample of $z\sim 6$ OI Absorption-Associated Galaxies with ALMA

Authors: Yunjing Wu, Zheng Cai, Jianan Li, Kristian Finlator, Marcel Neeleman, J. Xavier Prochaska, Bjorn H. C. Emonts, Shiwu Zhang, Feige Wang, Jinyi Yang, Ran Wang, Xiaohui Fan, Dandan Xu, Emmet Golden-Marx, Laura C. Keating, Joseph F. Hennawi

Abstract: We report the first statistical analyses of [CII] and dust continuum observations in six strong OI absorber fields at the end of the reionization epoch obtained by the Atacama Large Millimeter/Submillimeter Array (ALMA). Combined with one [CII] emitter reported in Wu et al. (2021), we detect one OI-associated [CII] emitter in six fields. At redshifts of OI-absorbers in non-detection fields, no emi… ▽ More We report the first statistical analyses of [CII] and dust continuum observations in six strong OI absorber fields at the end of the reionization epoch obtained by the Atacama Large Millimeter/Submillimeter Array (ALMA). Combined with one [CII] emitter reported in Wu et al. (2021), we detect one OI-associated [CII] emitter in six fields. At redshifts of OI-absorbers in non-detection fields, no emitters are brighter than our detection limit within impact parameters of 50 kpc and velocity offsets between $\pm200\ {\rm km\ s^{-1}}$. The averaged [CII]-detection upper limit is $< 0.06$ Jy ${\rm km\ s^{-1}}$ (3$σ$), corresponding to the [CII] luminosity of $L_{\rm [CII]} <5.8\times 10^7\ L_{\odot}$ and the [CII]-based star formation rate of ${\rm SFR_{\rm [CII]}} < 5.5$ $M_\odot$ yr$^{-1}$. Cosmological simulations suggest that only $\sim10^{-2.5}$ [CII] emitters around [OI] absorbers have comparable SFR to our detection limit. Although the detection in one out of six fields is reported, an order of magnitude number excess of emitters obtained from our ALMA observations supports that the contribution of massive galaxies that caused the metal enrichment cannot be ignored. Further, we also found 14 tentative galaxy candidates with S/N of $\approx4.3$ at large impact parameters ($>50$ kpc) and having larger outflow velocities within $\pm 600$ km s$^{-1}$. If these detections are confirmed in the future, then the mechanism of pushing metals at larger distances with higher velocities needs to be further explored from the theoretical side. △ Less

Submitted 8 November, 2023; v1 submitted 5 October, 2023; originally announced October 2023.

Comments: Accepted for publication in ApJ. Main text 10 pages, 5 figures

A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE): JWST Discovers an Overdensity around a Metal Absorption-selected Galaxy at $z\sim5.5$

Authors: Yunjing Wu, Feige Wang, Zheng Cai, Xiaohui Fan, Kristian Finlator, Jinyi Yang, Joseph F. Hennawi, Fengwu Sun, Jaclyn B. Champagne, Xiaojing Lin, Zihao Li, Zuyi Chen, Eduardo Bañados, George D. Becker, Sarah E. I. Bosman, Gstavo Bruzual, Stephane Charlot, Hsiao-Wen Chen, Jacopo Chevallard, Anna-Christina Eilers, Emanuele Paolo Farina, Xiangyu Jin, Hyunsung D. Jun, Koki Kakiichi, Mingyu Li , et al. (5 additional authors not shown)

Abstract: The launch of ${\it JWST}$ opens a new window for studying the connection between metal-line absorbers and galaxies at the end of the Epoch of Reionization (EoR). Previous studies have detected absorber-galaxy pairs in limited quantities through ground-based observations. To enhance our understanding of the relationship between absorbers and their host galaxies at $z>5$, we utilized the NIRCam Wid… ▽ More The launch of ${\it JWST}$ opens a new window for studying the connection between metal-line absorbers and galaxies at the end of the Epoch of Reionization (EoR). Previous studies have detected absorber-galaxy pairs in limited quantities through ground-based observations. To enhance our understanding of the relationship between absorbers and their host galaxies at $z>5$, we utilized the NIRCam Wide Field Slitless Spectroscopy (WFSS) to search for absorber-associated galaxies by detecting their rest-frame optical emission lines (e.g., [OIII] + H$β$). We report the discovery of a MgII-associated galaxy at $z=5.428$ using data from the ${\it JWST}$ ASPIRE program. The MgII absorber is detected on the spectrum of quasar J0305--3150 with a rest-frame equivalent width of 0.74$\mathring{A}$. The associated galaxy has an [OIII] luminosity of $10^{42.5}\ {\rm erg\ s^{-1}}$ with an impact parameter of 24.9 proper kiloparsecs (pkpc). The joint ${\it HST}$-${\it JWST}$ spectral energy distribution (SED) implies a stellar mass and star-formation rate of ${\rm M_* \approx 10^{8.8}}$ ${\rm M_{\odot}}$, ${\rm SFR}\approx 10\ {\rm M_{\odot}\ yr^{-1}}$. Its [OIII] equivalent width and stellar mass are typical of [OIII] emitters at this redshift. Furthermore, connecting the outflow starting time to the SED-derived stellar age, the outflow velocity of this galaxy is $\sim300\ {\rm km\ s^{-1}}$, consistent with theoretical expectations. We identified six additional [OIII] emitters with impact parameters of up to $\sim300$ pkpc at similar redshifts ($|dv|<1000\ {\rm km\ s^{-1}}$). The observed number is consistent with that in cosmological simulations. This pilot study suggests that systematically investigating the absorber-galaxy connection within the ASPIRE program will provide insights into the metal-enrichment history in the early universe. △ Less

Submitted 8 November, 2023; v1 submitted 28 September, 2023; originally announced September 2023.

Comments: Accepted for publication in ApJL. Main text 8 pages, 4 figures. For more information of the JWST ASPIRE program please check https://aspire-quasar.github.io/index.html

CLEAR: Survey Overview, Data Analysis and Products

Authors: Raymond C. Simons, Casey Papovich, Ivelina G. Momcheva, Gabriel Brammer, Vicente Estrada-Carpenter, Steven L. Finkelstein, Catherine M. Gosmeyer, Jasleen Matharu, Jonathan R. Trump, Bren E. Backhaus, Yingjie Cheng, Nikko J. Cleri, Henry C. Ferguson, Kristian Finlator, Mauro Giavalisco, Zhiyuan Ji, Intae Jung, Jennifer M. Lotz, Rosalia O'Brien, Rosalind E. Skelton, Vithal Tilvi, Benjamin Weiner

Abstract: We present an overview of the CANDELS Lyman-a Emission At Reionization (CLEAR) survey. CLEAR is a 130 orbit program of the Hubble Space Telescope using the Wide Field Camera 3 (WFC3) IR G102 grism. CLEAR targets 12 pointings divided between the GOODS-N and GOODS-S fields of the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). Combined with existing spectroscopic data from other… ▽ More We present an overview of the CANDELS Lyman-a Emission At Reionization (CLEAR) survey. CLEAR is a 130 orbit program of the Hubble Space Telescope using the Wide Field Camera 3 (WFC3) IR G102 grism. CLEAR targets 12 pointings divided between the GOODS-N and GOODS-S fields of the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). Combined with existing spectroscopic data from other programs, the full CLEAR dataset includes spectroscopic imaging of these fields over 0.8-1.7 um. In this Paper, we describe the CLEAR survey, the survey strategy, the data acquisition, reduction, processing, and science products and catalogs released alongside this paper. The catalogs include emission line fluxes and redshifts derived from the combination of the photometry and grism spectroscopy for 6048 galaxies, primarily ranging from 0.2 < z < 3. We also provide an overview of CLEAR science goals and results. In conjunction with this Paper we provide links to electronic versions of the data products, including 1D + 2D extracted spectra and emission line maps. △ Less

Submitted 16 March, 2023; originally announced March 2023.

Comments: Accepted in ApJS. 29 pages, 14 figures. The data products (1D + 2D spectra, emission line maps) and catalogs (photometric and spectroscopic catalogs of GOODS-N and GOODS-S) are available at https://archive.stsci.edu/hlsp/clear

Brighter and More Massive Galaxies in the Vicinity of Lyman-alpha Nebulae

Authors: Natalie K. Wells, Moire K. M. Prescott, Kristian M. Finlator

Abstract: It has been well established in the local universe that galaxy properties differ based on the large-scale environment in which they reside. As luminous Lyman-alpha nebulae have been shown to trace overdense environments at z~2-3, comparing the properties of galaxies within Lyman-alpha nebulae systems to those in the field can provide insight into how and when locally-observed trends between galaxy… ▽ More It has been well established in the local universe that galaxy properties differ based on the large-scale environment in which they reside. As luminous Lyman-alpha nebulae have been shown to trace overdense environments at z~2-3, comparing the properties of galaxies within Lyman-alpha nebulae systems to those in the field can provide insight into how and when locally-observed trends between galaxy properties and environment emerged. Six Lyman-alpha nebulae were discovered at z~2.3 in a blind search of the GOODS-S extragalactic field, a region also covered by the 3D-HST spectroscopic survey. Utilizing 3D-HST data, we identified 86 galaxies in the vicinity of these nebulae and used statistical tests to compare their physical properties to galaxies elsewhere in the field. Galaxies lying within 320 proper kpc of a Lyman-alpha nebula are roughly half a magnitude brighter than those in the field, with higher stellar masses, higher star formation rates, and larger effective radii. Even when considering the effects of sample incompleteness, our study suggests that galaxies in overdensities at z~2.3 traced by Lyman-alpha nebulae are being influenced by their environment. Furthermore, Lyman-alpha nebula-associated galaxies lie on the same main sequence of star formation as field galaxies, but have a larger proportion of high-mass galaxies, consistent with the idea that galaxy evolution is accelerated in rich environments. Expanded surveys for Lyman-alpha nebulae in other deep extragalactic fields and galaxy spectroscopic follow-up with JWST will better constrain the demographics of Lyman-alpha nebula-associated galaxies. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Comments: 15 pages, 12 figures

Journal ref: Astrophys. J. 941 (2022) 180-190

The environments and hosts of metal absorption at z>5

Authors: Caitlin C. Doughty, Kristian M. Finlator

Abstract: A growing population of metal absorbers are observed at z>5, many showing strong evolution in incidence approaching the epoch of hydrogen reionization. Follow-up surveys examining fields around these metals have resulted in galaxy detections but the direct physical relationship between the detected galaxies and absorbers is unclear. Upcoming observations will illuminate this galaxy-absorber relati… ▽ More A growing population of metal absorbers are observed at z>5, many showing strong evolution in incidence approaching the epoch of hydrogen reionization. Follow-up surveys examining fields around these metals have resulted in galaxy detections but the direct physical relationship between the detected galaxies and absorbers is unclear. Upcoming observations will illuminate this galaxy-absorber relationship, but the theoretical framework for interpreting these observations is lacking. To inform future z>5 studies, we define the expected relationship between metals and galaxies using the Technicolor Dawn simulation to model metal absorption from z=5-7, encompassing the end of reionization. We find that metal absorber types and strengths are slightly better associated with their environment than with the traits of their host galaxies, as absorption system strengths are more strongly correlated with the local galaxy overdensity than the stellar mass of their host galaxy. For redshifts prior to the end of the epoch of reionization, strong high ionization transitions like C IV are more spatially correlated with brighter galaxies on scales of a few hundred proper kpc than are low ionization systems, due to the former's preference for environments with higher UVB amplitudes and those ions' relative rarity at z>6. Post-reionization, the galaxy counts near these high-ionization ions are reduced, and increase surrounding certain low-ionization ions due to a combination of their relative abundances and preferred denser gas phase. We conclude that galaxy-absorber relationships are expected to evolve rapidly such that high-ionization absorbers are better tracers of galaxies pre-reionization while low-ionization absorbers are better post-reionization. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: Accepted to MNRAS

ALMA reveals extended cool gas and hot ionized outflows in a typical star-forming galaxy at $z=7.13$

Authors: Hollis Akins, Seiji Fujimoto, Kristian Finlator, Darach Watson, Kirsten Knudsen, Johan Richard, Tom Bakx, Takuya Hashimoto, Akio Inoue, Hiroshi Matsuo, Michal Michalowski, Yoichi Tamura

Abstract: We present spatially-resolved morphological properties of [CII] 158 $μ$m, [OIII] 88 $μ$m, dust, and rest-frame ultraviolet (UV) continuum emission for A1689-zD1, a strongly lensed, sub-L* galaxy at $z=7.13$, by utilizing deep Atacama Large Millimeter/submillimeter Array (ALMA) and Hubble Space Telescope (HST) observations. While the [OIII] line and UV continuum are compact, the [CII] line is exten… ▽ More We present spatially-resolved morphological properties of [CII] 158 $μ$m, [OIII] 88 $μ$m, dust, and rest-frame ultraviolet (UV) continuum emission for A1689-zD1, a strongly lensed, sub-L* galaxy at $z=7.13$, by utilizing deep Atacama Large Millimeter/submillimeter Array (ALMA) and Hubble Space Telescope (HST) observations. While the [OIII] line and UV continuum are compact, the [CII] line is extended up to a radius of $r \sim 12$ kpc. Using multi-band rest-frame far-infrared (FIR) continuum data ranging from 52-400 $μ$m, we find an average dust temperature and emissivity index of $T_{\rm dust} = 41^{+17}_{-14}$ K and $β= 1.7^{+1.1}_{-0.7}$, respectively, across the galaxy. We find slight differences in the dust continuum profiles at different wavelengths, which may indicate that the dust temperature decreases with distance. We map the star-formation rate (SFR) via IR and UV luminosities and determine a total SFR of $37\pm 1~M_\odot~{\rm yr}^{-1}$ with an obscured fraction of $87\%$. While the [OIII] line is a good tracer of the SFR, the [CII] line shows deviation from the local $L_{\rm [CII]}$-SFR relations in the outskirts of the galaxy. Finally, we observe a clear difference in the line profile between [CII] and [OIII], with significant residuals ($\sim 5σ$) in the [OIII] line spectrum after subtracting a single Gaussian model. This suggests a possible origin of the extended [CII] structure from the cooling of hot ionized outflows. The extended [CII] and high-velocity [OIII] emission may both contribute in part to the high $L_{\rm [OIII]}$/$L_{\rm [CII]}$ ratios recently reported in $z>6$ galaxies. △ Less

Submitted 14 June, 2022; originally announced June 2022.

Comments: 22 pages, 14 figures, accepted to ApJ, presented at AAS240

CLEAR: The Evolution of Spatially Resolved Star Formation in Galaxies between $0.5\lesssim z \lesssim1.7$ using H$α$ Emission Line Maps

Authors: Jasleen Matharu, Casey Papovich, Raymond C. Simons, Ivelina Momcheva, Gabriel Brammer, Zhiyuan Ji, Bren E. Backhaus, Nikko J. Cleri, Vicente Estrada-Carpenter, Steven L. Finkelstein, Kristian Finlator, Mauro Giavalisco, Intae Jung, Adam Muzzin, Annalisa Pillepich, Jonathan R. Trump, Benjamin Weiner

Abstract: Using spatially resolved H-alpha emission line maps of star-forming galaxies, we study the evolution of gradients in galaxy assembly over a wide range in redshift ($0.5<z<1.7$). Our $z\sim0.5$ measurements come from deep Hubble Space Telescope WFC3 G102 grism spectroscopy obtained as part of the CANDELS Lyman-alpha Emission at Reionization (CLEAR) Experiment. For star-forming galaxies with Log… ▽ More Using spatially resolved H-alpha emission line maps of star-forming galaxies, we study the evolution of gradients in galaxy assembly over a wide range in redshift ($0.5<z<1.7$). Our $z\sim0.5$ measurements come from deep Hubble Space Telescope WFC3 G102 grism spectroscopy obtained as part of the CANDELS Lyman-alpha Emission at Reionization (CLEAR) Experiment. For star-forming galaxies with Log$(M_{*}/\mathrm{M}_{\odot})\geqslant8.96$, the mean H-alpha effective radius is $1.2\pm0.1$ times larger than that of the stellar continuum, implying inside-out growth via star formation. This measurement agrees within $1σ$ with those measured at $z\sim1$ and $z\sim1.7$ from the 3D-HST and KMOS-3D surveys respectively, implying no redshift evolution. However, we observe redshift evolution in the stellar mass surface density within 1 kiloparsec ($Σ_\mathrm{1kpc}$). Star-forming galaxies at $z\sim0.5$ with a stellar mass of Log$(M_{*}/\mathrm{M}_{\odot})=9.5$ have a ratio of $Σ_\mathrm{1kpc}$ in H-alpha relative to their stellar continuum that is lower by $(19\pm2)\%$ compared to $z\sim1$ galaxies. $Σ_{1\mathrm{kpc, H}α}$/$Σ_{1\mathrm{kpc,Cont}}$ decreases towards higher stellar masses. The majority of the redshift evolution in $Σ_{1\mathrm{kpc,H}α}$/$Σ_{1\mathrm{kpc,Cont}}$ versus stellar mass stems from the fact that Log($Σ_{1\mathrm{kpc, H}α}$) declines twice as much as Log($Σ_{1\mathrm{kpc, Cont}}$) from $z\sim 1$ to 0.5 (at a fixed stellar mass of Log$(M_{*}/\mathrm{M}_{\odot})=9.5$). By comparing our results to the TNG50 cosmological magneto-hydrodynamical simulation, we rule out dust as the driver of this evolution. Our results are consistent with inside-out quenching following in the wake of inside-out growth, the former of which drives the significant drop in $Σ_{1\mathrm{kpc, H}α}$ from $z\sim1$ to $z\sim0.5$. △ Less

Submitted 17 May, 2022; originally announced May 2022.

Comments: 25 pages, 12 figures, 4 tables, Submitted to ApJ

Assuming Ionization Equilibrium and the Impact on the Lyman-$\mathrmα$ Forest Power Spectrum during the End of Reionization at $8 \geq z \geq 5$

Authors: Samir Kušmić, Kristian Finlator, Laura Keating, Ezra Huscher

Abstract: We explore how the assumption of ionization equilibrium modulates the modeled intergalactic medium (IGM) at the end of the hydrogen Epoch of Reionization using the cosmological radiation hydrodynamic \textsc{Technicolor Dawn} simulation. In neutral and partially-ionized regions where the metagalactic ultraviolet background (UVB) is weak, the ionization timescale $t_\mathrm{ion}\equiv Γ^{-1}$ excee… ▽ More We explore how the assumption of ionization equilibrium modulates the modeled intergalactic medium (IGM) at the end of the hydrogen Epoch of Reionization using the cosmological radiation hydrodynamic \textsc{Technicolor Dawn} simulation. In neutral and partially-ionized regions where the metagalactic ultraviolet background (UVB) is weak, the ionization timescale $t_\mathrm{ion}\equiv Γ^{-1}$ exceeds the Hubble time. Assuming photoionization equilibrium in such regions artificially boosts the ionization rate, accelerating reionization. By contrast, the recombination time $t_\mathrm{rec} < t_\mathrm{ion}$ in photoionized regions, with the result that assuming photoionization equilibrium artificially increases the neutral hydrogen fraction. Using snapshots between $8 \geq z \geq 5$, we compare the predicted Lyman-$α$ forest flux power spectrum with and without the assumption of ionization equilibrium. Small scales ($k > 0.1$ rad s km$^{-1}$) exhibit reduced power from $7 \leq z \leq 5.5$ in the ionization equilibrium case while larger scales are unaffected. This occurs for the same reasons: ionization equilibrium artificially suppresses the neutral fraction in self-shielded gas and boosts ionizations in voids, suppressing small-scale fluctuations in the ionization field. When the volume-averaged neutral fraction drops below $10^{-4}$, the signature of non-equilibrium ionizations on the Lyman-$α$ forest (LAF) disappears. Comparing with recent observations indicates that these non-equilibrium effects are not yet observable in the LAF flux power spectrum. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Comments: This has been submitted and accepted to ApJ as of April 11, 2022

The evolution of the Si IV content in the Universe from the epoch of reionization to cosmic noon

Authors: V. D'Odorico, K. Finlator, S. Cristiani, G. Cupani, S. Perrotta, F. Calura, M. Cènturion, G. Becker, T. A. M. Berg, S. Lopez, S. Ellison, E. Pomante

Abstract: We investigate the abundance and distribution of metals in the high-redshift intergalactic medium and circum-galactic medium through the analysis of a sample of almost 600 SiIV absorption lines detected in high and intermediate resolution spectra of 147 quasars. The evolution of the number density of SiIV lines, the column density distribution function and the cosmic mass density are studied in th… ▽ More We investigate the abundance and distribution of metals in the high-redshift intergalactic medium and circum-galactic medium through the analysis of a sample of almost 600 SiIV absorption lines detected in high and intermediate resolution spectra of 147 quasars. The evolution of the number density of SiIV lines, the column density distribution function and the cosmic mass density are studied in the redshift interval 1.7 <= z <= 6.2 and for log N(SiIV) >= 12.5. All quantities show a rapid increase between z~6 and z< 5 and then an almost constant behaviour to z~2 in very good agreement with what is already observed for CIV absorption lines. The present results are challenging for numerical simulations: when simulations reproduce our SiIV results, they tend to underpredict the properties of CIV, and when the properties of CIV are reproduced, the number of strong SiIV lines (log N(SiIV) > 14) is overpredicted. △ Less

Submitted 24 February, 2022; originally announced February 2022.

Comments: 14 pages, 8 figures, accepted for publication in MNRAS

Using Multiple Emission Line Ratios to Constrain the Slope of the Dust Attenuation Law

Authors: Moire K. M. Prescott, Kristian M. Finlator, Nikko J. Cleri, Jonathan R. Trump, Casey Papovich

Abstract: We explore the possibility and practical limitations of using a three-line approach to measure both the slope and normalization of the dust attenuation law in individual galaxies. To do this, we focus on a sample of eleven galaxies with existing ground-based Balmer H$α$ and H$β$ measurements from slit spectra, plus space-based grism constraints on Paschen-$β$. When accounting for observational unc… ▽ More We explore the possibility and practical limitations of using a three-line approach to measure both the slope and normalization of the dust attenuation law in individual galaxies. To do this, we focus on a sample of eleven galaxies with existing ground-based Balmer H$α$ and H$β$ measurements from slit spectra, plus space-based grism constraints on Paschen-$β$. When accounting for observational uncertainties, we show that one galaxy has a well-constrained dust law slope and normalization in the range expected from theoretical arguments; this galaxy therefore provides an example of what may be possible in the future. However, most of the galaxies are best-fit by unusually steep or shallow slopes. We then explore whether additional astrophysical effects or observational biases could explain the elevated Paschen-$β$/H$α$ ratios driving these results. We find that galaxies with high Paschen-$β$/H$α$ ratios may be explained by slightly sub-unity covering fractions ($>$97%). Alternatively, differing slit losses for different lines can have a large impact on the results, emphasizing the importance of measuring all three lines with a consistent spectroscopic aperture. We conclude that while the three-line approach to constraining the shape of the dust attenuation law in individual galaxies is promising, deep observations and a consistent observational strategy will be required to minimize observational biases and to disentangle the astrophysically interesting effect of differing covering fractions. The James Webb Space Telescope will provide more sensitive measurements of Balmer and Paschen lines for galaxies at $z\approx0.3-2$, enabling uniform constraints on the optical-infrared dust attenuation law and its intrinsic variation. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: 7 figures, 3 tables, accepted to ApJ

CLEAR: Boosted Ly$α$ Transmission of the Intergalactic Medium in UV bright Galaxies

Authors: Intae Jung, Casey Papovich, Steven L. Finkelstein, Raymond C. Simons, Vicente Estrada-Carpenter, Bren E. Backhaus, Nikko J. Cleri, Kristian Finlator, Mauro Giavalisco, Zhiyuan Ji, Jasleen Matharu, Ivelina Momcheva, Amber N. Straughn, Jonathan R. Trump

Abstract: Reionization is an inhomogeneous process, thought to begin in small ionized bubbles of the intergalactic medium (IGM) around overdense regions of galaxies. Recent Lyman-alpha (Ly$α$) studies during the epoch of reionization show growing evidence that ionized bubbles formed earlier around brighter galaxies, suggesting higher IGM transmission of Ly$α$ from these galaxies. We investigate this problem… ▽ More Reionization is an inhomogeneous process, thought to begin in small ionized bubbles of the intergalactic medium (IGM) around overdense regions of galaxies. Recent Lyman-alpha (Ly$α$) studies during the epoch of reionization show growing evidence that ionized bubbles formed earlier around brighter galaxies, suggesting higher IGM transmission of Ly$α$ from these galaxies. We investigate this problem using IR slitless spectroscopy from the Hubble Space Telescope (HST) Wide-Field Camera 3 (WFC3) G102 grism observations of 148 galaxies selected via photometric redshifts at $6.0<z<8.2$. These galaxies have spectra extracted from the CANDELS Ly$α$ Emission at Reionization (CLEAR) survey. We combine the CLEAR data for 275 galaxies with the Keck/DEIMOS+MOSFIRE dataset from the Texas Spectroscopic Search for Ly$α$ Emission at the End of Reionization Survey. We then constrain the Ly$α$ equivalent-width (EW) distribution at $6.0<z<8.2$, which is described by an exponential form, $dN/d\text{EW}\propto\text{exp(-EW)}/W_0$, with the characteristic $e$-folding scale width ($W_0$). We confirm a significant drop of the Ly$α$ strength (or $W_0$) at $z>6$. Furthermore, we compare the redshift evolution of $W_0$ between galaxies at different UV luminosities. The UV-bright ($M_{\text{UV}}<-21$, or $L_{\text{UV}}>L^{*}$) galaxies show weaker evolution with a decrease of 0.4 ($\pm$0.2) dex in $W_0$ at $z>6$ while UV-faint ($M_{\text{UV}}>-21$, or $L_{\text{UV}}<L^{*}$) galaxies exhibit a significant drop by a factor of 0.7-0.8 ($\pm0.2$) dex in $W_0$ from $z<6$ to $z>6$. Our results add to the accumulating evidence that UV-bright galaxies exhibit boosted Ly$α$ transmission in the IGM, suggesting that reionization completes sooner in regions proximate to galaxies of higher UV luminosity. △ Less

Submitted 29 November, 2021; originally announced November 2021.

Comments: 21 pages, 9 figures, 2 tables, 1st revision submitted to ApJ

A [C II] 158$μ$m Emitter Associated with an OI Absorber at the End of the Reionization Epoch

Authors: Yunjing Wu, Zheng Cai, Marcel Neeleman, Kristian Finlator, Shiwu Zhang, J. Xavier Prochaska, Ran Wang, Bjorn H., C Emonts, Xiaohui Fan, Laura C. Keating, Feige Wang, Jinyi Yang, Joseph F. Hennawi, Junxian Wang

Abstract: The physical and chemical properties of the circumgalactic medium (CGM) at $z\gtrsim6$ have been studied successfully through the absorption in the spectra of background Quasi-Stellar Objects (QSOs). One of the most crucial questions is to investigate the nature and location of the source galaxies that give rise to these early metal absorbers. Theoretical models suggest that momentum-driven outflo… ▽ More The physical and chemical properties of the circumgalactic medium (CGM) at $z\gtrsim6$ have been studied successfully through the absorption in the spectra of background Quasi-Stellar Objects (QSOs). One of the most crucial questions is to investigate the nature and location of the source galaxies that give rise to these early metal absorbers. Theoretical models suggest that momentum-driven outflows from typical star-forming galaxies can eject metals into the CGM and the intergalactic medium (IGM) at z=5-6. Deep, dedicated surveys have searched for Ly$α$ emission associated with strong CIV absorbers at $z\approx 6$, but only a few Ly$α$ emitter candidates have been detected. Interpreting these detections is moreover ambiguous because Ly$α$ is a resonant line, raising the need for complementary techniques for detecting absorbers' host galaxies. Here, using Atacama Large Millimeter Array (ALMA), we report a [C II] 158$μ$m emitter associated with a strong low-ionization absorber, OI, at $z=5.978$. The projected impact parameter between OI and [C II] emitter is 20.0 kpc. The measured [C II] luminosity is $7.0\times 10^7$ solar luminosities. Further analysis indicates that strong OI absorbers may reside in the circumgalactic medium of massive halos one to two orders of magnitude more massive than expected values. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: Published in Nature Astronomy on Sep 27, 2021; 45 pages; authors' version. Link to the paper: https://www.nature.com/articles/s41550-021-01471-4

The effects of binary stars on galaxies and metal-enriched gas during reionization

Authors: Caitlin C. Doughty, Kristian M. Finlator

Abstract: Binary stars are abundant in nearby galaxies, but are typically unaccounted for in simulations of the high redshift Universe. Stellar population synthesis models that include the effects of binary evolution result in greater relative abundances of ionizing photons that could significantly affect the ambient ionizing background during the epoch of hydrogen reionization, additionally leading to diff… ▽ More Binary stars are abundant in nearby galaxies, but are typically unaccounted for in simulations of the high redshift Universe. Stellar population synthesis models that include the effects of binary evolution result in greater relative abundances of ionizing photons that could significantly affect the ambient ionizing background during the epoch of hydrogen reionization, additionally leading to differences in galaxy gas content and star formation. We use hydrodynamic cosmological simulations including in situ multifrequency radiative transfer to evaluate the effects of a high binary fraction in reionization-era galaxies on traits of the early intergalactic medium and the abundance of H I and He II ionizing photons. We further extend this to analyze the traits of enriched gas. In comparing metrics generated using a fiducial simulation assuming single stars with one incorporating a high binary fraction, we find that binary stars cause H I reionization to complete earlier and at an accelerated pace, while also increasing the abundances of high-ionization metals (C IV and Si IV) in simulated absorption spectra while reducing the abundance of low-ionization states (O I, Si II, and C II). However, through increased photoheating of galactic and circumgalactic gas, they simultaneously reduce the rate of star formation in low-mass galaxies, slowing the ongoing process of enrichment and suppressing their own ionizing background. This potentially contributes to a slower He II reionization process at $z\geq5$, and further indicates that self-regulation of galaxies could be underestimated when neglecting binary stellar evolution. △ Less

Submitted 20 May, 2021; originally announced May 2021.

Comments: 19 pages, 14 figures. Accepted to MNRAS

SED Analysis of 13 Spectroscopically Confirmed Galaxies at z$\simeq$6 to Constrain UV-Slope, Model Dust Attenuation and Escape Fractions

Authors: Junehyoung Jeon, Rogier Windhorst, Seth H. Cohen, Rolf A. Jansen, Brent M. Smith, Timothy Carleton, Eiichi Egami, Kristian Finlator, Linhua Jiang, Kartheik G. Iyer, Vihang Mehta

Abstract: The reionization of the hydrogen in the Universe is thought to have completed by redshift $z\simeq5.5-6$. To probe this era, galaxy observations in the Subaru Deep Field (SDF) have identified more than 100 galaxies at $z\simeq6$, many spectroscopically confirmed through follow-up observations. We model the spectral energy distributions (SEDs) of 13 SDF galaxies with the CIGALE and Dense Basis code… ▽ More The reionization of the hydrogen in the Universe is thought to have completed by redshift $z\simeq5.5-6$. To probe this era, galaxy observations in the Subaru Deep Field (SDF) have identified more than 100 galaxies at $z\simeq6$, many spectroscopically confirmed through follow-up observations. We model the spectral energy distributions (SEDs) of 13 SDF galaxies with the CIGALE and Dense Basis codes using available optical/IR data. Modeling deep IR photometry has the potential to constrain the galaxy's Lyman continuum (LyC) escape fraction (\fesc). We use the modeled nebular emission lines and find that the implied escape fractions ranges from \textbf{0 to 0.8 with a median of $\sim$0.35 for Dense Basis and $\sim$0.55 for CIGALE.} Significant uncertainties in the data exist, so that fitting results in a large range of \fesc\ for individual objects. The implied median \fesc-values may be high enough for galaxies to finish reionization by $z\sim6$. Furthermore, we find no strong trends between the UV-slope $β$ or \EBminV with model \fesc. If true, the lack of trends suggest that other factors besides nebular emission or dust extinction could have led to LyC escaping, such as the presence of holes in the ISM with sufficiently wide opening angles from outflows of supernovae and/or weak AGN, resulting in a range of implied \fesc-values depending on the viewing angle of each galaxy. The current \textit{HST, Spitzer} and ground-based photometric and model errors for the galaxies remain large, so IR spectroscopic observations with the \textit{James Webb Space Telescope} are needed to constrain this possibility. △ Less

Submitted 18 March, 2022; v1 submitted 11 November, 2020; originally announced November 2020.

Comments: 24 pages, 5 figures, 3 tables

Evolution of CIV Absorbers I. The Cosmic Incidence

Authors: Farhanul Hasan, Christopher W. Churchill, Bryson Stemock, Nigel L. Mathes, Nikole M. Nielsen, Kristian Finlator, Caitlin Doughty, Mark Croom, Glenn G. Kacprzak, Michael T. Murphy

Abstract: We present a large high-resolution study of the distribution and evolution of CIV absorbers, including the weakest population with equivalent widths $W_r<0.3$~Å. By searching 369 high-resolution, high signal-to-noise spectra of quasars at $1.1\leq z_{em} \leq5.3$ from Keck/HIRES and VLT/UVES, we find $1268$ CIV absorbers with $W_r \geq 0.05$~Å (our $\sim50\%$ completeness limit) at redshifts… ▽ More We present a large high-resolution study of the distribution and evolution of CIV absorbers, including the weakest population with equivalent widths $W_r<0.3$~Å. By searching 369 high-resolution, high signal-to-noise spectra of quasars at $1.1\leq z_{em} \leq5.3$ from Keck/HIRES and VLT/UVES, we find $1268$ CIV absorbers with $W_r \geq 0.05$~Å (our $\sim50\%$ completeness limit) at redshifts $1\leq z \leq4.75$. A Schechter function describes the observed equivalent width distribution with a transition from power-law to exponential decline at $W_r \gtrsim 0.5$~Å. The power-law slope $α$ rises by $\sim7\%$ and transition equivalent width $W_{\star}$ falls by $\sim\!20\%$ from $\langle z \rangle=1.7$ to $\langle z \rangle=3.6$. We find that the co-moving redshift path density, $dN/dX$, of $W_r \geq 0.05$~Å absorbers rises by $\sim1.8$ times from $z\simeq 4.0$ to $z\simeq 1.3$, while the $W_r \geq 0.6$~Å $dN/dX$ rises by a factor of $\sim8.5$. We quantify the observed evolution by a model in which $dN/dX$ decreases linearly with redshift. The model suggests that populations with larger $W_r$ thresholds evolve faster with redshift and appear later in the universe. The cosmological Technicolor Dawn simulations at $z=3-5$ over-produce the observed abundance of absorbers with $W_r<0.3$~Å, while yielding better agreement at higher $W_r$. Our empirical linear model successfully describes {CIV evolution in the simulations and the observed evolution of $W_r \geq 0.6$~Å CIV for the past $\sim12$ Gyr. Combining our measurements with the literature gives us a picture of CIV-absorbing structures becoming more numerous and/or larger in physical size over the last $\approx13$ Gyr of cosmic time ($z\sim6$ to $z\sim0$). △ Less

Submitted 2 October, 2020; v1 submitted 22 July, 2020; originally announced July 2020.

Comments: 27 pages, 13 figures, 7 tables; submitted to ApJ on July 22, 2020; revised September 28, 2020; accepted October 2, 2020

Luminous Lyman-alpha Emitters with Very Blue UV-continuum Slopes at Redshift 5.7 <= z <= 6.6

Authors: Linhua Jiang, Seth H. Cohen, Rogier A. Windhorst, Eiichi Egami, Kristian Finlator, Daniel Schaerer, Fengwu Sun

Abstract: We study six luminous Lyman-alpha emitters (LAEs) with very blue rest-frame UV continua at $5.7\le z \le 6.6$. These LAEs have previous HST and Spitzer IRAC observations. Combining our newly acquired HST images, we find that their UV-continuum slopes $β$ are in a range of $-3.4\le β\le -2.6$. Unlike previous, tentative detections of $β\simeq -3$ in photometrically selected, low-luminosity galaxies… ▽ More We study six luminous Lyman-alpha emitters (LAEs) with very blue rest-frame UV continua at $5.7\le z \le 6.6$. These LAEs have previous HST and Spitzer IRAC observations. Combining our newly acquired HST images, we find that their UV-continuum slopes $β$ are in a range of $-3.4\le β\le -2.6$. Unlike previous, tentative detections of $β\simeq -3$ in photometrically selected, low-luminosity galaxies, our LAEs are spectroscopically confirmed and luminous ($M_{\rm UV}<-20$ mag). We model their broadband spectral energy distributions (SEDs), and find that two $β\simeq-2.6\pm0.2$ galaxies can be well fitted with young and dust-free stellar populations. However, it becomes increasingly difficult to fit bluer galaxies. We explore further interpretations by including non-zero LyC escape fraction $f_{\rm esc}$, very low metallicities, and/or AGN contributions. Assuming $f_{\rm esc}\simeq0.2$, we achieve the bluest slopes $β\simeq-2.7$ when nebular emission is considered. This can nearly explain the SEDs of two galaxies with $β\simeq-2.8$ and --2.9 ($σ_β=0.15$). Larger $f_{\rm esc}$ values and very low metallicities are not favored by the strong nebular line emission (evidenced by the IRAC flux) or the observed (IRAC 1 - IRAC 2) color. Finally, we find that the $β\simeq-2.9$ galaxy can potentially be well explained by the combination of a very young population with a high $f_{\rm esc}$ ($\ge0.5$) and an old, dusty population. We are not able to produce two $β\simeq -3.4 \pm0.4$ galaxies. Future deep spectroscopic observations are needed to fully understand these galaxies. △ Less

Submitted 5 February, 2020; originally announced February 2020.

Comments: Published in ApJ on 2020 Feb 1; Authors' version (9 pages); See published version at https://iopscience.iop.org/article/10.3847/1538-4357/ab64ea

Journal ref: ApJ, 889, 90 (2020)

The Faint Host Galaxies of C IV Absorbers at z > 5

Authors: Kristian Finlator, Caitlin Doughty, Zheng Cai, Gonzalo Díaz

Abstract: We explore the expected galaxy environments of CIV absorbers at z>5 using the Technicolor Dawn simulations. These simulations reproduce the observed history of reionization, the z~6 galaxy stellar mass function, the Ly$α$ forest transmission at z>5, and the SiIV column density distribution (CDD) at z~5.5. Nonetheless, the CIV CDD remains underproduced. Comparison with observed CII/SiII equivalent… ▽ More We explore the expected galaxy environments of CIV absorbers at z>5 using the Technicolor Dawn simulations. These simulations reproduce the observed history of reionization, the z~6 galaxy stellar mass function, the Ly$α$ forest transmission at z>5, and the SiIV column density distribution (CDD) at z~5.5. Nonetheless, the CIV CDD remains underproduced. Comparison with observed CII/SiII equivalent width ratios and the CII line incidence suggests that a low carbon yield accounts for some, but not all, of the CIV discrepancy. Alternatively, a density-bounded escape scenario could harden the metagalactic ionizing background more dramatically even than binary stellar evolution, boosting the CIV CDD into near-agreement with observations. In this case galaxies ionize more efficiently and fewer are required to host a given high-ionization absorber. Absorbers' environments therefore constrain ionizing escape. Regardless of the escape scenario, galaxies correlate with CIV absorbers out to 300 proper kpc (pkpc). The correlation strengthens independently with galaxy luminosity and CIV column density. Around strong systems (log(N$_{\rm CIV}$/cm$^{-2}$)>14)), the overdensity of galaxies with M$_{\rm UV}$<-18 or log($L_{{\rm Ly}α}$/erg s$^{-1}$) > 41.9 declines from 200-300 within 100 pkpc to 40-60 within 250 pkpc. The previously-suggested association between strong CIV absorbers and Ly$α$ emitters at z>5 is not expected. It may arise if both populations inhabit large-scale voids, but for different reasons. Although most neighboring galaxies are too faint for HST, JWST will, with a single pointing, identify ~10 neighboring galaxies per strong CIV absorber at z > 5. Ground-based tests of these predictions are possible via deep surveys for Ly$α$ emission using integral field units. △ Less

Submitted 7 February, 2020; v1 submitted 10 January, 2020; originally announced January 2020.

Comments: 15 pages, 16 figures, resubmitted to MNRAS with improvements following first referee report

Weak Evolution of the Mass-metallicity Relation at Cosmic Dawn in the FirstLight Simulations

Authors: Ivanna Langan, Daniel Ceverino, Kristian Finlator

Abstract: Little is known about the mass-metallicity relation (MZR) in galaxies at cosmic dawn. Studying the first appearance of the MZR is one of the keys to understanding the formation and evolution of the first galaxies. In order to lay the groundwork for upcoming observational campaigns, we analyze 290 galaxies in halos spanning Mh = 10^9-10^11 Msun selected from the FirstLight (FL) cosmological zoom si… ▽ More Little is known about the mass-metallicity relation (MZR) in galaxies at cosmic dawn. Studying the first appearance of the MZR is one of the keys to understanding the formation and evolution of the first galaxies. In order to lay the groundwork for upcoming observational campaigns, we analyze 290 galaxies in halos spanning Mh = 10^9-10^11 Msun selected from the FirstLight (FL) cosmological zoom simulations to predict the MZR at z=5-8. Over this interval, the metallicity of FL galaxies with stellar mass Ms= 10^8 Msun declines by < 0.2 dex. This contrasts with the observed tendency for metallicities to increase at lower redshifts, and reflects weakly-evolving or even increasing gas fractions. We assess the use of the R3 strong-line diagnostic as a metallicity indicator, finding that it is informative for 12+log(O/H) < 8 but saturates to R3=3 at higher metallicities owing to a cancellation between enrichment and spectral softening. Nonetheless, campaigns with JWST should be able to detect a clear trend between R3 and stellar mass for Ms > 10^7.5 Msun. We caution that, at fixed metallicity, galaxies with higher specific star formation show higher R3 owing to their more intense radiation fields, indicating a potential for selection biases. △ Less

Submitted 5 March, 2020; v1 submitted 25 October, 2019; originally announced October 2019.

Comments: 6 pages, 4 figures, resubmitted to MNRAS. More discussion added, conclusions unchanged

Testing Galaxy Formation Simulations with Damped Lyman-$α$ Abundance and Metallicity Evolution

Authors: Sultan Hassan, Kristian Finlator, Romeel Davé, Christopher W. Churchill, J. Xavier Prochaska

Abstract: We examine the properties of damped Lyman-$α$ absorbers (DLAs) emerging from a single set of cosmological initial conditions in two state-of-the-art cosmological hydrodynamic simulations: {\sc Simba} and {\sc Technicolor Dawn}. The former includes star formation and black hole feedback treatments that yield a good match with low-redshift galaxy properties, while the latter uses multi-frequency rad… ▽ More We examine the properties of damped Lyman-$α$ absorbers (DLAs) emerging from a single set of cosmological initial conditions in two state-of-the-art cosmological hydrodynamic simulations: {\sc Simba} and {\sc Technicolor Dawn}. The former includes star formation and black hole feedback treatments that yield a good match with low-redshift galaxy properties, while the latter uses multi-frequency radiative transfer to model an inhomogeneous ultraviolet background (UVB) self-consistently and is calibrated to match the Thomson scattering optical depth, UVB amplitude, and Ly-$α$ forest mean transmission at $z>5$. Both simulations are in reasonable agreement with the measured stellar mass and star formation rate functions at $z\geq 3$, and both reproduce the observed neutral hydrogen cosmological mass density, $Ω_{\rm HI}(z)$. However, the DLA abundance and metallicity distribution are sensitive to the galactic outflows' feedback and the UVB amplitude. Adopting a strong UVB and/or slow outflows under-produces the observed DLA abundance, but yields broad agreement with the observed DLA metallicity distribution. By contrast, faster outflows eject metals to larger distances, yielding more metal-rich DLAs whose observational selection may be more sensitive to dust bias. The DLA metallicity distribution in models adopting an ${\rm H}_2$-regulated star formation recipe includes a tail extending to $[M/H] \ll -3$, lower than any DLA observed to date, owing to curtailed star formation in low-metallicity galaxies. Our results show that DLA observations play an imporant role in constraining key physical ingredients in galaxy formation models, complementing traditional ensemble statistics such as the stellar mass and star formation rate functions. △ Less

Submitted 9 January, 2020; v1 submitted 16 October, 2019; originally announced October 2019.

Comments: Accepted for publication in MNRAS

Evolution of Neutral Oxygen During the Epoch of Reionization and its Use in Estimating the Neutral Hydrogen Fraction

Authors: Caitlin Doughty, Kristian Finlator

Abstract: We use synthetic sightlines drawn through snapshots of the Technicolor Dawn simulations to explore how the statistics of neutral oxygen OI absorbers respond to hydrogen reionization. The ionization state of the circumgalactic medium (CGM) initially roughly tracks that of the intergalactic medium, but beginning at $z=8$ the CGM grows systematically more neutral owing to self-shielding. Weak absorbe… ▽ More We use synthetic sightlines drawn through snapshots of the Technicolor Dawn simulations to explore how the statistics of neutral oxygen OI absorbers respond to hydrogen reionization. The ionization state of the circumgalactic medium (CGM) initially roughly tracks that of the intergalactic medium, but beginning at $z=8$ the CGM grows systematically more neutral owing to self-shielding. Weak absorbers trace diffuse gas that lies farther from halos, hence they are ionized first, whereas stronger systems are less sensitive to reionization. The overall OI covering fraction decreases slowly with time owing to competition between ongoing enrichment and gradual encroachment of ionization fronts into increasingly overdense gas. While the declining covering fraction is partially offset by continued formation of new halos, the ionization of the diffuse gas causes the predicted line-of-sight incidence rate of OI absorbers to decline abruptly at the overlap epoch, in qualitative agreement with observations. In comparison to the recently-observed equivalent width (EW) distribution at $z\approx6$, the simulations underproduce systems with $EW \geq 0.1 \unicode{x212B}$, although they reproduce weaker systems with $EW \geq 0.05 \unicode{x212B}$. By $z\approx5$, the incidence of $EW < 0.1 \unicode{x212B}$ systems are overproduced, consistent with previous indications that the simulated ionizing background is too weak at $z<6$. The summed column densities of SiII and SiIV trace the total oxygen column, and hence the ratio of the OI and SiII + SiIV comoving mass densities traces the progress of reionization. This probe may prove particularly useful in the regime where $x_{HI} > 10\%$ △ Less

Submitted 22 August, 2019; originally announced August 2019.

Comments: 15 pages, 10 figures. Accepted to MNRAS

Conditions for Reionizing the Universe with A Low Galaxy Ionizing Photon Escape Fraction

Authors: Steven L. Finkelstein, Anson D'Aloisio, Jan-Pieter Paardekooper, Russell Ryan Jr., Peter Behroozi, Kristian Finlator, Rachael Livermore, Phoebe R. Upton Sanderbeck, Claudio Dalla Vecchia, Sadegh Khochfar

Abstract: We explore scenarios for reionizing the intergalactic medium with low galaxy ionizing photon escape fractions. We combine simulation-based halo-mass dependent escape fractions with an extrapolation of the observed galaxy rest-ultraviolet luminosity functions to solve for the reionization history from z=20 to z=4. We explore the posterior distributions for key unknown quantities, including the limi… ▽ More We explore scenarios for reionizing the intergalactic medium with low galaxy ionizing photon escape fractions. We combine simulation-based halo-mass dependent escape fractions with an extrapolation of the observed galaxy rest-ultraviolet luminosity functions to solve for the reionization history from z=20 to z=4. We explore the posterior distributions for key unknown quantities, including the limiting halo mass for star-formation, the ionizing photon production efficiency, and a potential contribution from active galactic nuclei (AGN). We marginalize over the allowable parameter space using a Markov Chain Monte Carlo method, finding a solution which satisfies the most model-independent constraints on reionization. Our fiducial model can match observational constraints with an average escape fraction of <5% throughout the bulk of the epoch of reionization if: i) galaxies form stars down to the atomic cooling limit before reionization and a photosuppression mass of log(M_h/Msol)~9 during/after reionization (-13<M_UV,lim<-11); ii) galaxies become more efficient producers of ionizing photons at higher redshifts and fainter magnitudes, and iii) there is a significant, but sub-dominant, contribution by AGN at z < 7. In this model the faintest galaxies (M_UV>-15) dominate the ionizing emissivity, leading to an earlier start to reionization and a smoother evolution of the ionized volume filling fraction than models which assume a single escape fraction at all redshifts and luminosities. The ionizing emissivity from this model is consistent with observations at z=4-5 (and below, when extrapolated), in contrast to some models which assume a single escape fraction. Our predicted ionized volume filling fraction at z=7 of Q_HII=78% (+\- 8%) is in ~1-2 sigma tension with observations of Lya emitters at z~7 and the damping wing analyses of the two known z>7 quasars, which prefer Q_HII,z=7~40-50%. △ Less

Submitted 2 May, 2019; v1 submitted 7 February, 2019; originally announced February 2019.

Comments: 45 pages, 21 figures, accepted for publication in the Astrophysical Journal

ARTIST: Fast radiative transfer for large-scale simulations of the epoch of reionisation

Authors: Margherita Molaro, Romeel Davé, Sultan Hassan, Mario G. Santos, Kristian Finlator

Abstract: We introduce the "Asymmetric Radiative Transfer In Shells Technique" (ARTIST), a new method for photon propagation on large scales that explicitly conserves photons, propagates photons at the speed of light, approximately accounts for photon directionality, and closely reproduces results of more detailed radiative transfer (RT) codes. Crucially, it is computationally fast enough to evolve the larg… ▽ More We introduce the "Asymmetric Radiative Transfer In Shells Technique" (ARTIST), a new method for photon propagation on large scales that explicitly conserves photons, propagates photons at the speed of light, approximately accounts for photon directionality, and closely reproduces results of more detailed radiative transfer (RT) codes. Crucially, it is computationally fast enough to evolve the large cosmological volumes required to predict the 21cm power spectrum on scales that will be probed by future experiments targeting the Epoch of Reionisation (EoR). Most semi-numerical models aimed at predicting the EoR 21cm signal make use of an excursion set formalism (ESF) approach, which achieves computational viability by compromising on photon conservation, constraining ionised regions to be spherical by construction, and not accounting for light-travel time. By implementing our RT method within the semi-numerical code SimFast21, we show that ARTIST predicts a significantly different evolution for the EoR ionisation field compared to the code's native ESF. In particular, ARTIST predicts a more gradual evolution of the volume-averaged ionisation fraction, and up to an order-of-magnitude difference in the ionisation power, depending on the physical parameters assumed. Its application to large-scale EoR simulations will therefore allow more physically-motivated constraints to be obtained for key EoR parameters, such as the escape fraction. △ Less

Submitted 7 January, 2020; v1 submitted 10 January, 2019; originally announced January 2019.

Comments: 18 pages, 16 figures, published by MNRAS (14 August 2019)

Journal ref: 2019MNRAS.489.5594M

Reionization in Technicolor

Authors: Kristian Finlator, Laura Keating, Benjamin D. Oppenheimer, Romeel Davé, Erik Zackrisson

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

Submitted 19 July, 2018; v1 submitted 30 April, 2018; originally announced May 2018.

Comments: 24 pages, 17 figures, accepted to MNRAS

Aligned metal absorbers and the ultraviolet background at the end of reionization

Authors: Caitlin Doughty, Kristian Finlator, Benjamin D. Oppenheimer, Romeel Davé, Erik Zackrisson

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

Submitted 18 January, 2018; originally announced January 2018.

Comments: 13 pages, 4 figures, accepted by MNRAS

Demographics of Star-forming Galaxies since $z\sim2.5$. I. The $UVJ$ Diagram in CANDELS

Authors: Jerome J. Fang, S. M. Faber, David C. Koo, Aldo Rodriguez-Puebla, Yicheng Guo, Guillermo Barro, Peter Behroozi, Gabriel Brammer, Zhu Chen, Avishai Dekel, Henry C. Ferguson, Eric Gawiser, Mauro Giavalisco, Jeyhan Kartaltepe, Dale D. Kocevski, Anton M. Koekemoer, Elizabeth J. McGrath, Daniel McIntosh, Jeffrey A. Newman, Camilla Pacifici, Viraj Pandya, Pablo G. Perez-Gonzalez, Joel R. Primack, Brett Salmon, Jonathan R. Trump , et al. (19 additional authors not shown)

Abstract: This is the first in a series of papers examining the demographics of star-forming galaxies at $0.2<z<2.5$ in CANDELS. We study 9,100 galaxies from GOODS-S and UDS having published values of redshifts, masses, star-formation rates (SFRs), and dust attenuation ($A_V$) derived from UV-optical SED fitting. In agreement with previous works, we find that the $UVJ$ colors of a galaxy are closely correla… ▽ More This is the first in a series of papers examining the demographics of star-forming galaxies at $0.2<z<2.5$ in CANDELS. We study 9,100 galaxies from GOODS-S and UDS having published values of redshifts, masses, star-formation rates (SFRs), and dust attenuation ($A_V$) derived from UV-optical SED fitting. In agreement with previous works, we find that the $UVJ$ colors of a galaxy are closely correlated with its specific star-formation rate (SSFR) and $A_V$. We define rotated $UVJ$ coordinate axes, termed $S_\mathrm{SED}$ and $C_\mathrm{SED}$, that are parallel and perpendicular to the star-forming sequence and derive a quantitative calibration that predicts SSFR from $C_\mathrm{SED}$ with an accuracy of ~0.2 dex. SFRs from UV-optical fitting and from UV+IR values based on Spitzer/MIPS 24 $μ\mathrm{m}$ agree well overall, but systematic differences of order 0.2 dex exist at high and low redshifts. A novel plotting scheme conveys the evolution of multiple galaxy properties simultaneously, and dust growth, as well as star-formation decline and quenching, exhibit "mass-accelerated evolution" ("downsizing"). A population of transition galaxies below the star-forming main sequence is identified. These objects are located between star-forming and quiescent galaxies in $UVJ$ space and have lower $A_V$ and smaller radii than galaxies on the main sequence. Their properties are consistent with their being in transit between the two regions. The relative numbers of quenched, transition, and star-forming galaxies are given as a function of mass and redshift. △ Less

Submitted 6 April, 2018; v1 submitted 15 October, 2017; originally announced October 2017.

Comments: 36 pages, 26 figures, ApJ accepted

Probing the Metal Enrichment of the Intergalactic Medium at $z=5-6$ Using the Hubble Space Telescope

Authors: Zheng Cai, Xiaohui Fan, Romeel Dave, Kristian Finlator, Ben Oppenheimer

Abstract: We test the galactic outflow model by probing associated galaxies of four strong intergalactic CIV absorbers at $z=5$--6 using the Hubble Space Telescope (HST) ACS ramp narrowband filters. The four strong CIV absorbers reside at $z=5.74$, $5.52$, $4.95$, and $4.87$, with column densities ranging from $N_{\rm{CIV}}=10^{13.8}$ cm$^{-2}$ to $10^{14.8}$ cm$^{-2}$. At $z=5.74$, we detect an i-dropout L… ▽ More We test the galactic outflow model by probing associated galaxies of four strong intergalactic CIV absorbers at $z=5$--6 using the Hubble Space Telescope (HST) ACS ramp narrowband filters. The four strong CIV absorbers reside at $z=5.74$, $5.52$, $4.95$, and $4.87$, with column densities ranging from $N_{\rm{CIV}}=10^{13.8}$ cm$^{-2}$ to $10^{14.8}$ cm$^{-2}$. At $z=5.74$, we detect an i-dropout Ly$α$ emitter (LAE) candidate with a projected impact parameter of 42 physical kpc from the CIV absorber. This LAE candidate has a Ly$α$-based star formation rate (SFR$_{\rm{Lyα}}$) of 2 $M_\odot$ yr$^{-1}$ and a UV-based SFR of 4 $M_\odot$ yr$^{-1}$. Although we cannot completely rule out that this $i$-dropout emitter may be an [OII] interloper, its measured properties are consistent with the CIV powering galaxy at $z=5.74$. For CIV absorbers at $z=4.95$ and $z=4.87$, although we detect two LAE candidates with impact parameters of 160 kpc and 200 kpc, such distances are larger than that predicted from the simulations. Therefore we treat them as non-detections. For the system at $z=5.52$, we do not detect LAE candidates, placing a 3-$σ$ upper limit of SFR$_{\rm{Lyα}}\approx 1.5\ M_\odot$ yr$^{-1}$. In summary, in these four cases, we only detect one plausible CIV source at $z=5.74$. Combining the modest SFR of the one detection and the three non-detections, our HST observations strongly support that smaller galaxies (SFR$_{\rm{Lyα}} \lesssim 2\ M_\odot$ yr$^{-1}$) are main sources of intergalactic CIV absorbers, and such small galaxies play a major role in the metal enrichment of the intergalactic medium at $z\gtrsim5$. △ Less

Submitted 28 September, 2017; originally announced September 2017.

Comments: Accepted for Publications in ApJL

Constraining the contribution of active galactic nuclei to reionisation

Authors: Sultan Hassan, Romeel Davé, Sourav Mitra, Kristian Finlator, Benedetta Ciardi, Mario G. Santos

Abstract: Recent results have suggested that active galactic nuclei (AGN) could provide enough photons to reionise the Universe. We assess the viability of this scenario using a semi-numerical framework for modeling reionisation, to which we add a quasar contribution by constructing a Quasar Halo Occupation Distribution (QHOD) based on Giallongo et al. observations. Assuming a constant QHOD, we find that an… ▽ More Recent results have suggested that active galactic nuclei (AGN) could provide enough photons to reionise the Universe. We assess the viability of this scenario using a semi-numerical framework for modeling reionisation, to which we add a quasar contribution by constructing a Quasar Halo Occupation Distribution (QHOD) based on Giallongo et al. observations. Assuming a constant QHOD, we find that an AGN-only model cannot simultaneously match observations of the optical depth $τ_e$, neutral fraction, and ionising emissivity. Such a model predicts $τ_e$ too low by $\sim 2σ$ relative to Planck constraints, and reionises the Universe at $z\lesssim 5$. Arbitrarily increasing the AGN emissivity to match these results yields a strong mismatch with the observed ionising emissivity at $z\sim 5$. If we instead assume a redshift-independent AGN luminosity function yielding an emissivity evolution like that assumed in Madau & Haardt model, then we can match $τ_e$ albeit with late reionisation, however such evolution is inconsistent with observations at $z\sim 4-6$ and poorly motivated physically. These results arise because AGN are more biased towards massive halos than typical reionising galaxies, resulting in stronger clustering and later formation times. AGN-dominated models produce larger ionising bubbles that are reflected in $\sim\times 2$ more 21cm power on all scales. A model with equal parts galaxies and AGN contribution is still (barely) consistent with observations, but could be distinguished using next-generation 21cm experiments HERA and SKA-low. We conclude that, even with recent claims of more faint AGN than previously thought, AGN are highly unlikely to dominate the ionising photon budget for reionisation. △ Less

Submitted 2 October, 2017; v1 submitted 15 May, 2017; originally announced May 2017.

Comments: 16 pages, 9 figures, matches the accepted version for publication in MNRAS, 2017

Journal ref: Hassan et al. 2018, MNRAS, 473, 227

Epoch of Reionisation 21cm Forecasting From MCMC-Constrained Semi-Numerical Models

Authors: Sultan Hassan, Romeel Dave, Kristian Finlator, Mario G. Santos

Abstract: The recent low value of Planck (2016) integrated optical depth to Thomson scattering suggests that the reionization occurred fairly suddenly, disfavoring extended reionization scenarios. This will have a significant impact on the 21cm power spectrum. Using a semi-numerical framework, we improve our model from Hassan et al. (2016) to include time-integrated ionisation and recombination effects, and… ▽ More The recent low value of Planck (2016) integrated optical depth to Thomson scattering suggests that the reionization occurred fairly suddenly, disfavoring extended reionization scenarios. This will have a significant impact on the 21cm power spectrum. Using a semi-numerical framework, we improve our model from Hassan et al. (2016) to include time-integrated ionisation and recombination effects, and find that this leads to more sudden reionisation. It also yields larger HII bubbles which leads to an order of magnitude more 21cm power on large scales, while suppressing the small scale ionization power. Local fluctuations in the neutral hydrogen density play the dominant role in boosting the 21cm power spectrum on large scales, while recombinations are subdominant. We use a Monte Carlo Markov Chain approach to constrain our model to observations of the star formation rate functions at z = 6,7,8 from Bouwens et al. (2015), the Planck (2016) optical depth measurements, and the Becker & Bolton (2013) ionising emissivity data at z~5. We then use this constrained model to perform 21cm forecasting for LOFAR, HERA, and SKA in order to determine how well such data can characterise the sources driving reionisation. We find that the 21cm power spectrum alone can somewhat constrain the halo mass dependence of ionising sources, the photon escape fraction and ionising amplitude, but combining the 21cm data with other current observations enables us to separately constrain all these parameters. Our framework illustrates how 21cm data can play a key role in understanding the sources and topology of reionisation as observations improve. △ Less

Submitted 24 March, 2017; v1 submitted 17 December, 2016; originally announced December 2016.

Comments: 20 pages, 16 figues, matches the accepted version for publication in MNRAS

Journal ref: Mon Not R Astron Soc (2017) 468 (1): 122-139

Gas Accretion and Galactic Chemical Evolution: Theory and Observations

Authors: Kristian Finlator

Abstract: This chapter reviews how galactic inflows influence galaxy metallicity. The goal is to discuss predictions from theoretical models, but particular emphasis is placed on the insights that result from using models to interpret observations. Even as the classical G-dwarf problem endures in the latest round of observational confirmation, a rich and tantalizing new phenomenology of relationships betwee… ▽ More This chapter reviews how galactic inflows influence galaxy metallicity. The goal is to discuss predictions from theoretical models, but particular emphasis is placed on the insights that result from using models to interpret observations. Even as the classical G-dwarf problem endures in the latest round of observational confirmation, a rich and tantalizing new phenomenology of relationships between $M_*$, $Z$, SFR, and gas fraction is emerging both in observations and in theoretical models. A consensus interpretation is emerging in which star-forming galaxies do most of their growing in a quiescent way that balances gas inflows and gas processing, and metal dilution with enrichment. Models that explicitly invoke this idea via equilibrium conditions can be used to infer inflow rates from observations, while models that do not assume equilibrium growth tend to recover it self-consistently. Mergers are an overall subdominant mechanism for delivering fresh gas to galaxies, but they trigger radial flows of previously-accreted gas that flatten radial gas-phase metallicity gradients and temporarily suppress central metallicities. Radial gradients are generically expected to be steep at early times and then flattened by mergers and enriched inflows of recycled gas at late times. However, further theoretical work is required in order to understand how to interpret observations. Likewise, more observational work is needed in order to understand how metallicity gradients evolve to high redshifts. △ Less

Submitted 2 December, 2016; originally announced December 2016.

Comments: Invited review to appear in Gas Accretion onto Galaxies, Astrophysics and Space Science Library, eds. A. J. Fox & R. Davé, to be published by Springer. 29 pages, 2 figures

The Minimum Halo Mass for Star Formation at z = 6 - 8

Authors: K. Finlator, M. K. M. Prescott, B. D. Oppenheimer, R. Davé, E. Zackrisson, R. C. Livermore, S. L. Finkelstein, R. Thompson, S. Huang

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

Submitted 20 September, 2016; originally announced September 2016.

Comments: 8 pages, 5 figures, re-submitted to MNRAS after incorporating referee's comments

The spectral evolution of the first Galaxies. III. Simulated James Webb Space Telescope spectra of reionization-epoch galaxies with Lyman continuum leakage

Authors: E. Zackrisson, C. Binggeli, K. Finlator, N. Y. Gnedin, J. -P. Paardekooper, I. Shimizu, A. K. Inoue, H. Jensen, G. Micheva, S. Khochfar, C. Dalla Vecchia

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

Submitted 12 January, 2017; v1 submitted 29 August, 2016; originally announced August 2016.

Comments: 12 pages, 8 figures, v.2: ApJ, accepted. Model grids are available from http://www.astro.uu.se/~ez/lycan/lycan.html

Equilibrium model prediction for the scatter in the star-forming main sequence

Authors: Sourav Mitra, Romeel Davé, Vimal Simha, Kristian Finlator

Abstract: The analytic "equilibrium model" for galaxy evolution using a mass balance equation is able to reproduce mean observed galaxy scaling relations between stellar mass, halo mass, star formation rate (SFR) and metallicity across the majority of cosmic time with a small number of parameters related to feedback. Here we aim to test this data-constrained model to quantify deviations from the mean relati… ▽ More The analytic "equilibrium model" for galaxy evolution using a mass balance equation is able to reproduce mean observed galaxy scaling relations between stellar mass, halo mass, star formation rate (SFR) and metallicity across the majority of cosmic time with a small number of parameters related to feedback. Here we aim to test this data-constrained model to quantify deviations from the mean relation between stellar mass and SFR, i.e. the star-forming galaxy main sequence (MS). We implement fluctuation in halo accretion rates parameterised from merger-based simulations, and quantify the intrinsic scatter introduced into the MS under the assumption that fluctuations in star formation follow baryonic inflow fluctuations. We predict the 1-sigma MS scatter to be ~ 0.2 - 0.25 dex over the stellar mass range 10^8 Mo to 10^11 Mo and a redshift range 0.5 < z < 3 for SFRs averaged over 100 Myr. The scatter increases modestly at z > 3, as well as by averaging over shorter timescales. The contribution from merger-induced star formation is generally small, around 5% today and 10 - 15% during the peak epoch of cosmic star formation. These results are generally consistent with available observations, suggesting that deviations from the MS primarily reflect stochasticity in the inflow rate owing to halo mergers. △ Less

Submitted 1 October, 2016; v1 submitted 23 June, 2016; originally announced June 2016.

Comments: 11 pages, MNRAS, accepted

The Soft, Fluctuating UVB at $z\sim6$ as Traced by C IV, SiIV, and CII

Authors: K. Finlator, B. D. Oppenheimer, R. Davé, E. Zackrisson, R. Thompson, S. Huang

Abstract: The sources that drove cosmological reionization left clues regarding their identity in the slope and inhomogeneity of the ultraviolet ionizing background (UVB): Bright quasars (QSOs) generate a hard UVB with predominantly large-scale fluctuations while Population II stars generate a softer one with smaller-scale fluctuations. Metal absorbers probe the UVB's slope because different ions are sensit… ▽ More The sources that drove cosmological reionization left clues regarding their identity in the slope and inhomogeneity of the ultraviolet ionizing background (UVB): Bright quasars (QSOs) generate a hard UVB with predominantly large-scale fluctuations while Population II stars generate a softer one with smaller-scale fluctuations. Metal absorbers probe the UVB's slope because different ions are sensitive to different energies. Likewise, they probe spatial fluctuations because they originate in regions where a galaxy-driven UVB is harder and more intense. We take a first step towards studying the reionization-epoch UVB's slope and inhomogeneity by comparing observations of 12 metal absorbers at $z\sim6$ versus predictions from a cosmological hydrodynamic simulation using three different UVBs: a soft, spatially-inhomogeneous "galaxies+QSOs" UVB; a homogeneous "galaxies+QSOs" UVB (Haardt & Madau 2012); and a QSOs-only model. All UVBs reproduce the observed column density distributions of CII, SiIV, and CIV reasonably well although high-column, high-ionization absorbers are underproduced, reflecting numerical limitations. With upper limits treated as detections, only a soft, fluctuating UVB reproduces both the observed SiIV/CIV and CII/CIV distributions. The QSOs-only UVB overpredicts both CIV/CII and CIV/SiIV, indicating that it is too hard. The Haardt & Madau (2012) UVB underpredicts CIV/SiIV, suggesting that it lacks amplifications near galaxies. Hence current observations prefer a soft, fluctuating UVB as expected from a predominantly Population II background although they cannot rule out a harder one. Future observations probing a factor of two deeper in metal column density will distinguish between the soft, fluctuating and QSOs-only UVBs. △ Less

Submitted 21 March, 2016; originally announced March 2016.

Comments: 13 pages, 10 figures, re-submitted to MNRAS after incorporating referee's suggestions

Physical Properties of Spectroscopically-Confirmed Galaxies at $z\ge6$. III. Stellar Populations from SED Modeling with Secure Ly$α$ Emission and Redshifts

Authors: Linhua Jiang, Kristian Finlator, Seth H. Cohen, Eiichi Egami, Rogier A. Windhorst, Xiaohui Fan, Romeel Dave, Nobunari Kashikawa, Matthew Mechtley, Masami Ouchi, Kazuhiro Shimasaku, Benjamin Clement

Abstract: We present a study of stellar populations in a sample of spectroscopically-confirmed Lyman-break galaxies (LBGs) and Ly$α$ emitters (LAEs) at $5.7<z<7$. These galaxies have deep optical and infrared images from Subaru, $HST$, and $Spitzer$/IRAC. We focus on a subset of 27 galaxies with IRAC detections, and characterize their stellar populations utilizing galaxy synthesis models based on the multi-… ▽ More We present a study of stellar populations in a sample of spectroscopically-confirmed Lyman-break galaxies (LBGs) and Ly$α$ emitters (LAEs) at $5.7<z<7$. These galaxies have deep optical and infrared images from Subaru, $HST$, and $Spitzer$/IRAC. We focus on a subset of 27 galaxies with IRAC detections, and characterize their stellar populations utilizing galaxy synthesis models based on the multi-band data and secure redshifts. By incorporating nebular emission estimated from the observed Ly$α$ flux, we are able to break the strong degeneracy of model spectra between young galaxies with prominent nebular emission and older galaxies with strong Balmer breaks. The results show that our galaxies cover a wide range of ages from several to a few hundred million years (Myr), and a wide range of stellar masses from $\sim10^8$ to $\sim10^{11}$ $M_{\odot}$. These galaxies can be roughly divided into an `old' subsample and a `young' subsample. The `old' subsample consists of galaxies older than 100 Myr, with stellar masses higher than $10^9$ $M_{\odot}$. The galaxies in the `young' subsample are younger than $\sim$30 Myr, with masses ranging between $\sim10^8$ and $\sim3\times10^9$ $M_{\odot}$. Both subsamples display a correlation between stellar mass and star-formation rate (SFR), but with very different normalizations. The average specific SFR (sSFR) of the `old' subsample is 3--4 Gyr$^{-1}$, consistent with previous studies of `normal' star-forming galaxies at $z\ge6$. The average sSFR of the `young' subsample is an order of magnitude higher, likely due to starburst activity. Our results also indicate little or no dust extinction in the majority of the galaxies, as already suggested by their steep rest-frame UV slopes. Finally, LAEs and LBGs with strong Ly$α$ emission are indistinguishable in terms of age, stellar mass, and SFR. △ Less

Submitted 4 November, 2015; originally announced November 2015.

Comments: 19 pages, 15 figures, accepted for publication in ApJ

Simulating the 21-cm signal from reionisation including non-linear ionisations and inhomogeneous recombinations

Authors: Sultan Hassan, Romeel Davé, Kristian Finlator, Mario G. Santos

Abstract: We explore the impact of incorporating physically motivated ionisation and recombination rates on the history and topology of cosmic reionisation, by incorporating inputs from small-volume hydrodynamic simulations into a semi-numerical code, SimFast21, that evolves reionisation on large scales. We employ radiative hydrodynamic simulations to parameterize the ionisation rate Rion and recombination… ▽ More We explore the impact of incorporating physically motivated ionisation and recombination rates on the history and topology of cosmic reionisation, by incorporating inputs from small-volume hydrodynamic simulations into a semi-numerical code, SimFast21, that evolves reionisation on large scales. We employ radiative hydrodynamic simulations to parameterize the ionisation rate Rion and recombination rate Rrec as functions of halo mass, overdensity and redshift. We find that Rion is super-linearly dependent on halo mass (Rion ~ Mh^1.41), in contrast to previous assumptions. We implement these scalings into SimFast21 to identify the ionized regions. We tune our models to be consistent with recent observations of the optical depth, ionizing emissivity, and neutral fraction by the end of reionisation. We require an average photon escape fraction fesc=0.04 within ~ 0.5 cMpc cells, independent of halo mass or redshift, to simultaneously match these data. We present predictions for the 21cm power spectrum, and show that it is converged with respect to simulation volume. We find that introducing superlinearly mass-dependent ionisations increases the duration of reionisation and boosts the small-scale 21cm power by ~ 2-3 at intermediate phases of reionisation. Introducing inhomogeneous recombinations reduces ionised bubble sizes and suppresses large-scale 21cm power by ~ 2-3. Moreover, gas clumping on sub-cell scales has a minimal effect on the 21cm power, indicating that robust predictions do not depend on the behaviour of kpc-scale structures. The superlinear ionisations significantly increase the median halo mass scale for ionising photon output to >10^10 Mo, giving greater hope for detecting most of ionising sources with next-generation facilities. These results highlight the importance of more accurately treating ionising sources and recombinations for modeling reionisation and its 21cm signal. △ Less

Submitted 5 February, 2016; v1 submitted 14 October, 2015; originally announced October 2015.

Comments: 20 pages, 14 figues, accepted to MNRAS

Merging galaxies produce outliers from the Fundamental Metallicity Relation

Authors: Asger Grønnow, Kristian Finlator, Lise Christensen

Abstract: From a large sample of $\approx 170,000$ local SDSS galaxies, we find that the Fundamental Metallicity Relation (FMR) has an overabundance of outliers, compared to what would be expected from a Gaussian distribution of residuals, with significantly lower metallicities than predicted from their stellar mass and star formation rate (SFR). This low-metallicity population has lower stellar masses, bim… ▽ More From a large sample of $\approx 170,000$ local SDSS galaxies, we find that the Fundamental Metallicity Relation (FMR) has an overabundance of outliers, compared to what would be expected from a Gaussian distribution of residuals, with significantly lower metallicities than predicted from their stellar mass and star formation rate (SFR). This low-metallicity population has lower stellar masses, bimodial specific SFRs with enhanced star formation within the aperture and smaller half-light radii than the general sample, and is hence a physically distinct population. We show that they are consistent with being galaxies that are merging or have recently merged with a satellite galaxy. In this scenario, low-metallicity gas flows in from large radii, diluting the metallicity of star-forming regions and enhancing the specific SFR until the inflowing gas is processed and the metallicity has recovered. We introduce a simple model in which mergers with a mass ratio larger than a minimum dilute the central galaxy's metallicity by an amount that is proportional to the stellar mass ratio for a constant time, and show that it provides an excellent fit to the distribution of FMR residuals. We find the dilution time-scale to be $τ=1.568_{-0.027}^{+0.029}$ Gyr, the average metallicity depression caused by a 1:1 merger to be $α=0.2480_{-0.0020}^{+0.0017}$ dex and the minimum mass ratio merger that can be discerned from the intrinsic Gaussian scatter in the FMR to be $ξ_\text{min}=0.2030_{-0.0095}^{+0.0127}$ (these are statistical errors only). From this we derive that the average metallicity depression caused by a merger with mass ratio between 1:5 and 1:1 is 0.114 dex. △ Less

Submitted 1 November, 2015; v1 submitted 1 June, 2015; originally announced June 2015.

Comments: 14 pages, 10 figures, published in MNRAS, updated to be essentially identical to the published version

Journal ref: 2015, MNRAS, 451, 4005

A Critical Assessment of Stellar Mass Measurement Methods

Authors: Bahram Mobasher, Tomas Dahlen, Henry C. Ferguson, Viviana Acquaviva, Guillermo Barro, Steven L. Finkelstein, Adriano Fontana, Ruth Gruetzbauch, Seth Johnson, Yu Lu, Casey J. Papovich, Janine Pforr, Mara Salvato, Rachel S. Somerville, Tommy Wiklind, Stijn Wuyts, Matthew L. N. Ashby, Eric Bell, Christopher J. Conselice, Mark E. Dickinson, Sandra M. Faber, Giovanni Fazio, Kristian Finlator, Audrey Galametz, Eric Gawiser , et al. (12 additional authors not shown)

Abstract: In this paper we perform a comprehensive study of the main sources of random and systematic errors in stellar mass measurement for galaxies using their Spectral Energy Distributions (SEDs). We use mock galaxy catalogs with simulated multi-waveband photometry (from U-band to mid-infrared) and known redshift, stellar mass, age and extinction for individual galaxies. Given different parameters affect… ▽ More In this paper we perform a comprehensive study of the main sources of random and systematic errors in stellar mass measurement for galaxies using their Spectral Energy Distributions (SEDs). We use mock galaxy catalogs with simulated multi-waveband photometry (from U-band to mid-infrared) and known redshift, stellar mass, age and extinction for individual galaxies. Given different parameters affecting stellar mass measurement (photometric S/N ratios, SED fitting errors, systematic effects, the inherent degeneracies and correlated errors), we formulated different simulated galaxy catalogs to quantify these effects individually. We studied the sensitivity of stellar mass estimates to the codes/methods used, population synthesis models, star formation histories, nebular emission line contributions, photometric uncertainties, extinction and age. For each simulated galaxy, the difference between the input stellar masses and those estimated using different simulation catalogs, $Δ\log(M)$, was calculated and used to identify the most fundamental parameters affecting stellar masses. We measured different components of the error budget, with the results listed as follows: (1). no significant bias was found among different codes/methods, with all having comparable scatter; (2). A source of error is found to be due to photometric uncertainties and low resolution in age and extinction grids; (3). The median of stellar masses among different methods provides a stable measure of the mass associated with any given galaxy; (4). The deviations in stellar mass strongly correlate with those in age, with a weaker correlation with extinction; (5). the scatter in the stellar masses due to free parameters are quantified, with the sensitivity of the stellar mass to both the population synthesis codes and inclusion of nebular emission lines studied. △ Less

Submitted 6 May, 2015; originally announced May 2015.

Comments: 33 pages, 20 Figures, Accepted for publication in Astrophysical Journal

The Reionization of Carbon

Authors: K. Finlator, R. Thompson, S. Huang, R. Davé, E. Zackrisson, B. D. Oppenheimer

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

Submitted 8 January, 2015; v1 submitted 15 December, 2014; originally announced December 2014.

Comments: 16 pages, 14 figures, accepted to MNRAS

Equilibrium Model Constraints on Baryon Cycling Across Cosmic Time

Authors: Sourav Mitra, Romeel Davé, Kristian Finlator

Abstract: Galaxies strongly self-regulate their growth via energetic feedback from stars, supernovae, and black holes, but these processes are among the least understood aspects of galaxy formation theory. We present an analytic galaxy evolution model that directly constrains such feedback processes from observed galaxy scaling relations. The equilibrium model, which is broadly valid for star-forming centra… ▽ More Galaxies strongly self-regulate their growth via energetic feedback from stars, supernovae, and black holes, but these processes are among the least understood aspects of galaxy formation theory. We present an analytic galaxy evolution model that directly constrains such feedback processes from observed galaxy scaling relations. The equilibrium model, which is broadly valid for star-forming central galaxies that dominate cosmic star formation, is based on the ansatz that galaxies live in a slowly-evolving equilibrium between inflows, outflows, and star formation. Using a Bayesian Monte Carlo Markov chain approach, we constrain our model to match observed galaxy scaling relations between stellar mass and halo mass, star formation rate, and metallicity from 0<z<2. A good fit (chi^2~1.6) is achieved with eight free parameters. We further show that constraining our model to any two of the three data sets also produces a fit to the third that is within reasonable systematic uncertainties. The resulting best-fit parameters that describe baryon cycling suggest galactic outflow scalings intermediate between energy and momentum-driven winds, a weak dependence of wind recycling time on mass, and a quenching mass scale that evolves modestly upwards with redshift. This model further predicts a stellar mass-star formation rate relation that is in good agreement with observations to z~6. Our results suggest that this simple analytic framework captures the basic physical processes required to model the mean evolution of stars and metals in galaxies, despite not incorporating many canonical ingredients of galaxy formation models such as merging or disk formation. △ Less

Submitted 25 June, 2015; v1 submitted 5 November, 2014; originally announced November 2014.

Comments: 19 pages, MNRAS, accepted

The Relation Between SFR and Stellar Mass for Galaxies at 3.5 $\le z\le$ 6.5 in CANDELS

Authors: Brett Salmon, Casey Papovich, Steven L. Finkelstein, Vithal Tilvi, Kristian Finlator, Peter Behroozi, Tomas Dahlen, Romeel Davé, Avishai Dekel, Mark Dickinson, Henry C. Ferguson, Mauro Giavalisco, James Long, Yu Lu, Bahram Mobasher, Naveen Reddy, Rachel S. Somerville, Risa H. Wechsler

Abstract: Distant star-forming galaxies show a correlation between their star formation rates (SFR) and stellar masses, and this has deep implications for galaxy formation. Here, we present a study on the evolution of the slope and scatter of the SFR-stellar mass relation for galaxies at $3.5\leq z\leq 6.5$ using multi-wavelength photometry in GOODS-S from the Cosmic Assembly Near-infrared Deep Extragalacti… ▽ More Distant star-forming galaxies show a correlation between their star formation rates (SFR) and stellar masses, and this has deep implications for galaxy formation. Here, we present a study on the evolution of the slope and scatter of the SFR-stellar mass relation for galaxies at $3.5\leq z\leq 6.5$ using multi-wavelength photometry in GOODS-S from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and Spitzer Extended Deep Survey. We describe an updated, Bayesian spectral-energy distribution fitting method that incorporates effects of nebular line emission, star formation histories that are constant or rising with time, and different dust attenuation prescriptions (starburst and Small Magellanic Cloud). From $z$=6.5 to $z$=3.5 star-forming galaxies in CANDELS follow a nearly unevolving correlation between stellar mass and SFR that follows SFR $\sim$ $M_\star^a$ with $a = 0.54 \pm 0.16$ at $z\sim 6$ and $0.70 \pm 0.21$ at $z\sim 4$. This evolution requires a star formation history that increases with decreasing redshift (on average, the SFRs of individual galaxies rise with time). The observed scatter in the SFR-stellar mass relation is tight, $σ(\log \mathrm{SFR}/\mathrm{M}_\odot$ yr$^{-1})< 0.3\ - $ 0.4 dex, for galaxies with $\log M_\star/\mathrm{M}_\odot > 9$ dex. Assuming that the SFR is tied to the net gas inflow rate (SFR $\sim$ $\dot{M}_\mathrm{gas}$), then the scatter in the gas inflow rate is also smaller than 0.3$-$0.4 dex for star-forming galaxies in these stellar mass and redshift ranges, at least when averaged over the timescale of star formation. We further show that the implied star formation history of objects selected on the basis of their co-moving number densities is consistent with the evolution in the SFR-stellar mass relation. △ Less

Submitted 13 January, 2015; v1 submitted 22 July, 2014; originally announced July 2014.

Comments: 31 pages, 24 figures, accepted for publication in ApJ

A Comparison between Semi-Analytic Model Predictions for the CANDELS Survey

Authors: Yu Lu, Risa H. Wechsler, Rachel S. Somerville, Darren Croton, Lauren Porter, Joel Primack, Peter S. Behroozi, Henry C. Ferguson, David C. Koo, Yicheng Guo, Mohammadtaher Safarzadeh, Kristian Finlator, Marco Castellano, Catherine E. White, Veronica Sommariva, Chris Moody

Abstract: We compare the predictions of three independently developed semi-analytic galaxy formation models that are being used to aid in the interpretation of results from the CANDELS survey. These models are each applied to the same set of halo merger trees extracted from the "Bolshoi" simulation and are carefully tuned to match the local galaxy stellar mass function using the powerful method of Bayesian… ▽ More We compare the predictions of three independently developed semi-analytic galaxy formation models that are being used to aid in the interpretation of results from the CANDELS survey. These models are each applied to the same set of halo merger trees extracted from the "Bolshoi" simulation and are carefully tuned to match the local galaxy stellar mass function using the powerful method of Bayesian Inference coupled with MCMC or by hand. The comparisons reveal that in spite of the significantly different parameterizations for star formation and feedback processes, the three models yield qualitatively similar predictions for the assembly histories of galaxy stellar mass and star formation over cosmic time. We show that the SAMs generally require strong outflows to suppress star formation in low-mass halos to match the present day stellar mass function. However, all of the models considered produce predictions for the star formation rates and metallicities of low-mass galaxies that are inconsistent with existing data and diverge between the models. We suggest that large differences in the metallicity relations and small differences in the stellar mass assembly histories of model galaxies stem from different assumptions for the outflow mass-loading factor. Importantly, while more accurate observational measurements for stellar mass, SFR and metallicity of galaxies at 1<z<5 will discriminate between models, the discrepancies between the models and existing data of these observables have already revealed challenging problems in understanding star formation and its feedback in galaxy formation. The three sets of models are being used to construct catalogs of mock galaxies on light cones that have the same geometry as the CANDELS survey, which should be particularly useful for quantifying the biases and uncertainties on measurements and inferences from the real observations. -ABRIDGED △ Less

Submitted 11 December, 2013; originally announced December 2013.

Comments: 29 pages, 14 figures, 1 table, submitted to ApJ

A Critical Assessment of Photometric Redshift Methods: A CANDELS Investigation

Authors: Tomas Dahlen, Bahram Mobasher, Sandra M. Faber, Henry C. Ferguson, Guillermo Barro, Steven L. Finkelstein, Kristian Finlator, Adriano Fontana, Ruth Gruetzbauch, Seth Johnson, Janine Pforr, Mara Salvato, Tommy Wiklind, Stijn Wuyts, Viviana Acquaviva, Mark E. Dickinson, Yicheng Guo, Jiasheng Huang, Kuang-Han Huang, Jeffrey A. Newman, Eric F. Bell, Christopher J. Conselice, Audrey Galametz, Eric Gawiser, Mauro Giavalisco , et al. (13 additional authors not shown)

Abstract: We present results from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) photometric redshift methods investigation. In this investigation, the results from eleven participants, each using a different combination of photometric redshift code, template spectral energy distributions (SEDs) and priors, are used to examine the properties of photometric redshifts applied to… ▽ More We present results from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) photometric redshift methods investigation. In this investigation, the results from eleven participants, each using a different combination of photometric redshift code, template spectral energy distributions (SEDs) and priors, are used to examine the properties of photometric redshifts applied to deep fields with broad-band multi-wavelength coverage. The photometry used includes U-band through mid-infrared filters and was derived using the TFIT method. Comparing the results, we find that there is no particular code or set of template SEDs that results in significantly better photometric redshifts compared to others. However, we find codes producing the lowest scatter and outlier fraction utilize a training sample to optimize photometric redshifts by adding zero-point offsets, template adjusting or adding extra smoothing errors. These results therefore stress the importance of the training procedure. We find a strong dependence of the photometric redshift accuracy on the signal-to-noise ratio of the photometry. On the other hand, we find a weak dependence of the photometric redshift scatter with redshift and galaxy color. We find that most photometric redshift codes quote redshift errors (e.g., 68% confidence intervals) that are too small compared to that expected from the spectroscopic control sample. We find that all codes show a statistically significant bias in the photometric redshifts. However, the bias is in all cases smaller than the scatter, the latter therefore dominates the errors. Finally, we find that combining results from multiple codes significantly decreases the photometric redshift scatter and outlier fraction. We discuss different ways of combining data to produce accurate photometric redshifts and error estimates. △ Less

Submitted 24 August, 2013; originally announced August 2013.

Comments: 26 pages, 15 figures, ApJ in press

The Host Halos of OI Absorbers in the Reionization Epoch

Authors: K. Finlator, J. A. Muñoz, B. D. Oppenheimer, S. Peng Oh, F. Özel, R. Davé

Abstract: We use a radiation hydrodynamic simulation of the hydrogen reionization epoch to study OI absorbers at z~6. The intergalactic medium (IGM) is reionized before it is enriched, hence OI absorption originates within dark matter halos. The predicted abundance of OI absorbers is in reasonable agreement with observations. At z=10, roughly 70% of sightlines through atomically-cooled halos encounter a vis… ▽ More We use a radiation hydrodynamic simulation of the hydrogen reionization epoch to study OI absorbers at z~6. The intergalactic medium (IGM) is reionized before it is enriched, hence OI absorption originates within dark matter halos. The predicted abundance of OI absorbers is in reasonable agreement with observations. At z=10, roughly 70% of sightlines through atomically-cooled halos encounter a visible (N_OI > 10^14 cm^-2) column. Reionization ionizes and removes gas from halos less massive than 10^8.4 M_0, but 20% of sightlines through more massive halos encounter visible columns even at z=5. The mass scale of absorber host halos is 10-100 times smaller than the halos of Lyman break galaxies and Lyman-alpha emitters, hence absorption probes the dominant ionizing sources more directly. OI absorbers have neutral hydrogen columns of 10^19-10^21 cm^-2, suggesting a close resemblance between objects selected in OI and HI absorption. Finally, the absorption in the foreground of the z=7.085 quasar ULASJ1120+0641 cannot originate in a dark matter halo because halo gas at the observed HI column density is enriched enough to violate the upper limits on the OI column. By contrast, gas at less than one third the cosmic mean density satisfies the constraints. Hence the foreground absorption likely originates in the IGM. △ Less

Submitted 30 August, 2013; v1 submitted 21 May, 2013; originally announced May 2013.

Comments: 21 pages, 15 figures, re-submitted to MNRAS after many improvements suggested by referee, results unchanged

The metallicity evolution of low mass galaxies: New constraints at intermediate redshift

Authors: Alaina Henry, Crystal L. Martin, Kristian Finlator, Alan Dressler

Abstract: We present abundance measurements from 26 emission-line selected galaxies at z~0.6-0.7. By reaching stellar masses as low as 10^8 M_{\sun}, these observations provide the first measurement of the intermediate redshift mass-metallicity (MZ) relation below 10^9 M_{\sun} For the portion of our sample above M > 10^9 M_{\sun} (8/26 galaxies), we find good agreement with previous measurements of the int… ▽ More We present abundance measurements from 26 emission-line selected galaxies at z~0.6-0.7. By reaching stellar masses as low as 10^8 M_{\sun}, these observations provide the first measurement of the intermediate redshift mass-metallicity (MZ) relation below 10^9 M_{\sun} For the portion of our sample above M > 10^9 M_{\sun} (8/26 galaxies), we find good agreement with previous measurements of the intermediate redshift MZ relation. Compared to the local relation, we measure an evolution that corresponds to a 0.12 dex decrease in oxygen abundances at intermediate redshifts. This result confirms the trend that metallicity evolution becomes more significant towards lower stellar masses, in keeping with a downsizing scenario where low mass galaxies evolve onto the local MZ relation at later cosmic times. We show that these galaxies follow the local fundamental metallicity relation, where objects with higher specific (mass-normalized) star formation rates (SFRs) have lower metallicities. Furthermore, we show that the galaxies in our sample lie on an extrapolation of the SFR-M_{*} relation (the star-forming main sequence). Leveraging the MZ relation and star-forming main sequence (and combining our data with higher mass measurements from the literature), we test models that assume an equilibrium between mass inflow, outflow and star formation. We find that outflows are required to describe the data. By comparing different outflow prescriptions, we show that momentum driven winds can describe the MZ relation; however, this model under-predicts the amount of star formation in low mass galaxies. This disagreement may indicate that preventive feedback from gas-heating has been overestimated, or it may signify a more fundamental deviation from the equilibrium assumption. △ Less

Submitted 15 April, 2013; originally announced April 2013.

Comments: Accepted for publication in ApJ

Physical Properties of Spectroscopically-Confirmed Galaxies at z >= 6. II. Morphology of the Rest-Frame UV Continuum and Lyman-alpha Emission

Authors: Linhua Jiang, Eiichi Egami, Xiaohui Fan, Rogier A. Windhorst, Seth H. Cohen, Romeel Dave, Kristian Finlator, Nobunari Kashikawa, Matthew Mechtley, Masami Ouchi, Kazuhiro Shimasaku

Abstract: We present a detailed structural and morphological study of a large sample of spectroscopically-confirmed galaxies at z >= 6, using deep HST near-IR broad-band images and Subaru optical narrow-band images. The galaxy sample consists of 51 Lyman-alpha emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, and 16 Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. These galaxies exhibit a wide range of rest-frame UV c… ▽ More We present a detailed structural and morphological study of a large sample of spectroscopically-confirmed galaxies at z >= 6, using deep HST near-IR broad-band images and Subaru optical narrow-band images. The galaxy sample consists of 51 Lyman-alpha emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, and 16 Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. These galaxies exhibit a wide range of rest-frame UV continuum morphology in the HST images, from compact features to multiple component systems. The fraction of merging/interacting galaxies reaches 40% ~ 50% at the brightest end of M_1500 <= -20.5 mag. The intrinsic half-light radii r_{hl,in}, after correction for PSF broadening, are roughly between r_{hl,in} ~ 0.05" (0.3 kpc) and 0.3" (1.7 kpc) at M_1500 <= -19.5 mag. The median r_{hl,in} value is 0.16" (~0.9 kpc). This is consistent with the sizes of bright LAEs and LBGs at z >= 6 in previous studies. In addition, more luminous galaxies tend to have larger sizes, exhibiting a weak size-luminosity relation r_{hl,in} \propto L^{0.14} at M_1500 <= -19.5 mag. The slope of 0.14 is significantly flatter than those in fainter LBG samples. We discuss the morphology of z >= 6 galaxies with nonparametric methods, including the CAS system and the Gini and M_20 parameters, and demonstrate their validity through simulations. We search for extended Lyman-alpha emission halos around LAEs at z ~ 5.7 and 6.5, by stacking a number of narrow-band images. We do not find evidence of extended halos predicted by cosmological simulations. Such Lyman-alpha halos, if they exist, could be weaker than predicted. Finally, we investigate any positional misalignment between UV continuum and Lyman-alpha emission in LAEs. While the two positions are generally consistent, several merging galaxies show significant positional differences. This is likely caused by a disturbed ISM distribution due to merging activity. △ Less

Submitted 5 July, 2013; v1 submitted 28 February, 2013; originally announced March 2013.

Comments: Accepted for publication in ApJ (22 pages, 13 figures, revised after referee report, galaxy simulations added)

Physical Properties of Spectroscopically-Confirmed Galaxies at z >= 6. I. Basic Characteristics of the Rest-Frame UV Continuum and Lyman-alpha emission

Authors: Linhua Jiang, Eiichi Egami, Matthew Mechtley, Xiaohui Fan, Seth H. Cohen, Rogier A. Windhorst, Romeel Dave, Kristian Finlator, Nobunari Kashikawa, Masami Ouchi, Kazuhiro Shimasaku

Abstract: We present deep HST near-IR and Spitzer mid-IR observations of a large sample of spectroscopically-confirmed galaxies at z >= 6. The sample consists of 51 Lyman-alpha emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, and 16 Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. The near-IR images were mostly obtained with WFC3 in the F125W and F160W bands, and the mid-IR images were obtained with IRAC in the 3.6um… ▽ More We present deep HST near-IR and Spitzer mid-IR observations of a large sample of spectroscopically-confirmed galaxies at z >= 6. The sample consists of 51 Lyman-alpha emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, and 16 Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. The near-IR images were mostly obtained with WFC3 in the F125W and F160W bands, and the mid-IR images were obtained with IRAC in the 3.6um and 4.5um bands. Our galaxies also have deep optical imaging data from Subaru Suprime-Cam. We utilize the multi-band data and secure redshifts to derive their rest-frame UV properties. These galaxies have steep UV continuum slopes roughly between beta ~ -1.5 and -3.5, with an average value of beta ~ -2.3, slightly steeper than the slopes of LBGs in previous studies. The slope shows little dependence on UV continuum luminosity except for a few of the brightest galaxies. We find a statistically significant excess of galaxies with slopes around beta ~ -3, suggesting the existence of very young stellar populations with extremely low metallicity and dust content. Our galaxies have moderately strong rest-frame Lyman-alpha equivalent width (EW) in a range of ~10 to ~200 Å. The star-formation rates are also moderate, from a few to a few tens solar masses per year. The LAEs and LBGs in this sample share many common properties, implying that LAEs represent a subset of LBGs with strong Lyman-alpha emission. Finally, the comparison of the UV luminosity functions between LAEs and LBGs suggests that there exists a substantial population of faint galaxies with weak Lyman-alpha emission (EW < 20 Å) that could be the dominant contribution to the total ionizing flux at z >= 6. △ Less

Submitted 3 June, 2013; v1 submitted 28 February, 2013; originally announced March 2013.

Comments: ApJ Accepted (31 pages, 18 figures)