.id) = [];
                            p2.interacts(p2.interacts==int)   = [];
                            p2.neigh_p(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
 