The gravity field maps of the satellite gravimetry missions GRACE (Gravity Recovery and Climate Experiment) and GRACE Follow-On are derived by means of precise orbit determination. The key observation is the biased inter-satellite range, which is measured primarily by a K-Band Ranging system (KBR) in GRACE and GRACE Follow-On. The GRACE Follow-On satellites are additionally equipped with a Laser Ranging Interferometer (LRI), which provides measurements with lower noise compared to the KBR. The biased range of KBR and LRI needs to be converted for gravity field recovery into an instantaneous range, i.e. the biased Euclidean distance between the satellites' center-of-mass at the same time. One contributor to the difference between measured and instantaneous range arises due to the non-zero travel time of electro-magnetic waves between the spacecraft. We revisit the calculation of the light time correction (LTC) from first principles considering general relativistic effects and state-of-the-art models of Earth's potential field. The novel analytical expressions for the LTC of KBR and LRI can circumvent numerical limitations of the classical approach. The dependency of the LTC on geopotential models and on the parameterization is studied, and afterwards the results are compared against the LTC provided in the official datasets of GRACE and GRACE Follow-On. It is shown that the new approach has a significantly lower noise, well below the instrument noise of current instruments, especially relevant for the LRI, and even if used with kinematic orbit products. This allows calculating the LTC accurate enough even for the next generation of gravimetric missions.

The Laser Ranging Interferometer (LRI) instrument on the Gravity Recovery and Climate Experiment (GRACE) Follow-On mission has provided the first laser interferometric range measurements between remote spacecraft, separated by approximately 220 km. Autonomous controls that lock the laser frequency to a cavity reference and establish the 5 degree of freedom two-way laser link between remote spacecraft succeeded on the first attempt. Active beam pointing based on differential wavefront sensing compensates spacecraft attitude fluctuations. The LRI has operated continuously without breaks in phase tracking for more than 50 days, and has shown biased range measurements similar to the primary ranging instrument based on microwaves, but with much less noise at a level of $1\,{\rm nm}/\sqrt{\rm Hz}$ at Fourier frequencies above 100 mHz.

The primary difficulty in measuring dynamical masses of galaxy clusters from galaxy data lies in the separation between true cluster members from interloping galaxies along the line of sight. We study the impact of membership contamination and incompleteness on cluster mass estimates obtained with 25 commonly used techniques applied to nearly 1000 mock clusters. We show that all methods overestimate or underestimate cluster masses when applied to contaminated or incomplete galaxy samples respectively. This appears to be the main source of the intrinsic scatter in the mass scaling relation. Applying corrections based on a prior knowledge of contamination and incompleteness can reduce the scatter to the level of shot noise expected for poorly sampled clusters. We establish an empirical model quantifying the effect of imperfect membership on cluster mass estimation and discuss its universal and method-dependent features. We find that both imperfect membership and the response of the mass estimators depend on cluster mass, effectively causing a flattening of the estimated - true mass relation. Imperfect membership thus alters cluster counts determined from spectroscopic surveys, hence the cosmological parameters that depend on such counts.

With the advent of wide-field cosmological surveys, we are approaching samples of hundreds of thousands of galaxy clusters. While such large numbers will help reduce statistical uncertainties, the control of systematics in cluster masses becomes ever more crucial. Here we examine the effects of an important source of systematic uncertainty in galaxy-based cluster mass estimation techniques: the presence of significant dynamical substructure. Dynamical substructure manifests as dynamically distinct subgroups in phase-space, indicating an 'unrelaxed' state. This issue affects around a quarter of clusters in a generally selected sample. We employ a set of mock clusters whose masses have been measured homogeneously with commonly-used galaxy-based mass estimation techniques (kinematic, richness, caustic, radial methods). We use these to study how the relation between observationally estimated and true cluster mass depends on the presence of substructure, as identified by various popular diagnostics. We find that the scatter for an ensemble of clusters does not increase dramatically for clusters with dynamical substructure. However, we find a systematic bias for all methods, such that clusters with significant substructure have higher measured masses than their relaxed counterparts. This bias depends on cluster mass: the most massive clusters are largely unaffected by the presence of significant substructure, but masses are significantly overestimated for lower mass clusters, by $\sim10\%$ at $10^{14}$ and $\geq20\%$ for $\leq10^{13.5}$. The use of cluster samples with different levels of substructure can, therefore, bias certain cosmological parameters up to a level comparable to the typical uncertainties in current cosmological studies.

