Studies of large-scale structures in the Universe, such as superstructures or cosmic voids, have been widely used to characterize the properties of the cosmic web through statistical analyses. On the other hand, the 2-point correlation function of large-scale tracers such as galaxies or halos provides a reliable statistical measure. However, this function applies to the spatial distribution of point-like objects, and therefore it is not appropriate for extended large structures which strongly depart from spherical symmetry. Here we present an analysis based on the standard correlation function formalism that can be applied to extended objects exhibiting arbitrary shapes. Following this approach, we compute the probability excess $\Xi$ of having spheres sharing parts of cosmic structures with respect to a realization corresponding to a distribution of the same structures in random positions. For this aim, we identify superstructures defined as Future Virialized Structures (FVSs) in semi-anaytic galaxies on the MPDL2 MultiDark simulation. We have also identified cosmic voids to provide a joint study of their relative distribution with respect to the superstructures. Our analysis suggests that $\Xi$ provides useful characterizations of the large scale distribution, as suggested from an analysis of sub-sets of the simulation. Even when superstructure properties may exhibit negligible variations across the sub-sets, $\Xi$ has the sensitivity to statistically distinguish sub-boxes that departs from the mean at larger scales. Thus, our methods can be applied in analysis of future surveys to provide characterizations of large-scale structure suitable to distinguish different theoretical scenarios.

The expanding complex pattern of filaments, walls and voids build the evolving cosmic web with material flowing from underdense onto high density regions. Here we explore the dynamical behaviour of voids and galaxies in void shells relative to neighboring overdense superstructures, using the Millenium Simulation and the main galaxy catalogue in Sloan Digital Sky Survey data. We define a correlation measure to estimate the tendency of voids to be located at a given distance from a superstructure. We find voids-in-clouds (S-types) preferentially located closer to superstructures than voids-in-voids (R-types) although we obtain that voids within $\sim40~\mathrm{Mpc}\,\mathrm{h}^{-1}$ of superstructures are infalling in a similar fashion independently of void type. Galaxies residing in void shells show infall towards the closest superstructure, along with the void global motion, with a differential velocity component depending on their relative position in the shell with respect to the direction to the superstructure. This effect is produced by void expansion and therefore is stronger for R-types. We also find that galaxies in void shells facing the superstrucure flow towards the overdensities faster than galaxies elsewere at the same relative distance to the superstructure. The results obtained for the simulation are also reproduced for the SDSS data with a linearized velocity field implementation.

We analyse the clustering of cosmic voids using a numerical simulation and the main galaxy sample from the Sloan Digital Sky Survey. We take into account the classification of voids into two types that resemble different evolutionary modes: those with a rising integrated density profile (void-in-void mode, or R-type) and voids with shells (void-in-cloud mode, or S-type). The results show that voids of the same type have stronger clustering than the full sample. We use the correlation analysis to define void clumps, associations with at least two voids separated by a distance of at most the mean void separation. In order to study the spatial configuration of void clumps, we compute the minimal spanning tree and analyse their multiplicity, maximum length and elongation parameter. We further study the dynamics of the smaller sphere that encloses all the voids in each clump. Although the global densities of void clumps are different according to their member-void types, the bulk motions of these spheres are remarkably lower than those of randomly placed spheres with the same radii distribution. In addition, the coherence of pairwise void motions does not strongly depend on whether voids belong to the same clump. Void clumps are useful to analyse the large-scale flows around voids, since voids embedded in large underdense regions are mostly in the void-in-void regime, were the expansion of the larger region produces the separation of voids. Similarly, voids around overdense regions form clumps that are in collapse, as reflected in the relative velocities of voids that are mostly approaching.

