We extend the random-walk model of Vitvitska et al. for predicting the spins of dark matter halos from their merger histories. Using updated merger rates, orbital parameter distributions, and N-body constraints we show that this model can accurately reproduce the distribution of spin parameters measured in N-body simulations when we include a weak correlation between the spins of halos and the angular momenta of infalling subhalos. We further show that this model is in approximate agreement with the correlation of the spin magnitude over time as determined from N-body simulations, while it slightly underpredicts the correlation in the direction of the spin vector measured from the same simulations. This model is useful for predicting spins from merger histories derived from non-N-body sources, thereby circumventing the need for very high resolution simulations to permit accurate measurements of spins. It may be particularly relevant to modeling systems which accumulate angular momentum from halos over time (such as galactic discs) - we show that this model makes small but significant changes in the distribution of galactic disc sizes computed using the Galacticus semi-analytic galaxy formation model.

We use a modified version of the Peak-Patch excursion set formalism to compute the mass and size distribution of QCD axion miniclusters from a fully non-Gaussian initial density field obtained from numerical simulations of axion string decay. We find strong agreement with N-Body simulations at significantly lower computational cost. We employ a spherical collapse model, and provide fitting functions for the modified barrier in the radiation era. The halo mass function at $z=99$ has a power-law distribution $M^{-0.68}$ for small masses and $M^{-0.35}$ for large masses within the range $10^{-15}\lesssim M\lesssim 10^{-9}M_{\odot}$, with all masses scaling as $(m_a/50\mu\mathrm{eV})^{-0.5}$. We construct merger trees to estimate the collapse redshift and concentration mass relation, $C(M)$, which is well described using analytical results from the initial power spectrum and linear growth. Using the calibrated analytic results to extrapolate to $z=0$, our method predicts a mean concentration $C\sim \mathcal{O}(\text{few})\times10^4$. The low computational cost of our method makes future investigation of the statistics of rare, dense miniclusters easy to achieve.

We present a same-level comparison of the most prominent inversion methods for the reconstruction of the matter density field in the quasi-linear regime from the Ly$\alpha$ forest flux. Moreover, we present a pathway for refining the reconstruction in the framework of numerical optimization. We apply this approach to construct a novel hybrid method. The methods which are used so far for matter reconstructions are the Richardson-Lucy algorithm, an iterative Gauss-Newton method and a statistical approach assuming a one-to-one correspondence between matter and flux. We study these methods for high spectral resolutions such that thermal broadening becomes relevant. The inversion methods are compared on synthetic data (generated with the lognormal approach) with respect to their performance, accuracy, their stability against noise, and their robustness against systematic uncertainties. We conclude that the iterative Gauss-Newton method offers the most accurate reconstruction, in particular at small S/N, but has also the largest numerical complexity and requires the strongest assumptions. The other two algorithms are faster, comparably precise at small noise-levels, and, in the case of the statistical approach, more robust against inaccurate assumptions on the thermal history of the intergalactic medium (IGM). We use these results to refine the statistical approach using regularization. Our new approach has low numerical complexity and makes few assumptions about the history of the IGM, and is shown to be the most accurate reconstruction at small S/N, even if the thermal history of the IGM is not known. Our code will be made publicly available under https://github.com/hmuellergoe/reglyman.

We report the discovery of 10-kpc scale [CII] 158um halos surrounding star-forming galaxies in the early Universe. We choose deep ALMA data of 18 galaxies each with a star-formation rate of ~ 10-70 Msun with no signature of AGN whose [CII] lines are individually detected at z=5.153-7.142, and conduct stacking of the [CII] lines and dust-continuum in the uv-visibility plane. The radial profiles of the surface brightnesses show a 10-kpc scale [CII] halo at the 9.2sigma level significantly extended more than the HST stellar continuum data by a factor of ~5 on the exponential-profile basis, as well as the dust continuum. We also compare the radial profiles of [CII] and Lya halos universally found in star-forming galaxies at this epoch, and find that the scale lengths agree within the 1sigma level. While two independent hydrodynamical zoom-in simulations match the dust and stellar continuum properties, the simulations cannot reproduce the extended [CII] line emission. The existence of the extended [CII] halo is the evidence of outflow remnants in the early galaxies and suggest that the outflows may be dominated by cold-mode outflows expelling the neutral gas.

