Researcher profile

Simeon Bird

Simeon Bird contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

16 published item(s)

preprint2024arXiv

PRIYA: A New Suite of Lyman-alpha Forest Simulations for Cosmology

We present the PRIYA suite of cosmological simulations, based on the code and hydrodynamic model of the ASTRID simulation, and designed for cosmological analyses of the Lyman-$α$ forest. Our simulation suite spans a $9$-dimensional parameter space, including $4$ cosmological parameters and $5$ astrophysical/thermal parameters. We have run $48$ low fidelity simulations with $1536^3$ particles in a $120$ Mpc/h box and $3$ high fidelity simulations with $3072^3$ particles in a $120$ Mpc/h box. All our simulations include a full physics model for galaxy formation, including supernova and AGN feedback, and thus also contain a realistic population of DLAs. We advance on earlier simulations suites by larger particle loads, by incorporating new physical models for patchy hydrogen and helium reionization, and by self-consistently incorporating a model for AGN feedback. We show that patchy helium reionization imprints an excess in the 1D flux power spectrum on large scales, which may allow future measurements of helium reionization bubble sizes. Simulation parameters are chosen based on a Latin hypercube design and a Gaussian process is used to interpolate to arbitrary parameter combinations. We build a multi-fidelity emulator for the 1D flux power spectrum and the mean IGM temperature. We show that our final interpolation error is $< 1\%$ and that our simulations produce a flux power spectrum converged at the percent level for $z=5.4$ - $2.2$. Our simulation suite will be used to interpret Lyman-$α$ forest 1D flux power spectra from SDSS and future DESI data releases.

preprint2022arXiv

A population of ultraviolet-dim protoclusters detected in absorption

Galaxy protoclusters, which will eventually grow into the massive clusters we see in the local universe, are usually traced by locating overdensities of galaxies. Large spectroscopic surveys of distant galaxies now exist, but their sensitivity depends mainly on a galaxy&#39;s star formation activity and dust content rather than its mass. Tracers of massive protoclusters that do not rely on their galaxy constituents are therefore needed. Here we report observations of Lyman-$α$ absorption in the spectra of a dense grid of background galaxies, which we use to locate a substantial number of candidate protoclusters at redshifts 2.2-2.8 via their intergalactic gas. We find that the structures producing the most absorption, most of which were previously unknown, contain surprisingly few galaxies compared to the dark matter content of their analogs in cosmological simulations. Nearly all are expected to be protoclusters, and we infer that half of their expected galaxy members are missing from our survey because they are unusually dim at rest-frame ultraviolet wavelengths. We attribute this to an unexpectedly strong and early influence of the protocluster environment on the evolution of these galaxies that reduced their star formation or increased their dust content.

preprint2022arXiv

Characterizing Protoclusters and Protogroups at z $\sim$ 2.5 Using Ly$α$ Tomography

Ly-$α$ tomography surveys have begun to produce three-dimensional (3D) maps of the intergalactic medium (IGM) opacity at $z \sim 2.5$ with Mpc resolution. These surveys provide an exciting new way to discover and characterize high-redshift overdensities, including the progenitors of today&#39;s massive groups and clusters of galaxies, known as protogroups and protoclusters. We use the IllustrisTNG-$300$ hydrodynamical simulation to build mock maps that realistically mimic those observed in the Ly-$α$ Tomographic IMACS Survey (LATIS). We introduce a novel method for delineating the boundaries of structures detected in 3D Ly-$α$ flux maps by applying the watershed algorithm. We provide estimators for the dark matter masses of these structures (at $z\sim2.5$), their descendant halo masses at $z=0$, and the corresponding uncertainties. We also investigate the completeness of this method for the detection of protogroups and protoclusters. Compared to earlier work, we apply and characterize our method over a wider mass range that extends to massive protogroups. We also assess the widely used fluctuating Gunn-Peterson approximation (FGPA) applied to dark-matter-only simulations; we conclude while it is adequate for estimating the Ly-$α$ absorption signal from moderate-to-massive protoclusters ($ \gtrsim 10^{14.2} M_{\odot}$), it artificially merges a minority of lower-mass structures with more massive neighbors. Our methods will be applied to current and future Ly-$α$~tomography surveys to create catalogs of overdensities and study environment-dependent galaxy evolution in the Cosmic~Noon~era.

