  We investigate the soliton dynamics of the electromagnetic wave propagating
in an inhomogeneous or deformed ferromagnet. The dynamics of magnetization and
the propagation of electromagnetic waves are governed by the
Landau-Lifshitz-Maxwell (LLM) equation, a certain coupling between the
Landau-Lifshitz and Maxwell's equations. In the framework of multiscale
analysis, we obtain the perturbed integral modified KdV (PIMKdV) equation.
Since the dynamic is governed by the nonlinear integro-differential equation,
we rely on numerical simulations to study the interaction of its mKdV solitons
with various types of inhomogeneities. Apart from simple one soliton
experiments with periodic or localised inhomogeneities, the numerical
simulations revealed an interesting dynamical scenario where the collision of
two solitons on a localised inhomogeneity create a bound state which then
produces either two separated solitons or a mKdV breather.