A tight relation between the [CII]158$\mu$m line luminosity and star formation rate is measured in local galaxies. At high redshift ($z>5$), though, a much larger scatter is observed, with a considerable (15-20\%) fraction of the outliers being [CII]-deficient. Moreover, the [CII] surface brightness ($\Sigma_{\rm CII}$) of these sources is systematically lower than expected from the local relation. To clarify the origin of such [CII]-deficiency we have developed an analytical model that fits local [CII] data, and has been validated against radiative transfer simulations performed with CLOUDY. The model predicts an overall increase of $\Sigma_{\rm CII}$ with the surface star formation rate ($\Sigma_*$). However, for $\Sigma_* > 1 M_\odot~{\rm yr}^{-1}~{\rm kpc}^{-2}$, $\Sigma_{\rm CII}$ saturates. We conclude that underluminous [CII] systems can result from a combination of three factors: (a) large upward deviations from the Kennicutt-Schmidt relation ($\kappa_s \gg 1$), parameterized by the "burstiness" parameter $\kappa_s$; (b) low metallicity; (c) low gas density, at least for the most extreme sources (e.g. CR7). Observations of [CII] emission alone cannot break the degeneracy among the above three parameters; this requires additional information coming from other emission lines (e.g. [OIII]88$\mu$m, CIII]1909A, CO lines). Simple formulae are given to interpret available data for low and high-$z$ galaxies.

Complex radiative transfer (RT) of the Lyman-$\alpha$ photons poses a theoretical challenge to galaxy surveys which infer the large-scale structure with Lyman-$\alpha$ emitters (LAEs). Guided by RT simulations, prior studies investigated the impact of RT on the large-scale LAE clustering, and claimed that RT induces a selection effect which results in an anisotropic distortion even in real space but in an otherwise negligible effect in redshift space. However, our previous study, which relies on a full radiative transfer code run on the Illustris simulations, shows that the anisotropic selection effect was drastically reduced with higher spatial resolution. Adopting the same simulation framework, we further study the impact of RT on the LAE clustering in redshift space. Since we measure LAE's radial position through a spectral peak of a Lyman-$\alpha$ emission, the frequency shift due to RT contaminates the redshift measurement and hence the inferred radial position in redshift space. We demonstrate that this additional RT offset suppresses the LAE clustering along the line of sight, which can be interpreted as a novel Fingers-of-God (FoG) effect. To assess the FoG effect, we develop a theoretical framework with particular emphasis on its connection with the underlying one-point velocity distribution which simultaneously takes into account the RT offset as well as the peculiar velocity that is commonly studied in the context of the Redshift Space Distortion (RSD). Although our findings strongly encourage a more careful RSD modeling in LAE surveys, we also seek a method to mitigate the additional FoG effect due to RT by making use of other information in a Lyman-$\alpha$ spectrum.

The intrinsic strength of the Ly$\alpha$ line in young, star-forming systems makes it a special tool for studying high-redshift galaxies. However, interpreting observations remains challenging due to the complex radiative transfer involved. Here, we combine state-of-the-art hydrodynamical simulations of 'Althaea', a prototypical Lyman Break Galaxy (LBG, stellar mass $M_{\star}$$\simeq$ $10^{10}{\rm M}_{\odot})$ at $z=7.2$, with detailed radiative transfer computations of dust/continuum, [CII] 158 $\mu$m, and Ly$\alpha$ to clarify the relation between the galaxy properties and its Ly$\alpha$ emission. Althaea exhibits low ($f_\alpha< 1\%$) Ly$\alpha$ escape fractions and Equivalent Widths, EW $\lesssim 6$ Angstrom for the simulated lines of sight, with a large scatter. The correlation between escape fraction and inclination is weak, as a result of the rather chaotic structure of high-redshift galaxies. Low $f_\alpha$ values persist even if we artificially remove neutral gas around star forming regions to mimick the presence of HII regions. The high attenuation is primarily caused by dust clumps co-located with young stellar clusters. We can turn Althaea into a Lyman Alpha Emitter (LAE) only if we artificially remove dust from the clumps, yielding EWs up to $22$ Angstrom. Our study suggests that the LBG-LAE duty-cycle required by recent clustering measurements poses the challenging problem of a dynamically changing dust attenuation. Finally, we find an anti-correlation between the magnitude of Ly$\alpha$-[CII] line velocity shift and Ly$\alpha$ luminosity.

