<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Sorry, Chitta, but this does not make much sense. k-point sampling is intimately related to periodicity: no periodicity, No k-point sampling. The dependence of any quantity of a non periodic system on the number of k point simply shows that the size of the supercell may be too small. Using k-point sampling may in this case speed up convergence towards the infinite-size limit, but in doing so one should know what one is doing and why, which this discussion shows not to be necessarily the case. The rule of thumb is: if your result depend too much on the number of k point, increase the size of the cell until it does not. The number of monolayers has nothing to do with the rate of convergence. The vacuum in between pairs of monolayers does. Regards — Stefano B<div class=""><br class=""><div><br class=""><blockquote type="cite" class=""><div class="">On 9 Oct 2018, at 06:34, Sai Pavan Chitta <<a href="mailto:chitta.sai@iitgn.ac.in" class="">chitta.sai@iitgn.ac.in</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="auto" class=""><div class="">Dear Sabike,</div><div dir="auto" class=""><br class=""></div><div dir="auto" class=""><div dir="auto" class="">In my experience of working with monolayers, I found that the number of k-points in the k-mesh in z direction aalso affects the convergence. The energy starts converging the moment I increase the number of k-points in z-direction. It also makes sense as you are using two monolayers in your model and thus requires more k-points in the z-direction. </div><div dir="auto" class=""><br class=""></div><div dir="auto" class="">Moreover, you can make the convergence threshold more stringent to about 1e-8 for a proper convergence.</div><div dir="auto" class=""><br class=""></div><div dir="auto" class="">Regards</div><div dir="auto" class=""><div data-smartmail="gmail_signature" dir="auto" class="">Chitta Sai Pavan<br class="">IIT Gandhinagar</div></div><br class=""><br class=""><div class="gmail_quote" dir="auto"><div dir="ltr" class="">On Tue, 9 Oct 2018, 08:18 Sabike Ghasemi, <<a href="mailto:s.ghasemi@du.ac.ir" class="">s.ghasemi@du.ac.ir</a>> wrote:<br class=""></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class=""><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt;" class=""><div class="">dear all<br class=""></div><div class="">i want to relax atomic coordinates and lattic constants for graphene bilayer. but it dosen't converge.<br class=""></div><div class="">input file is:<br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""> &CONTROL<br class=""> title = GRAPHENEBILAYER ,<br class=""> calculation = 'vc-relax' ,<br class=""> restart_mode = 'from_scratch' ,<br class=""> wf_collect = .true. ,<br class=""> outdir = './out' ,<br class=""> wfcdir = './out' ,<br class=""> pseudo_dir = './' ,<br class=""> prefix = 'GBI' ,<br class=""> verbosity = 'high' ,<br class=""> etot_conv_thr = 1.0D-4 ,<br class=""> forc_conv_thr = 1.0D-4 ,<br class=""> tstress = .true. ,<br class=""> tprnfor = .true. ,<br class=""> /<br class=""> &SYSTEM<br class=""> ibrav = 4,<br class=""> celldm(1) = 4.648726,<br class=""> celldm(3) = 9.512195384,<br class=""> nat = 4,<br class=""> ntyp = 1,<br class=""> ecutwfc = 50 ,<br class=""> ecutrho = 500 ,<br class=""> nosym = .false. ,<br class=""> occupations = 'smearing' ,<br class=""> degauss = 0.002 ,<br class=""> smearing = 'gaussian' ,<br class=""> /<br class=""> &ELECTRONS<br class=""> electron_maxstep = 200,<br class=""> conv_thr = 1.0D-6 ,<br class=""> mixing_beta = 0.3 ,<br class=""> diagonalization = 'david' ,<br class=""> /<br class=""> &IONS<br class=""> ion_dynamics = 'bfgs' ,<br class=""> /<br class=""> &CELL<br class=""> cell_dynamics = 'bfgs' ,<br class=""> cell_factor = 2.D0 ,<br class=""> cell_dofree = '2Dxy' ,<br class=""> /<br class="">ATOMIC_SPECIES<br class=""> C 12.00000 C.blyp-mt.UPF <br class="">ATOMIC_POSITIONS crystal <br class=""> C 0.333333330 0.666666660 0.572649573 <br class=""> C 1.000000000 0.000000000 0.572649753 <br class=""> C 0.666666660 0.333333330 0.427350427 <br class=""> C 1.000000000 0.000000000 0.427350427 <br class="">K_POINTS automatic <br class=""> 16 16 1 0 0 0 <br class=""><br class=""></div><div class=""><br class=""></div><div class="">best, <br class=""></div><div class="">sabike<br class=""></div></div></div>_______________________________________________<br class="">
users mailing list<br class="">
<a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noreferrer" class="">users@lists.quantum-espresso.org</a><br class="">
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer" target="_blank" class="">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div></div></div>
_______________________________________________<br class="">users mailing list<br class=""><a href="mailto:users@lists.quantum-espresso.org" class="">users@lists.quantum-espresso.org</a><br class="">https://lists.quantum-espresso.org/mailman/listinfo/users</div></blockquote></div><br class=""><div class="">
<span class="Apple-style-span" style="border-collapse: separate; font-variant-ligatures: normal; font-variant-position: normal; font-variant-numeric: normal; font-variant-alternates: normal; font-variant-east-asian: normal; line-height: normal; orphans: 2; text-indent: 0px; widows: 2; border-spacing: 0px;"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div class=""><div class=""><span class="Apple-style-span" style="border-collapse: separate; font-variant-ligatures: normal; font-variant-position: normal; font-variant-numeric: normal; font-variant-alternates: normal; font-variant-east-asian: normal; line-height: normal; orphans: 2; text-indent: 0px; widows: 2; border-spacing: 0px;"><span class="Apple-style-span" style="border-collapse: separate; font-variant-ligatures: normal; font-variant-position: normal; font-variant-numeric: normal; font-variant-alternates: normal; font-variant-east-asian: normal; line-height: normal; orphans: 2; text-indent: 0px; widows: 2; border-spacing: 0px;"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div class=""><span style="font-size: 12px;" class="">— </span></div><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div class=""><div style="margin: 0px;" class=""><font class="Apple-style-span"><span style="font-size: 12px;" class="">Stefano Baroni - Trieste — </span><a href="http://stefano.baroni.me" style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-stroke-width: 0px;" class="">http://stefano.baroni.me</a><span style="font-size: 12px;" class=""> </span></font></div></div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-stroke-width: 0px;" class=""><br class=""></div></div></div></div></div></div></div></span></span></div></div></div></span><br class="Apple-interchange-newline"><br class="Apple-interchange-newline">
</div>
<br class=""></div></body></html>