BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Simeon Bird

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

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAgr-qchep-phastro-ph.HEastro-ph.IMastro-ph.EPhep-exhep-thphysics.data-an

Catalog footprint

What is connected

38works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

38 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.

0 citations0 reviewsNo code linkedOpen access
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'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.

0 citations0 reviewsNo code linkedOpen access
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'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.

0 citations0 reviewsNo code linkedOpen access
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).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

The direct measurement of the Universe'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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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}}$).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cross-correlation between thermal Sunyaev-Zeldovich effect and the integrated Sachs-Wolfe effect

Large-angle fluctuations in the cosmic microwave background temperature induced by the integrated Sachs-Wolfe effect and Compton-y distortions from the thermal Sunyaev-Zeldovich (tSZ) effect are both due to line-of-sight density perturbations. Here we calculate the cross-correlation between these two signals. Measurement of this cross-correlation can be used to test the redshift distribution of the tSZ distortion, which has implications for the redshift at which astrophysical processes in clusters begin to operate. We also evaluate the detectability of a yT cross-correlation from exotic early-Universe sources in the presence of this late-time effect.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Determining the progenitors of merging black-hole binaries

We investigate a possible method for determining the progenitors of black hole (BH) mergers observed via their gravitational wave (GW) signal. We argue that measurements of the cross-correlation of the GW events with overlapping galaxy catalogs may provide an additional tool in determining if BH mergers trace the stellar mass of the Universe, as would be expected from mergers of the endpoints of stellar evolution. If on the other hand the BHs are of primordial origin, as has been recently suggested, their merging would be preferentially hosted by lower biased objects, and thus have a lower cross-correlation with luminous galaxies. Here we forecast the expected precision of the cross-correlation measurement for current and future GW detectors such as LIGO and the Einstein Telescope. We then predict how well these instruments can distinguish the model that identifies high-mass BH-BH mergers as the merger of primordial black holes that constitute the dark matter in the Universe from more traditional astrophysical sources.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Did LIGO detect dark matter?

We consider the possibility that the black-hole (BH) binary detected by LIGO may be a signature of dark matter. Interestingly enough, there remains a window for masses $20\,M_\odot \lesssim M_{\rm bh} \lesssim 100\, M_\odot$ where primordial black holes (PBHs) may constitute the dark matter. If two BHs in a galactic halo pass sufficiently close, they radiate enough energy in gravitational waves to become gravitationally bound. The bound BHs will rapidly spiral inward due to emission of gravitational radiation and ultimately merge. Uncertainties in the rate for such events arise from our imprecise knowledge of the phase-space structure of galactic halos on the smallest scales. Still, reasonable estimates span a range that overlaps the $2-53$ Gpc$^{-3}$ yr$^{-1}$ rate estimated from GW150914, thus raising the possibility that LIGO has detected PBH dark matter. PBH mergers are likely to be distributed spatially more like dark matter than luminous matter and have no optical nor neutrino counterparts. They may be distinguished from mergers of BHs from more traditional astrophysical sources through the observed mass spectrum, their high ellipticities, or their stochastic gravitational wave background. Next generation experiments will be invaluable in performing these tests.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Orbital eccentricities in primordial black holes binaries

It was recently suggested that the merger of $\sim30\,M_\odot$ primordial black holes (PBHs) may provide a significant number of events in gravitational-wave observatories over the next decade, if they make up an appreciable fraction of the dark matter. Here we show that measurement of the eccentricities of the inspiralling binary black holes can be used to distinguish these binaries from those produced by more traditional astrophysical mechanisms. These PBH binaries are formed on highly eccentric orbits and can then merge on timescales that in some cases are years or less, retaining some eccentricity in the last seconds before the merger. This is to be contrasted with massive-stellar-binary, globular-cluster, or other astrophysical origins for binary black holes (BBHs) in which the orbits have very effectively circularized by the time the BBH enters the observable LIGO window. Here we discuss the features of the gravitational-wave signals that indicate this eccentricity and forecast the sensitivity of LIGO and the Einstein Telescope to such effects. We show that if PBHs make up the dark matter, then roughly one event should have a detectable eccentricity given LIGO's expected sensitivity and observing time of six years. The Einstein Telescope should see $O(10)$ such events after ten years.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Simulating the Carbon Footprint of Galactic Halos

