  We have measured the specific heat of zincblende ZnS for several isotopic
compositions and over a broad temperature range (3 to 1100 K). We have compared
these results with calculations based on ab initio electronic band structures,
performed using both LDA and GGA exchange- correlation functionals. We have
compared the lattice dynamics obtained in this manner with experimental data
and have calculated the one-phonon and two-phonon densities of states. We have
also calculated mode Grueneisen parameters at a number of high symmetry points
of the Brillouin zone. The electronic part of our calculations has been used to
investigate the effect of the 3d core electrons of zinc on the spin-orbit
splitting of the top valence bands. The effect of these core electrons on the
band structure of the rock salt modification of ZnS is also discussed.
