scholarly journals Effect of separate initial conditions on the lyman-α forest in simulations

2021 ◽  
Vol 503 (2) ◽  
pp. 1668-1679
Author(s):  
M A Fernandez ◽  
Simeon Bird ◽  
Phoebe Upton Sanderbeck
Keyword(s):  
Power Spectrum ◽  
Transfer Functions ◽  
Thermal State ◽  
Initial Conditions ◽  
High Redshift ◽  
Doppler Width ◽  
Width Distribution ◽  
Cent Level ◽  
Species Specific ◽  
Simulation Parameters

ABSTRACT Using a set of high resolution simulations, we quantify the effect of species-specific initial transfer functions on probes of the intergalactic medium (IGM) via the Lyman-α forest. We focus on redshifts 2–6, after H i reionization. We explore the effect of these initial conditions on measures of the thermal state of the low density IGM: the curvature, Doppler width cutoff, and Doppler width distribution. We also examine the matter and flux power spectrum, and potential consequences for constraints on warm dark matter models. We find that the curvature statistic is at most affected at the $\approx 2{{\ \rm per\ cent}}$ level at z = 6. The Doppler width cutoff parameters are affected by $\approx 5{{\ \rm per\ cent}}$ for the intercept, and $\approx 8{{\ \rm per\ cent}}$ for the fit slope, though this is subdominant to sample variation. The Doppler width distribution shows a $\approx 30{{\ \rm per\ cent}}$ effect at z = 3, however the distribution is not fully converged with simulation box size and resolution. The flux power spectrum is at most affected by $\approx 5{{\ \rm per\ cent}}$ at high redshift and small scales. We discuss numerical convergence with simulation parameters.

scholarly journals Formation of supermassive black hole seeds in nuclear star clusters via gas accretion and runaway collisions

2021 ◽  
Author(s):  
Arpan Das ◽  
Dominik R G Schleicher ◽  
Nathan W C Leigh ◽  
Tjarda C N Boekholt
Keyword(s):  
Black Holes ◽  
Initial Conditions ◽  
High Redshift ◽  
Star Clusters ◽  
Gas Reservoir ◽  
Highly Sensitive ◽  
Stellar Collisions ◽  
Key Variables ◽  
Chaotic Nature ◽  
Gas Accretion

Abstract More than two hundred supermassive black holes (SMBHs) of masses ≳ 109 M⊙ have been discovered at z ≳ 6. One promising pathway for the formation of SMBHs is through the collapse of supermassive stars (SMSs) with masses ∼103 − 5 M⊙ into seed black holes which could grow upto few times 109 M⊙ SMBHs observed at z ∼ 7. In this paper, we explore how SMSs with masses ∼103 − 5 M⊙ could be formed via gas accretion and runaway stellar collisions in high-redshift, metal-poor nuclear star clusters (NSCs) using idealised N-body simulations. We explore physically motivated accretion scenarios, e.g. Bondi-Hoyle-Lyttleton accretion and Eddington accretion, as well as simplified scenarios such as constant accretions. While gas is present, the accretion timescale remains considerably shorter than the timescale for collisions with the most massive object (MMO). However, overall the timescale for collisions between any two stars in the cluster can become comparable or shorter than the accretion timescale, hence collisions still play a crucial role in determining the final mass of the SMSs. We find that the problem is highly sensitive to the initial conditions and our assumed recipe for the accretion, due to the highly chaotic nature of the problem. The key variables that determine the mass growth mechanism are the mass of the MMO and the gas reservoir that is available for the accretion. Depending on different conditions, SMSs of masses ∼103 − 5 M⊙ can form for all three accretion scenarios considered in this work.

scholarly journals The anisotropy of the power spectrum in periodic cosmological simulations