We compare simulations, including the Illustris simulations, to observations of CIV and CII absorption at z=2-4. These are the CIV column density distribution function in the column density range $10^{12} - 10^{15}$ cm$^{-2}$, the CIV equivalent width distribution at 0.1 - 2 Å, and the covering fractions and equivalent widths of CIV 1548 and CII 1337 around DLAs. In the context of the feedback models we investigate, all CIV observations favour the use of more energetic wind models, which are better able to enrich the gas surrounding halos. We propose two ways to achieve this; an increased wind velocity and an increase in wind thermal energy. However, even our most energetic wind models do not produce enough absorbers with CIV equivalent width > 0.6 Å, which in our simulations are associated with the most massive haloes. All simulations are in reasonable agreement with the CII covering fraction and equivalent widths around Damped Lyman-$α$ absorbers, although there is a moderate deficit in one bin 10 - 100 kpc from the DLA. Finally, we show that the CIV in our simulations is predominantly photoionized.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stochastic Gravitational-Wave Background due to Primordial Binary Black Hole Mergers

Recent Advanced LIGO detections of binary black hole mergers have prompted multiple studies investigating the possibility that the heavy GW150914 binary system was of primordial origin, and hence could be evidence for dark matter in the form of black holes. We compute the stochastic background arising from the incoherent superposition of such primordial binary black hole systems in the universe and compare it to the similar background spectrum due to binary black hole systems of stellar origin. We investigate the possibility of detecting this background with future gravitational wave detectors, and discuss the possibility of using the stochastic gravitational-wave background measurement to constrain the dark matter component in the form of black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Testing the effect of galactic feedback on the IGM at z ~ 6 with metal-line absorbers

We present models of low- and high-ionization metal-line absorbers (O I, C II, C IV and Mg II) during the end of the reionization epoch, at z ~ 6. Using four cosmological hydrodynamical simulations with different feedback schemes (including the Illustris and Sherwood simulations) and two different choices of hydro-solver, we investigate how the overall incidence rate and equivalent width distribution of metal-line absorbers varies with the galactic wind prescription. We find that the O I and C II absorbers are reasonably insensitive to the feedback scheme. All models, however, struggle to reproduce the observations of C IV and Mg II, which are probing down to lower overdensities than O I and C II at z ~ 6, suggesting that the metals in the simulations are not being transported out into the IGM efficiently enough. The situation is improved but not resolved if we choose a harder (but still reasonable) and/or (locally) increased UV background at z ~ 6.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Reproducing the Kinematics of Damped Lyman-alpha Systems

We examine the kinematic structure of Damped Lyman-alpha Systems (DLAs) in a series of cosmological hydrodynamic simulations using the AREPO code. We are able to match the distribution of velocity widths of associated low ionisation metal absorbers substantially better than earlier work. Our simulations produce a population of DLAs dominated by halos with virial velocities around 70 km/s, consistent with a picture of relatively small, faint objects. In addition, we reproduce the observed correlation between velocity width and metallicity and the equivalent width distribution of SiII. Some discrepancies of moderate statistical significance remain; too many of our spectra show absorption concentrated at the edge of the profile and there are slight differences in the exact shape of the velocity width distribution. We show that the improvement over previous work is mostly due to our strong feedback from star formation and our detailed modelling of the metal ionisation state.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The BlueTides Simulation: First Galaxies and Reionization

