We used the mark weighted correlation functions (MCFs), $W(s)$, to study the large scale structure of the Universe. We studied five types of MCFs with the weighting scheme $\rho^\alpha$, where $\rho$ is the local density, and $\alpha$ is taken as $-1,\ -0.5,\ 0,\ 0.5$, and 1. We found that different MCFs have very different amplitudes and scale-dependence. Some of the MCFs exhibit distinctive peaks and valleys that do not exist in the standard correlation functions. Their locations are robust against the redshifts and the background geometry, however it is unlikely that they can be used as ``standard rulers'' to probe the cosmic expansion history. Nonetheless we find that these features may be used to probe parameters related with the structure formation history, such as the values of $\sigma_8$ and the galaxy bias. Finally, after conducting a comprehensive analysis using the full shapes of the $W(s)$s and $W_{\Delta s}(\mu)$s, we found that, combining different types of MCFs can significantly improve the cosmological parameter constraints. Compared with using only the standard correlation function, the combinations of MCFs with $\alpha=0,\ 0.5,\ 1$ and $\alpha=0,\ -1,\ -0.5,\ 0.5,\ 1$ can improve the constraints on $\Omega_m$ and $w$ by $\approx30\%$ and $50\%$, respectively. We find highly significant evidence that MCFs can improve cosmological parameter constraints.

We introduce MANTRA, an annotated dataset of 4869 transient and 71207 non-transient object lightcurves built from the Catalina Real Time Transient Survey. We provide public access to this dataset as a plain text file to facilitate standardized quantitative comparison of astronomical transient event recognition algorithms. Some of the classes included in the dataset are: supernovae, cataclysmic variables, active galactic nuclei, high proper motion stars, blazars and flares. As an example of the tasks that can be performed on the dataset we experiment with multiple data pre-processing methods, feature selection techniques and popular machine learning algorithms (Support Vector Machines, Random Forests and Neural Networks). We assess quantitative performance in two classification tasks: binary (transient/non-transient) and eight-class classification. The best performing algorithm in both tasks is the Random Forest Classifier. It achieves an F1-score of 96.25% in the binary classification and 52.79% in the eight-class classification. For the eight-class classification, non-transients ( 96.83% ) is the class with the highest F1-score, while the lowest corresponds to high-proper-motion stars ( 16.79% ); for supernovae it achieves a value of 54.57% , close to the average across classes. The next release of MANTRA includes images and benchmarks with deep learning models.

We propose a light-weight deep convolutional neural network (CNN) to estimate the cosmological parameters from simulated 3-dimensional dark matter distributions with high accuracy. The training set is based on 465 realizations of a cubic box with a side length of $256\ h^{-1}\ \rm Mpc$, sampled with $128^3$ particles interpolated over a cubic grid of $128^3$ voxels. These volumes have cosmological parameters varying within the flat $\Lambda$CDM parameter space of $0.16 \leq \Omega_m \leq 0.46$ and $2.0 \leq 10^9 A_s \leq 2.3$. The neural network takes as an input cubes with $32^3$ voxels and has three convolution layers, three dense layers, together with some batch normalization and pooling layers. In the final predictions from the network we find a $2.5\%$ bias on the primordial amplitude $\sigma_8$ that can not easily be resolved by continued training. We correct this bias to obtain unprecedented accuracy in the cosmological parameter estimation with statistical uncertainties of $\delta \Omega_m$=0.0015 and $\delta \sigma_8$=0.0029, which are several times better than the results of previous CNN works. Compared with a 2-point analysis method using clustering region of 0-130 and 10-130 $h^{-1}$ Mpc, the CNN constraints are several times and an order of magnitude more precise, respectively. Finally, we conduct preliminary checks of the error-tolerance abilities of the neural network, and find that it exhibits robustness against smoothing, masking, random noise, global variation, rotation, reflection, and simulation resolution. Those effects are well understood in typical clustering analysis, but had not been tested before for the CNN approach. Our work shows that CNN can be more promising than people expected in deriving tight cosmological constraints from the cosmic large scale structure.

Supervised classification of temporal sequences of astronomical images into meaningful transient astrophysical phenomena has been considered a hard problem because it requires the intervention of human experts. The classifier uses the expert's knowledge to find heuristic features to process the images, for instance, by performing image subtraction or by extracting sparse information such as flux time series in the form of light curves. We present a successful deep learning approach that learns directly from imaging data. Our method models explicitly the spatio-temporal patterns with Deep Convolutional Neural Networks and Gated Recurrent Units. We train these deep neural networks using 1.3 million real astronomical images from the Catalina Real-Time Transient Survey to classify the sequences into five different types of astronomical transient classes. The TAO-Net (for Transient Astronomical Objects Network) architecture achieves on the five-type classification task an average F1-score of 54.58$\pm$13.32, almost nine points higher than the F1-score of 45.49 $\pm$ 13.75 from the random forest classification on light curves. The achievement TAO-Net opens the possibility to develop new deep-learning architectures for early transient detection. We make available the training dataset and trained models of TAO-Net to allow for future extensions of this work.

We present shape measurements of Milky Way-sized dark matter halos at redshift $z=0$ in a suite of 30 zoom simulations from the Auriga project. We compare the results in full magnetohydrodynamics against dark matter only simulations and find a strong influence of baryons in making dark matter haloes rounder at all radii compared to their dark matter only counterparts. At distances $\lesssim 30$ kpc, rounder dark matter distributions correlate with extended massive stellar discs and low core gas densities. We measure the alignment between the halo and the disc shapes at different radii and find a high degree of alignment at all radii for most of the galaxies. In some cases the alignment significantly changes as a function of radius implying that the halo shape twists; this effect correlates with recently formed bulges and is almost absent in the dark matter only simulations. In a comparison against observational constraints we find that $20\%$ of halos in our sample are consistent with observational results derived from the Pal 5 stream that favours an almost spherical shape. Including baryons is a required element to achieve this level of agreement. In contrast, none of the simulations (neither dark matter only nor with baryons) match the constraints derived from the Sagittarius stream that favour an oblate dark matter halo.

