We calculated spatial correlation functions of galaxies, $\xi(r)$, structure functions, $g(r)=1 +\xi(r)$, gradient functions, $\gamma(r)= d \log g(r)/ d \log r$, and fractal dimension functions, $D(r)= 3+\gamma(r)$, using dark matter particles of the biased $\Lambda$ cold dark matter (CDM) simulation, observed galaxies of the Sloan Digital Sky Survey (SDSS), and simulated galaxies of the Millennium and EAGLE simulations. We analysed how these functions describe fractal and biasing properties of the cosmic web. The correlation functions of the biased $\Lambda$CDM model samples at small distances (particle and galaxy separations), $r \le 2.25$~\Mpc, describe the distribution of matter inside dark matter (DM) halos. In real and simulated galaxy samples, only the brightest galaxies in clusters are visible, and the transition from clusters to filaments occurs at a distance $r \approx 0.8 - 1.5$~\Mpc. Real and simulated galaxies of low luminosity, $M_r \ge -19$, have almost identical correlation lengths and amplitudes, indicating that dwarf galaxies are satellites of brighter galaxies, and do not form a smooth population in voids. The combination of several physical processes (e.g. the formation of halos along the caustics of particle trajectories and the phase synchronisation of density perturbations on various scales) transforms the initial random density field to the current highly non-random density field. Galaxy formation is suppressed in voids, which increases the amplitudes of correlation functions and power spectra of galaxies, and increases the large-scale bias parameter. The combined evidence leads to the large-scale bias parameter of $L_\star$ galaxies the value $b_\star =1.85 \pm 0.15$. We find $r_0(L_\star) = 7.20 \pm 0.19$ for the correlation length of $L_\star$ galaxies.

We study how the indirect observation of dark matter substructures in the Milky Way, using recent stellar stream studies, translates into constraints for different dark matter models. Particularly, we use the measured number of dark subhalos in the mass range $10^7-10^9\,{\rm M_{\odot}}$ to constrain modifications of the subhalo mass function compared to the cold dark matter scenario. We obtain the lower bounds $m_{\rm WDM} > 3.2\,{\rm keV}$ and $m_{\rm FDM} > 5.2\times 10^{-21}\,{\rm eV}$ on the warm dark matter and fuzzy dark matter particle mass, respectively. When dark matter is coupled to a dark radiation bath, we find that kinetic decoupling must take place at temperatures higher than $T_{\rm kd} > 0.7\,{\rm keV}$. We also discuss future prospects of stellar stream observations.

We investigate the evolution of superclusters and supercluster cocoons (basins of attraction), and the influence of cosmological parameters to the evolution. We perform numerical simulations of the evolution of the cosmic web for different cosmological models: the LCDM model with a conventional value of the dark energy (DE) density, the open model OCDM with no DE, the standard SCDM model with no DE, and the Hyper-DE HCDM model with an enhanced DE density value. We find ensembles of superclusters of these models for five evolutionary stages, corresponding to the present epoch z = 0, and to redshifts z = 1, 3, 10, 30. We use diameters of the largest superclusters and the number of superclusters as percolation functions to describe properties of the ensemble of superclusters in the cosmic web. We analyse the size and mass distribution of superclusters in models and in real Sloan Digital Sky Survey (SDSS) based samples. In all models numbers and volumes of supercluster cocoons are independent on cosmological epochs. Supercluster masses increase with time, and geometrical sizes in comoving coordinates decrease with time, for all models. LCDM, OCDM and HCDM models have almost similar percolation parameters. This suggests that the essential parameter, which defines the evolution of superclusters, is the matter density. The DE density influences the growth of the amplitude of density perturbations, and the growth of masses of superclusters, albeit significantly less strongly. The HCDM model has the largest speed of the growth of the amplitude of density fluctuations, and the largest growth of supercluster masses during the evolution. Geometrical diameters and numbers of HCDM superclusters at high threshold densities are larger than for LCDM and OCDM superclusters. SCDM model has about two times more superclusters than other models; SCDM superclusters have smaller diameters and masses.

