  Linear Maxwell equations for transverse magnetic (TM) polarized fields
support single frequency surface plasmon polaritons (SPPs) localized at the
interface of a metal and a dielectric. Metals are typically dispersive, i.e.
the dielectric function depends on the frequency. We prove the bifurcation of
localized SPPs in dispersive media in the presence of a cubic nonlinearity and
provide an asymptotic expansion of the solution and the frequency. The problem
is reduced to a system of nonlinear differential equations in one spatial
dimension by assuming a plane wave dependence in the direction tangential to
the (flat) interfaces. The number of interfaces is arbitrary and the nonlinear
system is solved in a subspace of functions with the $H^1$-Sobolev regularity
in each material layer. The corresponding linear system is an operator pencil
in the frequency parameter due to the material dispersion. Because of the
TM-polarization the problem cannot be reduced to a scalar equation.
  The studied bifurcation occurs at a simple isolated eigenvalue of the pencil.
For geometries consisting of two or three homogeneous layers we provide
explicit conditions on the existence of eigenvalues and on their simpleness and
isolatedness. Real frequencies are shown to exist in the nonlinear setting in
the case of PT-symmetric materials. We also apply a finite difference numerical
method to the nonlinear system and compute bifurcating curves.
