To get random rotation with POP particles, add a point wrangle at SOP level with a more specific version of the following code-
v@w=rand(x-value,y-value,z-value);
Then in the dop network, add a POP-torque node and set the amount to zero, so that it takes the initial "w" value. "w" is angular velocity..all according to keikaku.