Matter Overdensities Research Articles

Gravitational stability of dark energy in galaxies and clusters of galaxies

We analyze the behavior of the scalar field as dark energy of the Universe in a static world of galaxies and clusters of galaxies. We find the analytical solutions of evolution equations of the density and velocity perturbations of dark matter and dark energy, which interact only gravitationally, along with the perturbations of metric in a static world with background Minkowski metric. It was shown that quintessential and phantom dark energy in the static world of galaxies and clusters of galaxies is gravitationally stable and can only oscillate by the influence of self-gravity. In the gravitational field of dark matter perturbations, it is able to condense monotonically, but the amplitude of density and velocity perturbations on all scales remains small. It was also illustrated that the “accretion” of phantom dark energy in the region of dark matter overdensities causes formation of dark energy underdensities-the regions with negative amplitude of density perturbations of dark energy.

General relativistic density perturbations

We investigate a general relativistic mechanism in which spikes generate matter overdensities in the early universe. When the cosmological fluid is tilted, the tilt provides another mechanism in generating matter inhomogeneities. We numerically investigate the effect of a sign change in the tilt, when there is a spike but the tilt does not change sign, and when the spike and the sign change in the tilt coincide. We find that the tilt plays the primary role in generating matter inhomogeneities, and it does so by creating both local overdensities and underdensities. We discuss of the physical implications of the work.

The scale-dependent signature of primordial non-Gaussianity in the large-scale structure of cosmic reionization

(ABRIDGED)The rise of cosmic structure depends upon the statistical distribution of initial density fluctuations generated by inflation. While the simplest models predict an almost perfectly Gaussian distribution, more-general models predict a level of primordial non-Gaussianity (PNG) that observations might yet be sensitive enough to detect. Recent Planck Collaboration measurements of the CMB temperature anisotropy bispectrum significantly tighten the observational limits, but they are still far from the PNG level predicted by the simplest models of inflation. Probing levels below CMB sensitivities will require other methods, such as searching for the statistical imprint of PNG on galactic halo clustering. During the epoch of reionization (EoR), the first stars and galaxies released radiation into the intergalactic medium (IGM) that created ionized patches whose large-scale geometry and evolution reflected the underlying abundance and large-scale clustering of the star-forming galaxies. This statistical connection between ionized patches in the IGM and galactic halos suggests that observing reionization may be another way to constrain PNG. We employ the linear perturbation theory of reionization and semi-analytic models based on the excursion-set formalism to model the effects of PNG on the EoR. We quantify the effects of PNG on the large-scale structure of reionization by deriving the ionized density bias, i.e. ratio of ionized atomic to total matter overdensities in Fourier space, at small wavenumber. Just as previous studies found that PNG creates a scale-dependent signature in the halo bias, so, too, we find a scale-dependent signature in the ionized density bias. Our results, which differ significantly from previous attempts in the literature to characterize this PNG signature, will be applied elsewhere to predict its observable consequences, e.g. in the cosmic 21cm background.

Perturbed recombination from dark matter annihilation

We show that dark matter annihilation around the time of recombination can lead to growing ionization fraction perturbations, that track the linear collapse of matter over-densities. This amplifies small scale cosmological perturbations to the free electron density by a significant amount compared to the usual acoustic oscillations. Electron density perturbations distort the CMB, inducing secondary non-gaussianity. We calculate the CMB bispectrum from recombination, that is marginally observable by Planck. Even though electron perturbations can be markedly boosted compared with the Standard Model prediction, the dark matter effect in the CMB bispectrum turns out to be small and will be difficult to disentangle from the standard model in the foreseeable future.

Symmetries and consistency relations in the large scale structure of the universe

We study the symmetries enjoyed by the Newtonian equations of motion of the non-relativistic dark matter fluid coupled to gravity which give rise to the phenomenon of gravitational instability. We also discuss some consistency relations involving the soft limit of the (n+1)-correlator functions of matter and galaxy overdensities.

Cosmological perturbations in theories with non-minimal coupling between curvature and matter

In this work, we examine how the presence of a non-minimal coupling between spacetime curvature and matter affects the evolution of cosmological perturbations on a homogeneous and isotropic Universe, and hence the formation of large-scale structure. This framework places constraints on the terms which arise due to the coupling with matter and, in particular, on the modified growth of matter density perturbations. We derive approximate analytical solutions for the evolution of matter overdensities during the matter dominated era and discuss the compatibility of the obtained results with the hypothesis that the late time acceleration of the Universe is driven by a non-minimal coupling.

