<div dir="ltr"><div class="" style="font-size: 12px; line-height: normal; color: rgb(69, 69, 69);"><span class="" style="line-height: normal;"><font face="Helvetica Neue">Dear</font></span><span class="" style="font-family: ".PingFang SC"; line-height: normal;">,</span><font face="Helvetica Neue"> all</font></div><div class="" style="font-size: 12px; line-height: normal; font-family: "Helvetica Neue"; color: rgb(69, 69, 69);"><br></div><div class="" style="font-size: 12px; line-height: normal; font-family: "Helvetica Neue"; color: rgb(69, 69, 69);">I did some calculations with DFT+U functional recently, they all obtained from the CP2K program. <br></div><div class="" style="font-size: 12px; line-height: normal; font-family: "Helvetica Neue"; color: rgb(69, 69, 69);"><br></div><div class="" style="font-size: 12px; line-height: normal; font-family: "Helvetica Neue"; color: rgb(69, 69, 69);">At first, I calculated the geometry optimization of the structre, and obtained the total energy directly from the last optimized structure (let's say E1). Then<span class="" style="line-height: normal; font-family: ".PingFang SC";">,</span> I ran a single point energy calculation with the same optimized structure and got the total energy E2. </div><div class="" style="font-size: 12px; line-height: normal; font-family: "Helvetica Neue"; color: rgb(69, 69, 69);"><br></div><div class="" style="font-size: 12px; line-height: normal; font-family: "Helvetica Neue"; color: rgb(69, 69, 69);">However, those two total energies (E1 & E2) I collected from the exact same structure have significant difference<span class="" style="line-height: normal; font-family: ".PingFang SC";">,</span> E2-E1 = 1.2 eV.  </div><div class="" style="font-size: 12px; line-height: normal; font-family: "Helvetica Neue"; color: rgb(69, 69, 69);"><br></div><div class="" style="font-size: 12px; line-height: normal; font-family: "Helvetica Neue"; color: rgb(69, 69, 69);">When checking their output results<span class="" style="line-height: normal; font-family: ".PingFang SC";">,</span>  I find that the total energy difference could be attributed to two parts, DFT+U energy and  Core Hamiltonian energy. The input and output files are attached below. </div><div class="" style="font-size: 12px; line-height: normal; font-family: "Helvetica Neue"; color: rgb(69, 69, 69);"><br></div><div class="" style="font-size: 12px; line-height: normal; color: rgb(69, 69, 69);"><span class="" style="line-height: normal;"><font face="Helvetica Neue">As far as I know</font></span><span class="" style="font-family: ".PingFang SC"; line-height: normal;">,</span><font face="Helvetica Neue"> the total energy calculated from the same settings and structures should be the same. I have no idea about these results</font><font face="Helvetica Neue">. Could you please help me check what is wrong with my calculation? Any comments will be greatly appreciated. </font></div><div class="" style="font-size: 12px; line-height: normal; color: rgb(69, 69, 69);"><font face="Helvetica Neue"><br></font></div><div class="" style="font-size: 12px; line-height: normal; color: rgb(69, 69, 69);"><font face="Helvetica Neue">Best Regards,</font></div><div class="" style="font-size: 12px; line-height: normal; color: rgb(69, 69, 69);"><font face="Helvetica Neue"><br></font></div><div class="" style="font-size: 12px; line-height: normal; color: rgb(69, 69, 69);"><font face="Helvetica Neue">Xiaohui</font></div><div class="" style="font-size: 12px; line-height: normal; color: rgb(69, 69, 69);"><font face="Helvetica Neue"><br></font></div><div class="" style="font-size: 12px; line-height: normal; color: rgb(69, 69, 69);"><font face="Helvetica Neue"><br></font></div><div class="" style="font-size: 12px; line-height: normal; color: rgb(69, 69, 69);"><font face="Helvetica Neue"><br></font></div></div>