In this paper, the basic framework of PIMC method is firstly described in detail. In particular, four techniques developed for PIMC simulations to deal with the Fermi sign problem are reviewed. EOS data of various WDM systems obtained with PIMC simulations are then summarized, with which the advantages and disadvantages of PIMC method are illustrated. Finally, some future prospects of PIMC method are discussed.