preprint2022arXiv

Improved selection of extremely red quasars with boxy CIV lines in BOSS

Extremely red quasars (ERQs) are an interesting sample of quasars in the Baryon Oscillation Spectroscopic Sample (BOSS) in the redshift range of $2.0 - 3.4$ and have extreme red colours of $i-W3\ge4.6$. Core ERQs have strong CIV emission lines with rest equivalent width of $\ge100$Å. Many core ERQs also have CIV line profiles with peculiar boxy shapes which distinguish them from normal blue quasars. We show, using a combination of kernel density estimation and local outlier factor analyses on a space of the $i-W3$ colour, CIV rest equivalent width and line kurtosis, that core ERQs likely represent a separate population rather than a smooth transition between normal blue quasars and the quasars in the tail of the colour-REW distribution. We apply our analyses to find new criteria for selecting ERQs in this 3D parameter space. Our final selection produces $133$ quasars, which are \emph{three} times more likely to have a visually verified CIV broad absorption line feature than the previous core ERQ sample. We further show that our newly selected sample are extreme objects in the intersection of the WISE AGN catalogue with the MILLIQUAS quasar catalogue in the colour-colour space of ($W1-W2$, $W2-W3$). This paper validates an improved selection method for red quasars which can be applied to future datasets such as the quasar catalogue from the Dark Energy Spectroscopic Instrument (DESI).

preprint2022arXiv

Snowmass2021: Vera C. Rubin Observatory as a Flagship Dark Matter Experiment

Establishing that Vera C. Rubin Observatory is a flagship dark matter experiment is an essential pathway toward understanding the physical nature of dark matter. In the past two decades, wide-field astronomical surveys and terrestrial laboratories have jointly created a phase transition in the ecosystem of dark matter models and probes. Going forward, any robust understanding of dark matter requires astronomical observations, which still provide the only empirical evidence for dark matter to date. We have a unique opportunity right now to create a dark matter experiment with Rubin Observatory Legacy Survey of Space and Time (LSST). This experiment will be a coordinated effort to perform dark matter research, and provide a large collaborative team of scientists with the necessary organizational and funding supports. This approach leverages existing investments in Rubin. Studies of dark matter with Rubin LSST will also guide the design of, and confirm the results from, other dark matter experiments. Supporting a collaborative team to carry out a dark matter experiment with Rubin LSST is the key to achieving the dark matter science goals that have already been identified as high priority by the high-energy physics and astronomy communities.

preprint2022arXiv

The ASTRID Simulation: Galaxy Formation and Reionization

We introduce the ASTRID simulation, a large-scale cosmological hydrodynamic simulation in a $250$ Mpc/h box with $2\times 5500^3$ particles. ASTRID contains a large number of high redshift galaxies, which can be compared to future survey data, and resolves galaxies in halos more massive than $2\times 10^9 M_\odot$. ASTRID has been run from $z=99$ to $z=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 summarise the implemented models, and the technical choices we took when developing the simulation code. We validate the model, showing good agreement with observed UV luminosity functions, galaxy stellar mass functions and specific star formation rates. 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$, halos with $M \sim 2\times 10^9 M_\odot$ which have been reionized have a star formation rate $1.5$ times greater than those which have not yet been reionized.

preprint2022arXiv

The ASTRID simulation: the evolution of Supermassive Black Holes

