Hello,
I want to calculate spin-spin coupling and it works with the spinspin keyword in the property section. However, if I specify the atom pairs to use, e.g. spinspin 1 3 4, to calclulate only certain coupling constants, it constantly crashes saying: ngai_ddot_patch:ngai_dot_patch:g_b indices out of range :: -996
sample input:
geometry
symmetry c1
C -0.00000 0.00000 0.00000
C 0.00000 1.52903 0.00000
H -1.01818 -0.39435 0.00322
H 0.51073 -0.39466 0.88063
H 0.50637 -0.39522 -0.88284
H -0.51442 1.92344 0.87863
H -0.50458 1.92372 -0.88415
H 1.01777 1.92418 0.00586
end
dft
odft
xc ssb-d
direct
iterations 200
grid xfine
end
basis spherical
* library cc-pvtz
end
property
spinspin 1 3 4
end
task dft property
|