Boiling and condensation are among the most critical and challenging topics in hydrodynamics since it involves several complex phenomena including turbulence effects, heat transfer, phase change, multiphase flow with large density ratio, etc. Particle method research community has paid attention to this topic during the past decades, but so far there have been no accurate, robust and unified particle-based boiling-condensation models for coping with such problems. The present paper aims at presenting a Lagrangian Incompressible-compressible particle method, for modeling a multiphase flow with boiling and condensation, by coupling the incompressible moving particle semi-implicit (MPS) method and the weakly compressible smoothed particle hydrodynamics (SPH) method. The coupling strategies developed by Lind et al. (JCP, 2016) are adopted. An improved potential-based surface tension model is developed. The mass transfer based on coupled gas-liquid-structure heat transfer is modeled in a straightforward manner by increasing and decreasing the gas particles’ mass, and injecting gas particles if needed. The gas volume expansion and contraction are considered by coupling the particle splitting and fusion technology. Two benchmark simulations are carried out to validate the accuracy and stability of the present MPS-SPH model. By comparing with other reference results, it is indicated that the present MPS-SPH model can give a basically-accurate prediction for the boiling and condensation phenomena (e.g. liquid-gas interface and pressures) and that the phase change model shows satisfactory performance to retain computational stability when steams bubble are growing and collapsing.