Our goal is to find the relation between the two-point correlation functions of the projected and spatial density fields of galaxies, and their influence to biasing, fractal and other geometrical properties of the cosmic web. Using spatial (3D) and projected (2D) density fields we calculate 3D and 2D correlation functions of galaxies, $\xi(r)$, structure functions, $g(r)=1 +\xi(r)$, and fractal functions, $\gamma(r)= d \log g(r)/ d \log r$, for a biased $\Lambda$ cold dark matter (CDM) simulation. We analyse how these functions describe biasing, fractal and other properties of the cosmic web. We compare the correlation functions of spatial and projected density fields as descriptors of the cosmic web. Dominant elements of the cosmic web are clusters and filaments, separated by voids filling most of the volume. In individual 2D sheets the positions of clusters and filament do not coincide. As a result, in projection clusters and filaments fill in 2D voids. This leads to the decrease of amplitudes of correlation functions (and power spectra) in projection. For this reason amplitudes of 2D correlation functions are lower than amplitudes of 3D correlation functions, the difference is the larger, the thicker are 2D samples. Spatial correlation functions of galaxies contain valuable information on geometrical properties of the cosmic web, not available in angular correlation functions. 2D correlation functions do not contain information on voids in 3D density field, thus 3D correlation functions cannot be calculated from 2D correlation functions.

Primordial black hole (PBH) dark matter (DM) nonlinear small-scale structure formation begins before the epoch of recombination due to large Poisson density fluctuations. Those small-scale effects still survive today, distinguishing physics of PBH DM structure formation from the one involving WIMP DM. We construct an analytic model for the small-scale PBH velocities that reproduces the velocity floor seen in numerical simulations, and investigate how these motions impact PBH accretion bounds at different redshifts. We find that the effect is small at the time of recombination, leaving the cosmic microwave background bounds on PBH abundance unchanged. However, already at $z=20$ the PBH internal motion significantly reduces their accretion due to the additional $1/v^6$ suppression, affecting the 21 cm bounds. Today the accretion bounds arising from dwarf galaxies or smaller PBH substructures are all reduced by the PBH velocity floor. We also investigate the feasibility for the PBH clusters to coherently accrete gas leading to a possible enhancement proportional to the cluster's occupation number but find this effect to be insignificant for PBH around $10 M_{\odot}$ or lighter. Those results should be reconsidered if the initial PBH distribution is not Poisson, for example, in the case of large initial PBH clustering.

We study the effect of dynamical friction on globular clusters and on the stars evaporated from the globular clusters (stellar streams) moving in a galactic halo. Due to dynamical friction, the position of a globular cluster (GC) as a stream progenitor starts to shift with respect to its original position in the reference frame of initial GC orbit. Therefore the stars that have evaporated at different times have different mean position with respect to the GC position. This shifting results in a certain asymmetry in stellar density distribution between the leading and trailing arms of the stream. The degree of the asymmetry depends on the characteristics of the environment in which the GC and the stream stars move. As GCs are located mainly in outer parts of a galaxy, this makes dynamical friction a unique probe to constrain the underlying dark matter spatial density and velocity distributions. For a GC NGC 3201 we compared our theoretical shift estimates with available observations. Due to large uncertainties in current observation data, we can only conclude that the derived estimates have the same order of magnitude.

Surprisingly high masses of the black holes inferred from the LIGO & Virgo gravitational wave measurements have lead to speculations that the observed mergers might be due to ${\cal O}(10) M_\odot$ primordial black holes (PBHs). Furthermore, it has been suggested that the whole amount of dark matter (DM) might be in that exotic form. We investigate constraints on the PBH DM using NuSTAR Galactic center (GC) X-ray data. We used a robust Monte Carlo approach in conjunction with a radiatively inefficient PBH accretion model with commonly accepted model parameters. Compared to previous studies we allowed for multiple forms of DM density profiles. Most importantly, our study includes treatment of the gas turbulence, which significantly modifies the relative velocity between PBHs and gas. We show that inclusion of the effects of gas turbulence and the uncertainties related to the DM density profile reduces significantly the gas accretion onto PBHs compared to the claimed values in previous papers. It is highly improbable to obtain accreting PBHs brighter than the NuSTAR point source limit using observationally determined gas velocities. As such, one can safely conclude that GC X-ray observations cannot rule out ${\cal O}(10) M_\odot$ PBH DM.