We perform a statistical study of the global motion of cosmic voids using both a numerical simulation and observational data. We analyse their relation to large-scale mass flows and the physical effects that drive those motions. We analyse the bulk motions of voids, defined by the mean velocity of haloes in the surrounding shells in the numerical simulation, and by galaxies in the Sloan Digital Sky Survey Data Release 7. We find void mean bulk velocities close to 400 km/s, comparable to those of haloes (~ 500-600 km/s), depending on void size and the large-scale environment. Statistically, small voids move faster than large ones, and voids in relatively higher density environments have higher bulk velocities than those placed in large underdense regions. Also, we analyze the mean mass density around voids finding, as expected, large-scale overdensities (underdensities) along (opposite to) the void motion direction, suggesting that void motions respond to a pull-push mechanism. This contrasts with massive cluster motions who are mainly governed by the pull of the large-scale overdense regions. Our analysis of void pairwise velocities shows how their relative motions are generated by large-scale density fluctuations. In agreement with linear theory, voids embedded in low (high) density regions mutually recede (attract) each other, providing the general mechanism to understand the bimodal behavior of void motions. In order to compare the theoretical results and the observations we have inferred void motions in the SDSS using linear theory, finding that the estimated observational void motions are in statisticalagreement with the results of the simulation. Regarding large-scale flows, our results suggest a scenario of galaxies and galaxy systems flowing away from void centers with the additional, and morerelevant, contribution of the void bulk motion to the total velocity.

We compute the bulk motions of cosmic voids, using a $\Lambda$CDM numerical simulation considering the mean velocities of the dark matter inside the void itself and that of the haloes in the surrounding shell. We find coincident values of these two measures in the range $\sim$ 300-400 km/s, not far from the expected mean peculiar velocities of groups and galaxy clusters. When analysing the distribution of the pairwise relative velocities of voids, we find a remarkable bimodal behaviour consistent with an excess of both systematically approaching and receding voids. We determine that the origin of this bimodality resides in the void large scale environment, since once voids are classified into void-in-void (R-type) or void-in-cloud (S-type), R-types are found mutually receding away, while S-types approach each other. The magnitude of these systematic relative velocities account for more than 100 km/s, reaching large coherence lengths of up to 200 h$^{-1}$ Mpc . We have used samples of voids from the Sloan Digital Sky Survey Data Release 7 (SDSS-DR7) and the peculiar velocity field inferred from linear theory, finding fully consistent results with the simulation predictions. Thus, their relative motion suggests a scenario of a sparkling Universe, with approaching and receding voids according to their local environment.

We study the dependence of the properties of group galaxies on the surrounding large-scale environment, using SDSS-DR7 data. Galaxies are ranked according to their luminosity within each group and classified morphologically by the S\'ersic index. We have considered samples of the host groups in superstructures of galaxies, and elsewhere. We find a significant dependence of the properties of late-type brightest group galaxies on the large-scale environment: they show statistically significant higher luminosities and stellar masses, redder u-r colours, lower star formation activity and longer star-formation time-scale when embedded in superstructures. By contrast, the properties of the early-type brightest group galaxies are remarkably similar regardless of the group global environment. The other group member galaxies exhibit only the local influence of the group they inhabit. Our analysis comprises tests against the dependence on the host group luminosity and we argue that group brightest member properties are not only determined by the host halo, but also by the large-scale structure which can influence the accretion process onto their late-type brightest galaxies.

Inspired on the well known dynamical dichotomy predicted in voids, where some underdense regions expand whereas others collapse due to overdense surrounding regions, we explored the interplay between the void inner dynamics and its large scale environment. The environment is classified depending on its density as in previous works. We analyse the dynamical properties of void-centered spherical shells at different void-centric distances depending on this classification. The above dynamical properties are given by the angular distribution of the radial velocity field, its smoothness, the field dependence on the tracer density and shape, and the field departures from linear theory. We found that the velocity field in expanding voids follows more closely the linear prediction, with a more smooth velocity field. However when using velocity tracers with large densities such deviations increase. Voids with sizes around $18\,h^{-1}\,Mpc$ are in a transition regime between regions with expansion overpredicted and underpredicted from linear theory. We also found that velocity smoothness increases as the void radius, indicating the laminar flow dominates the expansion of larger voids (more than $18\,h^{-1}\,Mpc$). The correlations observed suggest that nonlinear dynamics of the inner regions of voids could be dependent on the evolution of the surrounding structures. These also indicate possible scale couplings between the void inner expansion and the large scale regions where voids are embedded. These results shed some light to the origin of nonlinearities in voids, going beyond the fact that voids just quickly becomes nonlinear as they become emptier.