We present a Bayesian phase-space reconstruction of the cosmic large-scale matter density and velocity fields from the SDSS-III Baryon Oscillations Spectroscopic Survey Data Release 12 (BOSS DR12) CMASS galaxy clustering catalogue. We rely on a given $\Lambda$CDM cosmology, a mesh resolution in the range of 6-10 $h^{-1}$ Mpc, and a lognormal-Poisson model with a redshift dependent nonlinear bias. The bias parameters are derived from the data and a general renormalised perturbation theory approach. We use combined Gibbs and Hamiltonian sampling, implemented in the \textsc{argo} code, to iteratively reconstruct the dark matter density field and the coherent peculiar velocities of individual galaxies, correcting hereby for coherent redshift space distortions (RSD). Our tests relying on accurate $N$-body based mock galaxy catalogues, show unbiased real space power spectra of the nonlinear density field up to $k\sim0.2\, h$ Mpc$^{-1}$, and vanishing quadrupoles down to $r\sim20\,h^{-1}$ Mpc. We also demonstrate that the nonlinear cosmic web can be obtained from the tidal field tensor based on the Gaussian component of the reconstructed density field. We find that the reconstructed velocities have a statistical correlation coefficient compared to the true velocities of each individual lightcone mock galaxy of $r\sim0.68$ including about 10% of satellite galaxies with virial motions (about $r=0.75$ without satellites). The power spectra of the velocity divergence agree well with theoretical predictions up to $k\sim0.2\,h\,{\rm Mpc}^{-1}$. This work will be especially useful to improve, e.g. BAO reconstructions, kinematic Sunyaev-Zeldovich (kSZ), integrated Sachs-Wolfe (ISW) measurements, or environmental studies.

We study the signatures of reionization and ionizing properties of the early galaxies in the cross-correlations between the 21cm emission from the spin-flip transition of neutral hydrogen (HI) and the underlying galaxy population, in particular a sub-population of galaxies visible as Lyman Alpha Emitters (LAEs). With both observables simultaneously derived from a $z\simeq6.6$ hydrodynamical simulation (GADGET-2) snapshot post-processed with a radiative transfer code (pCRASH) and a dust model, we perform a parameter study and aim to constrain both the average intergalactic medium (IGM) ionization state ($1-\langle \chi_{HI} \rangle$) and the reionization topology (outside-in versus inside-out). We find that in our model LAEs occupy the densest and most-ionized regions resulting in a very strong anti-correlation between the LAEs and the 21cm emission. A 1000h SKA-LOW1 - Subaru Hyper Suprime Cam experiment can provide constraints on $\langle \chi_{HI} \rangle$, allowing us to distinguish between IGM ionization levels of 50%, 25%, 10% and fully ionized at scales $r<10$ comoving Mpc (assuming foreground avoidance for SKA). Our results support the inside-out reionization scenario where the densest knots (under-dense voids) are ionized first (last) for $\langle \chi_{HI} \rangle >= 0.1$. Further, 1000h SKA-LOW1 observations should be able to confirm the inside-out scenario by detecting a lower 21cm brightness temperature (by about 2-10 mK) in the densest regions ($> 2$ arcminute scales) hosting LAEs compared to lower-density regions devoid of them.

We calculate Lyman Alpha Emitter (LAE) angular correlation functions (ACFs) at $z \simeq 6.6$ and the fraction of lifetime (for the 100 Myrs preceding $z\simeq6.6$) galaxies spend as Lyman Break Galaxies (LBGs) or as LBGs with Lyman Alpha (Ly$\alpha$) emission using a model that combines SPH cosmological simulations (GADGET-2), dust attenuation and a radiative transfer code (pCRASH). The ACFs are a powerful tool that significantly narrows the 3D parameter space allowed by LAE Ly$\alpha$ and UV luminosity functions (LFs) alone. With this work, we simultaneously constrain the escape fraction of ionizing photons $f_{esc}=0.05-0.5$, the mean fraction of neutral hydrogen in the intergalactic medium (IGM) $\langle \chi_{HI} \rangle \leq 0.01$ and the dust-dependent ratio of the escape fractions of Ly$\alpha$ and UV continuum photons $f_{\alpha}/f_c=0.6-1.2$. Our results show that reionization has the largest impact on the amplitude of the ACFs, and its imprints are clearly distinguishable from those of $f_{esc}$ and $f_\alpha/f_c$. We also show that galaxies with a critical stellar mass of $M_* = 10^{8.5} (10^{9.5})M_{\odot}$ produce enough luminosity to stay visible as LBGs (LAEs). Finally, the fraction of time during the past 100 Myrs prior to z=6.6 a galaxy spends as a LBG or as a LBG with Ly$\alpha$ emission increases with the UV magnitude (and the stellar mass $M_*$): considering observed (dust and IGM attenuated) luminosities, the fraction of time a galaxy spends as a LBG (LAE) increases from 65% to 100% (0-100%) as $M_{UV}$ decreases from $M_{UV} = -18.0$ to $-23.5$ ($M_*$ increases from $10^8-10^{10.5} M_{\odot}$). Thus in our model the brightest (most massive) LBGs most often show Ly$\alpha$ emission.