2021 ◽  
Vol 503 (4) ◽  
pp. 5638-5645
Author(s):  
Gábor Rácz ◽  
István Szapudi ◽  
István Csabai ◽  
László Dobos
Keyword(s):  
Dark Matter ◽  
Power Spectrum ◽  
Large Scale ◽  
Cold Dark Matter ◽  
Initial Conditions ◽  
Power Spectra ◽  
Cosmological Simulations ◽  
Spherical Collapse ◽  
Lower Amplitude ◽  
Hydrodynamical Simulations

ABSTRACT The classical gravitational force on a torus is anisotropic and always lower than Newton’s 1/r2 law. We demonstrate the effects of periodicity in dark matter only N-body simulations of spherical collapse and standard Lambda cold dark matter (ΛCDM) initial conditions. Periodic boundary conditions cause an overall negative and anisotropic bias in cosmological simulations of cosmic structure formation. The lower amplitude of power spectra of small periodic simulations is a consequence of the missing large-scale modes and the equally important smaller periodic forces. The effect is most significant when the largest mildly non-linear scales are comparable to the linear size of the simulation box, as often is the case for high-resolution hydrodynamical simulations. Spherical collapse morphs into a shape similar to an octahedron. The anisotropic growth distorts the large-scale ΛCDM dark matter structures. We introduce the direction-dependent power spectrum invariant under the octahedral group of the simulation volume and show that the results break spherical symmetry.

scholarly journals EDGE: a new approach to suppressing numerical diffusion in adaptive mesh simulations of galaxy formation

2020 ◽  
Vol 501 (2) ◽  
pp. 1755-1765
Author(s):  
Andrew Pontzen ◽  
Martin P Rey ◽  
Corentin Cadiou ◽  
Oscar Agertz ◽  
Romain Teyssier ◽  
...  
Keyword(s):  
Star Formation ◽  
Galaxy Formation ◽  
Adaptive Mesh Refinement ◽  
Large Scale ◽  
Initial Conditions ◽  
Dwarf Galaxy ◽  
High Redshift ◽  
Morphological Structure ◽  
Adaptive Mesh ◽  
Numerical Diffusion

ABSTRACT We introduce a new method to mitigate numerical diffusion in adaptive mesh refinement (AMR) simulations of cosmological galaxy formation, and study its impact on a simulated dwarf galaxy as part of the ‘EDGE’ project. The target galaxy has a maximum circular velocity of $21\, \mathrm{km}\, \mathrm{s}^{-1}$ but evolves in a region that is moving at up to $90\, \mathrm{km}\, \mathrm{s}^{-1}$ relative to the hydrodynamic grid. In the absence of any mitigation, diffusion softens the filaments feeding our galaxy. As a result, gas is unphysically held in the circumgalactic medium around the galaxy for $320\, \mathrm{Myr}$, delaying the onset of star formation until cooling and collapse eventually triggers an initial starburst at z = 9. Using genetic modification, we produce ‘velocity-zeroed’ initial conditions in which the grid-relative streaming is strongly suppressed; by design, the change does not significantly modify the large-scale structure or dark matter accretion history. The resulting simulation recovers a more physical, gradual onset of star formation starting at z = 17. While the final stellar masses are nearly consistent ($4.8 \times 10^6\, \mathrm{M}_{\odot }$ and $4.4\times 10^6\, \mathrm{M}_{\odot }$ for unmodified and velocity-zeroed, respectively), the dynamical and morphological structure of the z = 0 dwarf galaxies are markedly different due to the contrasting histories. Our approach to diffusion suppression is suitable for any AMR zoom cosmological galaxy formation simulations, and is especially recommended for those of small galaxies at high redshift.

scholarly journals On the road to per cent accuracy – III. Non-linear reaction of the matter power spectrum to massive neutrinos

2019 ◽  
Vol 491 (3) ◽  
pp. 3101-3107 ◽  
Author(s):  
M Cataneo ◽  
J D Emberson ◽  
D Inman ◽  
J Harnois-Déraps ◽  
C Heymans
Keyword(s):  
Power Spectrum ◽  
Cold Dark Matter ◽  
New Physics ◽  
Model Reaction ◽  
Matter Power Spectrum ◽  
Non Linear ◽  
Cent Level ◽  
Massive Neutrinos ◽  
On The Road ◽  
Total Matter