We study the kinematical properties of galaxies in the Epoch of Reionization via the [CII] 158$\mu$m line emission. The line profile provides information on the kinematics as well as structural properties such as the presence of a disk and satellites. To understand how these properties are encoded in the line profile, first we develop analytical models from which we identify disk inclination and gas turbulent motions as the key parameters affecting the line profile. To gain further insights, we use "Althaea", a highly-resolved ($30\, \rm pc$) simulated prototypical Lyman Break Galaxy, in the redshift range $z = 6-7$, when the galaxy is in a very active assembling phase. Based on morphology, we select three main dynamical stages: I) Merger , II) Spiral Disk, and III) Disturbed Disk. We identify spectral signatures of merger events, spiral arms, and extra-planar flows in I), II), and III), respectively. We derive a generalised dynamical mass vs. [CII]-line FWHM relation. If precise information on the galaxy inclination is (not) available, the returned mass estimate is accurate within a factor $2$ ($4$). A Tully-Fisher relation is found for the observed high-$z$ galaxies, i.e. $L_{\rm[CII]}\propto (FWHM)^{1.80\pm 0.35}$ for which we provide a simple, physically-based interpretation. Finally, we perform mock ALMA simulations to check the detectability of [CII]. When seen face-on, Althaea is always detected at $> 5\sigma$; in the edge-on case it remains undetected because the larger intrinsic FWHM pushes the line peak flux below detection limit. This suggests that some of the reported non-detections might be due to inclination effects.

We study the formation and evolution of a sample of Lyman Break Galaxies in the Epoch of Reionization by using high-resolution ($\sim 10 \,{\rm pc}$), cosmological zoom-in simulations part of the SERRA suite. In SERRA, we follow the interstellar medium (ISM) thermo-chemical non-equilibrium evolution, and perform on-the-fly radiative transfer of the interstellar radiation field (ISRF). The simulation outputs are post-processed to compute the emission of far infrared lines ([CII], [NII], and [OIII]). At $z=8$, the most massive galaxy, `Freesia', has an age $t_\star \simeq 409\,{\rm Myr}$, stellar mass $M_{\star} \simeq 4.2\times 10^9 {\rm M}_{\odot}$, and a star formation rate ${\rm SFR} \simeq 11.5\,{\rm M}_{\odot}{\rm yr}^{-1}$, due to a recent burst. Freesia has two stellar components (A and B) separated by $\simeq 2.5\, {\rm kpc}$; other 11 galaxies are found within $56.9 \pm 21.6 \, {\rm kpc}$. The mean ISRF in the Habing band is $G = 7.9\, G_0$ and is spatially uniform; in contrast, the ionisation parameter is $U = 2^{+20}_{-2} \times 10^{-3}$, and has a patchy distribution peaked at the location of star-forming sites. The resulting ionising escape fraction from Freesia is $f_{\rm esc}\simeq 2\%$. While [CII] emission is extended (radius 1.54 kpc), [OIII] is concentrated in Freesia-A (0.85 kpc), where the ratio $\Sigma_{\rm [OIII]}/\Sigma_{\rm [CII]} \simeq 10$. As many high-$z$ galaxies, Freesia lies below the local [CII]-SFR relation. We show that this is the general consequence of a starburst phase (pushing the galaxy above the Kennicutt-Schmidt relation) which disrupts/photodissociates the emitting molecular clouds around star-forming sites. Metallicity has a sub-dominant impact on the amplitude of [CII]-SFR deviations.

The most massive and luminous galaxies in the Universe serve as powerful probes to study the formation of structure, the assembly of mass, and cosmology. However, their detailed formation and evolution is still barely understood. Here we extract a sample of massive mock galaxies from the semi-analytical model of galaxy formation (SAM) GALACTICUS from the MultiDark-Galaxies, by replicating the CMASS photometric selection from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). The comparison of the GALACTICUS CMASS-mock with BOSS-CMASS data allows us to explore different aspects of the massive galaxy population at $0.5<z<0.6$, including the galaxy-halo connection and the galaxy clustering. We find good agreement between our modelled galaxies and observations regarding the galaxy-halo connection, but our CMASS-mock over-estimates the clustering amplitude of the 2-point correlation function, due to a smaller number density compared to BOSS, a lack of blue objects, and a small intrinsic scatter in stellar mass at fixed halo mass of $<0.1$ dex. To alleviate this problem, we construct an alternative mock catalogue mimicking the CMASS colour-magnitude distribution by randomly down-sampling the SAM catalogue. This CMASS-mock reproduces the clustering of CMASS galaxies within 1$\sigma$ and shows some environmental dependency of star formation properties that could be connected to the quenching of star formation and the assembly bias.

Lyman-$\alpha$ emitters (LAEs) are a promising probe of the large-scale structure at high redshift, $z\gtrsim 2$. In particular, the Hobby-Eberly Telescope Dark Energy Experiment aims at observing LAEs at 1.9 $<z<$ 3.5 to measure the Baryon Acoustic Oscillation (BAO) scale and the Redshift-Space Distortion (RSD). However, Zheng et al. (2011) pointed out that the complicated radiative transfer (RT) of the resonant Lyman-$\alpha$ emission line generates an anisotropic selection bias in the LAE clustering on large scales, $s\gtrsim 10$ Mpc. This effect could potentially induce a systematic error in the BAO and RSD measurements. Also, Croft et al. (2016) claims an observational evidence of the effect in the Lyman-$\alpha$ intensity map, albeit statistically insignificant. We aim at quantifying the impact of the Lyman-$\alpha$ RT on the large-scale galaxy clustering in detail. For this purpose, we study the correlations between the large-scale environment and the ratio of an apparent Lyman-$\alpha$ luminosity to an intrinsic one, which we call the `observed fraction', at $2<z<6$. We apply our Lyman-$\alpha$ RT code by post-processing the full Illustris simulations. We simply assume that the intrinsic luminosity of the Lyman-$\alpha$ emission is proportional to the star formation rate of galaxies in Illustris, yielding a sufficiently large sample of LAEs to measure the anisotropic selection bias. We find little correlations between large-scale environment and the observed fraction induced by the RT, and hence a smaller anisotropic selection bias than what was claimed by Zheng et al. (2011). We argue that the anisotropy was overestimated in the previous work due to the insufficient spatial resolution: it is important to keep the resolution such that it resolves the high density region down to the scale of the interstellar medium, $\sim1$ physical kpc. (abridged)