The Dark Energy Spectroscopic Instrument (DESI) Bright Galaxy Survey (BGS) will be a survey of bright, low redshift galaxies, which is planned to cover an area of ~14,000 sq deg in 3 passes. Each pass will cover the survey area with ~2000 pointings, each of area ~8 sq deg. The BGS is currently proposed to consist of a bright high priority sample to an r-band magnitude limit r ~ 19.5, with a fainter low priority sample to r ~ 20. The geometry of the DESI fibre positioners in the focal plane of the telescope affects the completeness of the survey, and has a non-trivial impact on clustering measurements. Using a BGS mock catalogue, we show that completeness due to fibre assignment primarily depends on the surface density of galaxies. Completeness is high (>95%) in low density regions, but very low (<10%) in the centre of massive clusters. We apply the pair inverse probability (PIP) weighting correction to clustering measurements from a BGS mock which has been through the fibre assignment algorithm. This method is only unbiased if it is possible to observe every galaxy pair. To facilitate this, we randomly promote a small fraction of the fainter sample to be high priority, and dither the set of tile positions by a small angle. We show that inverse pair weighting combined with angular upweighting provides an unbiased correction to galaxy clustering measurements for the complete 3 pass survey, and also after 1 pass, which is highly incomplete.

Outflows and rotation are two ubiquitous kinematic features in the gas kinematics of galaxies. Here we introduce a semi-analytic model to quantify how rotating outflows impact the morphology of the Lyman-$\alpha$ emission line. The model is contrasted against Monte Carlo radiative transfer simulations of outflowing gas with additional solid body rotation. We explore a range of neutral Hydrogen optical depth of $10^5 \leq \tau_{\mathrm{H}} \leq10^7$, rotational velocity $0 \leq v_{\mathrm{rot}}/\mathrm{km\ s}^{-1} \leq 100$ and outflow velocity $0\leq v_{\mathrm{out}}/\mathrm{km\ s}^{-1} \leq 50$. We find three consequences of rotation. First, it introduces a dependency with viewing angle; second it broadens the line and third it increases the flux at the line's center. For all viewing angles, the semi-analytic model reproduces the radiative transfer results for the line width and flux change at the line's center within a $7\%$ and $50\%$ precision for an optical depth of $\tau_{\mathrm{H}}=10^5$, respectively, and within $2\%$ and $1\%$ for an optical depth of $\tau_{\mathrm{H}}=10^7$. Using this model we also show that the peaks of integrated spectra taken from opposite sides of an edge-on rotating gas distribution should have a separation of $\frac{1}{2}v_{\mathrm{rot}}$. The semi-analytic model presented here is a convenient tool to introduce rotational kinematics as a post-processing step of idealized Monte Carlo simulations; it provides a framework to interpret \lya spectra in systems where rotation is expected or directly measured through kinematic maps.

The Emission Line Galaxy survey made by the Dark Energy Spectroscopic Instrument (DESI) survey will be created from five passes of the instrument on the sky. On each pass, the constrained mobility of the ends of the fibres in the DESI focal plane means that the angular-distribution of targets that can be observed is limited. Thus, the clustering of samples constructed using a limited number of passes will be strongly affected by missing targets. In two recent papers, we showed how the effect of missing galaxies can be corrected when calculating the correlation function using a weighting scheme for pairs. Using mock galaxy catalogues we now show that this method provides an unbiased estimator of the true correlation function for the DESI survey after any number of passes. We use multiple mocks to determine the expected errors given one to four passes, compared to an idealised survey observing an equivalent number of randomly selected targets. On BAO scales, we find that the error is a factor 2 worse after one pass, but that after three or more passes, the errors are very similar. Thus we find that the fibre assignment strategy enforced by the design of DESI will not affect the cosmological measurements to be made by the survey, and can be removed as a potential risk for this experiment.

The $\beta$-skeleton is a mathematical method to construct graphs from a set of points that has been widely applied in the areas of image analysis, machine learning, visual perception, and pattern recognition. In this work, we apply the $\beta$-skeleton to study the cosmic web. We use this tool on observed and simulated data to identify the filamentary structures and characterize the statistical properties of the skeleton. In particular, we compare the $\beta$-skeletons built from SDSS-III galaxies to those obtained from MD-PATCHY mocks, and also to mocks directly built from the Big MultiDark $N$-body simulation. We find that the $\beta$-skeleton is able to reveal the underlying structures in observed and simulated samples without any parameter fine-tuning. A different degree of sparseness can be obtained by adjusting the value of $\beta$; in addition, the statistical properties of the length and direction of the skeleton connections show a clear dependence on redshift space distortions (RSDs), cosmological effects and galaxy bias. We also find that the $N$-body simulation accurately reproduces the RSD effect in the data, while the MD-PATCHY mocks appear to underestimate its magnitude. Our proof-of-concept study shows that the statistical properties of the $\beta$-skeleton can be used to probe cosmological parameters and galaxy evolution.

We use simulations to build an analytic probability distribution for the asphericity in the satellite distribution around Local Group (LG) type galaxies in the Lambda Cold Dark Matter (LCDM) paradigm. We use this distribution to estimate the atypicality of the satellite distributions in the LG even when the underlying simulations do not have enough systems fully resembling the LG in terms of its typical masses, separation and kinematics. We demonstrate the method using three different simulations (Illustris-1, Illustris-1-Dark and ELVIS) and a number of satellites ranging from 11 to 15. Detailed results differ greatly among the simulations suggesting a strong influence of the typical DM halo mass, the number of satellites and the simulated baryonic effects. However, there are three common trends. First, at most $2\%$ of the pairs are expected to have satellite distributions with the same asphericity as the LG; second, at most $80\%$ of the pairs have a halo with a satellite distribution as aspherical as in M31; and third, at most $4\%$ of the pairs have a halo with satellite distribution as planar as in the MW. These quantitative results place the LG at the level of a $3\sigma$ outlier in the LCDM paradigm. We suggest that understanding the reasons for this atypicality requires quantifying the asphericity probability distribution as a function of halo mass and large scale environment. The approach presented here can facilitate that kind of study and other comparisons between different numerical setups and choices to study satellites around LG pairs in simulations.

Star-forming Compact Dwarf Galaxies (CDGs) resemble the expected pristine conditions of the first galaxies in the Universe and are the best systems to test models on primordial galaxy formation and evolution. Here we report on one of such CDGs, Tololo 1214-277, which presents a broad, single peaked, highly symmetric Ly$\alpha$ emission line that had evaded theoretical interpretation so far. In this paper we reproduce for the first time these line features with two different physically motivated kinematic models: an interstellar medium composed by outflowing clumps with random motions and an homogeneous gaseous sphere undergoing solid body rotation. The multiphase model requires a clump velocity dispersion of $54.3\pm 0.6$ km s$^{-1}$ with outflows of $54.3\pm 5.1$ km s$^{-1}$, while the bulk rotation velocity is constrained to be $348^{+75}_{-48}$ km s$^{-1}$. We argue that the results from the multiphase model provide a correct interpretation of the data. In that case the clump velocity dispersion implies a dynamical mass of $2\times 10^{9}$ M$_{\odot}$, ten times its baryonic mass. If future kinematic maps of Tololo 1214-277 confirm the velocities suggested by the multiphase model, it would provide additional support to expect such kinematic state in primordial galaxies, opening the opportunity to use the models and methods presented in this paper to constrain the physics of star formation and feedback in the early generation of Ly-$\alpha$ emitting galaxies.

