For any delta_e do you get all zeros?<div>And if delta_e = 0?</div><div><br></div><div>About the settings, they look fine to me. </div><div>Why are you using a Fermi-Dirac smearing for GaAs, it is not a metal. </div><div>Anyway, this should not change anything with respect to the energy gap.</div><div>Have you considered to run the same calculation with multiple cells instead of k-points? <br><br></div><div>Regards</div><div>Marcella</div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, June 8, 2022 at 10:49:31 PM UTC+2 nwi...@berkeley.edu wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Marcella, I looked further and played with delta_e and there is something that seems incorrect to me.<div><br></div><div>If I let delta_e -> 0 then the density column is equal to 0 everywhere. It looks like this might be due to precision? The column only lists 3 decimal places.<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, June 6, 2022 at 1:02:40 PM UTC-7 Nicholas Winner wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Thank you Marcella. I am clear about the Fermi level setting now. <br><div><br></div><div>Any ideas as to why calculated gap is so large (first image)? The input file I have attached has pretty reasonable parameters I think:</div><div><ul><li>multigrain can accommodate all basis functions</li><li>precision settings are above average </li><li>TZVP basis sets</li><li>Tight MP k-point grid of 16x16x16</li><li>EPS_SCF of 1e-6</li></ul><div>Not sure what else could be causing such a large discrepancy without published PBE results. Is it still an issue of interpreting the DOS?</div></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, June 3, 2022 at 2:38:05 AM UTC-7 Marcella Iannuzzi wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><br></div><div><br></div>Hi, <div><br></div><div>The Fermi energy is not well defined for electronic structures with an energy gap. </div><div>Using the Fermi-Dirac distribution, CP2K sets the Fermi energy somewhere in the middle of the gap.</div><div><br></div><div><br></div><div>The histogram calculated by activating PRINT%DOS is just counting the number of eigenvalues falling into each histogram bit.</div><div>The partitioning of the histogram in bits is determined from the DELTA_E keyword, the default of this is 0.001.</div><div>The final value of the density (sum of states per bit) is then normalised by the number of NMO. </div><div>The column occupation reports the sum of the occupation numbers of the eigenstates with energy falling into the corresponding bit.</div><div>For small DELTA_E, there are bits with zero density and zero occupation number. </div><div>Setting DELTA_E to zero, the histogram is replaced by simply the list of eigenvalues. </div><div><br></div><div>When using k-point sampling the procedure is more or less the same, but each contribution to the density as well as each contribution to the occupation are multiplied by the weight of the corresponding k-point. </div><div><br></div><div>Regards</div><div>Marcella</div><div><br></div><div><br></div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Thursday, June 2, 2022 at 10:05:09 PM UTC+2 <a rel="nofollow">nwi...@berkeley.edu</a> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div>It is a static calculation, so iteration = 0. I redid the calculation with cell optimization. Didn't seem much of change. Here at the plots using the occupation first and the density second. I set fermi level to be 0, and it doesn't put it anywhere near the VBM like we expect. In the second case it looks like there is a gap, but its too big. <br></div><div><br></div><div>Any knowledge maybe about the definition of fermi level? Maybe that's where I'm tripping up.</div><img alt="Screen Shot 2022-06-02 at 1.00.02 PM.png" width="860px" height="578px" src="https://groups.google.com/group/cp2k/attach/409af07ceb3d2/Screen%20Shot%202022-06-02%20at%201.00.02%20PM.png?part=0.1&view=1"><div>Do you know anything about the fermi level definition?</div><div><img alt="Screen Shot 2022-06-02 at 12.59.51 PM.png" width="866px" height="578px" src="https://groups.google.com/group/cp2k/attach/409af07ceb3d2/Screen%20Shot%202022-06-02%20at%2012.59.51%20PM.png?part=0.2&view=1"><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Thursday, June 2, 2022 at 9:04:42 AM UTC-7 <a rel="nofollow">wave...@gmail.com</a> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi<div><br></div><div>It's not clear to me why you are looking at the first step. Did it converge in just one step?</div><div><br></div><div>At any rate, you'll notice a difference between &DOS and &PDOS. PDOS returns values per orbital (s, p, d, f). If I understand it correctly, the density is the DOS as a solid state physicist would understand it, and occupation gives you a fractional filling at that energy. Not quite the same thing.</div><div><br></div><div>Regards</div><div>Sam<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, June 1, 2022 at 7:39:20 PM UTC+2 <a rel="nofollow">nwi...@berkeley.edu</a> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">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" width="534px" height="350px" src="https://groups.google.com/group/cp2k/attach/3b325621f9218/Screen%20Shot%202022-06-01%20at%2010.32.10%20AM.png?part=0.1&view=1"><br></div></blockquote></div></blockquote></div></blockquote></div></blockquote></div></blockquote></div></blockquote></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/b029f99c-f044-498e-a763-39f92341fec0n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/b029f99c-f044-498e-a763-39f92341fec0n%40googlegroups.com</a>.<br />