We present a detailed 1.2 mm continuum and CS spectral line study of a large sample of 69 massive star forming regions in very early stages of evolution, most of them prior to building up an ultracompact HII region. The continuum data show a zoo of different morphologies and give detailed information on the spatial distributions, the masses, column densities and average densities of the whole sample. Fitting the radial intensity profiles shows that three parameters are needed to describe the spatial distribution of the sources: constant emission from the center out to a few arcsec radius followed by a first power law intensity distribution which steepens further outside into a second power law distribution. The mean inner power law intensity index mi (I~r^(-mi)) is 1.2 corresponding to density indices p (n~r^(-p)) of 1.6. In total the density distribution of our massive star formations sites seem to be not too different from their low-mass counterparts, but we show that setting tight constrains on the density indices is very difficult and subject to many possible errors. The local densities we derive from CS calculations are higher (up to one order of magnitude) than the mean densities we find via the mm-continuum. Such inhomogeneous density distribution reflects most likely the ubiquitous phenomenon of clumping and fragmentation in molecular clouds. Linewidth-mass relations show a departure from virial equilibrium in the stages of strongly collapsing cores.

15 pages, 13 jpeg-figures. Astrophysical Journal, in press