We found by studying the interaction energy of neon dimer that Pauli repulsions between valence electrons on different atoms were underestimated, which caused bonds between heteroatoms to be too strong, and hydrogen bonds to be too strong as well. We show in this section that this effect causes problems in other systems as well, including systems that were previously well-described by the old version of eFF.
Cyclic alkanes with more than three carbons are nonplanar, and the old eFF reproduced well the magnitude of key dihedral angles measuring the twist of cycloalkanes out of planarity. The new eFF does slightly worse, making cyclobutane, cyclopentane, and cyclohexane more planar than they should be (Table 5.4). This is likely a consequence of the too-small barrier to ethane rotation discussed previously (1.6 kcal/mol versus 3 kcal/mol exact). We find that other intermediate range steric interactions are underestimated as well (Table 5.5), including the gauche versus trans butane interaction (0.4 versus 0.9 kcal/mol B3LYP), and the 1,3-diaxial interaction (3.3 kcal/mol versus 5.9 kcal/mol B3LYP).
The old eFF was used to study the dynamics of hydrogen plasmas, but the new eFF is no longer useful for that purpose, because it improperly stabilizes triangular and tetrahedral clusters of hydrogen atoms (triangular is 28 kcal/mol more stable than + H, and tetrahedral is 28 kcal/mol more stable than ). This is a clear consequence of insufficient Pauli repulsion between s-like electrons.
We have been curious to study bonding in electron-rich solids, and were delighted to find that the new eFF could optimize a cluster of twelve boron atoms in an icosaheral arrangement (Figure 5.27). Such an icosahedral arrangement mirrors the units found in -boron solid, and it was not possible to obtain such a structure with the old eFF. With the new eFF, the boron-boron distances are nearly identical to those found with B3LYP (1.677 versus 1.673 B3LYP), and a slight shear distortion is also reproduced.
However, we also discover with the new eFF a boron-centered amorphous jumble of atoms 600 kcal/mol more stable than the icosahedral structure; such a structure does not have any special stability according to B3LYP. It is possible to destabilize this structure by adjusting other parameters in the eFF, and to obtain stable structures for , larger boron hydrides, and carboranes; but we would prefer to stay with a single set of parameters that treated all systems consistently.
Based on the above results, it is not possible to unequivocally recommend the new eFF over the old eFF, even though the new eFF describes with higher accuracy atoms, molecules with lone pairs, atom hydrides, and single and double bonds.