The cosmic web is one of the most striking features of the distribution of galaxies and dark matter on the largest scales in the Universe. It is composed of dense regions packed full of galaxies, long filamentary bridges, flattened sheets and vast low density voids. The study of the cosmic web has focused primarily on the identification of such features, and on understanding the environmental effects on galaxy formation and halo assembly. As such, a variety of different methods have been devised to classify the cosmic web -- depending on the data at hand, be it numerical simulations, large sky surveys or other. In this paper we bring twelve of these methods together and apply them to the same data set in order to understand how they compare. In general these cosmic web classifiers have been designed with different cosmological goals in mind, and to study different questions. Therefore one would not {\it a priori} expect agreement between different techniques however, many of these methods do converge on the identification of specific features. In this paper we study the agreements and disparities of the different methods. For example, each method finds that knots inhabit higher density regions than filaments, etc. and that voids have the lowest densities. For a given web environment, we find substantial overlap in the density range assigned by each web classification scheme. We also compare classifications on a halo-by-halo basis; for example, we find that 9 of 12 methods classify around a third of group-mass haloes (i.e. $M_{\rm halo}\sim10^{13.5}h^{-1}M_{\odot}$) as being in filaments. Lastly, so that any future cosmic web classification scheme can be compared to the 12 methods used here, we have made all the data used in this paper public.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

We use bootstrapping to estimate the bias of concentration estimates on N-body dark matter halos as a function of particle number. We find that algorithms based on the maximum radial velocity and radial particle binning tend to overestimate the concentration by 15%-20% for halos sampled with 200 particles and by 7% - 10% for halos sampled with 500 particles. To control this bias at low particle numbers we propose a new algorithm that estimates halo concentrations based on the integrated mass profile. The method uses the full particle information without any binning, making it reliable in cases when low numerical resolution becomes a limitation for other methods. This method reduces the bias to less than 3% for halos sampled with 200-500 particles. The velocity and density methods have to use halos with at least 4000 particles in order to keep the biases down to the same low level. We also show that the mass-concentration relationship could be shallower than expected once the biases of the different concentration measurements are taken into account. These results show that bootstrapping and the concentration estimates based on the integrated mass profile are valuable tools to probe the internal structure of dark matter halos in numerical simulations.

We revisit calculations of nebular hydrogen Lya and HeII 1640A line strengths for population III galaxies, undergoing continuous and bursts of star formation. We focus on initial mass functions (IMFs) motivated by recent theoretical studies, which generally span a lower range of stellar masses than earlier works. We also account for case-B departures and the stochastic sampling of the IMF. In agreement with previous works, we find that departures from case-B can enhance the Lya flux by a factor of a few, but we argue that this enhancement is driven mainly by collisional excitation and ionization, and not due to photoionization from the n = 2 state of atomic hydrogen. The increased sensitivity of the Lya flux to the high-energy end of the galaxy spectrum makes it more subject to stochastic sampling of the IMF. The latter introduces a dispersion in the predicted nebular line fluxes around the deterministic value by as much as a factor of ~4. In contrast, the stochastic sampling of the IMF has less impact on the emerging Lyman Werner (LW) photon flux. When case-B departures and stochasticity effects are combined, nebular line emission from population III galaxies can be up to one order of magnitude brighter than predicted by 'standard' calculations that do not include these effects. This enhances the prospects for detection with future facilities such as JWST and large, groundbased telescopes.

In this paper we study the impact of cosmic variance and observational uncertainties in constraining the mass and occupation fraction, $f_{\rm occ}$, of dark matter halos hosting Ly-$\alpha$ Emitting Galaxies (LAEs) at high redshift. To this end, we construct mock catalogs from an N-body simulation to match the typical size of observed fields at $z=3.1$ ($\sim 1 {\rm deg^2}$). In our model a dark matter halo with mass in the range $M_{\rm min}<M_{\mathrm h}<M_{\rm max}$ can only host one detectable LAE at most. We proceed to explore the parameter space determined by $M_{\rm min}$,$M_{\rm max}$ and $f_{\rm occ}$ with a Markov Chain Monte-Carlo algorithm using the angular correlation function (ACF) and the LAEs number density as observational constraints. We find that the preferred minimum and maximum masses in our model span a wide range $10^{10.0}h^{-1}{\rm{M_{\odot}}}\leq M_{\rm min} \leq 10^{11.1}h^{-1}{\rm{M_{\odot}}}$ , $10^{11.0}h^{-1}{\rm{M_{\odot}}}\leq M_{\rm max} \leq 10^{13.0}h^{-1}{\rm{M_{\odot}}}$; followed by a wide range in the occupation fraction $0.02\leq f_{\rm occ} \leq 0.30$. As a consequence the median mass, $M_{50}$, of all the consistent models has a large uncertainty $M_{50} = 3.16^{+9.34}_{-2.37}\times 10^{10}$$h^{-1}{\rm{M_{\odot}}}$. However, we find that the same individual models have a relatively tight $1\sigma$ scatter around the median mass $\Delta M_{1\sigma} = 0.55^{+0.11}_{-0.31}$ dex. We are also able to show that \focc\ is uniquely determined by $M_{\rm min}$, regardless of $M_{\rm max}$. We argue that upcoming large surveys covering at least $25$ deg$^{2}$ should be able to put tighter constraints on $M_{\rm min}$ and $f_{\rm occ}$ through the LAE number density distribution width constructed over several fields of $\sim 1$ deg$^{2}$.

