<div dir="ltr">Hi there,<div><br></div><div>I'm using QE 6.5 to run a relaxation simulation with FeSi(NaCl structure) supercell. I put 8 atoms in the cell and I used simple cubic for crystal symmetry. Here's the input script I used:</div><div><br></div><div>&control<br>calculation='vc-relax',<br>prefix='FeSi_vc-relax_fixed'<br>pseudo_dir='/home/kchen/PSEUDO/'<br>/<br>&system<br>ibrav=1,<br>celldm(1)=9.032891,<br>nat=8,ntyp=2,<br>ecutwfc=100,<br>ecutrho=1500,<br>occupations='fixed'<br>/<br>&electrons<br>mixing_beta=0.3,<br>electron_maxstep=300<br>/<br>&ions<br>/<br>&cell<br>cell_dynamics='bfgs'<br>/<br>ATOMIC_SPECIES<br>Fe 55.845 Fe.pbe-spn-rrkjus_psl.1.0.0.UPF<br>Si 28.0855 Si.pbe-n-rrkjus_psl.1.0.0.UPF<br>ATOMIC_POSITIONS (crystal)<br>Fe 0  0  0<br>Si  0.5  0  0<br>Fe 0.5  0.5  0<br>Si  0  0.5  0<br>Fe 0.5  0  0.5<br>Si  0  0  0.5<br>Fe 0  0.5  0.5<br>Si  0.5  0.5  0.5<br>K_POINTS (automatic)<br>6 6 6 1 1 1<br></div><div><br></div><div>I got the final Fe-Si bond length as 2.33 A and Heat of Formation as 0.04eV/atom. However, in the reference, <a href="https://journals.aps.org/prb/pdf/10.1103/PhysRevB.85.245102">https://journals.aps.org/prb/pdf/10.1103/PhysRevB.85.245102</a>, page 6, Table VI, the bond length is around 2.37A, which is close to my result but the heat of formation is 0.25 eV/atom. I don't know what I'm doing wrong. For potential energy of Fe and Si. I used 1 Fe atom in bcc crystal and 2 Si atoms in fcc crystal for their potential energies. You can find the script I used in the attachment. The potential energy for bcc Fe and dia Si I got are -247.789Ry/atom and -11.4201Ry/atom respectively. </div><div><br></div><div>I also tried SCF calculation with the bond length provided in the reference, but it did not help a lot.</div><div><br></div><div>Thank you for your help</div><div><br></div><div>Kefan</div><div><br></div><div>Ph.D. student</div><div>Department of Chemical and Materials Engineering</div><div>University of Nevada, Reno, USA</div><div><br></div></div>