The EDGES experiment has recently measured an anomalous global 21-cm spectrum due to hydrogen absorptions at redshifts of about $z\sim 17$. Model independently, the unusually low temperature of baryons probed by this observable sets strong constraints on any physical process that transfers energy into the baryonic environment at such redshifts. Here we make use of the 21-cm spectrum to derive bounds on the energy injection due to a possible population of ${\cal O}(1-100) M_\odot$ primordial black holes, which induce a wide spectrum of radiation during the accretion of the surrounding gas. After calculating the total radiative intensity of a primordial black hole population, we estimate the amount of heat and ionisations produced in the baryonic gas and compute the resulting thermal history of the Universe with a modified version of RECFAST code. Finally, by imposing that the temperature of the gas at $z\sim 17$ does not exceed the indications of EDGES, we constrain the possible abundance of primordial black holes. Depending on uncertainties related to the accretion model, we find that ${\cal O}(10) M_\odot$ primordial black holes can only contribute to a fraction $f_{\rm PBH}<(1-10^{-3})$ of the total dark matter abundance.

We study possible effects of a dark matter (DM) core on the maximum mass of a neutron star (NS), on the mass-radius relation and on the NS tidal deformability parameter $\Lambda$. We show that all these quantities would in general be reduced in the presence of a DM core. In particular, our calculations indicate that the presence of a DM core with a mass fraction $\sim 5\%$ could affect significantly the interpretation of these NS data as constraints on the nuclear equation of state (EOS), potentially excluding some EOS models on the basis of the measured mass of PSR J0348+0432, while allowing other EOS models to become consistent with the LIGO/Virgo upper limit on $\Lambda$. Specific scenarios for generating such DM cores are explored in an Appendix.

The recent observation of the 21-cm absorption signal by the EDGES experiment provides a new observational window into the dynamics of the young Universe. Based on this result, we constrain the properties of ALPs and excited dark matter via the energy injections into the gaseous medium that these models produce. In particular, we derive bounds that outperform the present constraints for energy injections in the 10.2 eV to 50 keV range by analysing the intensity of Lyman-$\alpha$ photons and the electron emissivity produced by ionising radiation. Our results also show that once the standard temperature of the soft radiation background is assumed, the explanations of the observed 3.5 keV photon signal within the considered dark matter model are excluded at 95\% confidence level.

Motivated by the recent detection of the gravitational wave signal emitted by a binary neutron star merger, we analyse the possible impact of dark matter on such signals. We show that dark matter cores in merging neutron stars may yield an observable supplementary peak in the gravitational wave power spectral density following the merger, which could be distinguished from the features produced by the neutron components.

We investigate the potential of large X-ray-selected AGN samples for detecting baryonic acoustic oscillations (BAO). Though AGN selection in X-ray band is very clean and efficient, it does not provide redshift information, and thus needs to be complemented with an optical follow-up. The main focus of this study is (i) to find the requirements needed for the quality of the optical follow-up and (ii) to formulate the optimal strategy of the X-ray survey, in order to detect the BAO. We demonstrate that redshift accuracy of $\sigma_0=10^{-2}$ at $z=1$ and the catastrophic failure rate of $f_{\rm fail}\lesssim 30\%$ are sufficient for a reliable detection of BAO in future X-ray surveys. Spectroscopic quality redshifts ($\sigma_0=10^{-3}$ and $f_{\rm fail}\sim 0$) will boost the confidence level of the BAO detection by a factor of $\sim 2$. For meaningful detection of BAO, X-ray surveys of moderate depth of $F_{\rm lim}\sim {\rm few~} 10^{-15}$ erg/s/cm$^2$ covering sky area from a few hundred to $\sim$ten thousand square degrees are required. The optimal strategy for the BAO detection does not necessarily require full sky coverage. For example, in a 1000-day long survey by an eROSITA type telescope, an optimal strategy would be to survey a sky area of $\sim 9000$ deg$^2$, yielding a $\sim 16\sigma$ BAO detection. A similar detection will be achieved by ATHENA+ or WFXT class telescopes in a survey with a duration of 100 days, covering a similar sky area. XMM-Newton can achieve a marginal BAO detection in a 100-day survey covering $\sim 400$ deg$^2$. These surveys would demand a moderate-to-high cost in terms the optical follow-ups, requiring determination of redshifts of $\sim 10^5$ (XMM-Newton) to $\sim 3\times 10^6$ objects (eROSITA, ATHENA+, and WFXT) in these sky areas.

