Elastic molecular collision models for Monte Carlo simulation of rarefied gas flow, such as the variable hard sphere (VHS), variable soft sphere (VSS), and variable sphere (VS) molecular models, were extended to quantum mechanical scattering, and then applied to the simulation of spherical expansion of He4. The profiles of number density, flow velocity, and parallel and perpendicular kinetic temperatures obtained with the VHS, VSS, and VS models were in good agreement with those yielded by direct calculation of quantum mechanical scattering by Koura [Phys. Fluids 11, 3174 (1999)].