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.

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.

Cosmic voids are becoming key players in testing the physics of our Universe. Here we concentrate on the abundances and the dynamics of voids as these are among the best candidates to provide information on cosmological parameters. Cai, Padilla \& Li (2014) use the abundance of voids to tell apart Hu \& Sawicki $f(R)$ models from General Relativity. An interesting result is that even though, as expected, voids in the dark matter field are emptier in $f(R)$ gravity due to the fifth force expelling away from the void centres, this result is reversed when haloes are used to find voids. The abundance of voids in this case becomes even lower in $f(R)$ compared to GR for large voids. Still, the differences are significant and this provides a way to tell apart these models. The velocity field differences between $f(R)$ and GR, on the other hand, are the same for halo voids and for dark matter voids. Paz et al. (2013), concentrate on the velocity profiles around voids. First they show the necessity of four parameters to describe the density profiles around voids given two distinct void populations, voids-in-voids and voids-in-clouds. This profile is used to predict peculiar velocities around voids, and the combination of the latter with void density profiles allows the construction of model void-galaxy cross-correlation functions with redshift space distortions. When these models are tuned to fit the measured correlation functions for voids and galaxies in the Sloan Digital Sky Survey, small voids are found to be of the void-in-cloud type, whereas larger ones are consistent with being void-in-void. This is a novel result that is obtained directly from redshift space data around voids. These profiles can be used to remove systematics on void-galaxy Alcock-Pacinsky tests coming from redshift-space distortions.

Using the redshift-space distortions of void-galaxy cross-correlation function we analyse the dynamics of voids embedded in different environments. We compute the void-galaxy cross-correlation function in the Sloan Digital Sky Survey (SDSS) in terms of distances taken along the line of sight and projected into the sky. We analyse the distortions on the cross-correlation isodensity levels and we find anisotropic isocontours consistent with expansion for large voids with smoothly rising density profiles and collapse for small voids with overdense shells surrounding them. Based on the linear approach of gravitational collapse theory we developed a parametric model of the void-galaxy redshift space cross-correlation function. We show that this model can be used to successfully recover the underlying velocity and density profiles of voids from redshift space samples. By applying this technique to real data, we confirm the twofold nature of void dynamics: large voids typically are in an expansion phase whereas small voids tend to be surrounded by overdense and collapsing regions. These results are obtained from the SDSS spectroscopic galaxy catalogue and also from semi-analytic mock galaxy catalogues, thus supporting the viability of the standard LCDM model to reproduce large scale structure and dynamics.

We perform a statistical study focused on void environments. We examine galaxy density profiles around voids in the SDSS, finding a correlation between void--centric distance to the shell of maximum density and void radius when a maximum in overdensity exists. We analyze voids with and without a surrounding over-dense shell in the SDSS. We find that small voids are more frequently surrounded by over-dense shells whereas the radial galaxy density profile of large voids tends to rise smoothly towards the mean galaxy density. We analyse the fraction of voids surrounded by overdense shells finding a continuous trend with void radius. The differences between voids with and without an overdense shell around them can be understood in terms of whether the voids are, on average, in the process of collapsing or continuing their expansion, respectively, in agreement with previous theoretical expectations. We use numerical simulations coupled to semi-analytic models of galaxy formation in order to test and interpret our results. The very good agreement between the mock catalog results and the observations provides additional support to the viability of a LCDM model to reproduce the large scale structure of the universe as defined by the void network, in a way which has not been analysed previously.

We study the relative fraction of low and high surface brightness galaxies (LSBGs and HSBGs) at void walls in the SDSS DR7. We focus on galaxies in equal local density environments. We assume that the host dark-matter halo mass (for which we use SDSS group masses) is a good indicator of local density. This analysis allows to examine the behavior of the abundance of LSBG and HSBG galaxies at a fixed local density and distinguish the large-scale environment defined by the void geometry. We compare galaxies in the field, and in the void walls; the latter are defined as the volume of void shells of radius equal to that of the void. We find a significant decrement, a factor $\sim 4$, of the relative fraction of blue, active star-forming LSBGs in equal mass groups at the void walls and the field. This decrement is consistent with an increase of the fraction of blue, active star-forming HSBGs. By contrast, red LSBGs and HSBGs show negligible changes. We argue that these results are consistent with a scenario where LSBGs with blue colors and strong star formation activity at the void walls are fueled by gas from the expanding void regions. This process could lead to LSBG to HSBG transformations.