We present the evolution of black holes (BHs) and their relationship with their host galaxies in Astrid, a large-volume cosmological hydrodynamical simulation with box size 250 $h^{-1} \rm Mpc$ containing $2\times5500^3$ particles evolved to z=3. Astrid statistically models BH gas accretion and AGN feedback to their environments, applies a power-law distribution for BH seed mass $M_{\rm sd}$, uses a dynamical friction model for BH dynamics and executes a physical treatment of BH mergers. The BH population is broadly consistent with empirical constraints on the BH mass function, the bright end of the luminosity functions, and the time evolution of BH mass and accretion rate density. The BH mass and accretion exhibit a tight correlation with host stellar mass and star formation rate. We trace BHs seeded before z>10 down to z=3, finding that BHs carry virtually no imprint of the initial $M_{\rm sd}$ except those with the smallest $M_{\rm sd}$, where less than 50\% of them have doubled in mass. Gas accretion is the dominant channel for BH growth compared to BH mergers. With dynamical friction, Astrid predicts a significant delay for BH mergers after the first encounter of a BH pair, with a typical elapse time of about 200 Myrs. There are in total $4.5 \times 10^5$ BH mergers in Astrid at z>3, $\sim 10^3$ of which have X-ray detectable EM counterparts: a bright kpc scale dual AGN with $L_X>10^{43}$ erg/s. BHs with $M_{\rm BH} \sim 10^{7-8} M_{\odot}$ experience the most frequent mergers. Galaxies that host BH mergers are unbiased tracers of the overall $M_{\rm BH} - M_{*}$ relation. Massive ($>10^{11} M_{\odot}$) galaxies have a high occupation number (>10) of BHs, and hence host the majority of BH mergers.

preprint2022arXiv

The low redshift Lyman-$α$ Forest as a constraint for models of AGN feedback

We study the sensitivity of the $z=0.1$ Lyman-$α$ Forest observables, such as the column density distribution function (CDD), flux PDF, flux power spectrum, and line width distribution, to sub-grid models of active galactic nuclei (AGN) feedback using the Illustris and IllustrisTNG (TNG) cosmological simulations. The two simulations share an identical Ultraviolet Background (UVB) prescription and similar cosmological parameters, but TNG features an entirely reworked AGN feedback model. Due to changes in the AGN radio mode model, the original Illustris simulations have a factor of 2-3 fewer Lyman-$α$ absorbers than TNG at column densities $N_{\rm HI}< 10^{15.5}$ cm$^{-2}$. We compare the simulated forest statistics to UV data from the Cosmic Origins Spectrograph (COS) and find that neither simulation can reproduce the slope of the absorber distribution. Both Illustris and TNG also produce significantly smaller line width distributions than observed in the COS data. We show that TNG is in much better agreement with the observed $z=0.1$ flux power spectrum than Illustris. We explore which statistics can disentangle the effects of AGN feedback from alternative UVB models by rescaling the UVB of Illustris to produce a CDD match to TNG. While this UVB rescaling is degenerate with the effect of AGN feedback on the CDD, the amplitude and shape of the flux PDF and 1D flux power spectrum change in a way distinct from a scaling of the UVB. Our study suggests that the $z=0.1$ Lyman-$α$ forest observables can be used as a diagnostic of AGN feedback models.

preprint2021arXiv

Densified pupil spectrograph as high-precision radial velocimetry: From direct measurement of the Universe&#39;s expansion history to characterization of nearby habitable planet candidates

The direct measurement of the Universe&#39;s expansion history and the search for terrestrial planets in habitable zones around solar-type stars require extremely high-precision radial velocity measures over a decade. This study proposes an approach for enabling high-precision radial velocity measurements from space. The concept presents a combination of a high-dispersion densified pupil spectrograph and a novel telescope line-of-sight monitor. The precision of the radial velocity measurements is determined by combining the spectrophotometric accuracy and the quality of the absorption lines in the recorded spectrum. Therefore, a highly dispersive densified pupil spectrograph proposed to perform stable spectroscopy can be utilized for high-precision radial velocity measures. A concept involving the telescope line-of-sight monitor is developed to minimize the change of the telescope line-of-sight over a decade. This monitor allows the precise measurement of a long-term telescope drift without any significant impact on the Airy disk when the densified pupil spectra are recorded. We analytically derive the uncertainty of the radial velocity measurements, which is caused by the residual offset of the line-of-sights at two epochs. We find that the error could be reduced down to approximately 1 $cm/s$, and the precision will be limited by another factor (e.g., wavelength calibration uncertainty). A combination of the high precision spectrophotometry and the high spectral resolving power could open a new path toward the characterization of nearby non-transiting habitable planet candidates orbiting late-type stars. We present two simple and compact high-dispersed densified pupil spectrograph designs for the cosmology and exoplanet sciences.

