Hi everyone, <div><br></div><div>I know that DOS and points aren't what people first think of when using cp2k, but the functionality does seem to be there in some limited way... </div><div><br></div><div>Anyway, I'm wondering how to properly use it, because my test case doesn't look right. I'm running a point calculation on GaAs with a 16x16x16 Monkhorst-Pack pack grid and then printing the total DOS. I then take the columns from the DOS file (head attached) and plot energy (column 1) against occupation (column 3) using the reported Fermi energy from the output file of ~0.14 Ha. This results in a DOS where the Fermi energy does not line up with the VBM. What's more, if I visually inspect this DOS, the band gap of GaAs is >1eV. While that is closer to the experimental value, the reported band gap using GGA PBE is closer to 0.3eV.  </div><div><br></div><div>Questions:</div><div>(1) What is the distinction between columns 2 and 3 in the DOS file? They give similar, but not equivalent looking dos plots.</div><div>(2) Is the Fermi energy printed in the output file the same as the Fermi energy of the DOS, or is it defined different?</div><div><br></div><div>I've attached my input file for reference. Thanks in advance.</div><div><p> # DOS at iteration step i = 0</p>
<p>    Energy[a.u.]       Density     Occupation</p>
<p>    -0.40682067         0.0009         0.0957</p>
<p>    -0.40582067         0.0015         0.1543</p>
<p>    -0.40482067         0.0000         0.0000</p>
<p>    -0.40382067         0.0000         0.0000</p>
<p>    -0.40282067         0.0034         0.3516</p>
<p>    -0.40182067         0.0026         0.2656</p>
<p>    -0.40082067         0.0038         0.3828</p>
<p>    -0.39982067         0.0000         0.0000</p></div><div><br></div><div><img alt="Screen Shot 2022-06-01 at 10.32.10 AM.png" data-iml="275456" width="534px" height="350px" src="cid:a178e7d2-0b9d-4fe8-b8a0-118454057bf0"><br></div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/3c7ccd12-f111-40ac-8c31-b950d71ee41bn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/3c7ccd12-f111-40ac-8c31-b950d71ee41bn%40googlegroups.com</a>.<br />