The nature of dark matter, dark energy and large-scale gravity pose some of the most pressing questions in cosmology today. These fundamental questions require highly precise measurements, and a number of wide-field spectroscopic survey instruments are being designed to meet this requirement. A key component in these experiments is the development of a simulation tool to forecast science performance, define requirement flow-downs, optimize implementation, demonstrate feasibility, and prepare for exploitation. We present SPOKES (SPectrOscopic KEn Simulation), an end-to-end simulation facility for spectroscopic cosmological surveys designed to address this challenge. SPOKES is based on an integrated infrastructure, modular function organization, coherent data handling and fast data access. These key features allow reproducibility of pipeline runs, enable ease of use and provide flexibility to update functions within the pipeline. The cyclic nature of the pipeline offers the possibility to make the science output an efficient measure for design optimization and feasibility testing. We present the architecture, first science, and computational performance results of the simulation pipeline. The framework is general, but for the benchmark tests, we use the Dark Energy Spectrometer (DESpec), one of the early concepts for the upcoming project, the Dark Energy Spectroscopic Instrument (DESI). We discuss how the SPOKES framework enables a rigorous process to optimize and exploit spectroscopic survey experiments in order to derive high-precision cosmological measurements optimally.

We propose a method based on the redshift dependence of the Alcock-Paczynski (AP) test to measure the expansion history of the Universe. It uses the isotropy of the galaxy density gradient field to constrain cosmological parameters. If the density parameter $\Omega_m$ or the dark energy equation of state $w$ are incorrectly chosen, the gradient field appears to be anisotropic with the degree of anisotropy varying with redshift. We use this effect to constrain the cosmological parameters governing the expansion history of the Universe. Although redshift-space distortions (RSD) induced by galaxy peculiar velocities also produce anisotropies in the gradient field, these effects are close to uniform in magnitude over a large range of redshift. This makes the redshift variation of the gradient field anisotropy relatively insensitive to the RSD. By testing the method on mock surveys drawn from the Horizon Run 3 cosmological N-body simulations, we demonstrate that the cosmological parameters can be estimated without bias. Our method is complementary to the baryon acoustic oscillation or topology methods as it depends on $D_AH$, the product of the angular diameter distance and the Hubble parameter.

We present a new method to find voids in cosmological simulations based on the tidal and the velocity shear tensors definitions of the cosmic web. We use the fractional anisotropy (FA) computed from the eigenvalues of each web scheme as a void tracer. We identify voids using a watershed transform based on the local minima of the FA field without making any assumption on the shape or structure of the voids. We test the method on the Bolshoi simulation and report on the abundance and radial averaged profiles for the density, velocity and fractional anisotropy. We find that voids in the velocity shear web are smaller than voids in the tidal web, with a particular overabundance of very small voids in the inner region of filaments/sheets. We classify voids as subcompensated/overcompansated depending on the absence/presence of an overdense matter ridge in their density profile, finding that close to $65\%$ and $35\%$ of the total population are classified into each category, respectively. Finally, we find evidence for the existence of universal profiles from the radially averaged profiles for density, velocity and fractional anisotropy. This requires that the radial coordinate is normalized to the effective radius of each void. Put together, all these results show that the FA is a reliable tracer for voids, which can be used in complementarity to other existing methods and tracers.

We use the Bolshoi Simulation to find the most probable location of the Local Group (LG) in the cosmic web. Our LG simulacra are pairs of halos with isolation and kinematic properties consistent with observations. The cosmic web is defined using a tidal tensor approach. We find that the LG's preferred location is regions with a dark matter overdensity close to the cosmic average. This makes filaments and sheets the preferred environment. We also find a strong alignment between the LG and the cosmic web. The orbital angular momentum is preferentially perpendicular to the smallest tidal eigenvector, while the vector connecting the two halos is strongly aligned along the smallest tidal eigenvector and perpendicular to the largest tidal eigenvector; the pair lies and moves along filaments and sheets. We do not find any evidence for an alignment between the spin of each halo in the pair and the cosmic web.

We explore the characteristics of the cosmic web around Local Group(LG) like pairs using a cosmological simulation in the $\Lambda$CDM cosmology. We use the Hessian of the gravitational potential to classify regions on scales of $\sim 2$ Mpc as a peak, sheet, filament or void. The sample of LG counterparts is represented by two samples of halo pairs. The first is a general sample composed by pairs with similar masses and isolation criteria as observed for the LG. The second is a subset with additional observed kinematic constraints such as relative pair velocity and separation. We find that the pairs in the LG sample with all constraints are: (i) Preferentially located in filaments and sheets, (ii) Located in in a narrow range of local overdensity $0<\delta<2$, web ellipticity $0.1<e<1.0$ and prolateness $-0.4<p<0.4$. (iii) Strongly aligned with the cosmic web. The alignments are such that the pair orbital angular momentum tends to be perpendicular to the smallest tidal eigenvector, $\hat{e}_3$, which lies along the filament direction or the sheet plane. A stronger alignment is present for the vector linking the two halos with the vector $\hat{e}_3$. Additionally, we fail to find a strong correlation of the spin of each halo in the pair with the cosmic web. All these trends are expected to a great extent from the selection on the LG total mass on the general sample. Applied to the observed LG, there is a potential conflict between the alignments of the different planes of satellites and the numerical evidence for satellite accretion along filaments; the direction defined by $\hat{e}_3$. This highlights the relevance of achieving a precise characterization of the place of the LG in the cosmic web in the cosmological context provided by $\Lambda$CDM.

We analyzed the Einstein radius, $\theta_E$, in our sample of SL2S galaxy groups, and compared it with $R_A$ (the distance from the arcs to the center of the lens), using three different approaches: 1.- the velocity dispersion obtained from weak lensing assuming a Singular Isothermal Sphere profile ($\theta_{E,I}$), 2.- a strong lensing analytical method ($\theta_{E,II}$) combined with a velocity dispersion-concentration relation derived from numerical simulations designed to mimic our group sample, 3.- strong lensing modeling ($\theta_{E,III}$) of eleven groups (with four new models presented in this work) using HST and CFHT images. Finally, $R_A$ was analyzed as a function of redshift $z$ to investigate possible correlations with L, N, and the richness-to-luminosity ratio (N/L). We found a correlation between $\theta_{E}$ and $R_A$, but with large scatter. We estimate $\theta_{E,I}$ = (2.2 $\pm$ 0.9) + (0.7 $\pm$ 0.2)$R_A$, $\theta_{E,II}$ = (0.4 $\pm$ 1.5) + (1.1 $\pm$ 0.4)$R_A$, and $\theta_{E,III}$ = (0.4 $\pm$ 1.5) + (0.9 $\pm$ 0.3)$R_A$ for each method respectively. We found a weak evidence of anti-correlation between $R_A$ and $z$, with Log$R_A$ = (0.58$\pm$0.06) - (0.04$\pm$0.1)$z$, suggesting a possible evolution of the Einstein radius with $z$, as reported previously by other authors. Our results also show that $R_A$ is correlated with L and N (more luminous and richer groups have greater $R_A$), and a possible correlation between $R_A$ and the N/L ratio. Our analysis indicates that $R_A$ is correlated with $\theta_E$ in our sample, making $R_A$ useful to characterize properties like L and N (and possible N/L) in galaxy groups. Additionally, we present evidence suggesting that the Einstein radius evolves with $z$.