The recently claimed anomaly in the measurement of the 21 cm hydrogen absorption signal by EDGES at $z\sim 17$, if cosmological, requires the existence of new physics. The possible attempts to resolve the anomaly rely on either (i) cooling the hydrogen gas via new dark matter-hydrogen interactions or (ii) modifying the soft photon background beyond the standard CMB one, as possibly suggested also by the ARCADE~2 excess. We argue that solutions belonging to the first class are generally in tension with cosmological dark matter probes once simple dark sector models are considered. Therefore, we propose soft photon emission by light dark matter as a natural solution to the 21 cm anomaly, studying a few realizations of this scenario. We find that the signal singles out a photophilic dark matter candidate characterised by an enhanced collective decay mechanism, such as axion mini-clusters.

We study the surface brightness fluctuations of the cosmic X-ray background (CXB) using Chandra data of XBOOTES. After masking out resolved sources we compute the power spectrum of fluctuations of the unresolved CXB for angular scales from ~2 arcsec to ~3 deg. The non-trivial large-scale structure (LSS) signal dominates over the shot-noise of unresolved point sources at all scales above ~1 arcmin and is produced mainly by the intracluster medium (ICM) of unresolved clusters and groups of galaxies, as shown in our previous publication. The shot-noise-subtracted power spectrum of CXB fluctuations has a power-law shape with the slope of $\Gamma = 0.96 \pm 0.06$. Its energy spectrum is well described by the redshifted emission spectrum of optically-thin plasma with the best-fit temperature of $T \approx 1.3$ keV and the best-fit redshift of $z \approx 0.40$. They are in good agreement with theoretical expectations based on the X-ray luminosity function and scaling relations of clusters. From these values we estimate the typical mass and luminosity of the objects responsible for CXB fluctuations, $M_{500} \sim 10^{13.6}\,{\rm M}_{\odot}/h$ and $L_{0.5-2.0\,{\rm keV}} \sim 10^{42.5}$ erg/s. On the other hand, the flux-weighted mean temperature and redshift of resolved clusters are $T \approx 2.4$ keV and $z \approx 0.23$, confirming that fluctuations of unresolved CXB are caused by cooler (i.e. less massive) and more distant clusters, as expected. We show that the power spectrum shape is sensitive to the ICM structure all the way to the outskirts, out to $\sim{\rm few}\times R_{500}$. We also look for possible contribution of the warm-hot intergalactic medium (WHIM) to the observed CXB fluctuations. Our results underline the significant diagnostics potential of the CXB fluctuation analysis in studying the ICM structure in clusters.

