<html><head><meta http-equiv="Content-Type" content="text/html; charset=UTF-8" /></head><body style='font-size: 10pt; font-family: Verdana,Geneva,sans-serif'>
<p>Dear Balint,</p>
<p>Thank you for your help.</p>
<p>Unfortunately, I still have problems.</p>
<p>I actually did some calculations using either "ReadInitialCharges = Yes" or "ReadInitialCharges = No", but surprisingly I got the same overlap matrix (whereas the Hamiltonian in "hamsqr.dat" was different in the two cases).<br /><br /></p>
<p>Obviously, I still have the problem that C^T SC is not an identity matrix, since S and C are the same.</p>
<p>If it may help, I report below the input for H<sub>2</sub>O calculation:</p>
<p>Geometry = GenFormat {<br />3  C<br /> O H<br />     1    1    0.00000000000E+00  -0.10000000000E+01   0.00000000000E+00<br />     2    2    0.00000000000E+00   0.00000000000E+00   0.78306400000E+00<br />     3    2    0.00000000000E+00   0.00000000000E+00  -0.78306400000E+00<br />}<br /><br />Driver = ConjugateGradient {<br />  MovedAtoms = 1:3<br />  MaxForceComponent = 1.000000000000000E-008  # Extremely small!<br />  MaxSteps = 100<br />  OutputPrefix = "geom.out"<br />}<br /><br />Hamiltonian = DFTB {<br />  SCC = Yes<br />  SCCTolerance = 1.0E-010  # Extremely small!<br />  MaxSCCIterations = 1000<br />  MaxAngularMomentum = {<br />    O = "p"<br />    H = "s"<br />  }<br />  Charge = 0.000000000000000E+000<br />  SpinPolarisation = {}<br />  Eigensolver = Standard {}<br />  SlaterKosterFiles = {<br />    O-O = "O-O.skf"<br />    O-H = "O-H.skf"<br />    H-O = "O-H.skf"<br />    H-H = "H-H.skf"<br />  }<br />  ReadInitialCharges = Yes<br />}<br /><br />Options = {<br />  WriteAutotestTag = Yes<br />  WriteHS = Yes<br />}<br /><br />Analysis = {<br />WriteEigenvectors=Yes<br />EigenvectorsAsTxt=Yes<br />}<br /><br /><br />ParserOptions = {<br />  ParserVersion = 4<br />}</p>
<p><br /></p>
<p>Do you have any idea what am I doing wrong?</p>
<p>Best regards,</p>
<p>Alessandro</p>
<p><br /></p>
<p>Il 2017-10-17 08:25 Bálint Aradi ha scritto:</p>
<blockquote type="cite" style="padding: 0 0.4em; border-left: #1010ff 2px solid; margin: 0"><!-- html ignored --><!-- head ignored --><!-- meta ignored -->
<div class="pre" style="margin: 0; padding: 0; font-family: monospace">Dear Alessandro,<br /><br />
<blockquote type="cite" style="padding: 0 0.4em; border-left: #1010ff 2px solid; margin: 0">Let me clarify: if I have the following file "eigenvec.out" (H_2 O<br /> calculation):</blockquote>
<br />
<blockquote type="cite" style="padding: 0 0.4em; border-left: #1010ff 2px solid; margin: 0">Are the eigenvector listed by columns? I mean, is the eigenvector matrix<br /> as following?</blockquote>
<br /> Yes, the first column represents the components of the respective<br /> eigenvector.<br /><br />
<blockquote type="cite" style="padding: 0 0.4em; border-left: #1010ff 2px solid; margin: 0">I am puzzled about this point because I tried to compute C^T SC (using<br /> the overlap matrix in the file "oversqr.dat" as S) but I did not obtain<br /> an identity matrix, as I expected.</blockquote>
<br /> Please note, that if you switch on the writeHS option, there program<br /> just builds H and S and writes them and stops. There are no SCC-cycles<br /> made to reach convergence. Therefore, if you obtained your eigenvectors<br /> from an SCC-calculation, you should make sure, that your read in the<br /> charges.bin of that calculation (ReadInitialCharges = Yes), when asking<br /> DFTB+ to write H and S. Then, the overlap you obtain should be<br /> compatible with the eigenvectors and give the expected norm.<br /><br />   Best regards,<br /><br />   Bálint</div>
<br />
<div class="pre" style="margin: 0; padding: 0; font-family: monospace">_______________________________________________<br /> DFTB-Plus-User mailing list<br /><a href="mailto:DFTB-Plus-User@mailman.zfn.uni-bremen.de">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">https://mailman.zfn.uni-bremen.de/cgi-bin/mailman/listinfo/dftb-plus-user</a></div>
</blockquote>
<div>-- <br />
<div class="pre" style="margin: 0; padding: 0; font-family: monospace">
<blockquote type="cite" style="padding: 0 0.4em; border-left: #1010ff 2px solid; margin: 0"><span style="color: #808080; font-size: 10pt;">Alessandro Landi, Ph.D. Student</span></blockquote>
<blockquote type="cite" style="padding: 0 0.4em; border-left: #1010ff 2px solid; margin: 0"><span style="color: #808080; font-size: 10pt;">Università degli Studi di Salerno </span><br /><span style="color: #808080; font-size: 10pt;"> Dipartimento di Chimica e Biologia "Adolfo Zambelli"</span><br /><span style="color: #808080; font-size: 10pt;">Via Giovanni Paolo II, 132 - 84084 Fisciano (SA), Italy</span><br /><span style="color: #808080; font-size: 10pt;"> Phone number +39 089969390</span></blockquote>
</div>
</div>
</body></html>