This article is the second in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our aim is to quantify the scatter, systematic bias and completeness of cluster masses derived from a diverse set of 25 galaxy-based methods using two contrasting mock galaxy catalogues based on a sophisticated halo occupation model and a semi-analytic model. Analysing 968 clusters, we find a wide range in the RMS errors in log M200c delivered by the different methods (0.18 to 1.08 dex, i.e., a factor of ~1.5 to 12), with abundance matching and richness methods providing the best results, irrespective of the input model assumptions. In addition, certain methods produce a significant number of catastrophic cases where the mass is under- or over-estimated by a factor greater than 10. Given the steeply falling high-mass end of the cluster mass function, we recommend that richness or abundance matching-based methods are used in conjunction with these methods as a sanity check for studies selecting high mass clusters. We see a stronger correlation of the recovered to input number of galaxies for both catalogues in comparison with the group/cluster mass, however, this does not guarantee that the correct member galaxies are being selected. We do not observe significantly higher scatter for either mock galaxy catalogues. Our results have implications for cosmological analyses that utilise the masses, richnesses, or abundances of clusters, which have different uncertainties when different methods are used.

We present a Bayesian reconstruction algorithm to generate unbiased samples of the underlying dark matter field from halo catalogues. Our new contribution consists of implementing a non-Poisson likelihood including a deterministic non-linear and scale-dependent bias. In particular we present the Hamiltonian equations of motions for the negative binomial (NB) probability distribution function. This permits us to efficiently sample the posterior distribution function of density fields given a sample of galaxies using the Hamiltonian Monte Carlo technique implemented in the Argo code. We have tested our algorithm with the Bolshoi $N$-body simulation at redshift $z = 0$, inferring the underlying dark matter density field from sub-samples of the halo catalogue with biases smaller and larger than one. Our method shows that we can draw closely unbiased samples (compatible within 1-$\sigma$) from the posterior distribution up to scales of about $k$~1 h/Mpc in terms of power-spectra and cell-to-cell correlations. We find that a Poisson likelihood yields reconstructions with power spectra deviating more than 10% at $k$=0.2 h/Mpc. Our reconstruction algorithm is especially suited for emission line galaxy data for which a complex non-linear stochastic biasing treatment beyond Poissonity becomes indispensable.

Recently the Chiral Cosmological Model (CCM) coupled to cold dark matter (CDM) has been investigated as $\sigma$CDM model to study the observed accelerated expansion of the Universe. Dark sector fields (as Dark Energy content) coupled to cosmic dust were considered as the source of Einstein gravity in Friedmann-Robertson-Walker (FRW) cosmology. Such model had a beginning at the matter-dominated era. The purposes of our present investigation are two folds: to extend <<life>> of the $\sigma$CDM for earlier times to radiation-dominated era and to take into account variation of the exponential potential via variation of the interaction parameter $\lambda $. We use Markov Chain Monte Carlo (MCMC) procedure to investigate possible values of initial conditions constrained by the measured amount of the dark matter, dark energy and radiation component today. Our analysis includes dark energy contribution to critical density, the ratio of the kinetic and potential energies, deceleration parameter, effective equation of state and evolution of DE equation of state with variation of coupling constant $\lambda $. A comparison with the $\Lambda$CDM model was performed. A new feature of the model is the existence of some values of potential coupling constant, leading to a $\sigma$CDM solution without transit into accelerated expansion epoch.

We present a method to produce mock galaxy catalogues with efficient perturbation theory schemes, which match the number density, power spectra and bispectra in real and in redshift space from N-body simulations. The essential contribution of this work is the way in which we constrain the bias parameters in the PATCHY-code. In addition of aiming at reproducing the two-point statistics, we seek the set of bias parameters, which constrain the univariate halo probability distribution function (PDF) encoding higher-order correlation functions. We demonstrate that halo catalogues based on the same underlying dark matter field with a fix halo number density, and accurately matching the power spectrum (within 2%), can lead to very different bispectra depending on the adopted halo bias model. A model ignoring the shape of the halo PDF can lead to deviations up to factors of 2. The catalogues obtained additionally constraining the shape of the halo PDF can significantly lower the discrepancy in the three-point statistics, yielding closely unbiased bispectra both in real and in redshift space; which are in general compatible with those corresponding to an N-body simulation within 10% (deviating at most up to 20%). Our calculations show that the constant linear bias of ~2 for Luminous Red Galaxy (LRG) like galaxies seen in the power spectrum, mainly comes from sampling halos in high density peaks, choosing a high density threshold rather than from a factor multiplying the dark matter density field. Our method contributes towards an efficient modelling of the halo/galaxy distribution required to estimate uncertainties in the clustering measurements from galaxy redshift surveys. We have also demonstrated that it represents a powerful tool to test various bias models.

We study the statistical inference of the cosmological dark matter density field from non-Gaussian, non-linear and non-Poisson biased distributed tracers. We have implemented a Bayesian posterior sampling computer-code solving this problem and tested it with mock data based on N-body simulations.

We build a physical model for high-redshift Lyman Alpha emitters (LAEs) by coupling state of the art cosmological simulations (GADGET-2) with a dust model and a radiative transfer code (pCRASH). We post-process the cosmological simulation with pCRASH using five different values of the escape fraction of hydrogen ionizing photons (f_esc=0.05,0.25,0.5,0.75,0.95) until reionization is complete, i.e. the average neutral hydrogen fraction drops to <X_HI>~10^-4. Then, the only free-parameter left to match model results to the observed Lya and UV luminosity functions of LAEs at z~6.6 is the relative escape of Lyman Alpha (Lya) and continuum photons from the galactic environment (f_alpha/f_c). We find a three-dimensional degeneracy such that the theoretical model can be reconciled with observations for an IGM Lya transmission <T_alpha>_LAE~38-50% (which translates to <X_HI>~0.5-10^-4 for Gaussian emission lines), f_esc~0.05-0.50 and f_alpha/f_c~0.6-1.8.