ABSTRACT We analytically model the non-linear effects induced by massive neutrinos on the total matter power spectrum using the halo model reaction framework of Cataneo et al. In this approach, the halo model is used to determine the relative change to the matter power spectrum caused by new physics beyond the concordance cosmology. Using standard fitting functions for the halo abundance and the halo mass–concentration relation, the total matter power spectrum in the presence of massive neutrinos is predicted to per cent-level accuracy, out to $k=10 \,{ h}\,{\rm Mpc}^{-1}$. We find that refining the prescriptions for the halo properties using N-body simulations improves the recovered accuracy to better than 1 per cent. This paper serves as another demonstration for how the halo model reaction framework, in combination with a single suite of standard Λ cold dark matter (ΛCDM) simulations, can recover per cent-level accurate predictions for beyond ΛCDM matter power spectra, well into the non-linear regime.

scholarly journals Probing the thermal state of the intergalactic medium at z > 5 with the transmission spikes in high-resolution  Ly α forest spectra

2020 ◽  
Vol 494 (4) ◽  
pp. 5091-5109 ◽  
Author(s):  
Prakash Gaikwad ◽  
Michael Rauch ◽  
Martin G Haehnelt ◽  
Ewald Puchwein ◽  
James S Bolton ◽  
...  
Keyword(s):  
High Resolution ◽  
Radiative Transfer ◽  
Intergalactic Medium ◽  
Thermal State ◽  
Neutral Hydrogen ◽  
High Spectral Resolution ◽  
Profile Fitting ◽  
Density Relation ◽  
Width Distribution ◽  
Hydrodynamical Simulations

ABSTRACT We compare a sample of five high-resolution, high S/N  Ly α forest spectra of bright 6 < z < ∼6.5 QSOs aimed at spectrally resolving the last remaining transmission spikes at z > 5 with those obtained from mock absorption spectra from the Sherwoodand Sherwood–Relics simulation suites of hydrodynamical simulations of the intergalactic medium (IGM). We use a profile-fitting procedure for the inverted transmitted flux, 1 − F, similar to the widely used Voigt profile fitting of the transmitted flux F at lower redshifts, to characterize the transmission spikes that probe predominately underdense regions of the IGM. We are able to reproduce the width and height distributions of the transmission spikes, both with optically thin simulations of the post-reionization Universe using a homogeneous UV background and full radiative transfer simulations of a late reionization model. We find that the width of the fitted components of the simulated transmission spikes is very sensitive to the instantaneous temperature of the reionized IGM. The internal structures of the spikes are more prominent in low temperature models of the IGM. The width distribution of the observed transmission spikes, which require high spectral resolution (≤ 8  km s−1) to be resolved, is reproduced for optically thin simulations with a temperature at mean density of T0 = (11 000 ± 1600, 10 500 ± 2100, 12 000 ± 2200) K at z = (5.4, 5.6, 5.8). This is weakly dependent on the slope of the temperature-density relation, which is favoured to be moderately steeper than isothermal. In the inhomogeneous, late reionization, full radiative transfer simulations where islands of neutral hydrogen persist to z ∼ 5.3, the width distribution of the observed transmission spikes is consistent with the range of T0 caused by spatial fluctuations in the temperature–density relation.

scholarly journals ETHOS – an effective parametrization and classification for structure formation: the non-linear regime at z ≳ 5

2020 ◽  
Vol 498 (3) ◽  
pp. 3403-3419
Author(s):  
Sebastian Bohr ◽  
Jesús Zavala ◽  
Francis-Yan Cyr-Racine ◽  
Mark Vogelsberger ◽  
Torsten Bringmann ◽  
...  
Keyword(s):  
Power Spectrum ◽  
Structure Formation ◽  
Parameter Space ◽  
Broad Class ◽  
High Redshift ◽  
Acoustic Oscillations ◽  
Effective Parameters ◽  
Matter Power Spectrum ◽  
Wide Range ◽  
Non Linear

