jensengroup/xyz2mol

Purpose of 1.3 multiplier in checking covalent bonds

ljmartin opened this issue · 2 comments

Thanks for the great script.

Would someone please explain the purpose of the GetRCovalent covalent radius being multiplied by 1.3?
ie:

Rcov_i = pt.GetRcovalent(a_i.GetAtomicNum())*1.30

thanks

Hi Martin,

When you want to test of there is a covalent bond between to pairs of atoms you use the covalent radii of the atoms and check if the distance between the pair of atoms is shorter. It turns out that using only the covalent radii (from GetRCovalent) will not find bonds that are slightly longer then their equilibrium distance. So you add this fudge factor of 1.3 to make sure you also get those.

Great, thanks!