We present results of radiative transfer calculations to measure the impact of gas bulk rotation on the morphology of the Lyman $\alpha$ emission line in distant galaxies. We model a galaxy as a sphere with an homogeneous mixture of dust and hydrogen at a constant temperature. These spheres undergo solid-body rotation with maximum velocities in the range 0-300 \kms and neutral hydrogen optical depths in the range $\tau_{\rm H}=10^{5}-10^{7}$. We consider two types of source distributions in the sphere: central and homogeneous. Our main result is that rotation introduces a dependence of the line morphology with viewing angle and rotational velocity. Observations with a line of sight parallel to the rotation axis yield line morphologies similar to the static case. For lines of sight perpendicular to the rotation axis both the intensity at the line center and the line width increase with rotational velocity. Along the same line of sight, the line becomes single peaked at rotational velocities close to half the line width in the static case. Notably, we find that rotation does not induce any spatial anisotropy in the integrated line flux, the escape fraction or the average number of scatterings. This is because Lyman {\alpha} scattering through a rotating solid-body proceeds identical as in the static case. The only difference is the doppler shift from the different regions in the sphere that move with respect to the observer. This allows us to derive an analytic approximation for the viewing-angle dependence of the emerging spectrum, as a function of rotational velocity.

Although commonly assumed to be spherical, dark matter halos are predicted to be non-spherical by N-body simulations and their asphericity has a potential impact on the systematic uncertainties in dark matter searches. The evaluation of these uncertainties is the main aim of this work, where we study the impact of aspherical dark matter density distributions in Milky-Way-like halos on direct and indirect searches. Using data from the large N-body cosmological simulation Bolshoi, we perform a statistical analysis and quantify the systematic uncertainties on the determination of local dark matter density and the so-called $J$ factors for dark matter annihilations and decays from the galactic center. We find that, due to our ignorance about the extent of the non-sphericity of the Milky Way dark matter halo, systematic uncertainties can be as large as 35%, within the 95% most probable region, for a spherically averaged value for the local density of 0.3-0.4 GeV/cm$^3$. Similarly, systematic uncertainties on the $J$ factors evaluated around the galactic center can be as large as 10% and 15%, within the 95% most probable region, for dark matter annihilations and decays, respectively.

We study the alignment of dark matter haloes with the cosmic web characterized by the tidal and velocity shear fields. We focus on the alignment of their shape, angular momentum and peculiar velocities. We use a cosmological N-body simulation that allows to study dark matter halos spanning almost five orders of magnitude in mass ($10^{9}$-$10^{14}$) $h^{-1}$$M_{\odot}$ and spatial scales of $(0.5$-$1.0)$ $h^{-1}$ Mpc to define the cosmic web. We find that the halo shape presents the strongest alignment along the smallest tidal eigenvector, e.g. along filaments and walls, with a signal that gets stronger as the halo mass increases. In the case of the velocity shear field only massive halos $>10^{12}$ $h^{-1}$$M_{\odot}$ tend to have their shapes aligned along the largest tidal eigenvector; that is, perpendicular to filaments and walls. For the angular momentum we find alignment signals only for halos more massive than $10^{12}$ $h^{-1}$$M_{\odot}$ both in the tidal and velocity shear webs where the preferences are for it to be parallel to the middle eigenvector; perpendicular to filaments and parallel to walls. Finally, the peculiar velocities show a strong alignment along the smallest tidal eigenvector for all halo masses; halos move along filaments and walls. In the velocity shear the same alignment is present but weaker and only for haloes less massive than $10^{12}$ $h^{-1}$$M_{\odot}$. Our results clearly show that the two different algorithms we used to define the cosmic web describe different physical aspects of non-linear collapse and should be used in a complementary way to understand the effect of the cosmic web on galaxy evolution.

We present Simulated Annealing fiber-to-target allocation simulations for the proposed DESI and 4MOST massively multiplexed spectroscopic surveys, and for both Poisson and realistically clustered mock target samples. We simulate both Echidna and theta-phi actuator designs, including the restrictions caused by the physical actuator characteristics during repositioning. For DESI, with theta-phi actuators, used in 5 passes over the sky for a mock ELG/LRG/QSO sample, with matched fiber and target densities, a total target allocation yield of 89.3% was achieved, but only 83.7% for the high-priority Ly-alpha QSOs. If Echidna actuators are used with the same pitch and number of passes, the yield increases by 5.7% and 16% respectively. Echidna also allows a factor-of-two increase in the number of close Ly-alpha QSO pairs that can be observed. Echidna spine tilt causes a variable loss of throughput, with average loss being the same as the loss at the rms tilt. With a natural tilt minimization scheme, we find an rms tilt always close to 0.58 x maximum. There is an additional but much smaller defocus loss, equivalent to an average defocus of 30microns. These tilt losses offset the gains in yield for Echidna, but because the survey strategy is driven by the higher priority targets, a clear survey speed advantage remains. For 4MOST, high and low latitude sample mock catalogs were supplied by the 4MOST team, and allocations were carried out with the proposed Echidna-based positioner geometry. At high latitudes, the resulting target completeness was 85.3% for LR targets and 78.9% for HR targets. At low latitude, the target completeness was 93.9% for LR targets and 71.2% for HR targets.

We estimate the expected distribution of displacements between the two dominant dark matter (DM) peaks (DM-DM displacements) and between DM and gaseous baryon peak (DM-gas displacements) in dark matter halos with masses larger than $10^{13}$ Msun/h. We use as a benchmark the observation of SL2S J08544-0121, which is the lowest mass system ($1.0\times 10^{14}$ Msun/h) observed so far featuring a bi-modal dark matter distribution with a dislocated gas component. We find that $(50 \pm 10)$% of the dark matter halos with circular velocities in the range 300 km/s to 700 km/s (groups) show DM-DM displacements equal or larger than $186 \pm 30$ kpc/h as observed in SL2S J08544-0121. For dark matter halos with circular velocities larger than 700 km/s (clusters) this fraction rises to 70 $\pm$ 10%. Using the same simulation we estimate the DM-gas displacements and find that 0.1 to 1.0% of the groups should present separations equal or larger than $87\pm 14$kpc/h corresponding to our observational benchmark; for clusters this fraction rises to (7 $\pm$ 3)%, consistent with previous studies of dark matter to baryon separations. Considering both constraints on the DM-DM and DM-gas displacements we find that the number density of groups similar to SL2S J08544-0121 is $\sim 6.0\times 10^{-7}$ Mpc$^{-3}$, three times larger than the estimated value for clusters. These results open up the possibility for a new statistical test of LCDM by looking for DM-gas displacements in low mass clusters and groups.