This paper is the first in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our primary aim is to test the performance of these cluster mass estimation techniques on a diverse set of models that will increase in complexity. We begin by providing participating methods with data from a simple model that delivers idealised clusters, enabling us to quantify the underlying scatter intrinsic to these mass estimation techniques. The mock catalogue is based on a Halo Occupation Distribution (HOD) model that assumes spherical Navarro, Frenk and White (NFW) haloes truncated at R_200, with no substructure nor colour segregation, and with isotropic, isothermal Maxwellian velocities. We find that, above 10^14 M_solar, recovered cluster masses are correlated with the true underlying cluster mass with an intrinsic scatter of typically a factor of two. Below 10^14 M_solar, the scatter rises as the number of member galaxies drops and rapidly approaches an order of magnitude. We find that richness-based methods deliver the lowest scatter, but it is not clear whether such accuracy may simply be the result of using an over-simplistic model to populate the galaxies in their haloes. Even when given the true cluster membership, large scatter is observed for the majority non-richness-based approaches, suggesting that mass reconstruction with a low number of dynamical tracers is inherently problematic.

We study for the first time a three-dimensional octahedron constellation for a space-based gravitational wave detector, which we call the Octahedral Gravitational Observatory (OGO). With six spacecraft the constellation is able to remove laser frequency noise and acceleration disturbances from the gravitational wave signal without needing LISA-like drag-free control, thereby simplifying the payloads and placing less stringent demands on the thrusters. We generalize LISA's time-delay interferometry to displacement-noise free interferometry (DFI) by deriving a set of generators for those combinations of the data streams that cancel laser and acceleration noise. However, the three-dimensional configuration makes orbit selection complicated. So far, only a halo orbit near the Lagrangian point L1 has been found to be stable enough, and this allows only short arms up to 1400 km. We derive the sensitivity curve of OGO with this arm length, resulting in a peak sensitivity of about $2\times10^{-23}\,\mathrm{Hz}^{-1/2}$ near 100 Hz. We compare this version of OGO to the present generation of ground-based detectors and to some future detectors. We also investigate the scientific potentials of such a detector, which include observing gravitational waves from compact binary coalescences, the stochastic background and pulsars as well as the possibility to test alternative theories of gravity. We find a mediocre performance level for this short-arm-length detector, between those of initial and advanced ground-based detectors. Thus, actually building a space-based detector of this specific configuration does not seem very efficient. However, when alternative orbits that allow for longer detector arms can be found, a detector with much improved science output could be constructed using the octahedron configuration and DFI solutions demonstrated in this paper. (abridged)

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using gravitational waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

We have investigated the distribution and evolution of ~3100 intergalactic HI absorbers with HI column densities log N(HI) = [12.75, 17.0] at 1.9 < z < 3.2, using 18 high resolution, high S/N quasar spectra obtained from the ESO VLT/UVES archive. We used two sets of Voigt profile fitting analysis, one including all the available high-order Lyman lines to obtain reliable HI column densities of saturated lines, and another using only the Ly-alpha lines. There is no significant difference between the results from the two fits. Combining our results with literature data, the mean number density at 0 < z < 4 is not well described by a single power law and strongly suggests that its evolution slows down at z < 1.5 at the high and low column density ranges. We also divided our entire HI absorbers at 1.9 < z < 3.2 into two samples, the unenriched forest and the CIV-enriched forest, depending on whether HI lines are associated with CIV at log N(CIV) > 12.2 within a given velocity range. The entire HI column density distribution function (CDDF) can be described as the combination of these two well-characterised populations which overlap at log N(HI) ~ 15. At log N(HI) < 15, the unenriched forest dominates, showing a similar power-law distribution to the entire forest. The CIV-enriched forest dominates at log N(HI) > 15, with its distribution function proportional to N(HI)^(-1.45). However, it starts to flatten out at lower N(HI), since the enriched forest fraction decreases with decreasing N(HI). The deviation from the power law at log N(HI) = [14, 17] shown in the CDDF for the entire HI sample is a result of combining two different HI populations with a different CDDF shape. The total HI mass density relative to the critical density is Omega(HI) ~ 1.6 x 10^(-6) h^(-1), where the enriched forest accounts for ~40% of Omega(HI).

