<div dir="ltr">Hi,<div><br></div><div>First of all, a big thanks to Marcella for fixing this in rev. 14153.</div><div><br></div><div>We just noticed that there still seems to be a bug in the title line of the cube files.<br></div><div>E.g. for a system with 56 electrons (see attachment), we find</div><div><br></div><div><div>head  -n2 PROJ-WFN_00028_1-1_0.cube</div><div>-Quickstep-</div><div>  WAVEFUNCTION           28 spin            1 i.e. HOMO -            0</div></div><div><br></div><div><br></div><div><div>head  -n2 PROJ-WFN_00029_1-1_0.cube</div><div>-Quickstep-</div><div>  WAVEFUNCTION           29 spin            1 i.e. LUMO +           56</div></div><div><br></div><div><br></div><div>I.e. the number counting for the empty states is incorrect.</div><div><br></div><div>The problem probably lies around line 1394 of qs_scf_post_gpw.F</div><div><br></div><div>  WRITE(title,*) "WAVEFUNCTION ",index_mo," spin ",ispin," i.e. LUMO + ",ifirst+ivector-2<br></div><div><br></div><div>There are quite a few index variables in this function (ifirst, ivector, index_mo, homo and lumo) and I am not quite sure where the mistake lies, but I guess for the author it will be easy to spot.</div><div><br></div><div>Please also confirm that (as I suspect) the error only concerns the "LUMO + .." index and that the index_mo is printed correctly.</div><div><br></div><div>Best,</div><div>Leopold</div><div><blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><div><div><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><p></p></blockquote></div></div></div></blockquote></div></div>