Recent observations constrained the tangential velocity of M31 with respect to the Milky Way (MW) to be v_tan<34.4 km/s and the radial velocity to be in the range v_rad=-109+/- 4.4 km/s (van der Marel et al. 2012). In this study we use a large volume high resolution N-body cosmological simulation (Bolshoi) together with three constrained simulations to statistically study this kinematics in the context of the LCDM. The comparison of the ensembles of simulated pairs with the observed LG at the 1-sigma level in the uncertainties has been done with respect to the radial and tangential velocities, the reduced orbital energy (e_tot), angular momentum (l_orb) and the dimensionless spin parameter, lambda. Our main results are: (i) the preferred radial and tangential velocities for pairs in LCDM are v_rad=-80+/-20 km/s, v_tan=50+/-10 km/s, (ii) pairs around that region are 3 to 13 times more common than pairs within the observational values, (iii) 15%to 24% of LG-like pairs in LCDM have energy and angular momentum consistent with observations while (iv) 9% to 13% of pairs in the same sample show similar values in the inferred dimensionless spin parameter. It follows that within current observational uncertainties the quasi-conserved quantities that characterize the orbit of the LG, i.e. e_tot, r_orb and lambda, do not challenge the standard LCDM model, but the model is in tension with regard to the actual values of the radial and tangential velocities. This might hint to a problem of the LCDM model to reproduce the observed LG.

The alignment of DM halos and the surrounding large scale structure (LSS) is examined in the context of the cosmic web. Halo spin, shape and the orbital angular momentum of subhaloes is investigated relative to the LSS using the eigenvectors of the velocity shear tensor evaluated on a grid with a scale of 1 Mpc/h, deep within the non-linear regime. Knots, filaments, sheets and voids are associated with regions that are collapsing along 3, 2, 1 or 0 principal directions simultaneously. Each halo is tagged with a web classification (i.e. knot halo, filament halo, etc) according to the nature of the collapse at the halo's position. The full distribution of shear eigenvalues is found to be substantially different from that tagged to haloes, indicating that the observed velocity shear is significantly biased. We find that larger mass haloes live in regions where the shear is more isotropic, namely the expansion or collapse is more spherical. A correlation is found between the halo's shape and the eigenvectors of the shear tensor, with the longest (shortest) axis of the halo's shape being aligned with the slowest (fastest) collapsing eigenvector. This correlation is web independent, suggesting that the velocity shear is a fundamental tracer of the halo alignment. A similar result is found for the alignment of halo spin with the cosmic web. It has been shown that high mass haloes exhibit a spin flip with respect to the LSS: we find the mass at which this spin flip occurs is web dependent and not universal as suggested previously. Although weaker than haloes, subhalo orbits too exhibit an alignment with the LSS, providing a possible insight into the highly correlated co-rotation of the Milky Way's satellite system. The present study suggests that the velocity shear tensor constitutes the natural framework for studying the directional properties of the non-linear LSS and of halos and galaxies.

Observations of distant galaxies play a key role in improving our understanding of the Epoch of Reionization (EoR). The observed Ly-alpha emission line strength - quantified by its restframe equivalent width (EW) - provides a valuable diagnostic of stellar populations and dust in galaxies during and after the EoR. In this paper we quantify the effects of star formation stochasticity on the predicted Ly-alpha EW in dwarf galaxies, using the publicly available code SLUG ('Stochastically Light Up Galaxies'). We compute the number of hydrogen ionizing photons, as well as flux in the Far UV for a set of models with star formation rates (SFR) in the range 10-3-1 Msol/yr. From these fluxes we compute the luminosity, L-alpha, and the EW of the Ly-alpha line. We find that stochasticity alone induces a broad distribution in L-alpha and EW at a fixed SFR, and that the widths of these distributions decrease with increasing SFR. We parameterize the EW probability density function (PDF) as an SFR-dependent double power law. We find that it is possible to have EW as low as ~EW0/4 and as high as ~3 times the EW0, where EW0 denotes the expected EW in the absence of stochasticity. We argue that stochasticity may therefore be important when linking drop-out and narrow-band selected galaxies, when identifying population III galaxies, and that it may help to explain the large EW (EW > 100 - 200 A) observed for a fraction of Ly-alpha- selected galaxies. Finally, we show that stochasticity can also affect the inferred escape fraction of ionizing photons from galaxies. In particular, we argue that stochasticity may simultaneously explain the observed anomalous ratios of the Lyman continuum flux density to the (non-ionizing) UV continuum density in so-called Lyman-Bump galaxies at z = 3.1, as well as the absence of such objects among a sample of z = 1.3 drop-out galaxies.

There are thousands of confirmed detections of star forming galaxies at high redshift (z > 4). These observations rely primarily on the detection of the spectral Lyman Break and the Lyman-alpha emission line. Theoretical modelling of these sources helps to interpret the observations in the framework of the standard cosmological paradigm. We present results from the High-z MareNostrum Project, aimed at constructing a panchromatic picture of the high redshift galaxy evolution that will improve our understanding of young star forming galaxies. Our simulation successfully reproduces the observational constraints from Lyman Break Galaxies and Lyman-alpha emitters at 5 < z < 7 . Based on this model we make predictions on the expected Far Infrared (FIR) emission that should be observed for LAEs. These predictions will help to settle down the question on the dust content of massive high-z galaxies, an issue that will be feasible to probe observationally with the Atacama Large Millimetre Array (ALMA).

A new approach for the classification of the cosmic web is presented. In extension of the previous work of Hahn et al. (2007) and Forero-Romero et al. (2009) the new algorithm is based on the analysis of the velocity shear tensor rather than the gravitational tidal tensor. The procedure consists of the construction of the the shear tensor at each (grid) point in space and the evaluation of its three eigenvectors. A given point is classified to be either a void, sheet, filament or a knot according to the number of eigenvalues above a certain threshold, 0, 1, 2, or 3 respectively. The threshold is treated as a free parameter that defines the web. The algorithm has been applied to a dark matter only, high resolution simulation of a box of side-length 64$h^{-1}$Mpc and N = $1024^3$ particles with the framework of the WMAP5/LCDM model. The resulting velocity based cosmic web resolves structures down to <0.1$h^{-1}$Mpc scales, as opposed to the ~1$h^{-1}$Mpc scale of the tidal based web. The under-dense regions are made of extended voids bisected by planar sheets, whose density is also below the mean. The over-dense regions are vastly dominated by the linear filaments and knots. The resolution achieved by the velocity based cosmic web provides a platform for studying the formation of halos and galaxies within the framework of the cosmic web.

