Hopefully fix LT with Ewald

SunnySuite · Dec 18, 2024 · 55cf954 · 55cf954
1 parent b659055
commit 55cf954
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/Spiral/LuttingerTisza.jl b/src/Spiral/LuttingerTisza.jl
@@ -25,7 +25,7 @@ function luttinger_tisza_exchange(sys::System; k, ϵ=0)
         A = Sunny.precompute_dipole_ewald_at_wavevector(sys.crystal, (1,1,1), k) * sys.ewald.μ0_μB²
         A = reshape(A, Na, Na)
         for i in 1:Na, j in 1:Na
-            J_k[:, i, :, j] += gs[i]' * A[i, j] * gs[j] / 2
+            J_k[:, i, :, j] += sys.gs[i]' * A[i, j] * sys.gs[j] / 2
         end
     end