We perform a statistical analysis of the peculiar velocity field around dark matter haloes in numerical simulations. We examine different properties of the infall of material onto haloes and its relation to central halo shapes and the shape of the large scale surrounding regions (LSSR). We find that the amplitude of the infall velocity field along the halo shape minor axis is larger than that along the major axis. This is consistent for general triaxial haloes, and for both prolate and oblate systems. We also report a strong anisotropy of the velocity field along the principal axes of the LSSR. The infall velocity field around dark matter haloes reaches a maximum value along the direction of the minor axis of the LSSR, whereas along the direction of its major axis, it exhibits the smallest velocities. We also analyse the dependence of the matter velocity field on the local environment. The amplitude of the infall velocity at high local density regions is larger than at low local density regions. The velocity field tends to be more laminar along the direction towards the minor axis of the LSSR, where the mean ratio between flow velocity and velocity dispersion is of order unity and nearly constant up to scales of 15 Mpc/h. We also detect anisotropies in the outflowing component of the velocity field, showing a maximum amplitude along the surrounding LSSR major axis.

We perform a statistical study on the distribution and dynamics of voids in the 2dFGRS. Our statistics are tested and calibrated using mock 2dFGRS catalogues. We analyse the dynamics around voids in mock and real 2dFGRS surveys. The void-galaxy cross-correlation redshift-space distortions show evidence of the continuing growth of voids, confirming another prediction of the hierarchical clustering scenario. A non-linear outflow model can be used to provide quantitative estimates of the outflow velocities around 2dFGRS voids. These results are consistent with maximum outflows of 110km/s, 210km/s and 270km/s for voids of <r_void>=7.5, 12.5 and 17.5h^-1Mpc, assuming a galaxy bias of b=1. As an application for future surveys, our study of the mock catalogues shows that direct measurements of the expansion of voids can be obtained using peculiar velocity data. We find that it would also be possible to detect differences in the velocity dispersion of galaxies in the directions parallel and perpendicular to the void walls.

We perform a statistical study of the characteristics of galaxies in voids and void walls in the SDSS and 2dFGRS catalogues. We investigate dependencies of the distribution of galaxy spectral types and colours as a function of the relative position to the void centres for different luminosity and local density ranges. We find a trend towards bluer, star forming galaxies in void walls beyond the local density dependence. These results indicate that luminosity and local density do not entirely determine the distribution of galaxy properties such as colours and spectral types, and point towards a large scale modulation of star formation. We argue that this effect is due to the lower accretion and merger history of galaxies arriving at void walls from the emptier inner void regions.

We study statistical properties of voids in the distribution of mass, dark-matter haloes and galaxies (B_J<-16) in a LCDM numerical simulation populated with galaxies using a semi-analytic galaxy formation model(GALFORM, Cole et al. 2000). We find that the properties of voids selected from GALFORM galaxies are compatible with those of voids identified from a population of haloes with mass M>10^11.5 M_sun/h, similar to the median halo mass, M_med=10^11.3 M_sun/h. We also find that the number density of galaxy- and halo-defined voids can be up to two orders of magnitude higher than mass-defined voids for large void radii. As expected, there are outflow velocities which show their maximum at larger void-centric distances for larger voids (well described by a linear relation). The void-centric distance where this maximum occurs, follows a suitable power law fit of the form, log(d_vmax)=(r_void/A)^B. At sufficiently large distances, we find mild infall motions onto the sub-dense regions. We find that a similar analysis in redshift space would make both outflows and infalls to appear with a lower amplitude. We also find that the velocity dispersion of galaxies and haloes is larger in the direction parallel to the void walls by ~10-20%. We analyse the void-tracer cross-correlation functions in real- and redshift-space as a function of separation, and also as a function of separations parallel and perpendicular to the line of sight. The distortion pattern observed in xi(sigma,pi) is that of an elongation along the line of sight which extends out to large separations. Positive xi contours evidence finger-of-god motions at the void walls. Elongations along the line of sight are roughly comparable between galaxy-, halo- and mass-defined voids. (Abridged)