We introduce the BlueTides simulation and report initial results for the luminosity functions of the first galaxies and AGN, and their contribution to reionization. BlueTides was run on the BlueWaters cluster at NCSA from $z=99$ to $z=8.0$ and includes 2$\times$7040$^3$ particles in a $400$Mpc/h per side box, making it the largest hydrodynamic simulation ever performed at high redshift. BlueTides includes a pressure-entropy formulation of smoothed particle hydrodynamics, gas cooling, star formation (including molecular hydrogen), black hole growth and models for stellar and AGN feedback processes. The star formation rate density in the simulation is a good match to current observational data at $z\sim 8-10$. We find good agreement between observations and the predicted galaxy luminosity function in the currently observable range $-18\le M_{\mathrm UV} \le -22.5$ with some dust extinction required to match the abundance of brighter objects. BlueTides implements a patchy reionization model that produces a fluctuating UV background. BlueTides predicts number counts for galaxies fainter than current observational limits which are consistent with extrapolating the faint end slope of the luminosity function with a power law index $α\sim -1.8$ at $z\sim 8$ and redshift dependence of $α\sim (1+z)^{-0.4}$. The AGN population has a luminosity function well fit by a power law with a slope $α\sim -2.4$ that compares favourably with the deepest CANDELS-Goods fields. We investigate how these luminosity functions affect the progress of reionization, and find that a high Lyman-$α$ escape fraction ($f_\mathrm{esc} \sim 0.5$) is required if galaxies dominate the ionising photon budget during reionization. Smaller galaxy escape fractions imply a large contribution from faint AGN (down to $M_\mathrm{UV}=-12$) which results in a rapid reionization, disfavoured by current observations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The formation of Milky Way-mass disk galaxies in the first 500 million years of a cold dark matter universe

Whether among the myriad tiny proto-galaxies there exists a population with similarities to present day galaxies is an open question. We show, using BlueTides, the first hydrodynamic simulation large enough to resolve the relevant scales, that the first massive galaxies to form are %in fact predicted to have extensive rotationally-supported disks. Although their morphology resembles in some ways Milky-way types seen at much lower redshifts, these high-redshift galaxies are smaller, denser, and richer in gas than their low redshift counterparts. From a kinematic analysis of a statistical sample of 216 galaxies at redshift $z=8-10$ we have found that disk galaxies make up 70\% of the population of galaxies with stellar mass $10^{10} M_\odot$ or greater. Cold Dark Matter cosmology therefore makes specific predictions for the population of large galaxies 500 million years after the Big Bang. We argue that wide-field satellite telescopes (e.g. WFIRST) will in the near future discover these first massive disk galaxies. The simplicity of their structure and formation history should make possible new tests of cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Impact of Galactic Feedback on the Circumgalactic Medium

Galactic feedback strongly affects the way galactic environments are enriched. We examine this connection by performing a suite of cosmological hydrodynamic simulations, exploring a range of parameters based on the galaxy formation model developed in Vogelsberger et al. 2013 (henceforth V13). We examine the effects of AGN feedback, wind mass loading, wind specific energy, and wind metal-loading on the properties of the circumgalactic medium (CGM) of galaxies with $M_\text{halo} > 10^{11} M_\odot$. Note that while the V13 model was tuned to match observations including the stellar mass function, no explicit tuning was done for the CGM. The wind energy per unit outflow mass has the most significant effect on the CGM enrichment. High energy winds launch metals far beyond the virial radius. AGN feedback also has a significant effect, but only at $z < 3$. We compare to high redshift HI and CIV observations. All our simulations produce the observed number of Damped Lyman-$α$ Absorbers. At lower column density, several of our simulations produce enough Lyman Limit Systems (LLS) $100$ kpc from the galaxy, but in all cases the LLS abundance falls off with distance faster than observations, with too few LLS at $200$ kpc. Further, in all models the CIV abundance drops off too sharply with distance, with too little CIV $100$-$200$ kpc from the galaxy. Higher energy wind models produce more extended CIV but also produce less stars, in tension with star-formation rate density observations. This highlights the fact that circumgalactic observations are a strong constraint on galactic feedback models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Constraints on ionising photon production from the large-scale Lyman-alpha forest

Recent work has shown that the z~2.5 Lyman-alpha forest on large scales encodes information about the galaxy and quasar populations that keep the intergalactic medium photoionized. We present the first forecasts for constraining the populations with data from current and next-generation surveys. At a minimum the forest should tell us whether galaxies or, conversely, quasars dominate the photon production. The number density and clustering strength of the ionising sources might be estimated to sub-10% precision with a DESI-like survey if degeneracies (e.g., with the photon mean-free-path, small-scale clustering power normalization and potentially other astrophysical effects) can be broken by prior information. We demonstrate that, when inhomogeneous ionisation is correctly handled, constraints on dark energy do not degrade.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Damped Lyman-alpha absorbers as a probe of stellar feedback

