<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Arial,Helvetica,sans-serif;" dir="ltr">
<p>Hello,</p>
<p><br>
</p>
<p>I am currently trying to run a basic band structure calculation for a periodic silicon 2-atom PUC with spin orbit coupling. My dftb_in.hsd file reads as:</p>
<p></p>
<div><br>
Geometry = GenFormat {<br>
    2  F<br>
 Si<br>
    1 1    0 0 0<br>
    2 1    0.25 0.25 0.25<br>
    0 0 0<br>
    2.715               2.715                0<br>
    0                  2.715                2.715<br>
    2.715                0                   2.715<br>
}<br>
<br>
Driver = ConjugateGradient {<br>
  MovedAtoms = Si<br>
  MaxForceComponent = 1E-4<br>
  MaxSteps = 100<br>
  OutputPrefix = "geom.out"<br>
}<br>
<br>
Hamiltonian = DFTB {<br>
  SCC = Yes<br>
  SCCTolerance = 1e-5<br>
  SlaterKosterFiles {<br>
    Si-Si = "Si-Si.skf"<br>
  }<br>
  MaxAngularMomentum {<br>
    Si = "d"<br>
  }<br>
SpinOrbit = {<br>
Si [eV] = {0 0.044 0}<br>
}<br>
KPointsAndWeights = {  # 2x2x2 MP-scheme<br>
0.25 0.25 0.25   1.0<br>
0.25 0.25 -0.25   1.0<br>
0.25 -0.25 0.25   1.0<br>
0.25 -0.25 -0.25   1.0<br>
}<br>
}<br>
<br>
Regardless of the spin coupling constants I enter (with the exception of non-zero values for the s orbital), I obtain the following error after 10 geometry steps:<br>
<br>
<span>-> Failure in diagonalisation routine ZHEEVR, unable to complete Cholesky factorization of B     46</span><br>
<br>
Any ideas on what I can do to fix this? For reference, I am using the si-band parameter set.<br>
Changing to the pbc parameter set, geometry convergence DOES occur. However, after sampling along the k-lines L - G - X and using dp_bands, I obtain the error:<br>
<br>
<div>  File "/apps/dftbp/1.3.1/lib/python2.7/site-packages/dptools/bandout.py", line 69, in fromfile<br>
    nkpt = len(eigvalarrays) / nspin<br>
ZeroDivisionError: integer division or modulo by zero<br>
</div>
<br>
I assume this means I have to set nspin /= 0 somewhere, but I cannot figure out how to do so.<br>
<br>
Kind regards,<br>
Lachlan Oberg<br>
<br>
</div>
<p></p>
</div>
</body>
</html>