                  % Check for existing interaction
                    if (any(p1.neigh_pid == p2.id))
                        % Get interaction object
                        int = findobj(p1.interacts,'elem2',p2);
                        
                        % Compute separation between elements
                        int.kinemat = int.kinemat.setRelPos(p1,p2);
                        
                        % Check if separation is greater than cutoff distance
                        % Assumption: 