<div dir="ltr">Dear Marco,<br><br>just check Eqs. 5 and 6 of Dudarev et al. PRB 57(3), 1505 (1998). Please, note the Kronecker delta in Eq 6.<br><br>Best regards,<br><br>Matthias<br><br><br>On Monday, 9 February 2015 19:02:24 UTC+1, Marco  wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hello,<br><br>Sorry about this simple question but I was hoping someone could explain to me how the potential is calculated in the DFT+U Mulliken method.<br>In lines 1379-1383 of the dft_plus_u.f90 file, the potential dE[U]/dP(i,j) is defined as:<br><br>IF (isgf == jsgf) THEN<br>  v_block(isgf,isgf) = u_minus_j*(0.5_dp - fspin*q_matrix(i,i))<br>ELSE<br>  v_block(isgf,jsgf) = -u_minus_j*fspin*q_matrix(j,i)<br>END IF<br><br>Why is the v_block(isgf,jsgf) not defined as:<br>v_block(isgf,jsgf) = u_minus_j*(0.5_dp - fspin*q_matrix(j,i))<br><br>Best regards,<br>Marco<br></div></blockquote></div>