Fluctuations of the surface brightness of cosmic X-ray background (CXB) carry unique information about faint and low luminosity source populations, which is inaccessible for conventional large-scale structure (LSS) studies based on resolved sources. We used Chandra data of the XBOOTES field ($\sim9\,\mathrm{deg^2}$) to conduct the most accurate measurement to date of the power spectrum of fluctuations of the unresolved CXB on the angular scales of $\sim3\,$arcsec $-$ $\sim17\,$arcmin. We find that at sub-arcmin angular scales, the power spectrum is consistent with the AGN shot noise, without much need for any significant contribution from their one-halo term. This is consistent with the theoretical expectation that low-luminosity AGN reside alone in their dark matter halos. However, at larger angular scales we detect a significant LSS signal above the AGN shot noise. Its power spectrum, obtained after subtracting the AGN shot noise, follows a power law with the slope of $-0.8\pm0.1$ and its amplitude is much larger than what can be plausibly explained by the two-halo term of AGN. We demonstrate that the detected LSS signal is produced by unresolved clusters and groups of galaxies. For the flux limit of the XBOOTES survey, their flux-weighted mean redshift equals $\left<z\right>\sim0.3$, and the mean temperature of their intracluster medium (ICM), $\left<T\right>\approx 1.4$ keV, corresponds to the mass of $M_{500} \sim 10^{13.5}\,\mathrm{M}_\odot$. The power spectrum of CXB fluctuations carries information about the redshift distribution of these objects and the spatial structure of their ICM on the linear scales of up to $\sim$Mpc, i.e. of the order of the virial radius.

Theories where the Planck scale is dynamically generated from dimensionless interactions provide predictive inflationary potentials and super-Planckian field variations. We first study the minimal single-field realisation in the low-energy effective field theory limit, finding the predictions $n_s \approx 0.96$ for the spectral index and $r \approx 0.13$ for the tensor-to-scalar ratio, which can be reduced down to $\approx 0.04$ in presence of large couplings. Next we consider agravity as a dimensionless quantum gravity theory finding a multi-field inflation that converges towards an attractor trajectory that predicts $n_s\approx 0.96$ and $0.003<r<0.13$, interpolating between the quadratic and Starobinsky inflation. These theories relate the smallness of the weak scale to the smallness of inflationary perturbations: both arise naturally because of small couplings, implying a reheating temperature of $10^{7-9}$ GeV. A measurement of $r$ by Keck/Bicep3 would give us information on quantum gravity in the dimensionless scenario.

Our goal is to find a minimalistic model that describes the luminosity function and large-scale clustering bias of the X-ray-selected AGN in the general framework of the concordance LCDM model. We assume that a simple population-averaged scaling relation between the AGN X-ray luminosity L_X and the host dark matter halo mass M_h exists. With such a relation, the AGN X-ray luminosity function can be computed from the halo mass function. Using the concordance LCDM halo mass function for the latter, we obtain the M_h-L_X relation required to match the redshift-dependent AGN X-ray luminosity function known from X-ray observations. We find that with a simple power-law-scaling M_h \propto L_X^\Gamma(z), our model can successfully reproduce the observed X-ray luminosity function. Furthermore, we automatically obtain predictions for the large-scale AGN clustering amplitudes and their dependence on the luminosity and redshift, which seem to be compatible with AGN clustering measurements. Our model also includes the redshift-dependent AGN duty cycle which peaks at the redshift z ~ 1, and its peak value is consistent with unity, suggesting that on average there is no more than one AGN per dark matter halo. For a typical X-ray-selected AGN at z ~ 1, our best-fit M_h-L_X scaling implies low Eddington ratio L_X/L_Edd ~ 10^{-4}-10^{-3} (2-10 keV band, no bolometric correction applied) and correspondingly large mass growth e-folding times, suggesting that the typical X-ray AGN are dominantly fueled via relatively inefficient 'hot-halo' accretion mode.

We investigate the effect of the window function on the multipole power spectrum in two different ways. First, we consider the convolved power spectrum including the window effect, which is obtained by following the familiar (FKP) method developed by Feldman, Kaiser and Peacock. We show how the convolved multipole power spectrum is related to the original power spectrum, using the multipole moments of the window function. Second, we investigate the deconvolved power spectrum, which is obtained by using the Fourier deconvolution theorem. In the second approach, we measure the multipole power spectrum deconvolved from the window effect. We demonstrate how to deal with the window effect in these two approaches, applying them to the Sloan Digital Sky Survey (SDSS) luminous red galaxy (LRG) sample.

