Variational Monte Carlo calculation of dynamic multipole polarizabilities and van der Waals coefficients of the PsH system