Hi, Flo<div><br></div><div>  Besides the issue on the sign of determinant, I have another concern on the et_coupling subroutine. As shown in the attached file, the one-electron operator coupling between two Slater determinants is the summation of the corresponding couplings between the orbitals weighted by the inverse of overlap matrix. The derivation in details is given in "Methods of Molecular Quantum Mechanics, Second Edition. (R. McWeeny), Page 61-66".</div>
<div> </div><div>  Apparently, there is an index swap between the orbital coupling and overlap matrix. Therefore, in "et_coupling.F",</div><div><br></div><div>                b(i)=b(i)+rest_MO(2)%matrix%local_data(k,j)*inverse_mat%local_data(k,j)*tttt</div>
<div><br></div><div>should be changed to</div><div> </div><div>                 b(i)=b(i)+rest_MO(2)%matrix%local_data(k,j)*inverse_mat%local_data(j,k).</div><div><br></div><div>Similar modification ought to be applied on a(i) too. </div>
<div>                  </div><div>Could you please check this issue too ?</div><div><br></div><div>Thanks.</div><div><br></div><div>Hanning</div><div><br><br><div class="gmail_quote">On Wed, Nov 24, 2010 at 11:08 AM, flo <span dir="ltr"><<a href="mailto:fsch...@pci.uzh.ch">fsch...@pci.uzh.ch</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Hi Hanning,<br>
<br>
I agree, that it should matter there. But indeed I wrote something<br>
wrong in my first reply. It should give the correct sign of the<br>
determinant. As your test indicates it does not, that is definitely a<br>
bug. I will have a look into this as soon as I got a little time.<br>
<br>
Thank you<br>
<div><div></div><div class="h5">Flo<br>
<br>
--<br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>.<br>
To unsubscribe from this group, send email to <a href="mailto:cp2k%2Bun...@googlegroups.com">cp2k+uns...@googlegroups.com</a>.<br>
For more options, visit this group at <a href="http://groups.google.com/group/cp2k?hl=en" target="_blank">http://groups.google.com/group/cp2k?hl=en</a>.<br>
<br>
</div></div></blockquote></div><br></div>