Dear Or,
You can approximate a diffuse reflecting surface by using the material plug-in named "Specular_Standard_Material" and setting the parameter m_sigmaSlope equal to 785.3981, which is Pi divided by 4 in mrad, and by setting the parameter m_reflectivity to your preferred value between 0 and 1.
I have tested qualitatively this approximation using the attached Tonatiuh file to simulate the reflective behaviour of a flat square diffuse reflecting surface when subjected to incident solar radiation, and as it can bee seen from the attached Tonatiuh scereen capture the rays reflected by the surface seem to be scattered in all directions.
If you could analyse quantitatively this preliminary qualitative result to convince yourself that the approach I am suggesting to simulate a diffuse reflecting surface is correct and share the results with all of us it will be great. It will be useful to test if the behaviour of the surface is independent of the direction of the sunlight as I expect it to be.
Best regards,