We propose a connection between the visualisation of cosmic matter and structure formation in the Cartesian tradition and that used by contemporary astrophysics. More precisely, we identify cosmological simulations of large scale structure in the Universe with the system of vortices in Descartes physics. This connection operates at different levels of the images: their representational purpose; the theoretical systems behind their use; and, finally, their function and materiality as visual productions. A skilled use of image analysis is necessary to stress the continuities and peculiarities between different epochs and disciplines.

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 present theoretical results for the expected fraction of Lyman Break Galaxies (LBGs) to be detected as strong Lyman alpha emitters (LAEs) in the redshift range 5 < z < 7. We base our analysis on the 2-billion particle SPH simulation MareNostrum High-z Universe. We approximate galaxies as static dusty slabs with an additional clumpy dust distribution affecting stellar populations younger than 25 Myr. The model for the Lyman alpha escape fraction is based on the results of our Monte-Carlo radiative transfer code (CLARA) for a slab configuration. We also fix the transmission of Lyman alpha photons through the intergalactic medium to a constant value of 50% at all redshifts. From the results of this model we calculate xLya, the fraction of Lyman Break Galaxies with Lyman alpha equivalent width (EW) larger than 50 Angstrom. We find a remarkable agreement with observational data at 4.5 < z < 6. For bright (-22 < MUV < -20.5) and faint (-20.5 < MUV < -18.5) galaxies our model predicts xLya = 0.02 \pm 0.01 and xLya = 0.47 \pm 0.01 while observers report xLya = 0.08 \pm 0.02 and xLya = 0.47 \pm 0.16, respectively. Additional evolution of the extinction model at redshift z \sim 7, that decreases the intensity of transmitted Lyman alpha radiation by a factor of f_T = 0.4 as to match the LAE luminosity function at z \sim 6.5, naturally provides a good match for the recently reported xLya fractions at z > 6.3. Exploring different toy models for the Lyman alpha escape fraction, we show that a decreasing Lyman alpha escape fraction with increasing UV galaxy luminosity is a key element in our model to explain the of larger xLya fractions for fainter LBGs.

We present the implementation of a halo based method for the reconstruction of the cosmic mass density field. The method employs the mass density distribution of dark matter haloes and its environments computed from cosmological N-body simulations and convolves it with a halo catalog to reconstruct the dark matter density field determined by the distribution of haloes. We applied the method to the group catalog of Yang etal (2007) built from the SDSS Data Release 4. As result we obtain reconstructions of the cosmic mass density field that are independent on any explicit assumption of bias. We describe in detail the implementation of the method, present a detailed characterization of the reconstructed density field (mean mass density distribution, correlation function and counts in cells) and the results of the classification of large scale environments (filaments, voids, peaks and sheets) in our reconstruction. Applications of the method include morphological studies of the galaxy population on large scales and the realization of constrained simulations.

We make detailed theoretical predictions for the assembly properties of the Local Group (LG) in the standard LambdaCDM cosmological model. We use three cosmological N-body dark matter simulations from the CLUES project, which are designed to reproduce the main dynamical features of the matter distribution down to the scale of a few Mpc around the LG. Additionally, we use the results of an unconstrained simulation with a sixty times larger volume to calibrate the influence of cosmic variance. We characterize the Mass Aggregation History (MAH) for each halo by three characteristic times, the formation, assembly and last major merger times. A major merger is defined by a minimal mass ratio of 10:1. We find that the three LGs share a similar MAH with formation and last major merger epochs placed on average \approx 10 - 12 Gyr ago. Between 12% and 17% of the halos in the mass range 5 x 10^11 Msol/h < M_h < 5 x 10^12 Msol/h have a similar MAH. In a set of pairs of halos within the same mass range, a fraction of 1% to 3% share similar formation properties as both halos in the simulated LG. An unsolved question posed by our results is the dynamical origin of the MAH of the LGs. The isolation criteria commonly used to define LG-like halos in unconstrained simulations do not narrow down the halo population into a set with quiet MAHs, nor does a further constraint to reside in a low density environment. The quiet MAH of the LGs provides a favorable environment for the formation of disk galaxies like the Milky Way and M31. The timing for the beginning of the last major merger in the Milky Way dark matter halo matches with the gas rich merger origin for the thick component in the galactic disk. Our results support the view that the specific large and mid scale environment around the Local Group play a critical role in shaping its MAH and hence its baryonic structure at present.

Using CLARA (Code for Lyman Alpha Radiation Analysis) we constrain the escape fraction of Lyman-Alpha radiation in galaxies in the redshift range 5<z<7, based on the MareNostrum High-z Universe, a SPH cosmological simulation with more than 2 billion particles. We approximate Lyman-Alpha Emitters (LAEs) as dusty gaseous slabs with Lyman-Alpha radiation sources homogeneously mixed in the gas. Escape fractions for such a configuration and for different gas and dust contents are calculated using our newly developed radiative transfer code CLARA. The results are applied to the MareNostrum High-z Universe numerical galaxies. The model shows a weak redshift evolution and good agreement with estimations of the escape fraction as a function of reddening from observations at z \sim 2.2 and z \sim 3. We extend the slab model by including additional dust in a clumpy component in order to reproduce the UV con- tinuum luminosity function and UV colours at redshifts z>~5. The LAE Luminosity Function (LF) based on the extended clumpy model reproduces broadly the bright end of the LF derived from observations at z \sim 5 and z \sim 6. At z \sim 7 our model over-predicts the LF by roughly a factor of four, presumably because the effects of the neutral intergalactic medium are not taken into account. The remaining tension between the observed and simulated faint end of the LF, both in the UV-continuum and Lyman-Alpha at redshifts z \sim 5 and z \sim 6 points towards an overabundance of simulated LAEs hosted in haloes of masses 1.0x10^10h-1Msol < Mh < 4.0x10^10h-1Msol. Given the difficulties in explaining the observed overabundance by dust absorption, a probable origin of the mismatch are the high star formation rates in the simulated haloes around the quoted mass range. A more efficient supernova feedback should be able to regulate the star formation process in the shallow potential wells of these haloes.

