<div dir="ltr"><div class="gmail_default" style="font-family:tahoma,sans-serif;color:#000000"><div style="font-family:arial,helvetica,sans-serif;font-size:16px"><div style="background-color:rgb(250,250,250)">Hello all,<br><br>I have gone through all the archived discussion about cohesive energy calculation.<br><br>I am trying to calculate the cohesive energy of Ti which is <span style="color:rgb(34,34,34);font-family:arial,sans-serif;background-color:rgb(255,255,255)"> 4.85eV. But I am getting 5.23 eV which is closer. But I have found from a previous discussion here <span style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif">( </span><span id="gmail-m_8759704390713340179OBJ_PREFIX_DWT93_com_zimbra_url" style="color:rgb(51,102,153)"><a href="https://www.mail-archive.com/users@lists.quantum-espresso.org/msg11410.html" rel="nofollow noopener noreferrer" target="_blank" style="color:rgb(51,102,153);text-decoration-line:none;font-family:arial,helvetica,sans-serif">https://www.mail-archive.com/users@lists.quantum-espresso.org/msg11410.html ) </a></span></span><span style="font-size:12pt;color:rgb(34,34,34);font-family:arial,sans-serif">that, </span><span style="font-size:12pt">the final magnetization should be 2 Bohr mag/ cell at the end of calculation for Ti and smearing contribution of energy should be 0.0Ry. </span><br><div><div>But I am getting to get 3.83 Bohr mag/cell and a significant smearing contribution = -0.03295688 Ry.<br></div><div><br></div></div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;background-color:rgb(255,255,255)">Anyone help me to sort this. what should I change in the input?<br><br><br></span><div>here is my input,</div><br><div><div>&control</div><div>   calculation      = 'scf'</div><div>   restart_mode     = 'from_scratch'</div><div>   pseudo_dir       = '/work/shared/s-tih/pseudo/'</div><div>   prefix           = 'Tifree_deg1_50_8k_1'</div><div>   wf_collect       = .true.</div><div>   tstress          = .true.</div><div>   tprnfor          = .true.</div><div>   forc_conv_thr    = 1.0d-6</div><div>   verbosity        = 'high'</div><div> /</div><div>&system</div><div>   ibrav            = 1</div><div>   celldm(1)     = 20</div><div>   nat              = 1</div><div>   ntyp            = 1</div><div>   ecutwfc       = 50      //( I have performed convergence study. But since degauss is reduced to 0.01, i increased ecut) </div><div>   ecutrho       = 400</div><div>   nspin          = 2</div><div>   starting_magnetization(1) = 1        // I think, the final magnetisation should be 2 bohr mag/ cell at the end of calculation.</div><div>   nosym            = .true                    /// I hope this is enough to break the symmetry</div><div>   nbnd             = 100</div><div>   occupations      = 'smearing',  smearing='mp', degauss=0.01D0</div><div> /</div><div>&electrons</div><div>   diagonalization  = 'cg'</div><div>   mixing_beta      = 0.7d00</div><div>   conv_thr         = 1.0d-8     // I used even higher convergence since smearing is reduced.</div><div>/</div><div>ATOMIC_SPECIES</div><div>  Ti  47.8670  Ti.pbe-spn-rrkjus_psl.1.0.0.UPF</div><br><br><div>ATOMIC_POSITIONS (crystal)</div><div>   Ti     0.5 0.5 0.5</div><br><div>K_POINTS {automatic}</div><div>  8  8  5  0  0  0</div><br></div></div><div style="background-color:rgb(250,250,250)"><br>I don't really know this is the correct way to post a reply. All input is welcome. <br></div><div style="background-color:rgb(250,250,250)"><br></div></div></div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div style="font-size:12.8px"><span style="font-family:"Times New Roman";font-size:medium">Thanks & Regards,</span><br></div><div style="font-size:12.8px"><b style="font-family:"Times New Roman";font-size:medium">------------------------------------------------------------------------------------------------------</b></div><div style="line-height:22px"><div style="line-height:22px"><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:16px">Ajmalghan MUTHALI</span></font></div><div style="line-height:22px"><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:16px"><br></span></font></div><div style="line-height:22px"><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:16px">Post doctorate researcher</span></font></div><div style="line-height:22px"><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:16px">Laboratoire ICB </span></font></div><div style="line-height:22px"><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:16px">UMR 6303 CNRS-Université de Bourgogne </span></font></div><div style="line-height:22px"><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:16px">9 Avenue Alain Savary, BP 47870 </span></font></div><div style="line-height:22px"><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:16px">F-21078 DIJON Cedex, France </span></font></div><div style="line-height:22px"><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:16px">Tel: +33-(0)7.69.28.19.91 </span></font></div></div><div style="line-height:22px"><span style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif;font-size:16px">Email : </span><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:16px"><a href="mailto:ajmalghan.muthali@u-bourgogne.fr" target="_blank">ajmalghan.muthali@u-bourgogne.fr</a></span></font></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>