Multi-hump bright solitons in a Schrödinger–mKdV system

We consider the problem of energy transport in a Davydov model along an anharmonic crystal medium obeying quartic longitudinal interactions corresponding to rigid interacting particles. The Zabusky and Kruskal unidirectional continuum limit of the original discrete equations reduces, in the long wave approximation, to a coupled system between the linear Schrödinger (LS) equation and the modified Korteweg–de Vries (mKdV) equation. Single- and two-hump bright soliton solutions for this LS–mKdV system are predicted to exist by variational means and numerically confirmed. The one-hump bright solitons are found to be the anharmonic supersonic analogue of the Davydov's solitons while the two-hump (in both components) bright solitons are found to be a novel type of soliton consisting of a two-soliton solution of mKdV trapped by the wave function associated to the LS equation. This two-hump soliton solution, as a two component solution, represents a new class of polaron solution to be contrasted with the two-soliton interaction phenomena from soliton theory, as revealed by a variational approach and direct numerical results for the two-soliton solution.