Galaxy groups provide the means for a great diversity of studies that contribute to a better understanding of the structure of the universe on a large scale and allow the properties of galaxies to be linked to those of the host halos. However, the identification of galaxy systems is a challenging task and therefore it is necessary to improve the techniques involved as much as possible. In view of the large present and forthcoming galaxy catalogues, we propose, implement, and evaluate an algorithm that combines the two most popular techniques used to identify galaxy systems. The algorithm can be easily applied to any spectroscopic galaxy catalogue, but here we demonstrate its use on the Sloan Digital Sky Survey. Assuming that a galaxy group is a gravitationally bounded system that has at least one bright galaxy, we begin by identifying groups with a Friends-Of-Friends algorithm adapted to fit this definition. In view of the shortcomings of this method, particularly the lack of ability to identify low-number groups, and consequently the inability to study the occupation of halos throughout the mass range, we improve it by adding a halo-based procedure. To assess the performance, we construct a mock catalogue from a semi-analytical model to compare the groups identified using our method with those obtained from the simulation. The comparison of groups extracted using our method with those of a mock catalogue shows that the proposed algorithm provides excellent results. The modifications introduced to the Friends-Of-Friends algorithm in the first part of the procedure to fit the adopted group definition gave reliable groups. Furthermore, the incorporation of the halo-based method reduces the interlopers while more accurately reproducing the number of galaxies per group. As a specific application, we use the algorithm to extract groups from the Sloan Digital Sky Survey.

We present a new regular grid search algorithm for quick fixed-radius nearest-neighbor lookup developed in Python. This module indexes a set of k-dimensional points in a regular grid, with optional periodic conditions, providing a fast approach for nearest neighbors queries. In this first installment we provide three types of queries: $bubble$, $shell$ and the $nth-nearest$; as well as three different metrics of interest in astronomy: the $euclidean$ and two distance functions in spherical coordinates of varying precision, $haversine$ and $Vincenty$; and the possibility of providing a custom distance function. This package results particularly useful for large datasets where a brute-force search turns impractical.

We introduce an algorithm aimed to identify large-scale filaments founded on the conception that these structures are bridges of matter that connect high density peaks. Our method is based on two standard tools, the Minimal Spanning Tree (MST) and the Friends of Friends (FoF) algorithm. Briefly, the complete process consists of five stages. Initially, we use the FoF algorithm to select intermediate density regions to stave off underdense zones (voids). Next, we build the MST by restricting it only to the regions defined in the previous step. Then the tree is pruned according to the length of its branches keeping the most dominant and discarding the more tenuous ones. Finally, the filaments are individualised according to the mass of the ends and smoothed using a B-spline fitting routine. To assess the results of the filament finder we apply it to a cosmological simulation. By focusing our analysis on those filaments whose halos at the ends have large masses, we found that the radial density profile, at scales around $1\, h^{-1} \mathrm{Mpc}$, follow a power law function with index -2. Even though the method only relies on halo positions, it is capable to recover the expected velocity field in filamentary structures. Large infall velocities coming from low density environment approach perpendicularly to the filaments and diverges toward to the ends. By studying the transverse velocity dispersion, we estimate the dynamical linear density following Eisenstein et al. (1997), finding a good correspondence between that and the actual mass per unit length of the filaments.

The tidal torque theory (TTT) relates the origin and evolution of angular momentum with the environment in which dark matter (DM) haloes form. The deviations introduced by late non-linearities are commonly thought as noise in the model. In this work, we analyze a cosmological simulation looking for systematics on these deviations, finding that the classification of DM haloes according to their angular momentum growth results in samples with different internal alignment, spin parameter distribution and assembly history. Based on this classification, we obtain that low mass haloes are embedded in denser environments if they have acquired angular momentum below the TTT expectations (L haloes), whereas at high masses enhanced clustering is typically associated with higher angular momentum growths (W haloes). Additionally, we find that the low mass signal has a weak dependence on the direction, whereas the high mass signal is entirely due to the structure perpendicular to the angular momentum. Finally, we study the anisotropy of the matter distribution around haloes as a function of their mass. We find that the angular momentum direction of W (L) haloes remains statistically perpendicular (parallel) to the surrounding structure across the mass range $11<\mathrm{log}(M/h^{-1}\mathrm{M}_{\odot})<14$, whereas haloes following TTT show a "spin flip" mass consistent with previously reported values ($\sim 5 \times 10^{12}$ $h^{-1}\mathrm{M}_\odot$). Hence, whether the spin flip mass of the deviated samples is highly shifted or straightly undefined, our results indicate that is remarkably connected to the haloes angular momentum growth.

