The existence of shock waves in discrete systems such as chains of two-level atoms, Heisenberg ferromagnets, nonlinear Schrödinger chains, etc., has been investigated. In particular, in the work  the existence of solitons of bright and dark type in Salerno's equation have been considered. The stability regions of these solutions have been studies as a function of the parameters of the system, finding, quite surprisingly, that on the borders of these regions anomalous dispersion is possible leading to shock waves formation. In the papers [65,P12] the dynamics of these shocks have been characterized as a function of the background wave-number inside the Brillouin zone. It was shown that shock waves with rectangular wave front followed by a train of solitons and by background radiation, can exist. In  this analysis was extended to chains of two-level atoms in the presence of exchange and dipole-dipole interaction describing Frenkel excitons. In particular it was found that the exchange interaction helps the formation of the shocks (both of bright and dark type) while the dipole-dipole interaction works against it, making unstable the background field. The possibility of shock waves in anisotropic Heisenberg quantum chains was investigated in . The analysis has been done by using the coherent state representation and the stationary phase approximation to derive quasiclassical equation of motion. It was shown that in this approximation values of the anisotropy parameter exist for which anomalous dispersion and shock wave formation is possible. From a physical point of view these shock waves separate different magnetization regions and therefore are relevant for the formation magnetic domains. This analysis has been extended to the case of the pure quantum regime in .
In  a mechanical model for a class of "multiple" sine-Gordon equations was introduced and the quasi-soliton solution of the double sine-Gordon system investigated. In particular, in the papers [19, 20] the anharmonicity effects in the small oscillations around the 4p-kink solutions of the double sine-Gordon system were considered. In the papers [4,11, 14] it has been investigated, both from an analytical and numerical point of view, the effect of adiabatic perturbations on the sine-Gordon system.