Search for: All records

Abstract We present an analysis searching for dual active galactic nuclei (AGN) among 62 high-redshift (2.5 <z< 3.5) X-ray sources selected from the X-UDS, AEGIS-XD, CDF-S, and COSMOS-Legacy Chandra surveys. We aim to quantify the frequency of dual AGN in the high-redshift Universe, which holds implications for black hole merger timescales and low-frequency gravitational wave detection rates. We analyze each X-ray source using BAYMAX, an analysis tool that calculates the Bayes factor for whether a given archival Chandra AGN is more likely a single or dual point source. We find no strong evidence for dual AGN in any individual source in our sample. We increase our sensitivity to search for dual AGN across the sample by comparing our measured distribution of Bayes factors to that expected from a sample composed entirely of single point sources and find no evidence for dual AGN in the sample distribution. Although our analysis utilizes one of the largest Chandra catalogs of high-zX-ray point sources available to study, the findings remain limited by the modest number of sources observed at the highest spatial resolution with Chandra and the typical count rates of the detected sources. Our nondetection allows us to place an upper limit on the X-ray dual AGN fraction at 2.5 <z< 3.5 of 4.8% at the 95% confidence level. Expanding substantially on these results at X-ray wavelengths will require future surveys spanning larger sky areas and extending to fainter fluxes than has been possible with Chandra. We illustrate the potential of the AXIS mission concept in this regard. 

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 midpointz_mid, duration Δ_z, and asymmetryA_z. We further control the ionizing sources and radiation through the minimum halo massM_hand 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⁻¹) with matched reionization parameters. With a suite of (2 Gpc/h)³simulations using AMBER, we first constrain the redshift midpoint 6.0 <z_mid< 8.9 using the Planck 2018 Thomson optical depth result (95% CL). Then, assumingz_mid= 8, we find that the amplitude of $D_{\mathit{\ell }=3000}^{\mathrm{pkSZ}}$scales 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% lower $D_{\mathit{\ell }=3000}^{\mathrm{pkSZ}}$and a flatter slope in the $D_{\mathit{\ell }=3000}^{\mathrm{pkSZ}}-{\mathrm{\Delta }}_z$scaling relation, thereby affecting the constraints on Δ_zatℓ= 3000. Allowingz_midandλ_mfpto vary simultaneously, we get spectra consistent with the SPT result (95% CL) up to Δ_z= 12.8 (butA_z> 8 is needed to ensure the end of reionization beforez= 5.5). We show that constraints on the asymmetry require ∼0.1μk²measurement accuracy at multipoles other thanℓ= 3000. Finally, we find that the amplitude and shape of the kSZ spectrum are only weakly sensitive toM_hunder a fixed reionization history and radiation mean free path. 

Abstract We present CAMELS-ASTRID, the third suite of hydrodynamical simulations in the Cosmology and Astrophysics with MachinE Learning (CAMELS) project, along with new simulation sets that extend the model parameter space based on the previous frameworks of CAMELS-TNG and CAMELS-SIMBA, to provide broader training sets and testing grounds for machine-learning algorithms designed for cosmological studies. CAMELS-ASTRID employs the galaxy formation model following the ASTRID simulation and contains 2124 hydrodynamic simulation runs that vary three cosmological parameters (Ω_m,σ₈, Ω_b) and four parameters controlling stellar and active galactic nucleus (AGN) feedback. Compared to the existing TNG and SIMBA simulation suites in CAMELS, the fiducial model of ASTRID features the mildest AGN feedback and predicts the least baryonic effect on the matter power spectrum. The training set of ASTRID covers a broader variation in the galaxy populations and the baryonic impact on the matter power spectrum compared to its TNG and SIMBA counterparts, which can make machine-learning models trained on the ASTRID suite exhibit better extrapolation performance when tested on other hydrodynamic simulation sets. We also introduce extension simulation sets in CAMELS that widely explore 28 parameters in the TNG and SIMBA models, demonstrating the enormity of the overall galaxy formation model parameter space and the complex nonlinear interplay between cosmology and astrophysical processes. With the new simulation suites, we show that building robust machine-learning models favors training and testing on the largest possible diversity of galaxy formation models. We also demonstrate that it is possible to train accurate neural networks to infer cosmological parameters using the high-dimensional TNG-SB28 simulation set. 

