Hello,
I believe this discontinuity is expected. The component of the electric field normal to the boundary between two materials is discontinuous, according to the electromagnetic boundary conditions. The normal component of the electric displacement field is continuous. I'm not sure why this discontinuity wouldn't appear in the other software.