The four-year X-ray all-sky survey (eRASS) of the eROSITA telescope aboard the Spektrum-Roentgen-Gamma satellite will detect ~3 million active galactic nuclei (AGN) with a median redshift of z~1 and a typical luminosity of L_(0.5-2.0 keV) ~ 10^(44) erg/s. We show that this unprecedented AGN sample, complemented with redshift information, will supply us with outstanding opportunities for large-scale structure research. For the first time, detailed redshift- and luminosity-resolved studies of the bias factor for X-ray selected AGN will become possible. The eRASS AGN sample will not only improve the redshift- and luminosity-resolution of these studies, but will also expand their luminosity range beyond L_(0.5-2.0 keV) ~ 10^(44) erg/s, thus enabling a direct comparison of the clustering properties of luminous X-ray AGN and optical quasars. These studies will dramatically improve our understanding of the AGN environment, triggering mechanisms, the growth of supermassive black holes and their co-evolution with dark matter halos. The eRASS AGN sample will become a powerful cosmological probe. It will enable detecting baryonic acoustic oscillations (BAOs) for the first time with X-ray selected AGN. With the data from the entire extragalactic sky, BAO will be detected at a >~10sigma confidence level in the full redshift range and with ~8sigma confidence in the 0.8 < z < 2.0 range, which is currently not covered by any existing BAO surveys. To exploit the full potential of the eRASS AGN sample, photometric and spectroscopic surveys of large areas and a sufficient depth will be needed.

AGN are known to account for a major portion, if not all, of the cosmic X-ray background radiation. The dominant sharp spectral feature in their spectra is the 6.4 keV fluorescent line of iron, which may contribute to as much as ~ 5-10 % of the CXB spectral intensity at ~ 2-6 keV. Owing to cosmological redshift, the line photons detected at the energy E carry information about objects located at the redshift z=6.4/E-1. In particular, imprinted in their angular fluctuations is the information about the large-scale structure at redshift z. This opens the possibility of performing the Fe K_alpha line tomography of the cosmic large-scale structure. We show that detection of the tomographic signal at ~100 sigma confidence requires an all-sky survey by an instrument with an effective area of ~10 m^2 and field of view of ~1 deg^2. The signal is strongest for objects located at the redshift z~1 and at the angular scales corresponding to l ~ 100-300, therefore an optimal detection can be achieved with an instrument having a rather modest angular resolution of ~ 0.1-0.5 deg. For such an instrument, the CCD-type energy resolution of ~ 100-200 eV FWHM is entirely sufficient for the optimal separation of the signals coming from different redshifts. The gain in the signal strength that could potentially be achieved with energy resolution comparable to the line width is nullified by the photon counting and AGN discreteness noise. Among the currently planned and proposed missions, these requirements are best satisfied by LOFT, even though it was proposed for an entirely different purpose. Among others, clear detection should be achieved by WFXT (~ 20-35 sigma) and ATHENA (~ 10-20 sigma). eROSITA, in the course of its 4 years all-sky survey, will detect the tomographic signal only marginally.

We provide ingredients and recipes for computing signals of TeV-scale Dark Matter annihilations and decays in the Galaxy and beyond. For each DM channel, we present the energy spectra of electrons and positrons, antiprotons, antideuterons, gamma rays, neutrinos and antineutrinos e, mu, tau at production, computed by high-statistics simulations. We estimate the Monte Carlo uncertainty by comparing the results yielded by the Pythia and Herwig event generators. We then provide the propagation functions for charged particles in the Galaxy, for several DM distribution profiles and sets of propagation parameters. Propagation of electrons and positrons is performed with an improved semi-analytic method that takes into account position-dependent energy losses in the Milky Way. Using such propagation functions, we compute the energy spectra of electrons and positrons, antiprotons and antideuterons at the location of the Earth. We then present the gamma ray fluxes, both from prompt emission and from Inverse Compton scattering in the galactic halo. Finally, we provide the spectra of extragalactic gamma rays. All results are available in numerical form and ready to be consumed.

