I have added a dihedral potentialof form U(phi)=1/2*[K1(1-cos(phi))+K2(1-cos(2phi))+K3(1-cos(3phi))]