We study the existence and stability of multibreathers in Klein-Gordon chains with interactions that are not restricted to nearest neighbors. We provide a general framework where such long range effects can be taken into consideration for arbitrarily
varying (as a function of the node distance) linear couplings between arbitrary sets of neighbors in the chain. By examining special case examples such as three-site breathers with next-nearest-neighbors, we find {it crucial} modifications to the nearest-neighbor picture of one-dimensional oscillators being excited either in- or anti-phase. Configurations with nontrivial phase profiles, arise, as well as spontaneous symmetry breaking (pitchfork) bifurcations, when these states emerge from (or collide with) the ones with standard (0 or $pi$) phase difference profiles. Similar bifurcations, both of the supercritical and of the subcritical type emerge when examining four-site breathers with either next-nearest-neighbor or even interactions with the three-nearest one-dimensional neighbors. The latter setting can be thought of as a prototype for the two-dimensional building block, namely a square of lattice nodes, which is also examined. Our analytical predictions are found to be in very good agreement with numerical results.
We study the existence and stability of multisite discrete breathers in two prototypical non-square Klein-Gordon lattices, namely a honeycomb and a hexagonal one. In the honeycomb case we consider six-site configurations and find that for soft potent
ial and positive coupling the out-of-phase breather configuration and the charge-two vortex breather are linearly stable, while the in-phase and charge-one vortex states are unstable. In the hexagonal lattice, we first consider three-site configurations. In the case of soft potential and positive coupling, the in-phase configuration is unstable and the charge-one vortex is linearly stable. The out-of-phase configuration here is found to always be linearly unstable. We then turn to six-site configurations in the hexagonal lattice. The stability results in this case are the same as in the six-site configurations in the honeycomb lattice. For all configurations in both lattices, the stability results are reversed in the setting of either hard potential or negative coupling. The study is complemented by numerical simulations which are in very good agreement with the theoretical predictions. Since neither the form of the on-site potential nor the sign of the coupling parameter involved have been prescribed, this description can accommodate inverse-dispersive systems (e.g., supporting backward waves) such as transverse dust-lattice oscillations in dusty plasma (Debye) crystals or analogous modes in molecular chains.