According to the modern cosmological paradigm galaxies and galaxy systems form from tiny density perturbations generated during the very early phase of the evolution of the Universe. Using numerical simulations we study the evolution of phases of density perturbations of different scales to understand the formation and evolution of the cosmic web. We apply the wavelet analysis to follow the evolution of high-density regions (clusters and superclusters) of the cosmic web. We show that the positions of maxima and minima of density waves (their spatial phases) almost do not change during the evolution of the structure. Positions of extrema of density perturbations are the more stable, the larger is the wavelength of perturbations. Combining observational and simulation data we conclude that the skeleton of the cosmic web was present already in an early stage of structure evolution.

Our goal is to see how density waves of different scale combine to form voids between galaxy systems of various scale. We perform numerical simulations of structure formation in cubes of size 100 and 256 Mpc/h, with resolutions 256^3 and 512^3 particles and cells. To understand the role of density perturbations of various scale we cut power spectra at scales from 8 to 128 Mpc/h, using in all series identical initial random realisations. We find that small haloes and short filaments form all over the simulation box, if perturbations only up to scale 8 Mpc/h are present. The phenomenon of large multi-scale voids in the cosmic web requires the presence of an extended spectrum of primordial density perturbations. The void phenomenon is due to the action of two processes: the synchronisation of density perturbations of medium and large scales, and the suppression of galaxy formation in low-density regions by the combined action of negative sections of medium- and large-scale density perturbations, so that their densities are less than the mean density, and thus during the evolution their densities decrease.

We follow the evolution of galaxy systems in numerical simulation. Our goal is to understand the role of density perturbations of various scales in the formation and evolution of the cosmic web. We perform numerical simulations with the full power spectrum of perturbations, and with spectrum cut at long wavelengths. Additionally, we have one model, where we cut the intermediate waves. We analyze the density field and study the void sizes and density field clusters in different models. Our analysis shows that the fine structure (groups and clusters of galaxies) is created by small-scale density perturbations of scale $\leq 8$ \Mpc. Filaments of galaxies and clusters are created by perturbations of intermediate scale from $\sim 8$ to $\sim 32$ \Mpc, superclusters of galaxies by larger perturbations. We conclude that the scale of the pattern of the cosmic web is determined by density perturbations of scale up to $\sim 100$ \Mpc. Larger perturbations do not change the pattern of the web, but modulate the richness of galaxy systems, and make voids emptier. The stop of the increase of the scale of the pattern of the cosmic web with increasing scale of density perturbations can probably be explained as the freezing of the web at redshift $z\simeq 0.7$.