ABSTRACT We propose two effective parameters that fully characterize galactic-scale structure formation at high redshifts (z ≳ 5) for a variety of dark matter (DM) models that have a primordial cutoff in the matter power spectrum. Our description is within the recently proposed ETHOS framework and includes standard thermal warm DM (WDM) and models with dark acoustic oscillations (DAOs). To define and explore this parameter space, we use high-redshift zoom-in simulations that cover a wide range of non-linear scales from those where DM should behave as CDM (k ∼ 10 h Mpc−1), down to those characterized by the onset of galaxy formation (k ∼ 500 h Mpc−1). We show that the two physically motivated parameters hpeak and kpeak, the amplitude and scale of the first DAO peak, respectively, are sufficient to parametrize the linear matter power spectrum and classify the DM models as belonging to effective non-linear structure formation regions. These are defined by their relative departure from cold DM (kpeak → ∞) and WDM (hpeak = 0) according to the non-linear matter power spectrum and halo mass function. We identify a region where the DAOs still leave a distinct signature from WDM down to z = 5, while a large part of the DAO parameter space is shown to be degenerate with WDM. Our framework can then be used to seamlessly connect a broad class of particle DM models to their structure formation properties at high redshift without the need of additional N-body simulations.

scholarly journals A Model for the Sequence of Elliptical Galaxies

1987 ◽  
Vol 117 ◽  
pp. 367-367
Author(s):  
Rosemary F. G. Wyse ◽  
Bernard J. T. Jones
Keyword(s):  
Dark Matter ◽  
Power Spectrum ◽  
Galaxy Formation ◽  
Cold Dark Matter ◽  
Initial Conditions ◽  
Elliptical Galaxy ◽  
Elliptical Galaxies ◽  
Density Fluctuations ◽  
Analytic Approximation ◽  
Non Linear

We present a simple model for the formation of elliptical galaxies, based on a binary clustering hierarchy of dark matter, the chemical enrichment of the gas at each level being controlled by supernovae. The initial conditions for the non-linear phases of galaxy formation are set by the post-recombination power spectrum of density fluctuations. We investigate two models for this power spectrum - the first is a straightforward power law, |δk|2 ∝ kn, and the second is Peeble's analytic approximation to the emergent spectrum in a universe dominated by cold dark matter. The normalisation is chosen such that on some scale, say M ∼ 1012M⊙, the objects that condense out have properties - radius and velocity dispersion - resembling ‘typical’ galaxies. There is some ambiguity in this due to the poorly determined mass-to-light ratio of a typical elliptical galaxy — we look at two normalisations, σ1D ∼ 350kms−1 and σ1D ∼ 140kms−1. The choice determines which of Compton cooling or hydrogen cooling is more important during the galaxy formation period. The non-linear behaviour of the perturbations is treated by the homogeneous sphere approximation.

scholarly journals Comparing approximate methods for mock catalogues and covariance matrices II: power spectrum multipoles

2019 ◽  
Vol 485 (2) ◽  
pp. 2806-2824 ◽  
Author(s):  
Linda Blot ◽  
Martin Crocce ◽  
Emiliano Sefusatti ◽  
Martha Lippich ◽  
Ariel G Sánchez ◽  
...  
Keyword(s):  
Growth Rate ◽  
Distribution Function ◽  
Power Spectrum ◽  
Probability Distribution Function ◽  
Initial Conditions ◽  
Covariance Matrices ◽  
Full Density ◽  
Model Parameters ◽  
Approximate Methods ◽  
Parameter Constraints