The recent discovery of dusty galaxies well into the Epoch of Reionization (redshift $z>6$) poses challenging questions about the properties of the interstellar medium in these pristine systems. By combining state-of-the-art hydrodynamic and dust radiative transfer simulations, we address these questions focusing on the recently discovered dusty galaxy A2744_YD4 ($z=8.38$, Laporte et al. 2017}). We show that we can reproduce the observed spectral energy distribution (SED) only using different physical values with respect to the inferred ones by Laporte et al(2017), i.e. a star formation rate of $\mathrm{SFR} = 78\;\rm M_\odot \rm yr^{-1}$, a factor $\approx 4$ higher than deduced from simple Spectral Energy Distribution fitting. In this case we find: (a) dust attenuation (corresponding to $\tau_V=1.4$) is consistent with a Milky Way extinction curve; (b) the dust-to-metal ratio is low, $f_\mathrm{d} \sim 0.08$, implying that early dust formation is rather inefficient; (c) the luminosity-weighted dust temperature is high, $T_d=91\pm 23\, \rm K$, as a result of the intense ($\approx 100\times$ MW) interstellar radiation field; (d) due to the high $T_d$, the ALMA Band 7 detection can be explained by a limited dust mass, $M_d=1.6\times 10^6 $M$_\odot$. Finally, the high dust temperatures might solve the puzzling low infrared excess recently deduced for high-$z$ galaxies from the IRX-$\beta$ relation.

We present the public release of the MultiDark-Galaxies: three distinct galaxy catalogues derived from one of the Planck cosmology MultiDark simulations (i.e. MDPL2, with a volume of (1 Gpc/$h$)$^{3}$ and mass resolution of $1.5 \times 10^{9} M_{\odot}/h$) by applying the semi-analytic models GALACTICUS, SAG, and SAGE to it. We compare the three models and their conformity with observational data for a selection of fundamental properties of galaxies like stellar mass function, star formation rate, cold gas fractions, and metallicities - noting that they sometimes perform differently reflecting model designs and calibrations. We have further selected galaxy subsamples of the catalogues by number densities in stellar mass, cold gas mass, and star formation rate in order to study the clustering statistics of galaxies. We show that despite different treatment of orphan galaxies, i.e. galaxies that lost their dark-matter host halo due to the finite mass resolution of the N-body simulation or tidal stripping, the clustering signal is comparable, and reproduces the observations in all three models - in particular when selecting samples based upon stellar mass. Our catalogues provide a powerful tool to study galaxy formation within a volume comparable to those probed by on-going and future photometric and redshift surveys. All model data consisting of a range of galaxy properties - including broad-band SDSS magnitudes - are publicly available.