We estimate the expected distribution of displacements between the dark matter and gas cores in simulated clusters. We use the MareNostrum Universe, one of the largest non radiative, SPH Lambda CDM cosmological simulations. We find that projected 2-D displacements between dark matter and gas, equal or larger than the observed in the Bullet Cluster, are expected in 1% to 2% of the clusters with masses larger than 10^{14} Msun. The 2-D displacement distribution is roughly the same between redshifts 0<z<0.5 when multiplied by a factor of (1+z)^{-1/2}. We conclude that the separations between dark matter and gas as observed in the bullet cluster can be easily found in a Lambda CDM universe. Furthermore we find that the displacement distribution is not very sensitive to the normalization of the power spectrum. Upcoming surveys could extend the measurements of these displacements between dark matter and gas into large samples of hundreds of clusters, providing a potential test for Lambda CDM.

We discuss the rest-frame UV emission between 5< z < 7 from the MareNostrum High-z Universe, a SPH simulation done with more than 2 billion particles. Cosmological simulations of galaxy formation generally overpredict the UV restframe luminosity function at high redshift, both at the bright and faint ends. In this Letter we explore a dust attenuation model where a larger extinction is applied to star populations younger than a given age, mimicking the effect of a clumpy interstellar medium. We show that this scenario fits reasonably well both the UV luminosity functions and the UV-continuum slopes derived from observations. The model assumes a large obscuration for stars younger than 25 Myr from the gas clouds where they should be embedded at their formation time. We find that the optical depth in these clouds should be between 30 and 100 times larger than the mean optical depth for the homogeneous part of the interstellar medium. These values are one order of magnitude larger than those estimated in local galaxies. Therefore, we conclude that LambdaCDM predictions for the high-z UV emission can accommodate the current observations if we consider a dust extinction model based on the assumption of a clumpy environment at high redshift.

A dynamical classification of the cosmic web is proposed. The large scale environment is classified into four web types: voids, sheets, filaments and knots. The classification is based on the evaluation of the deformation tensor, i.e. the Hessian of the gravitational potential, on a grid. The classification is based on counting the number of eigenvalues above a certain threshold, lambda_th at each grid point, where the case of zero, one, two or three such eigenvalues corresponds to void, sheet, filament or a knot grid point. The collection of neighboring grid points, friends-of-friends, of the same web attribute constitutes voids, sheets, filaments and knots as web objects. A simple dynamical consideration suggests that lambda_th should be approximately unity, upon an appropriate scaling of the deformation tensor. The algorithm has been applied and tested against a suite of (dark matter only) cosmological N-body simulations. In particular, the dependence of the volume and mass filling fractions on lambda_th and on the resolution has been calculated for the four web types. Also, the percolation properties of voids and filaments have been studied. Our main findings are: (a) Already at lambda_th = 0.1 the resulting web classification reproduces the visual impression of the cosmic web. (b) Between 0.2 < lambda_th < 0.4, a system of percolated voids coexists with a net of interconected filaments. This suggests a reasonable choice for lambda_th as the parameter that defines the cosmic web. (c) The dynamical nature of the suggested classification provides a robust framework for incorporating environmental information into galaxy formation models, and in particular the semi-analytical ones.

We introduce the contour process to describe the geometrical properties of merger trees. The contour process produces a one-dimensional object, the contour walk, which is a translation of the merger tree. We portray the contour walk through its length and action. The length is proportional to to the number of progenitors in the tree, and the action can be interpreted as a proxy of the mean length of a branch in a merger tree. We obtain the contour walk for merger trees extracted from the public database of the Millennium Run and also for merger trees constructed with a public Monte-Carlo code which implements a Markovian algorithm. The trees correspond to halos of final masses between 10^{11} h^{-1} M_sol and 10^{14} h^{-1} M_sol. We study how the length and action of the walks evolve with the mass of the final halo. In all the cases, except for the action measured from Markovian trees, we find a transitional scale around 3 \times 10^{12} h^{-1} M_sol. As a general trend the length and action measured from the Markovian trees show a large scatter in comparison with the case of the Millennium Run trees.

We propose a general framework to scrutinize the performance of semi-analytic codes of galaxy formation. The approach is based on the analysis of the outputs from the model after a series of perturbations in the input parameters controlling the baryonic physics. The perturbations are chosen in a way that they do not change the results in the luminosity function or mass function of the galaxy population. We apply this approach on a particular semi-analytic model called GalICS. We chose to perturb the parameters controlling the efficiency of star formation and the efficiency of supernova feedback. We keep track of the baryonic and observable properties of the central galaxies in a sample of dark matter halos with masses ranging from 10^{10} M_sol to 10^{13} M_sol. We find very different responses depending on the halo mass. For small dark matter halos its central galaxy responds in a highly predictable way to small perturbation in the star formation and feedback efficiency. For massive dark matter halos, minor perturbations in the input parameters can induce large fluctuations on the properties of its central galaxy, at least $\sim 0.1$ in (B-V) color or $\sim 0.5$ mag in U or r filter, in a seemingly random fashion. We quantify this behavior through an objective scalar function we call predictability. We argue that finding the origin of this behavior needs additional information from other approximations and different semi-analytic codes. Furthermore, the implementation of an scalar objective function, such as the predictability, opens the door to quantitative benchmarking of semi-analytic codes based on its numerical performance.

We present the Lensed Mock Map Facility (LeMoMaF), a tool designed to perform mock weak lensing measurements on numerically simulated chunks of the universe. Coupling N-body simulations to a semi-analytical model of galaxy formation, LeMoMaF can create realistic lensed images and mock catalogues of galaxies, at wavelengths ranging from the UV to the submm. To demonstrate the power of such a tool we compute predictions of the source-lens clustering effect on the convergence statistics, and quantify the impact of weak lensing on galaxy counts in two different filters. We find that the source-lens clustering effect skews the probability density function of the convergence towards low values, with an intensity which strongly depends on the redshift distribution of galaxies. On the other hand, the degree of enhancement or depletion in galaxy counts due to weak lensing is independent of the source-lens clustering effect. We discuss the impact on the two-points shear statistics to be measured by future missions like SNAP and LSST. The source-lens clustering effect would bias the estimation of sigma_8 from two point statistics by 2% -5%. We conclude that accurate photometric redshifts for individual galaxies are necessary in order to quantify and isolate the source-lens clustering effect.