<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Dear Singh,<div><br></div><div> i) in the hexagonal cell the forces vanish by symmetry, that's why they are perfectly zero. In the rectangular cell this is no more the case. I think the problem arises because you did not specify the celldm(1) parameter with enough significant digits, perhaps, hence the structure does not correspond exactly to the optimal one.</div><div> Then keep in mind that the cell is rectangular, thus in principle you should use a different number of k-points in the two directions.</div><div><br></div><div> ii) I do not see the point of moving to larger supercells if you are in trouble with the small one. No reason to try such a high wfc and density cutoff either, in my opinion.</div><div><br></div><div> iii) In the hexagonal cell the forces look good to me, for practical purposes they are sufficiently small to say that you are at the theoretical energy minimum. The energy is much closer to the one of the primitive cell because the k-point sampling is correct this time.</div><div><br></div><div> About the 72 atom cell, well, first are you sure about the structure? If so, you need to specify your input file, otherwise if will be difficult to judge. </div><div><br></div><div><br></div><div>HTH</div><div><br></div><div><br></div><div>GS</div><div><br><div><div>Il giorno 22/nov/2010, alle ore 04.08, DHRUV SINGH ha scritto:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Dear PWSCF users,<br><br>I use a hexagonal unit cell with 2 basis atoms and lattice parameter (C-C bond length) of 2.6868 a.u. - this value was arrived at by energy minimization over a range of lattice parameters near this value. With a 16x16x1 kpoint grid it is seen that the energy per atom converges to -11.395192 Ry and the forces on the 2 atoms are zero (atleast until a good precision). The ultrasoft pseudopotential with the RRKJ modification is used for C and the results obtained are in excellent agreement with Marzari and Mounet PRB 71, 205214, 2005. <br>The essentials of the input file and the forces from the 2 atom unit cell are shown below<br>--------------------------------------------------------------------------------<br> &system<br> ibrav=0,<br> celldm(1)=2.6868,<br> nat=2,<br> ntyp=1,<br> ecutwfc=40.0,<br> ecutrho=480.0,<br> degauss=0.03,<br> nbnd=8<br> /<br> &electrons<br> mixing_beta = 0.5<br> conv_thr = 1.0d-8<br> /<br>ATOMIC_SPECIES<br>C 12.0 C.pbe-rrkjus.UPF<br><br>CELL_PARAMETERS {hexagonal}<br> 1.5 0.866025404 0.0<br> 1.5 -0.866025404 0.0<br> 0.0 0.0 20<br>ATOMIC_POSITIONS {crystal}<br> C 0.0000000 0.0000000 0.0000000<br> C 1/3.0 1/3.0 0.0000000<br><br>atom 1 type 1 force = -0.00000000 -0.00000000 0.00000000<br>atom 2 type 1 force = 0.00000000 0.00000000 0.00000000<br>--------------------------------------------------------------------------------<br><br>Using this lattice parameter I then try to do supercell calculations with rectangular and hexagonal supercells and I am facing the following problems in the supercell scf calculations :<br><br>i) With a 4 atom rectangular supercell and the same lattice parameter of 2.6868 au, atomic coordinates and input shown below, I find that the energy/atom coverges to -11.395182 Ry(very close to the 2 atom unit cell) for a 8x8x1 k point grid but there seem to be non-negligible net forces on all 4 atoms as shown below, <br>--------------------------------------------------------------------------------<br><br>&system<br> ibrav=0,<br> celldm(1)=2.6868,<br> nat=4,<br> ntyp=1,<br> ecutwfc=40.0,<br> ecutrho=480.0,<br> degauss=0.03,<br> nbnd=16<br> /<br> &electrons<br> mixing_beta = 0.5<br> conv_thr = 1.0d-8<br> /<br><br>CELL_PARAMETERS {cubic}<br> 3.0 0.0 0.0<br> 0.0 1.732050808 0.0<br> 0.0 0.0 14<br>ATOMIC_POSITIONS {alat}<br>C 0 0 0<br>C 1 0 0<br>C 1.5 0.866025404 0<br>C 2.5 0.866025404 0<br>--------------------------------------------------------------------------------<br><br> atom 1 type 1 force = 0.00310476 0.00000000 0.00000000<br> atom 2 type 1 force = -0.00310476 0.00000000 0.00000000<br> atom 3 type 1 force = 0.00310476 0.00000000 0.00000000<br> atom 4 type 1 force = -0.00310476 0.00000000 0.00000000<br><br>ii) With a 48 atom rectangular supercell and 2x2x1 kpoint grid I find energy/atom of -11.395152 Ry (again very close to the 2 atom unit cell) but non negligible forces: Forces on atom 1 are shown below<br>atom 1 type 1 force = 0.00274639 0.00000000 0.00000000<br>The story is the same as I keep increasing the supercell size ... for 60 atoms the forces on atom 1 are<br>atom 1 type 1 force = -0.00395744 0.00000000 0.00000000<br>The forces do not change even if I increase ecutwfc to 80.0 and ecutrho to 800.0. They also do not change if I increase the kpoint grid resolution<br><br>iii) Hexagonal supercells: I get slightly better answers with a hexagonal supercell in terms of forces. For example a 32 atom hexagonal supercell of graphene gives me an energy/atom of 11.395194 Ry for a 4x4x1 kpoint grid. The force on 1 of the atoms is <br>atom 2 type 1 force = 0.00001419 0.00000887 0.00000000<br>and the magnitude is similar for most other atoms. Once again I find that increasing ecutrho or ecutwfc doesnt matter. <br><br>However a more serious problem is seen with the hexagonal supercell. I find that scf calculations do not converge for a 72 atom supercell, with the following trend: <br><br>estimated scf accuracy < 7.42866516 Ry<br> estimated scf accuracy < 2.47702692 Ry<br> estimated scf accuracy < 0.05551173 Ry<br> estimated scf accuracy < 0.02263262 Ry<br> estimated scf accuracy < 0.01534456 Ry<br> estimated scf accuracy < 0.00799064 Ry<br> estimated scf accuracy < 0.00499060 Ry<br> estimated scf accuracy < 0.00272207 Ry<br> estimated scf accuracy < 0.00292201 Ry<br> estimated scf accuracy < 0.00255637 Ry<br> estimated scf accuracy < 0.00240574 Ry<br> estimated scf accuracy < 0.00239154 Ry<br> estimated scf accuracy < 0.00227859 Ry<br> estimated scf accuracy < 0.00214441 Ry<br> estimated scf accuracy < 0.00221359 Ry<br><br>I also tried changing the parameters ecutrho to 600.0; nband to 400 (from 72*4); decreasing degauss to 0.03 and trying mixing_mode='local-TF' as suggested but that does not help convergence one bit.<br><br>Please let me know if there is something else I could try which would help with these issues - decreasing the forces on each atom and helping convergence in the hexagonal supercell calculations.<br><br>Thank you very much<br>Dhruv Singh<br>Graduate Student, Purdue University ME<br><br>_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>http://www.democritos.it/mailman/listinfo/pw_forum<br></div></blockquote></div><br><div>
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div><span class="Apple-style-span" style="color: rgb(126, 126, 126); font-size: 16px; font-style: italic; "><br class="Apple-interchange-newline">§ Gabriele Sclauzero, EPFL SB ITP CSEA</span></div><div><font class="Apple-style-span" color="#7E7E7E"><i> PH H2 462, Station 3, CH-1015 Lausanne</i></font></div></span>
</div>
<br></div></body></html>