THE PSEUDO-EVOLUTION OF HALO MASS

A dark matter halo is commonly defined as a spherical overdensity of matter with respect to a reference density, such as the critical density or the mean matter density of the Universe. Such definitions can lead to a spurious pseudo-evolution of halo mass simply due to redshift evolution of the reference density, even if its physical density profile remains constant over time. We estimate the amount of such pseudo-evolution of mass between z=1 to 0 for halos identified in a large N-body simulation, and show that it accounts for almost the entire mass evolution of the majority of halos with M200 of about 1E12 solar masses and can be a significant fraction of the apparent mass growth even for cluster-sized halos. We estimate the magnitude of the pseudo-evolution assuming that halo density profiles remain static in physical coordinates, and show that this simple model predicts the pseudo-evolution of halos identified in numerical simulations to good accuracy, albeit with significant scatter. We discuss the impact of pseudo-evolution on the evolution of the halo mass function and show that the non-evolution of the low-mass end of the halo mass function is the result of a fortuitous cancellation between pseudo-evolution and the absorption of small halos into larger hosts. We also show that the evolution of the low mass end of the concentration-mass relation observed in simulations is almost entirely due to the pseudo-evolution of mass. Finally, we discuss the implications of our results for the interpretation of the evolution of various scaling relations between the observable properties of galaxies and galaxy clusters and their halo masses.

Second-order matter fluctuations via higher order galaxy correlators

