<div dir="ltr">Dear Developers,<div><br></div><div>Would it be possible to you to explain why there is a need for the loop listed below while calculating nonscc derivatives? </div><div><br></div><div> !! Add contribution to the force from atom 1 onto atom 2f using the</div><div>          !! symmetry in the blocks, and also the skew symmetry in the</div><div>          !! derivatives</div><div>          do ii = 1, 3</div><div>            deriv(ii,iAtom2f) = deriv(ii,iAtom2f) &</div><div>                &- sum(sqrDMTmp(1:nOrb2,1:nOrb1) &</div><div>                &* 2.0_dp*hPrimeTmp(1:nOrb2,1:nOrb1,ii)) &</div><div>                &+ sum(sqrEDMTmp(1:nOrb2,1:nOrb1)&</div><div>                &* 2.0_dp*sPrimeTmp(1:nOrb2,1:nOrb1,ii))</div><div>          end do</div><div><br></div><div>I read the comments above the loop. But, it is not clear to me from these comments that why the first loop in the derivative_nonSCC routine, i.e.,</div><div><br></div><div><div>          do ii = 1, 3</div><div>            deriv(ii,iAtom1) = deriv(ii,iAtom1) &</div><div>                &+ sum(sqrDMTmp(1:nOrb2,1:nOrb1)&</div><div>                &* 2.0_dp*hPrimeTmp(1:nOrb2,1:nOrb1,ii)) &</div><div>                &- sum(sqrEDMTmp(1:nOrb2,1:nOrb1)&</div><div>                &* 2.0_dp*sPrimeTmp(1:nOrb2,1:nOrb1,ii))</div><div>          end do</div></div><div><br></div><div>is not sufficient to describe the force on iAtom1? </div><div><br></div><div>Based on my understanding, the force on each atom should be related to the summation (of rho_uv * H_uv and rhoE_uv * S_uv) on all the orbitals related to the atom of interest and the orbitals corresponding to all the other neighbouring atoms and I could clearly see that this information is already considered in the first loop.</div><div><br></div><div>Kindly, help me to understand.</div><div><br></div><div>Thanks for your help,</div><div>Sincerely,<br>Sharma.</div><div><br></div><div>P. S.: May I also know, why derivative_nonSCC is separately made public despite its availability through the derivative_shift interface?</div><div><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><span style="color:rgb(0,0,0)"><br></span></div><div dir="ltr"><span style="color:rgb(0,0,0)">-------------------------------------------------------------------</span><br></div><div dir="ltr"><font color="#000000"><font size="2" face="arial, helvetica, sans-serif">Sharma</font></font></div><div dir="ltr"><font color="#000000" face="arial, helvetica, sans-serif"><a href="http://www.chem.rochester.edu/groups/huo/people/" target="_blank">http://www.chem.rochester.edu/groups/huo/people/</a></font></div><div dir="ltr"><br></div></div></div></div></div></div></div></div></div>
</div></div>