We examine the abundance, clustering and metallicity of Damped Lyman-alpha Absorbers (DLAs) in a suite of hydrodynamic cosmological simulations using the moving mesh code AREPO. We incorporate models of supernova and AGN feedback, as well as molecular hydrogen formation. We compare our simulations to the column density distribution function at $z=3$, the total DLA abundance at $z=2-4$, the measured DLA bias at $z=2.3$ and the DLA metallicity distribution at $z=2-4$. Our preferred models produce populations of DLAs in good agreement with most of these observations. The exception is the DLA abundance at $z < 3$, which we show requires stronger feedback in $10^{11-12} \, h^{-1} M_\odot$ mass halos. While the DLA population probes a wide range of halo masses, we find the cross-section is dominated by halos of mass $10^{10} - 10^{11} \, h^{-1} M_\odot$ and virial velocities $50 - 100 \;\mathrm{km/s}$. The simulated DLA population has a linear theory bias of $1.7$, whereas the observations require $2.17 \pm 0.2$. We show that non-linear growth increases the bias in our simulations to $2.3$ at $k=1\; \mathrm{Mpc/}h$, the smallest scale observed. The scale-dependence of the bias is, however, very different in the simulations compared against the observations. We show that, of the observations we consider, the DLA abundance and column density function provide the strongest constraints on the feedback model.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Massive Neutrinos and the Non-linear Matter Power Spectrum

We perform an extensive suite of N-body simulations of the matter power spectrum, incorporating massive neutrinos in the range M = 0.15-0.6 eV, probing the non-linear regime at scales k < 10 hMpc-1 at z < 3. We extend the widely used HALOFIT approximation to account for the effect of massive neutrinos on the power spectrum. In the strongly non-linear regime HALOFIT systematically over-predicts the suppression due to the free-streaming of the neutrinos. The maximal discrepancy occurs at k ~ 1 hMpc-1, and is at the level of 10% of the total suppression. Most published constraints on neutrino masses based on HALOFIT are not affected, as they rely on data probing the matter power spectrum in the linear or mildly non-linear regime. However, predictions for future galaxy, Lyman-alpha forest and weak lensing surveys extending to more non-linear scales will benefit from the improved approximation to the non-linear matter power spectrum we provide. Our approximation reproduces the induced neutrino suppression over the targeted scales and redshifts significantly better. We test its robustness with regard to changing cosmological parameters and a variety of modelling effects.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Properties of galaxies reproduced by a hydrodynamic simulation

Previous simulations of the growth of cosmic structures have broadly reproduced the 'cosmic web' of galaxies that we see in the Universe, but failed to create a mixed population of elliptical and spiral galaxies due to numerical inaccuracies and incomplete physical models. Moreover, because of computational constraints, they were unable to track the small scale evolution of gas and stars to the present epoch within a representative portion of the Universe. Here we report a simulation that starts 12 million years after the Big Bang, and traces 13 billion years of cosmic evolution with 12 billion resolution elements in a volume of $(106.5\,{\rm Mpc})^3$. It yields a reasonable population of ellipticals and spirals, reproduces the distribution of galaxies in clusters and statistics of hydrogen on large scales, and at the same time the metal and hydrogen content of galaxies on small scales.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Non-linear evolution of the cosmic neutrino background

