<p>Dear QE users:</p><p>I try to use QE-6.4.1 to simulate <img id="x_x_x__x0000_i1025" style="width: 60.74pt; height: 18pt; cursor: pointer;" src="cid:2609D217@08AB9A01.018E945D" crossorigin="use-credentials" data-outlook-trace="F:3|T:3" data-custom="AAMkADYyNGE5ZDVlLWI0M2ItNGNkYi04NWUzLTgwMDBiYmY4MWE1YgBGAAAAAADT1egYenlfTqNtjOrIuhcjBwBEgTiP%2FR2ETam0o5puADWEAAAAAAEJAABEgTiP%2FR2ETam0o5puADWEAALgunFEAAABEgAQAIdmOwWosi9FpfzF8EHcwG8%3D" naturalheight="24" naturalwidth="81" data-imagetype="AttachmentByCid" originalsrc="cid:image002.png@01D55F54.3DEAB060"><font face="Calibri"> as single molecule. The system is built in large super cell, and <em>total_charge=-4</em>  is added to input file leading to charged system. When I do the convergence test in term of the super cell size (20 Angs - 30 Angs ), the variation for total energy is small: -1038.6784 Ry for 20 Angs and 1038.7388 Ry for 30 Angs. But the band energy and Fermi energy (5.23 ev to 3.08 ev) shift visibly. The problem is that the band gap between band 104 (occupied) and 105 (unoccupied) become smaller with the increase of size (about 0.81 ev change to 0.11 ev). It seems the band gap is gradually disappear, which confuses me so much. The input file is attached with this email. It will be very appreciated if you can solve this for me. Thank you.</font></p><p><font face="Calibri"><br></font></p><p><font face="Calibri">Best wishes</font></p><p><font face="Calibri">Yang Zhou</font></p><p><font face="Calibri"><br></font></p><p><font face="Calibri">------</font></p><p><font face="Calibri">Yang Zhou</font></p><p><font face="Calibri">PhD student</font></p><p><font face="Calibri">University of Leeds</font></p><p><font face="Calibri">United Kingdom</font></p><div></div>