Voids in gravitational instability scenarios – I. Global density and velocity fields in an Einstein–de Sitter universe

The first results of an extensive study of the structure and dynamics of underdense regions in gravitational instability scenarios are presented. Instead of adopting spherically symmetric voids with some idealized initial density and velocity profile, underdense regions of a given size and depth, embedded in an initial density fluctuation field, are generated. In order to accomplish this in a consistent way, these initial conditions are set up by means of Bertschinger's constrained random field code. The generated particle samples of 643 particles in a box of side 100 Mpc are followed into the non-linear regime by Bertschinger's PM N- body code. In this way we address the dependence of the structure and kinematics of the void both on the initial depth of the void and on the fluctuation field in which it is embedded. In particular, this study provides some understanding of how far fluctuations on small scales modify the dynamics of the large-scale void, and especially of how far the properties of small structures inside the void are affected by the global properties of the void. One of the conspicuous features of the initial density fields inside protovoids appears to be the existence of a ‘void hierarchy’, with small voids embedded in larger voids. The survival of this hierarchy during the riot evolution of the void depends critically on the initial depth as well as on the clustering scenario involved. As well as presenting a qualitative discussion of the structure of underdense regions in initial density fields in different scenarios, and the results of simulations of the ensuing non-linear evolution, we concentrate in particular on a comparison of the global density and velocity fields in voids with predictions from linear theory as well as from the spherical outflow model. The relation between the initial linear depth, the resulting non-linear depth and the excess expansion velocities in voids is addressed. In addition, we find that, while near its centre a void becomes more and more spherical, the shape of its boundary is influenced to a large extent by the structures surrounding the void and therefore is generally more irregular. In this first study we concentrate on single voids in Einstein–de Sitter universes. The underdense regions considered are linear 1 σ0, 2 σ0 and 3 σ0 dips in fields that are Gaussian-smoothed on a scale of RG = 10 h–1 Mpc, approximately half the size of the Bootes void. These regions are studied in terms of the Cold Dark Matter and Hot Dark Matter scenarios as well as in terms of the scale-free scenarios P(k) ∝ k0, k–1 and k–2. The Hubble constant is taken to be H0 = 100 h km s–1 Mpc–1.