where CN A and CN B are coordination numbers of the atoms A and B in the system of interest. The k3 parameter is fixed and set to 4. As we will show, the large k3 value contributes to the presence of unphysical local minima on the PES. Details on how the k3 parameters were optimized are given below, resulting in the new expression for the D3S scheme: