tance
            if (addVL && separ < this.verlet_dist)
                p1.verlet_p(end+1) = p2;
            end
            
            % Check if interaction exists
            % Assumption: cutoff ratio applies to maximum radius
            if (separ >= this.cutoff * max(p1.radius,p2.radius))
                return;
            end
            
            % Create new interaction object by copying base object
            int = copy(this.b_interact);
            
            % Set handles to interectin