We explore the effects of environment on star formation in late-type galaxies by studying the dependence of the radial profiles of specific star formation rate (sSFR) on environment and the stellar mass, using a sample of 275 late-type galaxies drawn from the CALIFA survey. We consider three different discrete environments: field galaxies, galaxies in pairs, and galaxies in groups, with stellar masses $9\le \log(M_{\star}/M_{\odot}) \le 12$, and compare their sSFR profiles across the environments. Our results suggest that the stellar mass is the main factor determining the sSFR profiles of late-type galaxies; the influence of AGNs and bars are secondary. We find that the relative size of the bulge plays a key role in depressing star formation towards the center of late-type galaxies. The group environment determines clear differences in the sSFR profiles of galaxies. We find evidence of an outside-in action upon galaxies with stellar masses $9\le \log(M_{\star}/M_{\odot}) \le 10$ in groups. We find a much stronger suppression of star formation in the inner regions of massive galaxies in groups, which may be an indication of a different merger history.

We study the mass-richness relation using galaxy catalogues and images from the Sloan Digital Sky Survey. We use two independent methods, in the first one, we calibrate the scaling relation with weak-lensing mass estimates. In the second procedure we apply a background subtraction technique to derive the probability distribution, $P(M \mid N)$, that groups with $N$-members have a virialized halo mass $M$. Lensing masses are derived in different richness bins for two galaxy systems sets: the maxBCG catalogue and a catalogue based on a group finder algorithm developed by Yang et al. MaxBCG results are used to test the lensing methodology. The lensing mass-richness relation for the Yang et al. group sample shows a good agreement with $P(M \mid N)$ obtained independently with a straightforward procedure.

Merging galaxy systems provides observational evidence of the existence of dark matter and constraints on its properties. Therefore, statistical uniform samples of merging systems would be a powerful tool for several studies. In this work we presents a new methodology for merging systems identification and the results of its application to galaxy redshift surveys. We use as starting point a mock catalogue of galaxy systems, identified using traditional FoF algorithms, which experienced a major merger as indicated by its merger tree. Applying machine learning techniques in this training sample, and using several features computed from the observable properties of galaxy members, it is possible to select galaxy groups with a high probability of have been experienced a major merger. Next we apply clustering techniques on galaxy members in order to reconstruct the properties of the haloes involved in such merger. This methodology provides a highly reliable sample of merging systems with low contamination and precise recovered properties. We apply our techniques in samples of galaxy systems obtained from SDSS-DR7, WINGS and HeCS. Our results recover previously known merging systems and provide several new candidates. We present its measured properties and discuss future analysis on current and forthcoming samples.

Context. Halo occupation distribution (HOD) is a powerful statistic that allows the study of several aspects of the matter distribution in the Universe, such as evaluating semi-analytic models of galaxy formation or imposing constraints on cosmological models. Consequently, it is important to have a reliable method for estimating this statistic, taking full advantage of the available information on current and future galaxy surveys. Aims. The main goal of this project is to combine photometric and spectroscopic information using a discount method of background galaxies in order to extend the range of absolute magnitudes and to increase the upper limit of masses in which the HOD is estimated. We also evaluate the proposed method and apply it to estimating the HOD on the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) galaxy survey. Methods. We propose the background subtraction technique to mel information provided by spectroscopic galaxy groups and photometric survey of galaxies. To evaluate the feasibility of the method, we implement the proposed technique on a mock catalogue built from a semi-analytic model of galaxy formation. Furthermore, we apply the method to the SDSS DR7 using a galaxy group catalogue taken from spectroscopic version and the corresponding photometric galaxy survey. Results. We demonstrated the validity of the method using the mock catalogue.We applied this technique to obtain the SDSS DR7 HOD in absolute magnitudes ranging from $M=-21.5$ to $M=-16.0$ and masses up to $\simeq 10^{15} M_{\odot}$ throughout this range. On the brighter extreme, we found that our results are in excellent agreement with those obtained in previous works.

In the present work, we extend the classic halo model for the large-scale matter distribution including a triaxial model for the halo profiles and their alignments. In particular, we derive general expressions for the halo-matter cross correlation function. In addition, by numerical integration, we obtain instances of the cross-correlation function depending on the directions given by halo shape axes. These functions are called anisotropic cross-correlations. With the aim of comparing our theoretical results with the simulations, we compute averaged anisotropic correlations in cones with their symmetry axis along each shape direction of the centre halo. From these comparisons we characterise and quantify the alignment of dark matter haloes on the $\Lambda$CDM context by means of the presented anisotropic halo model. As our model requires multidimensional integral computation we implement a Monte Carlo method on GPU hardware which allows us to increase the precision of the results whereas it improves the performance of the computation.

