We study the dynamics of the collision of two solitary waves for the 2 and 3-dimensional Zakharov–Kuznetsov equation, a high-dimensional non-integrable version of the Korteweg-de Vries equation that appears as an asymptotic model in plasma physics. We describe the evolution of the solution behaving as a sum of 2-solitary waves of nearly equal speeds at time $$t=-\infty $$
t
=
-
∞
up to time $$t=+\infty $$
t
=
+
∞
. We show that this solution behaves as the sum of two modulated solitary waves and an error term which is small in $$H^1$$
H
1
for all time $$t \in {\mathbb {R}}$$
t
∈
R
. Finally, we also prove the stability of this solution for large times around the collision. The proofs are a non-trivial extension of the ones of Martel and Merle for the quartic generalized Korteweg-de Vries equation to higher dimensions. First, despite the non-explicit nature of the solitary wave, we construct an approximate solution in an intrinsic way by canceling the error to the equation only in the natural directions of scaling and translation. Then, to control the difference between a solution and the approximate solution, we use a modified energy functional and a refined modulation estimate in the transverse variable. Moreover, we rely on the Hamiltonian structure of the ODE governing the distance between the waves, which cannot be approximated by explicit solutions, to close the bootstrap estimates on the parameters. We hope that the techniques introduced here are robust and will prove useful in studying the collision phenomena for other focusing non-linear dispersive equations with non-explicit solitary waves.