Large-Large overlapping in pair_colloid

Hello,

I'm facing a problem in pair_colloid.cpp
I designed a small testcase containing just 2 particles, radius 0.01µm and varying the initial distance between them. They are underlaying DLVO potential by hybrid overlay colloid and yukawa/colloid
For larger initial distances the simulation is working. At a distance of 2µm I get an unspecific error (seems like division by zero). For distances smaller than 2µm I get the error: "Overlapping large/large in pair colloid" (line 187 of pair_colloid.cpp).

Following the equations inside pair_colloid.cpp this error message occurs when distance = radius1+radius2
In the initial state I have radius of 0.01µm and distance of for example 1µm -- despite the error message occurs.

Has anyone a hint, how to fix it?

The files attached are from a CFDEM case.

Thank you in advance
Henrik