The 4MOST consortium is currently halfway through a Conceptual Design study for ESO with the aim to develop a wide-field (>3 square degree, goal >5 square degree), high-multiplex (>1500 fibres, goal 3000 fibres) spectroscopic survey facility for an ESO 4m-class telescope (VISTA). 4MOST will run permanently on the telescope to perform a 5 year public survey yielding more than 20 million spectra at resolution R~5000 ({\lambda}=390-1000 nm) and more than 2 million spectra at R~20,000 (395-456.5 nm & 587-673 nm). The 4MOST design is especially intended to complement three key all-sky, space-based observatories of prime European interest: Gaia, eROSITA and Euclid. Initial design and performance estimates for the wide-field corrector concepts are presented. We consider two fibre positioner concepts, a well-known Phi-Theta system and a new R-Theta concept with a large patrol area. The spectrographs are fixed configuration two-arm spectrographs, with dedicated spectrographs for the high- and low-resolution. A full facility simulator is being developed to guide trade-off decisions regarding the optimal field-of-view, number of fibres needed, and the relative fraction of high-to-low resolution fibres. Mock catalogues with template spectra from seven Design Reference Surveys are simulated to verify the science requirements of 4MOST. The 4MOST consortium aims to deliver the full 4MOST facility by the end of 2018 and start delivering high-level data products for both consortium and ESO community targets a year later with yearly increments.

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.

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.

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

We analyse the proximity zone of the intergalactic matter around high-redshift quasars in a cosmological environment. In a box of 64 h-1 Mpc base length we employ dark matter only simulations. For estimating the hydrogen temperature and density distribution we use the effective equation of state. Hydrogen is assumed to be in photoionisation equilibrium with a model background flux which is fit to recent observations of the mean optical depth and transmission flux statistics. At redshifts z = 3, 4, and 4.8, we select model quasar positions at the centre of the 20 most massive halos and 100 less massive halos identified in the simulation. From each assumed quasar position we cast 100 random lines of sight for two box length including the changes in the ionisation fractions by the QSO flux field and derive mock Ly{\alpha} spectra. The proximity effect describes the dependence of the mean normalised optical depth {\xi} = {\tau}eff, QSO/{\tau}eff, Ly{\alpha} as a function of the ratio of the ionisation rate by the QSO and the background field, {\omega} = {\Gamma}QSO/{\Gamma}UVB, i.e. the profile {\xi} = (1 + {\omega}/a)-0.5, where a strength parameter a is introduced. The strength parameter measures the deviation from the theoretical background model and is used to quantify any influence of the environmental density field. We reproduce an unbiased measurement of the proximity effect which is not affected by the host halo mass. The scatter between different lines of sight and different quasar host positions increases with decreasing redshift. Around the host halos, we find only a slight average overdensity in the proximity zone at comoving radii of 1 < rc < 10h-1 Mpc. However, a clear power-law correlation of the strength parameter with the average overdensity in rc is found, showing an overestimation of the ionising background in overdense regions and an underestimation in underdense regions.

We study the influence of the environment on the evolution of galaxies by investigating the luminosity function (LF) of galaxies of different morphological types and colours at different environmental density levels. We construct the LFs separately for galaxies of different morphology (spiral and elliptical) and of different colours (red and blue) using data from the Sloan Digital Sky Survey (SDSS), correcting the luminosities for the intrinsic absorption. We use the global luminosity density field to define different environments, and analyse the environmental dependence of galaxy morphology and colour. The smoothed bootstrap method is used to calculate confidence regions of the derived luminosity functions. We find a strong environmental dependency for the LF of elliptical galaxies. The LF of spiral galaxies is almost environment independent, suggesting that spiral galaxy formation mechanisms are similar in different environments. Absorption by the intrinsic dust influences the bright-end of the LF of spiral galaxies. After attenuation correction, the brightest spiral galaxies are still about 0.5 mag less luminous than the brightest elliptical galaxies, except in the least dense environment, where spiral galaxies dominate the LF at every luminosity. Despite the extent of the SDSS survey, the influence of single rich superclusters is present in the galactic LF of the densest environment.

We present the new parallel version (pCRASH2) of the cosmological radiative transfer code CRASH2 for distributed memory supercomputing facilities. The code is based on a static domain decomposition strategy inspired by geometric dilution of photons in the optical thin case that ensures a favourable performance speed-up with increasing number of computational cores. Linear speed-up is ensured as long as the number of radiation sources is equal to the number of computational cores or larger. The propagation of rays is segmented and rays are only propagated through one sub-domain per time step to guarantee an optimal balance between communication and computation. We have extensively checked pCRASH2 with a standardised set of test cases to validate the parallelisation scheme. The parallel version of CRASH2 can easily handle the propagation of radiation from a large number of sources and is ready for the extension of the ionisation network to species other than hydrogen and helium.