preprint2021arXiv

Massive Black Hole Mergers with Orbital Information: Predictions from the ASTRID Simulation

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 a MBH seed population between $3\times 10^4M_{\odot}/h$ and $3\times 10^5M_{\odot}/h$ and a sub-grid dynamical friction (DF) model to follow the MBH dynamics down to $1.5\;\text{ckpc}/h$. We calculate 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{Myrs}$) due to DF. The binary hardening time is an order of magnitude longer than the DF time, especially for the seed-mass binaries ($M_\text{BH}<2M_\text{seed}$). As a result, only $\lesssim20\%$ 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 $>10^9\,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 $\sim 7$ higher than the prediction using the circular orbit assumption. The LISA events are expected at a similar rate, and comprise $\gtrsim 60\%$ seed-seed mergers, $\sim 30\%$ involving only one seed-mass MBH, and $\sim 10\%$ mergers of non-seed MBHs.

preprint2020arXiv

Automated Measurement of Quasar Redshift with a Gaussian Process

We develop an automated technique to measure quasar redshifts in the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey (SDSS). Our technique is an extension of an earlier Gaussian process method for detecting damped Lyman-alpha absorbers (DLAs) in quasar spectra with known redshifts. We apply this technique to a subsample of SDSS DR12 with BAL quasars removed and redshift larger than 2.15. We show that we are broadly competitive to existing quasar redshift estimators, disagreeing with the PCA redshift by more than 0.5 in only 0.38% of spectra. Our method produces a probabilistic density function for the quasar redshift, allowing quasar redshift uncertainty to be propagated to downstream users. We apply this method to detecting DLAs, accounting in a Bayesian fashion for redshift uncertainty. Compared to our earlier method with a known quasar redshift, we have a moderate decrease in our ability to detect DLAs, predominantly in the noisiest spectra. The area under curve drops from 0.96 to 0.91. Our code is publicly available.

preprint2020arXiv

Cosmic Filaments from Cosmic Strings

Cosmic strings are generically predicted in many extensions of the Standard Model of particle physics. We propose a new avenue for detecting cosmic strings through their effect on the filamentary structure in the cosmic web. Using cosmological simulations of the density wake from a cosmic string, we examine a variety of filament structure probes. We show that the largest effect of the cosmic string is an overdensity in the filament distribution around the string wake. The signal from the overdensity is stronger at higher redshift, and more robust with a wider field. We analyze the spatial distribution of filaments from a publicly available catalog of filaments built from SDSS galaxies. With existing data, we find no evidence for the presence of a cosmic string wake with string tension parameter $Gμ$ above $5\times 10^{-6}$. However, we project WFIRST will be able to detect a signal from such a wake at the $99\%$ confidence level at redshift $z=2$, with significantly higher confidence and the possibility of probing lower tensions ($Gμ\sim 10^{-6}$), at $z=10$. The sensitivity of this method is not competitive with constraints derived from the CMB. However, it provides an independent discovery channel at low redshift, which could be a smoking-gun in scenarios where the CMB bound can be weakened.

preprint2020arXiv

Detecting Multiple DLAs per Spectrum in SDSS DR12 with Gaussian Processes

