with D 0 being the maximum value of the given descriptor and k being a constant. We assume a bond to be broken if the change in the descriptor is smaller than a cutoff threshold γ cleave . This is done by investigating the first derivative of D f it with respect to r. The usage of γ cleave usually yields reasonable results for values 0.1% to 1%.