We study the CO line luminosity ($L_{\rm CO}$), the shape of the CO Spectral Line Energy Distribution (SLED), and the value of the CO-to-$\rm H_2$ conversion factor in galaxies in the Epoch of Reionization (EoR). To this aim, we construct a model that simultaneously takes into account the radiative transfer and the clumpy structure of giant molecular clouds (GMCs) where the CO lines are excited. We then use it to post-process state-of-the-art zoomed, high resolution ($30\, \rm{pc}$), cosmological simulation of a main-sequence ($M_{*}\approx10^{10}\, \rm{M_{\odot}}$, $SFR\approx 100\,\rm{M_{\odot}\, yr^{-1}}$) galaxy, "Alth{\ae}a", at $z\approx6$. We find that the CO emission traces the inner molecular disk ($r\approx 0.5 \,\rm{kpc}$) of Alth{\ae}a with the peak of the CO surface brightness co-located with that of the [CII] 158$\rm \mu m$ emission. Its $L_{\rm CO(1-0)}=10^{4.85}\, \rm{L_{\odot}}$ is comparable to that observed in local galaxies with similar stellar mass. The high ($\Sigma_{gas} \approx 220\, \rm M_{\odot}\, pc^{-2}$) gas surface density in Alth{\ae}a, its large Mach number (\mach$\approx 30$), and the warm kinetic temperature ($T_{k}\approx 45 \, \rm K$) of GMCs yield a CO SLED peaked at the CO(7-6) transition, i.e. at relatively high-$J$, and a CO-to-$\rm H_2$ conversion factor $\alpha_{\rm CO}\approx 1.5 \, \rm M_{\odot} \rm (K\, km\, s^{-1}\, pc^2)^{-1} $ lower than that of the Milky Way. The ALMA observing time required to detect (resolve) at 5$\sigma$ the CO(7-6) line from galaxies similar to Alth{\ae}a is $\approx13$ h ($\approx 38$ h).

Major mergers are considered to be a significant source of turbulence in clusters. We performed a numerical simulation of a major merger event using nested-grid initial conditions, adaptive mesh refinement, radiative cooling of primordial gas, and a homogeneous ultraviolet background. By calculating the microscopic viscosity on the basis of various theoretical assumptions and estimating the Kolmogorov length from the turbulent dissipation rate computed with a subgrid-scale model, we are able to demonstrate that most of the warm-hot intergalactic medium can sustain a fully turbulent state only if the magnetic suppression of the viscosity is considerable. Accepting this as premise, it turns out that ratios of turbulent and thermal quantities change only little in the course of the merger. This confirms the tight correlations between the mean thermal and non-thermal energy content for large samples of clusters in earlier studies, which can be interpreted as second self-similarity on top of the self-similarity for different halo masses. Another long-standing question is how and to which extent turbulence contributes to the support of the gas against gravity. From a global perspective, the ratio of turbulent and thermal pressures is significant for the clusters in our simulation. On the other hand, a local measure is provided by the compression rate, i.e. the growth rate of the divergence of the flow. Particularly for the intracluster medium, we find that the dominant contribution against gravity comes from thermal pressure, while compressible turbulence effectively counteracts the support. For this reason it appears to be too simplistic to consider turbulence merely as an effective enhancement of thermal energy.

In the context of structure formation with ultralight axion dark matter, we offer an alternative explanation for the mass relation of solitonic cores and their host halos observed in numerical simulations. Our argument is based entirely on the mass gain that occurs during major mergers of binary cores and largely independent of the initial core-halo mass relation assigned to hosts that have just collapsed. We find a relation between the halo mass $M_h$ and corresponding core mass $M_c$, $M_c\propto M_h^{2\beta-1}$, where $(1-\beta)$ is the core mass loss fraction. Following the evolution of core masses in stochastic merger trees, we find empirical evidence for our model. Our results are useful for statistically modeling the effects of dark matter cores on the properties of galaxies and their substructures in axion dark matter cosmologies.