The degeneracy between geometric (Alcock-Paczynski) and dynamic (redshift space) distortions in the pattern of the galaxy distribution has been a long standing problem in the study of the large scale structure of the universe. We examine the possibility of lifting this degeneracy and constraining cosmological parameters by using the Baryon Acoustic Oscillation (BAO) scale as a feature of known physical size, the sound horizon r_s ~ 150 Mpc. We callibrate this scale with the equivalent feature in the Cosmic Microwave Background (CMB). First, we construct a toy model of a power spectrum which includes the BAO as well as geometric and dynamic distortions. By adding a prior onto the sound horizon of ~1% we show, using a Fisher matrix analysis, that error ellipses for line of sight and tangential distortion parameters shrink by a factor of two for a 20(Gpc/h)^3 `DESpec/BigBOSS'-like galaxy survey including shot noise. This improvement is even more marked in smaller surveys. We also carry out a Monte Carlo Nested Sampling analysis on our parameter space. We find that Monte Carlo and Fisher methods can agree reasonably well for surveys with large volume but differ greatly for small volume surveys.

(Context) We calculate constraints from current and future cosmic microwave background (CMB) measurements on annihilating dark matter (DM) with masses below the electroweak scale: m_{DM}=5-100 GeV. In particular, we focus our attention on the lower end of this mass range, as DM particles with masses m_{DM} ~ 10 GeV have been recently claimed to be consistent with the CoGeNT and DAMA/LIBRA results, while also providing viable DM candidates to explain the measurements of Fermi and WMAP haze. (Aims) We study the model (in)dependence of CMB spectrum on particle physics DM models, large scale structure formation and cosmological uncertainties. We attempt to find a simple and practical recipe for estimating current and future CMB bounds on a broad class of DM annihilation models. (Results) We show that in the studied DM mass range the CMB signal of DM annihilations is independent of the details of large scale structure formation, distribution and profile of DM halos and other cosmological uncertainties. All particle physics models of DM annihilation can be described with only one parameter, the fraction of energy carried away by neutrinos in DM annihilation. As the main result we provide a simple and rather generic fitting formula for calculating CMB constraints on the annihilation cross section of light WIMPs. We show that thermal relic DM in the CoGeNT, DAMA favoured mass range is in a serious conflict with present CMB data for the annihilation channels with few neutrinos, and will definitely be tested by the Planck mission for all possible DM annihilation channels. Also, our findings strongly disfavor the claim that thermal relic DM annihilations with m_{DM} ~ 10 GeV and $<sigma_av> ~ 9x10^{-25} cm^{3}s^{-1} could be a cause for Fermi and WMAP haze.

We study in detail the model by Isidori and Kamenik that is claimed to explain the top quark forward-backward asymmetry at Tevatron, provide GeV-scale dark matter (DM), and possibly improve the agreement between data and theory in Tevatron W+jj events. We compute the DM thermal relic density, the spin-independent DM-nucleon scattering cross section, and the cosmic microwave background constraints on both Dirac and Majorana neutralino DM in the parameter space that explains the top asymmetry. A stable light neutralino is not allowed unless the local DM density is 3-4 times smaller than expected, in which case Dirac DM with mass around 3 GeV may be possible, to be tested by the Planck mission. The model predicts a too broad excess in the dijet distribution and a strong modification of the missing E_T distribution in W+jj events.

We analyze the recently published Fermi-LAT diffuse gamma-ray measurements in the context of leptonically annihilating or decaying dark matter (DM) with the aim to explain simultaneously the isotropic diffuse gamma-ray and the PAMELA, Fermi and HESS (PFH) anomalous $e^\pm$ data. Five different DM annihilation/decay channels $2e$, $2\mu$, $2\tau$, $4e$, or $4\mu$ (the latter two via an intermediate light particle $\phi$) are generated with PYTHIA. We calculate both the Galactic and extragalactic prompt and inverse Compton (IC) contributions to the resulting gamma-ray spectra. To find the Galactic IC spectra we use the interstellar radiation field model from the latest release of GALPROP. For the extragalactic signal we show that the amplitude of the prompt gamma-emission is very sensitive to the assumed model for the extragalactic background light. For our Galaxy we use the Einasto, NFW and Isothermal DM density profiles and include the effects of DM substructure assuming a simple subhalo model. Our calculations show that for the annihilating DM the extragalactic gamma-ray signal can dominate only if rather extreme power-law concentration-mass relation $C(M)$ is used, while more realistic $C(M)$ relations make the extragalactic component comparable or subdominant to the Galactic signal. For the decaying DM the Galactic signal always exceeds the extragalactic one. In the case of annihilating DM the PFH favored parameters can be ruled out only if power-law $C(M)$ relation is assumed. For DM decaying into $2\mu$ or $4\mu$ the PFH favored DM parameters are not in conflict with the Fermi gamma-ray data. We find that, due to the (almost) featureless Galactic IC spectrum and the DM halo substructure, annihilating DM may give a good simultaneous fit to the isotropic diffuse gamma-ray and to the PFH $e^\pm$ data without being in clear conflict with the other Fermi-LAT gamma-ray measurements.

The luminosity function of galaxies of the Early Data Release of the Sloan Digital Sky Survey (SDSS) and the Las Campanas Redshift Survey (LCRS) is calculated. The luminosity function depends on redshift, density of the environment and is different for the Norther and Southern slice of SDSS. Luminosity functions is used to derive the number and luminosity density fields of galaxies of the SDSS and LCRS surveys with a grid size of 1 h^{-1} Mpc for flat cosmological models with \Omega_m=0.3 and \Omega_\Lambda =0.7. We find that the luminosity function depends on the density of the environment: in high-density regions brightest galaxies are more luminous than in low-density regions by a factor up to 5.