The beam propagation method assumes that there are no reflections and that the beam travels along t he propagation axis. Surface plasmons usually involve large reflections, so I think you would need to use a FDTD approach (OptiFDTD).
Changing the cladding material from air to a larger dielectric constant will also generate reflections, so I think this application is better suited for OptiFDTD.
Regards