We present a revised version of our automated technique using Gaussian processes (GPs) to detect Damped Lyman-$α$ absorbers (DLAs) along quasar (QSO) sightlines. The main improvement is to allow our Gaussian process pipeline to detect multiple DLAs along a single sightline. Our DLA detections are regularised by an improved model for the absorption from the Lyman-$α$ forest which improves performance at high redshift. We also introduce a model for unresolved sub-DLAs which reduces mis-classifications of absorbers without detectable damping wings. We compare our results to those of two different large-scale DLA catalogues and provide a catalogue of the processed results of our Gaussian process pipeline using 158 825 Lyman-$α$ spectra from SDSS data release 12. We present updated estimates for the statistical properties of DLAs, including the column density distribution function (CDDF), line density ($dN/dX$), and neutral hydrogen density ($Ω_{\textrm{DLA}}$).

preprint2020arXiv

Inhomogeneous HeII reionization in Hydrodynamic Simulations

The reionization of the second electron of helium shapes the physical state of intergalactic gas at redshifts between 2 < z < 5. Because performing full in situ radiative transfer in hydrodynamic simulations is computationally expensive for large volumes, the physics of HeII reionization is often approximated by a uniform UV background model that does not capture the patchiness of reionization. We have devised a model that implements the effects of HeII reionization using semi-analytic calculations of the thermal state of intergalactic gas-- a way to bypass a full radiative transfer simulation while still realizing the physics of HeII reionization that affects observables such as the Lyman alpha forest. Here we present a publicly-available code that flexibly models inhomogeneous HeII reionization in simulations at a negligible computational cost. Because many of the parameters of HeII reionization are uncertain, our model is customizable from a set of free parameters. We show results from this code in MP-Gadget, where this model is implemented. We demonstrate the resulting temperature evolution and temperature-density relation of intergalactic gas-- consistent with recent measurements and previous radiative transfer simulations. We show that the impact of HeII reionization gives rise to subtle signatures in the one-dimensional statistics of the Lyman alpha forest at the level of several per cent, in agreement with previous findings. The flexible nature of these simulations are ideal for studies of HeII reionization and future observations of the HeII Lyman alpha forest.

preprint2020arXiv

Massive neutrinos and degeneracies in Lyman-alpha forest simulations

Using a suite of hydrodynamical simulations with cold dark matter, baryons, and neutrinos, we present a detailed study of the effect of massive neutrinos on the 1-D and 3-D flux power spectra of the Lyman-$α$ (Ly$α$) forest. The presence of massive neutrinos in cosmology induces a scale- and time-dependent suppression of structure formation that is strongest on small scales. Measuring this suppression is a key method for inferring neutrino masses from cosmological data, and is one of the main goals of ongoing and future surveys like eBOSS, DES, LSST, Euclid or DESI. The clustering in the Ly$α$ forest traces the quasi-linear power at late times and on small scales. In combination with observations of the cosmic microwave background, the forest therefore provides some of the tightest constraints on the sum of the neutrino masses. However there is a well-known degeneracy between $Σm_ν$ and the amplitude of perturbations in the linear matter power spectrum. We study the corresponding degeneracy in the 1-D flux power spectrum of the Ly$α$ forest, and for the first time also study this degeneracy in the 3-D flux power spectrum. We show that the non-linear effects of massive neutrinos on the Ly$α$ forest, beyond the effect of linear power amplitude suppression, are negligible, and this degeneracy persists in the Ly$α$ forest observables to a high precision. We discuss the implications of this degeneracy for choosing parametrisations of the Ly$α$ forest for cosmological analysis.

preprint2020arXiv

More accurate simulations with separate initial conditions for baryons and dark matter

We revisit techniques for performing cosmological simulations with both baryons and cold dark matter when each fluid has different initial conditions, as is the case at the end of the radiation era. Most simulations do not reproduce the linear prediction for the difference between the cold dark matter and baryon perturbations. We show that this is due to the common use of offset regular grids when setting up the particle initial conditions. The correct behaviour can be obtained without any loss of simulation resolution by using a Lagrangian glass for the baryon particles. We further show that the difference between cold dark matter and baryons may affect predictions for the Lyman-alpha forest flux power spectrum at the 5% level, potentially impacting current cosmological constraints.