We study the proximity effect in the Ly-a forest around high redshift quasars as a function of redshift and environment employing a set of 3D radiative transfer simulations. The analysis is based on dark matter only simulations at redshifts 3, 4, and 4.9 and, adopting an effective equation of state for the baryonic matter, we infer the HI densities and temperatures in the cosmological box. The UV background (UVB) and additional QSO radiation with Lyman limit flux of L_{\nu LL} = 1e31 and 1e32 erg / Hz s are implemented with a radiative transfer code until an equilibrium configuration is reached. We analyse mock spectra originating at the QSO in the most massive halo, in a random filament and in a void. The proximity effect is studied using flux transmission statistics, in particular with the normalised optical depth. Beyond a radius of r > 1 Mpc / h from the quasar, we measure a transmission profile consistent with geometric dilution of the QSO ionising radiation. A departure from geometric dilution is only seen, when strong absorbers intervene the line-of-sight. The cosmic density distribution around the QSO causes a large scatter in the normalised optical depth. The scatter decreases with increasing redshift and increasing QSO luminosity. The mean proximity effect provides an average signal that is biased by random large scale density enhancements on scales up to r \approx 15 Mpc / h. The distribution of the proximity effect strength provides a measure of the proximity effect along individual lines of sight. It shows a clear maximum almost without an environmental bias. Therefore it can be used for an unbiased estimate of the UVB. Differing spectral energy distributions between the QSO and the UVB modify the profile which can be reasonably well corrected analytically. A few Lyman limit systems have been identified that prevent the detection of the proximity effect due to shadowing.

The measurement of the scale of the baryon acoustic oscillations (BAO) in the galaxy power spectrum as a function of redshift is a promising method to constrain the equation-of-state parameter of the dark energy w. To measure the scale of the BAO precisely, a substantial volume of space must be surveyed. We test whether light-cone effects are important and whether the scaling relations used to compensate for an incorrect reference cosmology are in this case sufficiently accurate. We investigate the degeneracies in the cosmological parameters and the benefits of using the two-dimensional anisotropic power spectrum. Finally, we estimate the uncertainty with which w can be measured by proposed surveys at redshifts of about z=3 and z=1, respectively. In the simulated survey we find that light-cone effects are small and that the simple scaling relations used to correct for the cosmological distortion work fairly well even for large survey volumes. The analysis of the two-dimensional anisotropic power spectra enables an independent determination to be made of the apparent scale of the BAO, perpendicular and parallel to the line of sight. This is essential for two-parameter w-models, such as the redshift-dependent dark energy model w=w_0+(1-a)w_a. Using Planck priors for the matter and baryon density and Delta(H_0)=5% for the Hubble constant, we estimate that the BAO measurements of future surveys around z=3 and z=1 will be able to constrain, independently of other cosmological probes, a constant w to ~ 12% and ~ 11% (68% c.l.), respectively.

We extract groups of galaxies from the SDSS Data Release 5 with the purpose of studying the supercluster-void network and environmental properties of groups therein. We use a modified friends-of-friends (FoF) method with adopted variable linking length in transverse and radial direction to eliminate selection effects and to find reliably as many groups as possible to track the supercluster network. We take into account various selection effects due to the use of a magnitude limited sample. To determine linking length scaling we study the luminosity-density relation in observed groups. We follow the changes in group sizes and mean galaxy number densities within groups when shifting nearby groups to larger distances. As a result we show that the linking length should be a slowly growing function with distance. Our final sample contains 17143 groups in the equatorial, and 33219 groups in the northern part of the DR5 survey with membership N_g >= 2. The group catalogue is available at our web-site http://www.obs.ee/~erik/index.html. Due to a narrow magnitude window in the SDSS the group catalogue based on this survey has been obtained by moderately growing linking length scaling law up to redshift z = 0.12. Above this redshift the scaling law turns down. In the redshift range z=0.12 - 0.2 only the cores are detected. Along with applying weights when calculating luminosities it is possible to use groups for determination of the large-scale luminosity-density field.

We study the spatial distribution of loose groups from the Las Campanas Redshift Survey, comparing it with the supercluster-void network delineated by rich clusters of galaxies. We use density fields and the friends-of-friends algorithm to identify the members of superclusters of Abell clusters among the Las Campanas loose groups. We find that systems of loose groups tend to be oriented perpendicularly to the line-of-sight, and discuss possible reasons for that. We show that loose groups in richer systems (superclusters of Abell clusters) are themselves also richer and more massive than groups in systems without Abell clusters. Our results indicate that superclusters, as high density environments, have a major role in the formation and evolution of galaxy systems.

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.

We have studied the properties of Las Campanas Loose Groups in the neighbourhood of rich (Abell, APM and X-ray) clusters of galaxies. These loose groups show strong evidence of segregation measured in terms of the group richness and the group velocity dispersion: loose groups in the neighbourhood of a rich cluster are typically 2.5 times more massive and 1.6 times more luminous than groups on average, and these loose groups havevelocity dispersions 1.3 times larger than groups on average. This is evidence that the large-scale gravitational field causing the formation of rich clusters enhances the evolution of neighbouring poor systems, a phenomenon recently established in numerical simulations of group and cluster formation.

We have determined the mass function of loose groups of galaxies in the Las Campanas Redshift Survey. Loose groups of galaxies in the LCRS range in mass from $M \sim 10^{12} {\rm M}_{\sun}$ to $10^{15} {\rm M}_{\sun}$. We find that the sample is almost complete for masses in the interval $5\cdot 10^{13}-8\cdot 10^{14} {\rm M}_{\sun}$. Comparison of the observed mass function with theoretical mass functions obtained from N-body simulations shows good agreement with a CDM model with the parameters $\Omega_m = 0.3$, $\Omega_{\Lambda} = 0.7$ and the amplitude of perturbations about $\sigma_8=0.78-0.87$. For smaller masses the mass function of LCRS loose groups flattens out, differing considerably from the group mass function found by Girardi and Giuricin (2000) and from mass functions obtained by numerical simulations.

