From 6cf634cbf67ff00de4b813905397eaaa33eaee9d Mon Sep 17 00:00:00 2001
From: Kostas Vilkelis <kostasvilkelis@gmail.com>
Date: Wed, 31 May 2023 12:16:18 +0200
Subject: [PATCH] fix normalisation bug
---
codes/hf.py | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/codes/hf.py b/codes/hf.py
index b3ab4c6..d6f9ef9 100644
--- a/codes/hf.py
+++ b/codes/hf.py
@@ -51,7 +51,7 @@ def compute_mf(vals, vecs, filling, H_int):
E_F = utils.get_fermi_energy(vals, filling)
F = mean_field_F(vals, vecs, E_F=E_F)
rho = np.diag(np.average(F, axis=tuple([i for i in range(dim)])))
- exchange_mf = convolution(F, H_int) * nk ** (-2)
+ exchange_mf = convolution(F, H_int) * nk ** (-dim)
direct_mf = np.diag(np.einsum("i,ij->j", rho, H0_int))
return direct_mf - exchange_mf
--
GitLab