
<html><p>Hello everyone,</p><p>I am currently studying a supercell composed of 510 silicon atoms, including a divacancy. I am conducting a relaxation calculation using the PBE pseudopotential. Since silicon is tetravalent, each divacancy causes three neighboring silicon atoms to lose a bond. At a total charge of 0, I anticipate two distinct equilibrium structures.</p><ol><li>The three silicon atoms form two bonds, resulting in two short distances (~3 Å) and one longer distance (~3.5 Å) between them.</li><li>The three silicon atoms form only one bond, leading to one short distance and two longer distances between them.</li></ol><p>To achieve these two structures, I initialize the relaxation calculation by adjusting the distances between the three atoms to resemble the expected configurations. To expedite the calculation, I perform the relaxation in two stages: the first with force convergence (forc_conv) set to <span class="katex katex-mathml"><math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mn>1</mn><msup><mn>0</mn><mrow><mo>−</mo><mn>3</mn></mrow></msup></mrow><annotation encoding="application/x-tex">10^{-3}</annotation></semantics></math></span><span class="katex katex-html base mord" aria-hidden="true">10</span><span class="katex katex-html base mord msupsub vlist-t vlist-r vlist sizing reset-size6 size3 mtight" aria-hidden="true">−3</span> and a wavefunction cutoff energy (ecutwfc) of 10 Ry, and the second with forc_conv set to <span class="katex katex-mathml"><math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mn>1</mn><msup><mn>0</mn><mrow><mo>−</mo><mn>4</mn></mrow></msup></mrow><annotation encoding="application/x-tex">10^{-4}</annotation></semantics></math></span><span class="katex katex-html base mord" aria-hidden="true">10</span><span class="katex katex-html base mord msupsub vlist-t vlist-r vlist sizing reset-size6 size3 mtight" aria-hidden="true">−4</span> and ecutwfc of 25 Ry.</p><p>However, despite trying various approaches, I consistently observe three long distances of approximately 3.4 Å at the end of the second relaxation step, even though the structure appears correct after the first step.</p><p>Has anyone encountered a similar issue? Any suggestions on how to resolve this problem would be greatly appreciated.</p><p>Thank you.</p><p>Best regards,</p></html>