We present a study of a comparison of spin distributions of subhaloes found associated with a host halo. The subhaloes are found within two cosmological simulation families of Milky Way-like galaxies, namely the Aquarius and GHALO simulations. These two simulations use different gravity codes and cosmologies. We employ ten different substructure finders, which span a wide range of methodologies from simple overdensity in configuration space to full 6-d phase space analysis of particles.We subject the results to a common post-processing pipeline to analyse the results in a consistent manner, recovering the dimensionless spin parameter. We find that spin distribution is an excellent indicator of how well the removal of background particles (unbinding) has been carried out. We also find that the spin distribution decreases for substructure the nearer they are to the host halo's, and that the value of the spin parameter rises with enclosed mass towards the edge of the substructure. Finally subhaloes are less rotationally supported than field haloes, with the peak of the spin distribution having a lower spin parameter.

We present a detailed comparison of the substructure properties of a single Milky Way sized dark matter halo from the Aquarius suite at five different resolutions, as identified by a variety of different (sub-)halo finders for simulations of cosmic structure formation. These finders span a wide range of techniques and methodologies to extract and quantify substructures within a larger non-homogeneous background density (e.g. a host halo). This includes real-space, phase-space, velocity-space and time- space based finders, as well as finders employing a Voronoi tessellation, friends-of-friends techniques, or refined meshes as the starting point for locating substructure.A common post-processing pipeline was used to uniformly analyse the particle lists provided by each finder. We extract quantitative and comparable measures for the subhaloes, primarily focusing on mass and the peak of the rotation curve for this particular study. We find that all of the finders agree extremely well on the presence and location of substructure and even for properties relating to the inner part part of the subhalo (e.g. the maximum value of the rotation curve). For properties that rely on particles near the outer edge of the subhalo the agreement is at around the 20 per cent level. We find that basic properties (mass, maximum circular velocity) of a subhalo can be reliably recovered if the subhalo contains more than 100 particles although its presence can be reliably inferred for a lower particle number limit of 20. We finally note that the logarithmic slope of the subhalo cumulative number count is remarkably consistent and <1 for all the finders that reached high resolution. If correct, this would indicate that the larger and more massive, respectively, substructures are the most dynamically interesting and that higher levels of the (sub-)subhalo hierarchy become progressively less important.

In this paper we analyse the alignment of galaxy groups with the surrounding large scale structure traced by spectroscopic galaxies from the Sloan Digital Sky Survey Data Release 7. We characterise these alignments by means of an extension of the classical two-point cross-correlation function, developed by Paz et al. 2008 (arXiv:0804.4477, MNRAS 389 1127). We find a strong alignment signal between the projected major axis of group shapes and the surrounding galaxy distribution up to scales of 30 Mpc/h. This observed anisotropy signal becomes larger as the galaxy group mass increases, in excellent agreement with the corresponding predicted alignment obtained from mock catalogues and LCDM cosmological simulations. These measurements provide new direct evidence of the adequacy of the gravitational instability picture to describe the large-scale structure formation of our Universe.

The two-point correlation function is computed for galaxies and groups of galaxies selected using 3-dimensional information from the Updated Zwicky Galaxy Catalog - (UZC). The redshift space distortion of the correlation function $\xi(\sigma,\pi)$ in the directions parallel and perpendicular to the line of sight, induced by pairwise group peculiar velocities is evaluated. Two methods are used to characterize the pairwise velocity field of groups and galaxies. The best fitting parameters are $\gamma=-1.89 \pm 0.17$ and $r_0=9.7 \pm 4.5$ $h^{-1}$ Mpc for groups, and $\gamma=-2.00 \pm 0.03$ and $r_0=5.29 \pm 0.21$ $h^{-1}$ Mpc for galaxies. We find $\beta_{gx}=0.51 \pm 0.15$ for galaxies, while for groups we obtain a noisy estimate $\beta < 1.5$.

Our statistical studies comprise the mean relative velocity and velocity correlation for pairs of galaxies as a function of separation, and a bulk flow analysis determined on spheres of $10 h^{-1}$ $Mpc$ radius. In order to test the statistical analysis we use COBE normalized CDM numerical simulations with different density parameters and cosmological constant and take into account the Tully-Fisher (TF) scatter and a possible TF zero-point offset, as well as variations in the results of the simulations arising from different observer positions. We compare the results of the mock catalogs with samples of spiral galaxies taken from the Mark III catalog.