<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=utf-8"><meta name=Generator content="Microsoft Word 14 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
        {font-family:"Segoe UI";
        panose-1:2 11 5 2 4 2 4 2 2 3;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
span.StileMessaggioDiPostaElettronica17
        {mso-style-type:personal-reply;
        font-family:"Calibri","sans-serif";
        color:#1F497D;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-family:"Calibri","sans-serif";
        mso-fareast-language:EN-US;}
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:70.85pt 2.0cm 2.0cm 2.0cm;}
div.WordSection1
        {page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=IT link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Dear Roozbeh,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>the internal order of the p orbitals is not |px>, |py>, |pz> but |py> |pz> |px>. <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>This follows more naturally from the  order of angular momentum eigenstates |-1>, |0>, |1>  transformed into Tesseral wavefunctions. I think Ben can comment more on this, probably.<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Hence, the SK rotations agree with 1954 paper. <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>For hh(1,1) use ll**2+mm**2+nn**2=1, <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>hh(1,1) = Eyy = mm**2*sk(1) + (1 - mm**2)*sk(2)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>hh(2,2) = Ezz  = nn**2*sk(1) + (1 - nn**2)*sk(2)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>hh(3,3) = Exx  = ll**2*sk(1) + (1 - ll**2)*sk(2)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>hh(2,1) = Eyz = mm*nn* (sk(1)-sk(2))<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>hh(3,1) = Exy = mm*ll* (sk(1)-sk(2))<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>hh(3,2) = Exz = nn*ll* (sk(1)-sk(2))<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Also d-orbitals follow a similar ordering starting from |-2>, |-1>, |0>, |1>, |2>, <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>hence |dxy>, |dyz>, |dz^2-r^2>, |dxz>, |dx^2-y^2>    <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Alessandro<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Segoe UI","sans-serif"'>Da:</span></b><span style='font-size:10.0pt;font-family:"Segoe UI","sans-serif"'> DFTB-Plus-User [mailto:dftb-plus-user-bounces@mailman.zfn.uni-bremen.de] <b>Per conto di </b>Roozbeh Anvari<br><b>Inviato:</b> domenica 29 settembre 2019 02:03<br><b>A:</b> User list for DFTB+ related questions<br><b>Oggetto:</b> Re: [DFTB-Plus-User] Format of the Slater-Koster files<o:p></o:p></span></p><p class=MsoNormal><o:p> </o:p></p><div><p class=MsoNormal>   Hi,   <o:p></o:p></p><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>   I was wondering about how DFTB+ rotates H and S,  for example for a p-p interaction : ( reference:: DFTB+%SK.F90%subroutine pp(hh, ll, mm, nn, sk) ) : <o:p></o:p></p></div><div><p class=MsoNormal>    "<o:p></o:p></p></div><p class=MsoNormal>    hh(1,1) = (1.0_dp-nn**2-ll**2)*sk(1)+(nn**2+ll**2)*sk(2)<br>    hh(2,1) = nn*mm*sk(1)-nn*mm*sk(2)<br>    hh(3,1) = ll*mm*sk(1)-ll*mm*sk(2)<br>    hh(2,2) = nn**2*sk(1)+(1.0_dp-nn**2)*sk(2)<br>    hh(3,2) = nn*ll*sk(1)-nn*ll*sk(2)<br>    hh(3,3) = ll**2*sk(1)+(1.0_dp-ll**2)*sk(2)<o:p></o:p></p><div><p class=MsoNormal>    "<o:p></o:p></p></div><div><p class=MsoNormal>   while, <span style='font-size:10.0pt;font-family:"Arial","sans-serif"'>(Slater, J. C., & Koster, G. F. (1954) . <i>Physical Review</i>, <i>94</i>(6), 1498.) : </span> <o:p></o:p></p></div><div><p class=MsoNormal>      Exx = ll**2 *(pp_sigma) - (1-ll**2) *(pp_pi)   <o:p></o:p></p></div><div><p class=MsoNormal>      Exy = ll*mm *(pp_sigma) - ll*mm *(pp_pi)<o:p></o:p></p></div><div><p class=MsoNormal>      Exz = ll*nn *(pp_sigma) - ll*nn *(pp_pi)<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>here, could you please confirm the relative directions of the directional cosines and also sampled energy integrals in the .SKf files ?  i.e. ( (ll , mm , nn) == (dx, dy , dz ) &  (hh(1,1) , hh(1,2),  hh(1,3) == (Ex,x,Ex,y Ex,z)  ? ) <o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>2.<o:p></o:p></p></div><div><p class=MsoNormal>how accurate would be the calculated band structure without taking into account the potential shift  in H = H0 + shift ?  the E-K I have calculated using just H0/S from .SKF looks very similar to <o:p></o:p></p></div><div><p class=MsoNormal>    that given by dftb+ and dft, but band widths are significantly narrower, <o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Best regards <o:p></o:p></p></div><div><p class=MsoNormal>Roozbeh<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div></div><p class=MsoNormal><o:p> </o:p></p><div><div><p class=MsoNormal>On Thu, Sep 19, 2019 at 2:02 PM Roozbeh Anvari <<a href="mailto:roozbeh.anvari@gmail.com">roozbeh.anvari@gmail.com</a>> wrote:<o:p></o:p></p></div><blockquote style='border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm'><div><p class=MsoNormal>Dear Balint, <o:p></o:p></p><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>thank you for your help,<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Best regards,<o:p></o:p></p></div><div><p class=MsoNormal>Roozbeh<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div></div></div><p class=MsoNormal><o:p> </o:p></p><div><div><p class=MsoNormal>On Thu, Sep 19, 2019 at 5:08 AM Bálint Aradi <<a href="mailto:aradi@uni-bremen.de" target="_blank">aradi@uni-bremen.de</a>> wrote:<o:p></o:p></p></div><blockquote style='border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm'><p class=MsoNormal>Dear Roozbeh,<br><br>> I wonder why the overlap values are so large, are not the RAW values<br>> normalised ? <br><br>The first few lines are only placeholders (historical heritage...).<br>Aftert that, the overlap integrals are less than one. For homonuclear<br>cases, identical orbitals (e.g. H-H  Sss0) give values close to 1 (makes<br>sense), while different orbitals or heteronuclear cases yield values one<br>order of magnitude lower.<br><br>> To check this, I am searching where skHamCont and skOverCont  are read<br>> from  .SKF but have not found it yet, <br><br>The reading of skf and construction of the containers happens in<br>multiple steps at multiple places:<br><br><a href="https://github.com/dftbplus/dftbplus/blob/master/prog/dftb+/lib_type/oldskdata.F90" target="_blank">https://github.com/dftbplus/dftbplus/blob/master/prog/dftb+/lib_type/oldskdata.F90</a><br><br><a href="https://github.com/dftbplus/dftbplus/blob/1a31db05a60ad37ce9013b66dca4388a94b72817/prog/dftb%2B/lib_dftbplus/parser.F90#L2526" target="_blank">https://github.com/dftbplus/dftbplus/blob/1a31db05a60ad37ce9013b66dca4388a94b72817/prog/dftb%2B/lib_dftbplus/parser.F90#L2526</a><br><br><a href="https://github.com/dftbplus/dftbplus/blob/master/prog/dftb+/lib_dftb/slakocont.F90" target="_blank">https://github.com/dftbplus/dftbplus/blob/master/prog/dftb+/lib_dftb/slakocont.F90</a><br><br><br>I hope this helps.<br><br>  Best regards,<br><br>  Bálint<br><br><br>-- <br>Dr. Bálint Aradi<br>Bremen Center for Computational Materials Science, University of Bremen<br><a href="http://www.bccms.uni-bremen.de/cms/people/b-aradi/" target="_blank">http://www.bccms.uni-bremen.de/cms/people/b-aradi/</a><br><br><br>_______________________________________________<br>DFTB-Plus-User mailing list<br><a href="mailto:DFTB-Plus-User@mailman.zfn.uni-bremen.de" target="_blank">DFTB-Plus-User@mailman.zfn.uni-bremen.de</a><br><a href="https://mailman.zfn.uni-bremen.de/cgi-bin/mailman/listinfo/dftb-plus-user" target="_blank">https://mailman.zfn.uni-bremen.de/cgi-bin/mailman/listinfo/dftb-plus-user</a><o:p></o:p></p></blockquote></div></blockquote></div></div></body></html>