Commit 396f96d4 authored by Henning Janssen's avatar Henning Janssen

AMF double counting still not right

with the previous version for e.g. Holmium with DFT+U lowest lying states are spin-down
should make no difference for Hubbard 1
parent 4d81f0b5
......@@ -145,7 +145,7 @@ CONTAINS
! set diagonal terms and correct for non-spin-polarised case
!
DO ispin = 1,jspins
v_diag(ispin) = - u_htr * ( rho_tot - eta(ispin) ) + j_htr * ( rho_sig(ispin) - eta(ispin) )
v_diag(ispin) = - u_htr * ( rho_tot - 0.5*eta(0) ) + j_htr * ( rho_sig(ispin) - eta(ispin) )
DO m = -l,l
DO mp = -l,l
vs_mmp(m,mp,i_u,ispin) = vs_mmp(m,mp,i_u,ispin) * spin_deg
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment