Motivated by the recent detection of an enhanced clustering signal along the major axis of haloes in N-body simulations, we derive a formula for the anisotropic density distribution around haloes and voids on large scales. Our model, which assumes linear theory and that the formation and orientation of nonlinear structures are strongly correlated with the Lagrangian shear, is in good agreement with measurements. We also show that the measured amplitude is inconsistent with a model in which the alignment is produced by the initial inertia rather than shear tensor.

We build a model for the density and integrated Sachs-Wolfe (ISW) profile of supervoid and supercluster structures. Our model assumes that fluctuations evolve linearly from an initial Gaussian random field. We find these assumptions capable of describing N-body simulations and simulated ISW maps remarkably well on large scales. We construct an ISW map based on locations of superstructures identified previously in the SDSS Luminous Red Galaxy sample. A matched filter analysis of the cosmic microwave background confirms a signal at the $3.2-\sigma$ confidence level and estimates the radius of the underlying structures to be $55 \pm 28 ~h^{-1}$Mpc. The amplitude of the signal, however, is $2-\sigma$ higher than $\Lambda$CDM predictions.

We measure the matter probability distribution function (PDF) via counts in cells in a volume limited subsample of the Sloan Digital Sky Survey Luminous Red Galaxy Catalog on scales from $30 h^{-1}$Mpc to $150 h^{-1}$Mpc and estimate the linear Integrated Sachs--Wolfe effect produced by supervoids and superclusters in the tail of the PDF. We characterize the PDF by the variance, $S_3$, and $S_4$, and study in simulations the systematic effects due to finite volume, survey shape and redshift distortion. We compare our measurement to the prediction of $\Lambda$CDM with linear bias and find a good agreement. We use the moments to approximate the tail of the PDF with analytic functions. A simple Gaussian model for the superstructures appears to be consistent with the claim by Granett et al. that density fluctuations on $100 h^{-1}$Mpc scales produce hot and cold spots with $\Delta T \approx 10\mu K$ on the cosmic microwave background.