We investigate the non-linear evolution of the relic cosmic neutrino background by running large box-size, high resolution N-body simulations. Our set of simulations explore the properties of neutrinos in a reference $Λ$CDM model with total neutrino masses between 0.05-0.60 eV in cold dark matter haloes of mass $10^{11}-10^{15}$ $h^{-1}$M$_{\odot}$, over a redshift range $z=0-2$. We compute the halo mass function and show that it is reasonably well fitted by the Sheth-Tormen formula. More importantly, we focus on the CDM and neutrino properties of the density and peculiar velocity fields in the cosmological volume, inside and in the outskirts of virialized haloes. The dynamical state of the neutrino particles depends strongly on their momentum: whereas neutrinos in the low velocity tail behave similarly to CDM particles, neutrinos in the high velocity tail are not affected by the clustering of the underlying CDM component. We find that the neutrino (linear) unperturbed momentum distribution is modified and mass and redshift dependent deviations from the expected Fermi-Dirac distribution are in place both in the cosmological volume and inside haloes. The neutrino density profiles around virialized haloes have been carefully investigated and a simple fitting formula is provided. The neutrino profile, unlike the cold dark matter one, is found to be cored with core size and central density that depend on the neutrino mass, redshift and mass of the halo, for halos of masses larger than $\sim 10^{13.5}h^{-1}$M$_{\odot}$. For lower masses the neutrino profile is best fitted by a simple power-law relation in the range probed by the simulations. Our findings are particularly important in view of upcoming large-scale structure surveys, like Euclid, that are expected to probe the non-linear regime at the percent level with lensing and clustering observations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An efficient implementation of massive neutrinos in non-linear structure formation simulations

Massive neutrinos make up a fraction of the dark matter, but due to their large thermal velocities, cluster significantly less than cold dark matter (CDM) on small scales. An accurate theoretical modelling of their effect during the non-linear regime of structure formation is required in order to properly analyse current and upcoming high-precision large-scale structure data, and constrain the neutrino mass. Taking advantage of the fact that massive neutrinos remain linearly clustered up to late times, this paper treats the linear growth of neutrino overdensities in a non-linear CDM background. The evolution of the CDM component is obtained via N-body computations. The smooth neutrino component is evaluated from that background by solving the Boltzmann equation linearised with respect to the neutrino overdensity. CDM and neutrinos are simultaneously evolved in time, consistently accounting for their mutual gravitational influence. This method avoids the issue of shot-noise inherent to particle-based neutrino simulations, and, in contrast with standard Fourier-space methods, properly accounts for the non-linear potential wells in which the neutrinos evolve. Inside the most massive late-time clusters, where the escape velocity is larger than the neutrino thermal velocity, neutrinos can clump non-linearly, causing the method to formally break down. It is shown that this does not affect the total matter power spectrum, which can be very accurately computed on all relevant scales up to the present time.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constraints on Massive Neutrinos from the CFHTLS Angular Power Spectrum

We use the galaxy angular power spectrum at $z\sim0.5-1.2$ from the Canada-France-Hawaii-Telescope Legacy Survey Wide fields (CFHTLS-Wide) to constrain separately the total neutrino mass $\sum{m_ν}$ and the effective number of neutrino species $N_{\rm{eff}}$. This survey has recently benefited from an accurate calibration of the redshift distribution, allowing new measurements of the (non-linear) matter power spectrum in a unique range of scales and redshifts sensitive to neutrino free streaming. Our analysis makes use of a recent model for the effect of neutrinos on the weakly non-linear matter power spectrum derived from accurate N-body simulations. We show that CFHTLS, combined with WMAP7 and a prior on the Hubble constant provides an upper limit of $\sum{m_ν}<0.29\,$eV and $N_{\rm{eff}} =4.17^{+1.62}_{-1.26}$ (2$\,σ$ confidence levels). If we omit smaller scales which may be affected by non-linearities, these constraints become $\sum{m_ν}<0.41\,$eV and $N_{\rm{eff}} =3.98^{+2.02}_{-1.20}$ (2$\,σ$ confidence levels). Finally we show that the addition of other large scale structures probes can further improve these constraints, demonstrating that high redshift large volumes surveys such as CFHTLS are complementary to other cosmological probes of the neutrino mass.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Moving mesh cosmology: properties of neutral hydrogen in absorption

