From the Voigt profile fitting analysis of 183 intervening CIV systems at 1.7 < z < 3.3 in 23 high-quality UVES/VLT and HIRES/Keck QSO spectra, we find that a majority of CIV systems (~75%) display a well-characterised scaling relation between integrated column densities of HI and CIV with a negligible redshift evolution, when column densities of all the HI and CIV components are integrated within a given (-150, +150) km/sec range centred at the CIV flux minimum. The integrated CIV column density N(CIV, sys) increases with N(HI, sys) at log N(HI, sys) = 14.0--15.5 and log N(CIV, sys) = 11.8--14.0, then becomes almost independent of N(HI, sys) at log N(HI, sys) > 16, with a large scatter: at log N(HI, sys) = 14--22, log N(CIV, sys) = C1 / (log(NHI, sys) + C2) + C3, with C1 = -1.90+0.55, C2 = -14.11+0.19 and C3 = 14.76+0.17, respectively. The steep (flat) part is dominated by SiIV-free (SiIV-enriched) CIV systems. Extrapolating the N(HI, sys)-N(CIV, sys) relation implies that most absorbers with log N(HI) < 14 are virtually CIV-free. The N(HI, sys)-N(CIV, sys) relation does not hold for individual components, clumps or the integration velocity range less than +-100 km/sec. It is expected if the line-of-sight extent of CIV is smaller than HI and N(CIV, sys) decreases more rapidly than N(HI, sys) at the larger impact parameter, regardless of the location of the HI+CIV gas in the IGM filaments or in the intervening galactic halos.

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.

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

We present the online MultiDark Database -- a Virtual Observatory-oriented, relational database for hosting various cosmological simulations. The data is accessible via an SQL (Structured Query Language) query interface, which also allows users to directly pose scientific questions, as shown in a number of examples in this paper. Further examples for the usage of the database are given in its extensive online documentation (www.multidark.org). The database is based on the same technology as the Millennium Database, a fact that will greatly facilitate the usage of both suites of cosmological simulations. The first release of the MultiDark Database hosts two 8.6 billion particle cosmological N-body simulations: the Bolshoi (250/h Mpc simulation box, 1/h kpc resolution) and MultiDark Run1 simulation (MDR1, or BigBolshoi, 1000/h Mpc simulation box, 7/h kpc resolution). The extraction methods for halos/subhalos from the raw simulation data, and how this data is structured in the database are explained in this paper. With the first data release, users get full access to halo/subhalo catalogs, various profiles of the halos at redshifts z=0-15, and raw dark matter data for one time-step of the Bolshoi and four time-steps of the MultiDark simulation. Later releases will also include galaxy mock catalogs and additional merging trees for both simulations as well as new large volume simulations with high resolution. This project is further proof of the viability to store and present complex data using relational database technology. We encourage other simulators to publish their results in a similar manner.

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