NRPMD This is a imaginary sdsaads ssss def dHel(R): """ dv_ij/dR """ dHij = np.zeros((2,2,len(R))) ω,c = initModelParams() dHij[0,0,:] = c dHij[1,1,:] = -c return dHij $\alpha^2\sum$