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.