  We construct initial configurations for the scattering between kinks with
long-range tails. For this purpose, we exploit kink solutions in the presence
of Bogomol'nyi-Prasad-Sommerfield (BPS)-preserving impurities. This approach
offers a highly efficient method and effortless implementation with a
negligible computational cost. Our algorithm has a much smaller complexity than
the usual minimization method, becoming more than a hundred times faster in
some scenarios. Consequently, conducting kink-antikink simulations becomes
remarkably straightforward.
