<div dir="ltr"><div>Hi,</div><div><br></div><div>Have you checked the literature for band gap values obtained from DFT with the PBE functional for your compound? Often there is a mismatch between experiment and simulation.</div><div><br></div><div>Here are som pointers how you can improve your input<br></div><div><br></div><div> - check and converge the cutoff: https://www.cp2k.org/howto:converging_cutoff</div><div> - increase the accuracy: EPS_SCF 1.0e-6 or better 1.0e-7, EPS_DEFAULT 1.0e-12<br></div><div> - make sure your coordinates are ok. The forces on the atoms are large, this system is far from a minimum.</div><div><br></div><div>Best,</div><div>Fabian<br></div><br>On Monday, 25 May 2020 15:08:57 UTC+2, Ref  wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hello everyone,<div><br></div><div>I tried to calculate the DOS of YBCO, but in the end, on the resulting graph, the formed value of energy gap  is tens of times larger than the expected result (obtained from experimental data - it should be about 20 meV). I used standard potential GTH, functional PBE and  normal parameters. I performed geometry optimization for a primitive cell, and then multiplied it to size 3 x 3 x 3 in a graphics editor with some doping. (about 190 atoms are in the final system.)</div><div>What could be associated with such a strong discrepancy? What general recommendations could you give to improve the result?<br></div><div><br></div><div><br><div><div><div style="overflow:auto"><div style="max-height:10000px"><div dir="ltr">I've attached the input file that I've used. <div><br></div><div>Thanks in advance for your assistance!</div></div></div></div></div><div><div></div></div></div></div></div></blockquote></div>