In this paper we use the method of stochastic
integral due to Gaveau to construct the heat kernel for the
quaternionic Heisenberg groups, and then follow the line of
Coulhon et al. to deduce the uniformly boundedness of the Riesz
transforms on these nilponent Lie groups.