  We present analytical expressions for the Fourier analog of the CMB
three-point and four-point correlation functions, the spatial bispectrum and
trispectrum, of the Ostriker and Vishniac effect in the linear and mildly
nonlinear regime. Through this systematic study, we illustrate a technique to
tackle the calculation of such statistics making use of the effects of its
small-angle and vector-like properties through the Limber approximation.
Finally we discuss its configuration dependence and detectability in the
context of Gaussian theories for the currently favored flat Lambda CDM
cosmology.