We derive the halo mass function (HMF) for fuzzy dark matter (FDM) by solving the excursion set problem explicitly with a mass-dependent barrier function, which has not been done before. We find that compared to the naive approach of the Sheth--Tormen HMF for FDM, our approach has a higher cut off mass and the cut off mass changes less strongly with redshifts. Using merger trees constructed with a modified version of the Lacey & Cole formalism that accounts for suppressed small scale power and the scale-dependent growth of FDM haloes and the semi-analytic GALACTICUS code, we study the statistics of halo substructure including the effects from dynamical friction and tidal stripping. We find that if the dark matter is a mixture of cold dark matter (CDM) and FDM, there will be a suppression on the halo substructure on small scales which may be able to solve the Missing Satellites Problem faced by the pure CDM model. The suppression becomes stronger with increasing FDM fraction or decreasing FDM mass. Thus, it may be used to constrain the FDM model.

We compare the physical and morphological properties of z ~ 2 Lyman-alpha emitting galaxies (LAEs) identified in the HETDEX Pilot Survey and narrow band studies with those of z ~ 2 optical emission line galaxies (oELGs) identified via HST WFC3 infrared grism spectroscopy. Both sets of galaxies extend over the same range in stellar mass (7.5 < logM < 10.5), size (0.5 < R < 3.0 kpc), and star-formation rate (~1 < SFR < 100). Remarkably, a comparison of the most commonly used physical and morphological parameters -- stellar mass, half-light radius, UV slope, star formation rate, ellipticity, nearest neighbor distance, star formation surface density, specific star formation rate, [O III] luminosity, and [O III] equivalent width -- reveals no statistically significant differences between the populations. This suggests that the processes and conditions which regulate the escape of Ly-alpha from a z ~ 2 star-forming galaxy do not depend on these quantities. In particular, the lack of dependence on the UV slope suggests that Ly-alpha emission is not being significantly modulated by diffuse dust in the interstellar medium. We develop a simple model of Ly-alpha emission that connects LAEs to all high-redshift star forming galaxies where the escape of Ly-alpha depends on the sightline through the galaxy. Using this model, we find that mean solid angle for Ly-alpha escape is 2.4+/-0.8 steradians; this value is consistent with those calculated from other studies.

We present simulations of Lyman-Alpha radiation transfer in an isolated disk galaxy with a turbulence sub-grid model, multi-phase interstellar medium and detailed star formation modelling. We investigate the influence of inclination on the observed Lyman-Alpha properties for different snapshots. The Lyman-alpha spectrum, equivalent width distribution and escape fractions vary significantly with the detailed morphology of the disk, leading to variations from one snapshot to another. In particular, we find that supernova-driven cavities near star-forming regions in the simulation can dominate the transmitted Lyman-alpha fraction, suggesting a variability of LAEs on the timescales of the star formation activity.

We investigate the radiative transfer of Lyman alpha photons through simplified anisotropic gas distributions, which represent physically motivated extensions of the popular 'shell-models'. Our study is motivated by the notion that (i) shell models do not always reproduce observed Lyman alpha spectral line profiles, (ii) (typical) shell models do not allow for the escape of ionizing photons, and (iii) the observation & expectation that winds are more complex, anisotropic phenomena. We examine the influence of inclination on the Lyman alpha spectra, relative fluxes and escape fractions. We find the flux to be enhanced/suppressed by factors up to a few depending on the parameter range of the models, corresponding to a boost in equivalent width of the same amplitude if we neglect dust. In general, lower mean optical depths tend to reduce the impact of anisotropies as is expected. We find a correlation between an observed peak in the -- occasionally triple-peaked -- spectrum at the systemic velocity and the existence of a low optical depth cavity along the line of sight. This can be of importance in the search for ionizing photons leaking from high- redshift galaxies since these photons will also be able to escape through the cavity.

We investigate the correlations between the observed fraction of Lyman alpha emission from star-forming galaxies and the large scale structure by post-processing snapshots of a large, high resolution hydrodynamical simulation with a Lyman alpha radiative transfer code at redshifts z=4, 3, 2. We find correlations of the observed fraction with density, density gradient along the line of sight, velocity and velocity gradient along the line of sight, all within the same order of magnitude (tens of percent). Additionally, a correlation with the angular momentum of the dark matter halo is detected. In contrast to a previous study, we find no significant deformation of the 2-point correlation function due to selection effects from radiative transfer in the IGM within the limited statistics of the simulation volume.