separate into distinct solitary waves. They also constructed approximate solitary waves, using the resolution property, by truncating and iteratively cleaning a separated solitary wave as has been frequently done in numerical studies of other nonlinear dispersive wave equations. (Of course in this manner one does not have in general [*a priori*]{} knowledge of the speed $c_{s}$ or the value of $\tilde{\gamma}$ of the emerging solitary wave.) They used two such approximate solitary waves of different speeds to study their overtaking collision and observed that the interaction was not elastic, a fact indicating that the Benjamin equation is not integrable.

In [@CA], the authors considered solitary waves of the Benjamin equation and compared them to solitary waves of the full Euler equations for interfacial flows in the presence of surface tension when the parameters of the problem are close to the Benjamin equation regime of validity and also farther from it. The numerical scheme they used for approximating solitary waves of the Benjamin equation was based on a hybrid spatial discretization that employed fourth-order finite differences on a uniform grid for the derivatives, and the discrete Fourier transform for the nonlocal term. The resulting nonlinear system of equations was solved again 