We provide a formula for extracting the value of the rms of the linear matter fluctuations on a scale R directly from redshift surveys data. It allows to constrain the real-space amplitude of sigma_R without requiring any modeling of the nature and power spectrum of the matter distribution. Furthermore, the formalism is completely insensitive to the character of the bias function, namely its eventual scale or non-linear dependence. By contrasting measurements of sigma_R with predictions from linear perturbation theory, one can test for eventual departures from the standard description of gravity on large cosmological scales. The proposed estimator exploits the information contained in the 1-point moments and 2-point correlators of the matter and galaxy density fields, and it can be applied on cosmic scales where linear and semi-linear perturbative approximations of the evolution of matter overdensities offer a satisfactory description of the full underlying theory. We implement the test with N-body simulations to quantify potential systematics and successfully show that we are able to recover the present day value of sigma_8 `hidden' in the simulation. We also design a consistency check to gauge the soundness of the results inferred when the formalism is applied to real (as opposed to simulated) data. We expect that this approach will provide a sensitive probe of the clustering of matter when applied to future large redshift survey such as BigBOSS and EUCLID.

Spherical collapse in covariant Galileon theory

In this paper we study the evolution of a spherical matter overdensity in the context of the recently introduced Galileon field theory. Our analysis considers the complete covariant Lagrangian in four dimensions. This theory is composed by a potential and a standard kinetic term, a cubic kinetic term and two additional terms that include the coupling between the Galileon and the metric, to preserve the original properties of Galileons also in curved space-times. Here we extend previous studies, which considered both the quintessence and the cubic terms, by focussing on the role of the last two terms. The background evolution we consider is driven by a tracker solution. Studying scalar perturbations in the non-linear regime, we find constraints on the parameter of the model. We will show how the new terms contribute to the collapse phase and how they modify physical parameters, such as the linearized density contrast and the virial overdensity. The results show that the Galileon modifies substantially the dynamics of the collapse, thus making it possible to observationally constrain the parameters of this theory.

Generating Matter Inhomogeneities in General Relativity

In this Letter we discuss a natural general relativistic mechanism that causes inhomogeneities and hence generates matter perturbations in the early Universe. We concentrate on spikes, both incomplete spikes and recurring spikes, that naturally occur in the initial oscillatory regime of general cosmological models. In particular, we explicitly show that spikes occurring in a class of G2 models lead to inhomogeneities that, due to gravitational instability, leave small residual imprints on matter in the form of matter perturbations. The residual matter overdensities from recurring spikes are not local but form on surfaces. We discuss the potential physical consequences of the residual matter imprints and their possible effect on the subsequent formation of large-scale structure.

The shapes and alignments of dark matter halos

We presentmeasurements of the triaxial dark matter halo shapes and alignment correlation functions in theMillennium and Millennium-2 dark matter N-body simulations. These two simulations allow us to measure the distributions of halo shapes down to 10% of the virial radius over a halo mass range of 6 × 109–2 × 1014 h−1M⊙. We largely confirm previous results on the distributions of halo axis ratios asa function of halo mass, but we find that the median angle between halo major axes at different halo radiican vary by a factor of 2 between the Millennium-1 and 2 simulations because of the different massresolution. Thus, error in the shape determinations from limited resolution is potentially degenerate with the misalignment of halo inner and outer shapes used to constrain Brightest Cluster Galaxy alignments in previous works.We also present simplifying parameterizations for the 3-D halo-mass alignment correlation functions that arenecessary ingredients for triaxial halo models of large-scale structure and models of galaxyintrinsic alignments as contaminants for cosmic shear surveys. We measure strong alignments between halos of all masses and the surroundingdark matter overdensities out to several tens of h−1 Mpc, in agreement with observed shear-galaxy and cluster shape correlations. We use these measurements to forecast the contribution to the weak lensing signal around galaxy clusters from correlated mass along the line-of-sight. For prolate clusters with major axes aligned with the line-of-sight the fraction of the weak lensing signal from mass external to the cluster can be twice that predicted if the excess halo alignment correlation is assumed to be zero.

Gravity and large-scale nonlocal bias

The relationship between galaxy and matter overdensities, bias, is most often assumed to be local. This is however unstable under time evolution, we provide proofs under several sets of assumptions. In the simplest model galaxies are created locally and linearly biased at a single time, and subsequently move with the matter (no velocity bias) conserving their comoving number density (no merging). We show that, after this formation time, the bias becomes unavoidably non-local and non-linear at large scales. We identify the non-local gravitationally induced fields in which the galaxy overdensity can be expanded, showing that they can be constructed out of the invariants of the deformation tensor (Galileons). In addition, we show that this result persists if we include an arbitrary evolution of the comoving number density of tracers. We then include velocity bias, and show that new contributions appear, a dipole field being the signature at second order. We test these predictions by studying the dependence of halo overdensities in cells of fixed matter density: measurements in simulations show that departures from the mean bias relation are strongly correlated with the non-local gravitationally induced fields identified by our formalism. The effects on non-local bias seen in the simulations are most important for the most biased halos, as expected from our predictions. The non-locality seen in the simulations is not fully captured by assuming local bias in Lagrangian space. Accounting for these effects when modeling galaxy bias is essential for correctly describing the dependence on triangle shape of the galaxy bispectrum, and hence constraining cosmological parameters and primordial non-Gaussianity. We show that using our formalism we remove an important systematic in the determination of bias parameters from the galaxy bispectrum, particularly for luminous galaxies. (abridged)

Cosmological tests of general relativity: A principal component analysis

The next generation of weak lensing surveys will trace the evolution of matter perturbations and gravitational potentials from the matter dominated epoch until today. Along with constraining the dynamics of dark energy, they will probe the relations between matter overdensities, local curvature, and the Newtonian potential. We work with two functions of time and scale to account for any modifications of these relations in the linear regime from those in the LCDM model. We perform a Principal Component Analysis (PCA) to find the eigenmodes and eigenvalues of these functions for surveys like DES and LSST. This paper builds on and significantly extends the PCA analysis of Zhao et al. (2009) in several ways. In particular, we consider the impact of some of the systematic effects expected in weak lensing surveys. We also present the PCA in terms of other choices of the two functions needed to parameterize modified growth on linear scales, and discuss their merits. We analyze the degeneracy between the modified growth functions and other cosmological parameters, paying special attention to the effective equation of state w(z). Finally, we demonstrate the utility of the PCA as an efficient data compression stage which enables one to easily derive constraints on parameters of specific models without recalculating Fisher matrices from scratch.

Impact of primordial ultracompact minihaloes on the intergalactic medium and first structure formation

(Abridged) Ultracompact Minihaloes (UCMHs), which formed by dark matter accretion onto primordial black holes (PBHs) or initial dark matter overdensity produced by the primordial density perturbation, provide a new type of compact dark matter structure to ionize and heat the IGM after matter-radiation equality z_eq, which is much earlier than the formation of the first cosmological dark halo structure and later first stars. We show that dark matter annihilation density contributed by UCMHs can totally dominated over the homogenous dark matter annihilation background even for a tiny UCMH abundance, and provide a new gamma-ray background in the early Universe. The IGM ionization fraction x_ion and gas temperature T_m can be increased from the recombination residual and adiabatically cooling in the absence of energy injection to the highest value of x_ ion 0.1 and T_m ~ 5000 K at z>10 for the upper bound UCMH abundance constrained by the CMB optical depth. A small fraction of UCMHs are seeded by PBHs. The X-ray emission from gas accretion onto PBHs may totally dominated over dark matter annihilation and become the main cosmic ionization source, but the constraints of gas accretion rate and X-ray absorption by the baryon accumulation within the UCMHs and accretion feedback show that X-ray emission can only be a promising source much later than UCMH annihilation at z<z_m<1000, where z_m depends on the masses of PBHs, their host UCMHs, and the dark matter particles. Also, UCMH radiation including both annihilation and X-ray emission can significantly suppress the low mass first baryonic structure formation. The effects of UCMHs radiation on the baryonic structure evolution are quite small for the gas temperature after virialization, but more significant to enhance the gas chemical quantities such as the ionization fraction and molecular hydrogen abundance in the baryonic objects.

A SINFONI view of flies in the Spiderweb: a galaxy cluster in the making

The environment of the high-z radio galaxy PKS 1138-262 at z~2.2 is a prime example of a forming galaxy cluster. We use deep SINFONI data to perform a detailed study of the kinematics of the galaxies within 60 kpc of the radio core and we link this to the kinematics of the protocluster on the megaparsec scale. Identification of optical emission lines shows that 11 galaxies are at the redshift of the protocluster. The density of line emitters is more than an order of magnitude higher in the core of the protocluster than the larger scale environment. This implies a matter overdensity in the core of delta_m~70 which is similar to the outskirts of local galaxy clusters. The velocity distribution of the confirmed satellite galaxies shows a broad, double-peaked velocity structure with sigma=1360+/-206 km/s. A similar broad, double-peaked distribution was found in a previous study targeting the large scale protocluster structure, indicating that a common process is acting on both small and large scales. Including all spectroscopically confirmed protocluster galaxies, a velocity dispersion of 1013+/-87 km/s is found. We show that the protocluster has likely decoupled from the Hubble flow and is a dynamically evolved structure. Comparison to the Millenium simulation indicates that the protocluster velocity distribution is consistent with that of the most massive haloes at z~2, but we rule out that the protocluster is a fully virialized structure based on dynamical arguments and its X-ray luminosity. Comparison to merging haloes in the Millennium simulation shows that the structure as observed in and around the Spiderweb galaxy is best interpreted as being the result of a merger between two massive haloes. We propose that this merger can result in an increase in star formation and AGN activity in the protocluster core and is possibly an important stage in the evolution of massive cD galaxies.

THE NONLINEAR BIASING OF THE zCOSMOS GALAXIES UP TOz∼ 1 FROM THE 10k SAMPLE

We use the overdensity field reconstructed in the volume of the COSMOS area to study the nonlinear biasing of the zCOSMOS galaxies. The galaxy overdensity field is reconstructed using the current sample of ~8500 accurate zCOSMOS redshifts at I(AB)<22.5 out to z~1 on scales R from 8 to 12 Mpc/h. By comparing the probability distribution function (PDF) of galaxy density contrast delta_g to the lognormal approximation of the PDF of the mass density contrast delta, we obtain the mean biasing function b(delta,z,R) between the galaxy and matter overdensity field and its second moments b(hat) and b(tilde) up to z~1. Over the redshift interval 0.4<z<1 the conditional mean function <delta_g|delta> = b(delta,z,R) delta is of the following characteristic shape. The function vanishes in the most underdense regions and then sharply rises in a nonlinear way towards the mean densities. <delta_g|delta> is almost a linear tracer of the matter in the overdense regions, up to the most overdense regions in which it is nonlinear again and the local effective slope of <delta_g|delta> vs. delta is smaller than unity. The <delta_g|delta> function is evolving only slightly over the redshift interval 0.4<z<1. The linear biasing parameter increases from b(hat)=1.24+/-0.11 at z=0.4 to b(hat)=1.64+/-0.15 at z=1 for the M_B<-20-z sample of galaxies. b(hat) does not show any dependence on the smoothing scale from 8 to 12 Mpc/h, but increases with luminosity. The measured nonlinearity parameter b(tilde)/b(hat) is of the order of a few percent (but it can be consistent with 0) and it does not change with redshift, the smoothing scale or the luminosity. By matching the linear bias of galaxies to the halo bias, we infer that the M_B<-20-z galaxies reside in dark matter haloes with a characteristic mass of about 3-6 x 10^12 Msol, depending on the halo bias fit.

Evolution of spherical overdensity in thawing dark energy models

We study the general evolution of spherical over-densities for thawing class of dark energy models. We model dark energy with scalar fields having canonical as well as non-canonical kinetic energy. For non-canonical case, we consider models where the kinetic energy is of the Born-Infeld Form. We study various potentials like linear, inverse-square, exponential as well as PNGB-type. We also consider the case when dark energy is homogeneous as well as the case when it is inhomogeneous and virializes together with matter. Our study shows that models with linear potential in particular with Born-Infeld type kinetic term can have significant deviation from the $\Lambda$CDM model in terms of density contrast at the time of virialization. Although our approach is a simplified one to study the nonlinear evolution of matter overdensities inside the cluster and is not applicable to actual physical situation, it gives some interesting insights into the nonlinear clustering of matter in the presence of thawing class of dark energy models.

Cosmological perturbations inf(T)gravity

We investigate the cosmological perturbations in $f(T)$ gravity. Examining the pure gravitational perturbations in the scalar sector using a diagonal vierbein, we extract the corresponding dispersion relation, which provides a constraint on the $f(T)$ Ans\"atze that lead to a theory free of instabilities. Additionally, upon inclusion of the matter perturbations, we derive the fully perturbed equations of motion, and we study the growth of matter overdensities. We show that $f(T)$ gravity with $f(T)$ constant coincides with General Relativity, both at the background as well as at the first-order perturbation level. Applying our formalism to the power-law model we find that on large subhorizon scales ($\mathcal{O}(100\text{ }\text{ }\mathrm{Mpc})$ or larger), the evolution of matter overdensity will differ from $\ensuremath{\Lambda}\mathrm{CDM}$ cosmology. Finally, examining the linear perturbations of the vector and tensor sectors, we find that (for the standard choice of vierbein) $f(T)$ gravity is free of massive gravitons.

F(T) gravity mimicking dynamical dark energy. Background and perturbation analysis

We investigate f(T) cosmology in both the background, as well asin the perturbation level, and we present the general formalismfor reconstructing the equivalent one-parameter family of f(T)models for any given dynamical dark energy scenario. Despite thecompletely indistinguishable background behavior, theperturbations break this degeneracy and the growth histories ofall these models differ from one another. As an application wereconstruct the f(T) equivalent for quintessence, and we showthat the deviation of the matter overdensity evolution is strongfor small scales and weak for large scales, while it is negligiblefor large redshifts.

Minimizing the stochasticity of halos in large-scale structure surveys

In recent work (Seljak, Hamaus and Desjacques 2009) it was found that weighting central halo galaxies by halo mass can significantly suppress their stochasticity relative to the dark matter, well below the Poisson model expectation. In this paper we extend this study with the goal of finding the optimal mass-dependent halo weighting and use $N$-body simulations to perform a general analysis of halo stochasticity and its dependence on halo mass. We investigate the stochasticity matrix, defined as $C_{ij}\equiv<(\delta_i -b_i\delta_m)(\delta_j-b_j\delta_m)>$, where $\delta_m$ is the dark matter overdensity in Fourier space, $\delta_i$ the halo overdensity of the $i$-th halo mass bin and $b_i$ the halo bias. In contrast to the Poisson model predictions we detect nonvanishing correlations between different mass bins. We also find the diagonal terms to be sub-Poissonian for the highest-mass halos. The diagonalization of this matrix results in one large and one low eigenvalue, with the remaining eigenvalues close to the Poisson prediction $1/\bar{n}$, where $\bar{n}$ is the mean halo number density. The eigenmode with the lowest eigenvalue contains most of the information and the corresponding eigenvector provides an optimal weighting function to minimize the stochasticity between halos and dark matter. We find this optimal weighting function to match linear mass weighting at high masses, while at the low-mass end the weights approach a constant whose value depends on the low-mass cut in the halo mass function. Finally, we employ the halo model to derive the stochasticity matrix and the scale-dependent bias from an analytical perspective. It is remarkably successful in reproducing our numerical results and predicts that the stochasticity between halos and the dark matter can be reduced further when going to halo masses lower than we can resolve in current simulations.