We show that the amplitude of the second acoustic peak in the newly released BOOMERANG-98 and MAXIMA-I data is compatible with the standard primordial nucleosynthesis and with the locally broken-scale-invariant matter power spectrum suggested by recent measurements of the power spectrum in the range $20 - 200$ Mpc. If the slope of matter density perturbations on large scales is $n \approx 1$, the Hubble constant is $0.5< h < 0.75$, and r.m.s. mass fluctuations at 8 Mpc are $0.65 \le \sigma_8 \le 0.75$, then for a Universe approximately 14 Gyr old our best fit within the nucleosynthesis bound $\Omega_Bh^2= 0.019 \pm 0.0024$ requires $0.3 \le \Omega_m \le 0.5$. Cluster abundances further constraint the matter density to be $\Omega_m \approx 0.3$. The CMB data alone are not able to determine the detailed form of the matter power spectrum in the range $0.03<k<0.06$ \hmpc where deviations from the scale-invariant spectrum are expected to be most significant, but they do not contradict the existence of the previously claimed peak at $k\sim 0.05$ \hmpc, and a depression at $k\sim 0.035$ \hmpc.

We propose a $\Lambda$-inflation model which explains a large fraction of the COBE signal by cosmic gravitational waves. The primordial density perturbations fulfil both the contraints of large-scale microwave background and galaxy cluster normalization. The model is tested against the galaxy cluster power spectrum and the high-multipole angular CMB anisotropy.

The simulated matter distribution on large scales is studied using core-sampling, cluster analysis, inertia tensor analysis, and minimal spanning tree techniques. Seven simulations in large boxes for five cosmological models with COBE normalized CDM-like power spectra are studied. The wall-like Super Large Scale Structure with parameters similar to the observed one is found for the OCDM and LCDM models with Omega_m h = 0.3 & 0.245. In these simulations, the rich structure elements with a typical value for the largest extension of \~(30 - 50) h^-1 Mpc incorporate ~40% of matter with overdensity of about 10 above the mean. These rich elements are formed due to the anisotropic nonlinear compression of sheets with original size of ~(15 - 25) h^{-1} Mpc. They surround low-density regions with a typical diameter ~(50 - 70) h^{-1} Mpc. The statistical characteristics of these structures are found to be approximately consistent with observations and theoretical expectations. The cosmological models with higher matter density \Omega_m=1 in CDM with Harrison-Zeldovich or tilted power spectra cannot reproduce the characteristics of the observed galaxy distribution due to the very strong disruption of the rich structure elements. Another model with a broken scale invariant initial power spectrum (BCDM) shows not enough matter concentration in the rich structure elements.

The large supercluster structures obvious in recent galaxy redshift surveys are quantified using an one-dimensional cluster analysis (core sampling) and a three-dimensional cluster analysis based on the minimal spanning tree. The comparison with the LCRS reveals promising stable results. At a mean overdensity of about ten, the supercluster systems form huge wall-like structures comprising about 40% of all galaxies. The overdense clusters have a low mean transverse velocity dispersion of about 400 km/s, i.e. they look quite narrow in redshift space. We performed N-body simulations with large box sizes for six cosmological scenarios. The quantitative analysis shows that the observed structures can be understood best in low density models with $\Omega_m \le 0.5$ with or without a cosmological constant.

Using 3D hydrodynamical simulations of galaxy formation with supernova feedback and a multiphase medium, we derive theoretical relations analogous to the observed Tully-Fisher (TF) relations in various photometric bands. This paper examines the influence of self-regulation mechanisms including supernova feedback on galaxy luminosities and the TF relation in three cosmological scenarios (CDM, \Lambda CDM and BSI (broken scale invariance)). The galaxy catalogs derived from our hydrodynamical simulations lead to an acceptably small scatter in the theoretical TF relation amounting to \Delta M =0.2-0.4 in the I band, and increasing by 0.1 magnitude from the I-band to the B-band. Our results give strong evidence that the tightness of the TF relation cannot be attributed to supernova feedback alone. However, although eliminating supernova feedback hardly affects the scatter, it does influence the slope of the TF relation quite sensitively. With supernova feedback, L \propto V_c^{3-3.5} (depending on the strength of supernova feedback). Without it, L \propto V_c^{2} as predicted by the virial theorem with constant M/L. The luminosity functions in the B and K bands are quite sensitive to supernova feedback at the faint end studied here. We find that the faint end of the B-band luminosity function (-18 \leq M_B \leq -15) has a slope that is steeper than the Stromlo-APM estimate, but in rough agreement with the recent ESO Slice Project estimates.

