(p2.neigh_p==p1)        = [];
                            p2.neigh_pid(p2.neigh_pid==p1.id) = [];
                            
                            % Delete interaction object
                            delete(int);
                            rmv = true;
                        end
                    else
                        % Create new particle-particle interaction if needed
                        this.createInteractPP(drv,p1,p2,false);
                    end
                end
          