ABSTRACT We introduce the Astrid  simulation, a large-scale cosmological hydrodynamic simulation in a $$250 \, h^{-1}\mathrm{Mpc}$$ box with 2 × 55003 particles. Astrid contains a large number of high redshift galaxies, which can be compared to future survey data, and resolves galaxies in haloes more massive than $$2\times 10^9 \, \mathrm{M}_{\odot }$$. Astrid  has been run from z = 99 to 3. As a particular focus is modelling the high redshift Universe, it contains models for inhomogeneous hydrogen and helium reionization, baryon relative velocities and massive neutrinos, as well as supernova and AGN feedback. The black hole model includes mergers driven by dynamical friction rather than repositioning. We briefly summarize the implemented models, and the technical choices we took when developing the simulation code. We validate the model, showing good agreement with observed ultraviolet luminosity functions, galaxy stellar mass functions and specific star formation rates (SFRs). We show that the redshift at which a given galaxy underwent hydrogen reionization has a large effect on the halo gas fraction. Finally, at z = 6, haloes with $$M \sim 2\times 10^9 \, \mathrm{M}_{\odot }$$ which have been reionized have an SFR 1.5 times greater than those which have not yet been reionized. 

Abstract The statistics of galactic-scale quasar pairs can elucidate our understanding of the dynamical evolution of supermassive black hole (SMBH) pairs, the duty cycles of quasar activity in mergers, or even the nature of dark matter, but they have been challenging to measure at cosmic noon, the prime epoch of massive galaxy and SMBH formation. Here we measure a double quasar fraction of ∼6.2 ± 0.5 × 10⁻⁴integrated over ∼0.″3–3″ separations (projected physical separations of ∼3–30 kpc atz∼ 2) in luminous (L_bol> 10^45.8erg s⁻¹) unobscured quasars at 1.5 <z< 3.5 using Gaia EDR3-resolved pairs around SDSS DR16 quasars. The measurement was based on a sample of 60 Gaia-resolved double quasars (out of 487 Gaia pairs dominated by quasar+star superpositions) at these separations, corrected for pair completeness in Gaia, which we quantify as functions of pair separation, magnitude of the primary, and magnitude contrast. The double quasar fraction increases toward smaller separations by a factor of ∼5 over these scales. The division between physical quasar pairs and lensed quasars in our sample is currently unknown, requiring dedicated follow-up observations (in particular, deep, subarcsecond-resolution IR imaging for the closest pairs). Intriguingly, at this point, the observed pair statistics are in rough agreement with theoretical predictions both for the lensed quasar population in mock catalogs and for dual quasars in cosmological hydrodynamic simulations. Upcoming wide-field imaging/spectroscopic space missions such as Euclid, CSST, and Roman, combined with targeted follow-up observations, will conclusively measure the abundances and host galaxy properties of galactic-scale quasar pairs, offset AGNs, and subarcsecond lensed quasars across cosmic time. 

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 field $z_{\mathrm{re}}\left(\mathit{x}\right)$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 fraction ${\overline{x}}_{\mathrm{i}}\left(z\right)$. The code has the unique advantage of allowing users to directly specify the reionization history through the redshift midpointz_mid, duration Δ_z, and asymmetryA_zinput parameters. The reionization process is further controlled through the minimum halo massM_minfor galaxy formation and the radiation mean free pathl_mfpfor 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. 

ABSTRACT We examine massive black hole (MBH) mergers and their associated gravitational wave signals from the large-volume cosmological simulation Astrid . Astrid includes galaxy formation and black hole models recently updated with an MBH seed population between 3 × 104h−1M⊙ and 3 × 105h−1M⊙ and a sub-grid dynamical friction (DF) model to follow the MBH dynamics down to 1.5 ckpc h−1. We calculate the initial eccentricities of MBH orbits directly from the simulation at kpc-scales, and find orbital eccentricities above 0.7 for most MBH pairs before the numerical merger. After approximating unresolved evolution on scales below $${\sim 200\, \text{pc}}$$, we find that the in-simulation DF on large scales accounts for more than half of the total orbital decay time ($$\sim 500\, \text{Myr}$$) due to DF. The binary hardening time is an order of magnitude longer than the DF time, especially for the seed-mass binaries (MBH < 2Mseed). As a result, only $$\lesssim 20{{\rm per \,cent}}$$ of seed MBH pairs merge at z > 3 after considering both unresolved DF evolution and binary hardening. These z > 3 seed-mass mergers are hosted in a biased population of galaxies with the highest stellar masses of $$\gt 10^9\, {\rm M}_\odot$$. With the higher initial eccentricity prediction from Astrid , we estimate an expected merger rate of 0.3−0.7 per year from the z > 3 MBH population. This is a factor of ∼7 higher than the prediction using the circular orbit assumption. The Laser Interferometer Space Antenna events are expected at a similar rate, and comprise $$\gtrsim 60\,{\rm{per\,cent}}$$ seed-seed mergers, $$\sim 30\,{\rm{per\,cent}}$$ involving only one seed-mass MBH, and $$\sim 10\,{\rm{per\,cent}}$$ mergers of non-seed MBHs.