We present a study of the shape, size, and spatial orientation of superclusters of galaxies. Approximating superclusters by triaxial ellipsoids we show that superclusters are flattened, triaxial objects. We find that there are no spherical superclusters. The sizes of superclusters grow with their richness: the median semi-major axis of rich and poor superclusters (having at least 8 and less than 8 member clusters) is 42 and 31 Mpc, respectively. Similarly, the median semi-minor axis is 12 and 5 Mpc for rich and poor superclusters. The spatial orientation of superclusters, as determined from the axes of ellipsoids, is nearly random. We do not detect any preferable orientation of superclusters, neither with respect to the line of sight, nor relative to some other outstanding feature in the large scale structure, nor with respect to the directions of principal axes of adjacent superclusters.

The large scale gravitational potential distribution and its influence on the large-scale matter clustering is considered on the basis of six simulations. It is found that the mean separation between zero levels of the potential along random straight lines coincides with the theoretical expectations, but it scatters largely. A strong link of the initial potential and the structure evolution is shown. It is found that the under-dense and over-dense regions correlate with regions of positive and negative gravitational potential at large redshifts. The over-dense regions arise due to a slow matter flow into the negative potential regions, where more pronounced non-linear structures appear. Such regions are related to the formation of huge super-large scale structures seen in the galaxy distribution.

We calculate temperature anisotropies of the cosmic microwave background (CMB) for several initial power spectra of density perturbations with a built-in scale suggested by recent optical data on the spatial distribution of rich clusters of galaxies. Using cosmological models with different values of spectral index, baryon fraction, Hubble constant and cosmological constant, we compare the calculated radiation power spectrum with the CMB temperature anisotropies measured by the Saskatoon experiment. We show that spectra with a sharp peak at 120 h^{-1} Mpc are in agreement with the Saskatoon data. The combined evidence from cluster and CMB data favours the presence of a peak and a subsequent break in the initial matter power spectrum. Such feature is similar to the prediction of an inflationary model where an inflaton field is evolving through a kink in the potential.

We investigate properties of the correlation function of clusters of galaxies using geometrical models. On small scales the correlation function depends on the shape and the size of superclusters. On large scales it describes the geometry of the distribution of superclusters. If superclusters are distributed randomly then the correlation function on large scales is featureless. If superclusters and voids have a tendency to form a regular lattice then the correlation function on large scales has quasi-regularly spaced maxima and minima of decaying amplitude; i.e., it is oscillating. The period of oscillations is equal to the step size of the grid of the lattice. We calculate the power spectrum for our models and compare the geometrical information of the correlation function with other statistics. We find that geometric properties (the regularity of the distribution of clusters on large scales) are better quantified by the correlation function. We also analyse errors in the correlation function and the power spectrum by generating random realizations of models and finding the scatter of these realizations.

We use rich clusters of galaxies in the Northern and Southern Galactic hemispheres up to a redshift z=0.12 to determine the cluster correlation function. We show that superclusters of galaxies and voids between them form a moderately regular network. As a result the correlation function determined for clusters located in rich superclusters oscillates: it has a series of regularly spaced secondary maxima and minima. The scale of the supercluster-void network, determined from the period of oscillations, is P=115 Mpc. The amplitude of oscillations is about 3 times larger than the estimated error. We compare the observed cluster correlation function with similar functions derived for popular models of structure formation. We find that the production of the observed cluster correlation function in any model with a smooth transition of the power spectrum from a Harrison-Zeldovich regime with positive spectral index on long wavelengths to a negative spectral index on short wavelengths is highly unlikely. The power spectrum must have an extra peak located at the wavelength equal to the period of oscillations of the correlation function.

Using a new compilation of available data on galaxy clusters and superclusters we present evidence for a quasiregular three-dimensional network of rich superclusters and voids, with the regions of high density separated by about 120 Mpc. We calculate the power spectrum for clusters of galaxies; it has a peak on the wavelength equal to the step of the network; the excess in the amplitude of the spectrum over that of the cold dark matter model is by a factor of 1.4. The probability that the spectrum can be formed within the framework of the standard cosmogony is very small. If the cluster distribution reflects the distribution of all matter (luminous and dark), then there must exists some hithero unknown process that produces regular structure on large scales.

We have simulated the formation of large-scale structure arising from COBE-normalized spectra computed by convolving a primordial double-inflation perturbation spectrum with the CDM transfer function. Due to the broken scale invariance ('BSI') characterizing the primordial perturbation spectrum, this model has less small-scale power than the (COBE-normalized) standard CDM model. The particle-mesh code (with $512^3$ cells and $256^3$ particles) includes a model for thermodynamic evolution of baryons in addition to the usual gravitational dynamics of dark matter. It provides an estimate of the local gas temperature. In particular, our galaxy-finding procedure seeks peaks in the distribution of gas that has cooled. It exploits the fact that ``cold" particles trace visible matter better than average and thus provides a natural biasing mechanism. The basic picture of large-scale structure formation in the BSI model is the familiar hierarchical clustering scenario. We obtain particle in cell statistics, the galaxy correlation function, the cluster abundance and the cluster-cluster correlation function and statistics for large and small scale velocity fields. We also report here on a semi-quantitative study of the distribution of gas in different temperature ranges. Based on confrontation with observations and comparison with standard CDM, we conclude that the BSI scenario could represent a promising modification of the CDM picture capable of describing many details of large-scale structure formation.