F12 correct

This commit is contained in:
Anthony Scemama 2019-05-09 10:49:53 +02:00
parent d1efe67e5a
commit 5289d36b76
1 changed files with 3 additions and 4 deletions

View File

@ -28,10 +28,9 @@ let f12_integrals mo_basis =
and ijlk = F12.get_phys two_e_ints i j l k
in
if s' = Spin.other s then
0.375 *. ijkl +. 0.125 *. ijlk
(*
0.25 *. ijkl
*)
(* Minus sign because we swap spin variables
instead of orbital variables *)
0.375 *. ijkl -. 0.125 *. ijlk
else
0.25 *. (ijkl -. ijlk)
end