We analyse properties of galaxy groups and their dependence on the large scale environment as defined by superstructures. We find that group galaxy cross correlations depend only on group properties regardless the groups reside in superstructures. This indicates that the total galaxy density profile around groups is independent of the global environment. At a given global luminosity, a proxy to group total mass, groups have a larger stellar mass content by a factor 1.3, a relative excess independent of the group luminosity. Groups in superstructures have 40 per cent higher velocity dispersions and systematically larger minimal enclosing radii. We also find that the stellar population of galaxies in groups in superstructures is systematically older as infered from the galaxy spectra Dn4000 parameter. Although the galaxy number density profile of groups is independent of environment, the star formation rate and stellar mass profile of the groups residing in superstructures differs from groups elsewhere. For groups residing in superstructures, the combination of a larger stellar mass content and star formation rate produces a larger time scale for star formation regardless the distance to the group center. Our results provide evidence that groups in superstructures formed earlier than elsewhere, as expected in the assembly bias scenario.

We analyse the dependence of clustering properties of galaxies as a function of their large-scale environment. In order to characterize the environment on large scales, we use the catalogue of future virialized superstructures (FVS) by Luparello et al. and separate samples of luminous galaxies according to whether or not they belong to FVS. In order to avoid biases in the selection of galaxies, we have constructed different subsamples so that the distributions of luminosities and masses are comparable outside and within FVS. As expected, at large scales, there is a strong difference between the clustering of galaxies inside and outside FVS. However, this behaviour changes at scales r < 1 Mpc/h, where the correlations have similar amplitudes. The amplitude of the two-halo term of the correlation function for objects inside FVS does not depend on their mass, but rather on that of the FVS. This is confirmed by comparing this amplitude with that expected from extended Press-Schechter fits. In order to compare these observational results with current models for structure formation, we have performed a similar analysis using a semi-analytic implementation in a LCDM cosmological model. We find that the cross-correlation functions from the mock catalogue depend on the large-scale structures in a similar way to the observations. From our analysis, we conclude that the clustering of galaxies within the typical virialized regions of groups, mainly depends on the halo mass, irrespective of the large-scale environment.

We construct catalogues of present superstructures that, according to a LCDM scenario, will evolve into isolated, virialized structures in the future. We use a smoothed luminosity density map derived from galaxies in SDSS-DR7 data and separate high luminosity density peaks. The luminosity density map is obtained from a volume-limited sample of galaxies in the spectroscopic galaxy catalogue, within the SDSS-DR7 footprint area and in the redshift range 0.04 < z < 0.12. Other two samples are constructed for calibration and testing purposes, up to z = 0.10 and z = 0.15. The luminosity of each galaxy is spread using an Epanechnikov kernel of 8Mpc/h radius, and the map is constructed on a 1 Mpc/h cubic cells grid. Future virialized structures (FVS) are identified as regions with overdensity above a given threshold, calibrated using a LCDM numerical simulation, and the criteria presented by D\"unner et al. (2006). We assume a constant mass-to-luminosity ratio and impose the further condition of a minimum luminosity of 10^{12}Lsol. According to our calibrations with a numerical simulation, these criteria lead to a negligible contamination by less overdense (non FVS) superstructures.We present a catalogue of superstructures in the SDSS-DR7 area within redshift 0.04 < z < 0.12 and test the reliability of our method by studying different subsamples as well as a mock catalogue.We compute the luminosity and volume distributions of the superstructures finding that about 10% of the luminosity (mass) will end up in future virialized structures. The fraction of groups and X-ray clusters in these superstructures is higher for groups/clusters of higher mass, suggesting that future cluster mergers will involve the most massive systems. We also analyse known structures in the present Universe and compare with our catalogue of FVS.