We examine the distribution of neutral hydrogen in cosmological simulations carried out with the new moving-mesh code AREPO and compare it with the corresponding GADGET simulations based on the smoothed particle hydrodynamics (SPH) technique. The two codes use identical gravity solvers and baryonic physics implementations, but very different methods for solving the Euler equations, allowing us to assess how numerical effects associated with the hydro-solver impact the results of simulations. Here we focus on an analysis of the neutral gas, as detected in quasar absorption lines. We find that the high column density regime probed by Damped Lyman-alpha (DLA) and Lyman Limit Systems (LLS) exhibits significant differences between the codes. GADGET produces spurious artefacts in large halos in the form of gaseous clumps, boosting the LLS cross-section. Furthermore, it forms halos with denser central baryonic cores than AREPO, which leads to a substantially greater DLA cross-section from smaller halos. AREPO thus produces a significantly lower cumulative abundance of DLAs, which is intriguingly in much closer agreement with observations. The column density function, however, is not altered enough to significantly reduce the discrepancy with the observed value. For the low column density gas probed by the Lyman-alpha forest, the codes differ only at the level of a few percent, suggesting that this regime is quite well described by both methods, a fact that is reassuring for the many Lyman-alpha studies carried out with SPH thus far. While the residual differences are smaller than the errors on current Lyman-alpha forest data, we note that this will likely change for future precision experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Neutrino masses and cosmological parameters from a Euclid-like survey: Markov Chain Monte Carlo forecasts including theoretical errors

We present forecasts for the accuracy of determining the parameters of a minimal cosmological model and the total neutrino mass based on combined mock data for a future Euclid-like galaxy survey and Planck. We consider two different galaxy surveys: a spectroscopic redshift survey and a cosmic shear survey. We make use of the Monte Carlo Markov Chains (MCMC) technique and assume two sets of theoretical errors. The first error is meant to account for uncertainties in the modelling of the effect of neutrinos on the non-linear galaxy power spectrum and we assume this error to be fully correlated in Fourier space. The second error is meant to parametrize the overall residual uncertainties in modelling the non-linear galaxy power spectrum at small scales, and is conservatively assumed to be uncorrelated and to increase with the ratio of a given scale to the scale of non-linearity. It hence increases with wavenumber and decreases with redshift. With these two assumptions for the errors and assuming further conservatively that the uncorrelated error rises above 2% at k = 0.4 h/Mpc and z = 0.5, we find that a future Euclid-like cosmic shear/galaxy survey achieves a 1-sigma error on Mnu close to 32 meV/25 meV, sufficient for detecting the total neutrino mass with good significance. If the residual uncorrelated errors indeed rises rapidly towards smaller scales in the non-linear regime as we have assumed here then the data on non-linear scales does not increase the sensitivity to the total neutrino mass. Assuming instead a ten times smaller theoretical error with the same scale dependence, the error on the total neutrino mass decreases moderately from sigma(Mnu) = 18 meV to 14 meV when mildly non-linear scales with 0.1 h/Mpc < k < 0.6 h/Mpc are included in the analysis of the galaxy survey data.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Minimally Parametric Power Spectrum Reconstruction from the Lyman-alpha Forest

Current results from the Lyman alpha forest assume that the primordial power spectrum of density perturbations follows a simple power law form, with running. We present the first analysis of Lyman alpha data to study the effect of relaxing this strong assumption on primordial and astrophysical constraints. We perform a large suite of numerical simulations, using them to calibrate a minimally parametric framework for describing the power spectrum. Combined with cross-validation, a statistical technique which prevents over-fitting of the data, this framework allows us to reconstruct the power spectrum shape without strong prior assumptions. We find no evidence for deviation from scale-invariance; our analysis also shows that current Lyman alpha data do not have sufficient statistical power to robustly probe the shape of the power spectrum at these scales. In contrast, the ongoing Baryon Oscillation Sky Survey (BOSS) will be able to do so with high precision. Furthermore, this near-future data will be able to break degeneracies between the power spectrum shape and astrophysical parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Minimally Parametric Constraints on the Primordial Power Spectrum from Lyman-alpha

Current analyses of the Lyman-alpha forest assume that the primordial power spectrum of density perturbations obeys a simple power law, a strong theoretical assumption which should be tested. Employing a large suite of numerical simulations which drop this assumption, we reconstruct the shape of the primordial power spectrum using Lyman-alpha data from the Sloan Digital Sky Survey (SDSS). Our method combines a minimally parametric framework with cross-validation, a technique used to avoid over-fitting the data. Future work will involve predictions for the upcoming Baryon Oscillation Sky Survey (BOSS), which will provide new Lyman-alpha data with vastly decreased statistical errors.

0 citations0 reviewsNo code linkedOpen access