skip to main content


Title: AMBER: A Semi-numerical Abundance Matching Box for the Epoch of Reionization
Abstract

The Abundance Matching Box for the Epoch of Reionization (AMBER) is a semi-numerical code for modeling the cosmic dawn. The new algorithm is not based on the excursion set formalism for reionization, but takes the novel approach of calculating the reionization-redshift fieldzre(x)assuming that hydrogen gas encountering higher radiation intensity are photoionized earlier. Redshift values are assigned while matching the abundance of ionized mass according to a given mass-weighted ionization fractionx¯i(z). The code has the unique advantage of allowing users to directly specify the reionization history through the redshift midpointzmid, duration Δz, and asymmetryAzinput parameters. The reionization process is further controlled through the minimum halo massMminfor galaxy formation and the radiation mean free pathlmfpfor radiative transfer. We implement improved methods for constructing density, velocity, halo, and radiation fields, which are essential components for modeling reionization observables. We compare AMBER with two other semi-numerical methods and find that our code more accurately reproduces the results from radiation-hydrodynamic simulations. The parallelized code is over four orders of magnitude faster than radiative transfer simulations and will efficiently enable large-volume models, full-sky mock observations, and parameter-space studies. AMBER will be made publicly available to facilitate and transform studies of the Epoch of Reionization.

 
more » « less
Award ID(s):
2007390
NSF-PAR ID:
10364085
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
927
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 186
Size(s):
["Article No. 186"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Using the novel semi-numerical code for reionization AMBER, we model the patchy kinetic Sunyaev–Zel’dovich (kSZ) effect by directly specifying the reionization history with the redshift midpointzmid, duration Δz, and asymmetryAz. We further control the ionizing sources and radiation through the minimum halo massMhand the radiation mean free pathλmfp. AMBER reproduces the free-electron number density and the patchy kSZ power spectrum of radiation–hydrodynamic simulations at the target resolution (1 Mpch−1) with matched reionization parameters. With a suite of (2 Gpc/h)3simulations using AMBER, we first constrain the redshift midpoint 6.0 <zmid< 8.9 using the Planck 2018 Thomson optical depth result (95% CL). Then, assumingzmid= 8, we find that the amplitude ofD=3000pkSZscales linearly with the duration of reionization Δzand is consistent with the 1σupper limit from South Pole Telescope (SPT) results up to Δz< 5.1 (Δzencloses 5%–95% ionization). Moreover, a shorterλmfpcan lead to a ∼10% lowerD=3000pkSZand a flatter slope in theD=3000pkSZΔzscaling relation, thereby affecting the constraints on Δzat= 3000. Allowingzmidandλmfpto vary simultaneously, we get spectra consistent with the SPT result (95% CL) up to Δz= 12.8 (butAz> 8 is needed to ensure the end of reionization beforez= 5.5). We show that constraints on the asymmetry require ∼0.1μk2measurement accuracy at multipoles other than= 3000. Finally, we find that the amplitude and shape of the kSZ spectrum are only weakly sensitive toMhunder a fixed reionization history and radiation mean free path.

     
    more » « less
  2. Abstract

    Cosmic reionization was the last major phase transition of hydrogen from neutral to highly ionized in the intergalactic medium (IGM). Current observations show that the IGM is significantly neutral atz> 7 and largely ionized byz∼ 5.5. However, most methods to measure the IGM neutral fraction are highly model dependent and are limited to when the volume-averaged neutral fraction of the IGM is either relatively low (x¯HI103) or close to unity (x¯HI1). In particular, the neutral fraction evolution of the IGM at the critical redshift range ofz= 6–7 is poorly constrained. We present new constraints onx¯HIatz∼ 5.1–6.8 by analyzing deep optical spectra of 53 quasars at 5.73 <z< 7.09. We derive model-independent upper limits on the neutral hydrogen fraction based on the fraction of “dark” pixels identified in the Lyαand Lyβforests, without any assumptions on the IGM model or the intrinsic shape of the quasar continuum. They are the first model-independent constraints on the IGM neutral hydrogen fraction atz∼ 6.2–6.8 using quasar absorption measurements. Our results give upper limits ofx¯HI(z=6.3)<0.79±0.04(1σ),x¯HI(z=6.5)<0.87±0.03(1σ), andx¯HI(z=6.7)<0.940.09+0.06(1σ). The dark pixel fractions atz> 6.1 are consistent with the redshift evolution of the neutral fraction of the IGM derived from Planck 2018.

     
    more » « less
  3. Abstract

    The repeating fast radio burst FRB 20190520B is an anomaly of the FRB population thanks to its high dispersion measure (DM = 1205 pc cm−3) despite its low redshift ofzfrb= 0.241. This excess has been attributed to a large host contribution of DMhost≈ 900 pc cm−3, far larger than any other known FRB. In this paper, we describe spectroscopic observations of the FRB 20190520B field obtained as part of the FLIMFLAM survey, which yielded 701 galaxy redshifts in the field. We find multiple foreground galaxy groups and clusters, for which we then estimated halo masses by comparing their richness with numerical simulations. We discover two separateMhalo> 1014Mgalaxy clusters atz= 0.1867 and 0.2170 that are directly intersected by the FRB sight line within their characteristic halo radiusr200. Subtracting off their estimated DM contributions, as well that of the diffuse intergalactic medium, we estimate a host contribution ofDMhost=430220+140or280170+140pccm3(observed frame), depending on whether we assume that the halo gas extends tor200or 2 ×r200. This significantly smaller DMhost—no longer the largest known value—is now consistent with Hαemission measures of the host galaxy without invoking unusually high gas temperatures. Combined with the observed FRB scattering timescale, we estimate the turbulent fluctuation and geometric amplification factor of the scattering layer to beF˜G4.511(pc2km)1/3, suggesting that most of the gas is close to the FRB host. This result illustrates the importance of incorporating foreground data for FRB analyses both for understanding the nature of FRBs and to realize their potential as a cosmological probe.

     
    more » « less
  4. Abstract

    We present a multiwavelength analysis of the galaxy cluster SPT-CL J0607-4448 (SPT0607), which is one of the most distant clusters discovered by the South Pole Telescope atz= 1.4010 ± 0.0028. The high-redshift cluster shows clear signs of being relaxed with well-regulated feedback from the active galactic nucleus (AGN) in the brightest cluster galaxy (BCG). Using Chandra X-ray data, we construct thermodynamic profiles and determine the properties of the intracluster medium. The cool-core nature of the cluster is supported by a centrally peaked density profile and low central entropy (K0=189+11keV cm2), which we estimate assuming an isothermal temperature profile due to the limited spectral information given the distance to the cluster. Using the density profile and gas cooling time inferred from the X-ray data, we find a mass-cooling rateṀcool=10060+90Myr−1. From optical spectroscopy and photometry around the [Oii] emission line, we estimate that the BCG star formation rate isSFR[OII]=1.70.6+1.0Myr−1, roughly two orders of magnitude lower than the predicted mass-cooling rate. In addition, using ATCA radio data at 2.1 GHz, we measure a radio jet powerPcav=3.21.3+2.1×1044erg s−1, which is consistent with the X-ray cooling luminosity (Lcool=1.90.5+0.2×1044erg s−1withinrcool= 43 kpc). These findings suggest that SPT0607 is a relaxed, cool-core cluster with AGN-regulated cooling at an epoch shortly after cluster formation, implying that the balance between cooling and feedback can be reached quickly. We discuss the implications for these findings on the evolution of AGN feedback in galaxy clusters.

     
    more » « less
  5. Abstract

    Recently, the Hydrogen Epoch of Reionization Array (HERA) has produced the experiment’s first upper limits on the power spectrum of 21 cm fluctuations atz∼ 8 and 10. Here, we use several independent theoretical models to infer constraints on the intergalactic medium (IGM) and galaxies during the epoch of reionization from these limits. We find that the IGM must have been heated above the adiabatic-cooling threshold byz∼ 8, independent of uncertainties about IGM ionization and the radio background. Combining HERA limits with complementary observations constrains the spin temperature of thez∼ 8 neutral IGM to 27 KT¯S630 K (2.3 KT¯S640 K) at 68% (95%) confidence. They therefore also place a lower bound on X-ray heating, a previously unconstrained aspects of early galaxies. For example, if the cosmic microwave background dominates thez∼ 8 radio background, the new HERA limits imply that the first galaxies produced X-rays more efficiently than local ones. Thez∼ 10 limits require even earlier heating if dark-matter interactions cool the hydrogen gas. If an extra radio background is produced by galaxies, we rule out (at 95% confidence) the combination of high radio and low X-ray luminosities ofLr,ν/SFR > 4 × 1024W Hz−1M1yr andLX/SFR < 7.6 × 1039erg s−1M1yr. The new HERA upper limits neither support nor disfavor a cosmological interpretation of the recent Experiment to Detect the Global EOR Signature (EDGES) measurement. The framework described here provides a foundation for the interpretation of future HERA results.

     
    more » « less