We present a first-principles derivation of the variational equations
describing the dynamics of the interaction of a spatial soliton and a surface
plasmon polariton (SPP) propagating along a metal/dielectric interface. The
variational ansatz is based on the existence of solutions exhibiting
differentiated and spatially resolvable localized soliton and SPP components.
These states, referred to as soliplasmons, can be physically understood as
bound states of a soliton and a SPP. Their respective dispersion relations
permit the existence of a resonant interaction between them, as pointed out in
Ref.[1]. The existence of soliplasmon states and their interesting nonlinear
resonant behavior has been validated already by full-vector simulations of the
nonlinear Maxwell's equations, as reported in Ref.[2]. Here, we provide the
theoretical demonstration of the nonlinear resonator model previously
introduced in our previous work and analyze all the approximations needed to
obtain it. We also provide some extensions of the model to improve its
applicability