ABSTRACT We study the accuracy of several approximate methods for gravitational dynamics in terms of halo power spectrum multipoles and their estimated covariance matrix. We propagate the differences in covariances into parameter constraints related to growth rate of structure, Alcock–Paczynski distortions, and biasing. We consider seven methods in three broad categories: algorithms that solve for halo density evolution deterministically using Lagrangian trajectories (ICE–COLA, pinocchio, and peakpatch), methods that rely on halo assignment schemes on to dark matter overdensities calibrated with a target N-body run (halogen, patchy), and two standard assumptions about the full density probability distribution function (Gaussian and lognormal). We benchmark their performance against a set of three hundred N-body simulations, running similar sets of approximate simulations with matched initial conditions, for each method. We find that most methods reproduce the monopole to within $5{{\ \rm per\ cent}}$, while residuals for the quadrupole are sometimes larger and scale dependent. The variance of the multipoles is typically reproduced within $10{{\ \rm per\ cent}}$. Overall, we find that covariances built from approximate simulations yield errors on model parameters within $10{{\ \rm per\ cent}}$ of those from the N-body-based covariance.

scholarly journals Small-scale Substructure in Dark Matter Haloes: Where Does Galaxy Formation Come to an End?

2004 ◽  
Vol 220 ◽  
pp. 91-98 ◽  
Author(s):  
J. E. Taylor ◽  
J. Silk ◽  
A. Babul
Keyword(s):  
Dark Matter ◽  
Power Spectrum ◽  
Structure Formation ◽  
Galaxy Formation ◽  
Cold Dark Matter ◽  
High Redshift ◽  
Small Scale ◽  
Analytic Model ◽  
Model Predictions ◽  
Low Mass

Models of structure formation based on cold dark matter predict that most of the small dark matter haloes that first formed at high redshift would have merged into larger systems by the present epoch. Substructure in present-day haloes preserves the remains of these ancient systems, providing the only direct information we may ever have about the low-mass end of the power spectrum. We describe some recent attempts to model halo substructure down to very small masses, using a semi-analytic model of halo formation. We make a preliminary comparison between the model predictions, observations of substructure in lensed systems, and the properties of local satellite galaxies.

scholarly journals Power spectrum response of large-scale structure in 1D and in 3D: tests of prescriptions for post-collapse dynamics

2020 ◽  
Vol 499 (2) ◽  
pp. 1769-1787
Author(s):  
Anaëlle Halle ◽  
Takahiro Nishimichi ◽  
Atsushi Taruya ◽  
Stéphane Colombi ◽  
Francis Bernardeau
Keyword(s):  
Power Spectrum ◽  
Response Function ◽  
Large Scale ◽  
Initial Conditions ◽  
Large Scale Structure ◽  
Scale Structure ◽  
Adaptive Smoothing ◽  
Single Stream ◽  
Non Linear ◽  
Spectrum Response

ABSTRACT The power spectrum response function of the large-scale structure of the Universe describes how the evolved power spectrum is modified by a small change in initial power through non-linear mode coupling of gravitational evolution. It was previously found that the response function for the coupling from small to large scales is strongly suppressed in amplitude, especially at late times, compared to predictions from perturbation theory (PT) based on the single-stream approximation. One obvious explanation for this is that PT fails to describe the dynamics beyond shell crossing. We test this idea by comparing measurements in N-body simulations to prescriptions based on PT but augmented with adaptive smoothing to account for the formation of non-linear structures of various sizes in the multistream regime. We first start with one-dimensional (1D) cosmology, where the Zel’dovich approximation provides the exact solution in the single-stream regime. Similarly to the three-dimensional (3D) case, the response function of the large-scale modes exhibits a strong suppression in amplitude at small scales that cannot be explained by the Zel’dovich solution alone. However, by performing adaptive smoothing of initial conditions to identify haloes of different sizes and solving approximately post-collapse dynamics in the three-stream regime, agreement between theory and simulations drastically improves. We extend our analyses to the 3D case using the pinocchio algorithm, in which similar adaptive smoothing is implemented on the Lagrangian PT fields to identify haloes and is combined with a spherical halo prescription to account for post-collapse dynamics. Again, a suppression is found in the coupling between small- and large-scale modes and the agreement with simulations is improved.

Sign in / Sign up

Export Citation Format

Share Document