The breakdown phase in the plasma focus (PF) discharge is modelled by the Monte Carlo-particle in cell method. A surface ionization wave is derived. This numerical result is qualitatively confirmed in an experiment, performed by a dedicated model device with a PF-like geometry.