<div dir="ltr"><div><div>I have done my vc-relax successfully but while performing scf calculation I am getting an error: convergence is not achieved in 100 steps.Please help me in this issue.I have made ecutwfc=60.How to resolve this error.How can I converge my process.<br><br></div>regards,<br></div>Anindya Bose<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Mar 6, 2017 at 10:34 PM, Ari P Seitsonen <span dir="ltr"><<a href="mailto:Ari.P.Seitsonen@iki.fi" target="_blank">Ari.P.Seitsonen@iki.fi</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
Dear Anindya,<br>
<br>
  Please provide your affiliation - to respect the other subscribers of the forum (please go through the instructions on how to post to the forum).<br>
<br>
  Adding to the other items mentioned in the previous responses, I would like to point out that apparently you want to relax the lateral lattice constant of the cell, but without further options also the vertical lattice constant is optimised, and probably the vacuum is either reducing or increasing (it never goes to zero, due to numerical precision if nothing else; also the DFT+D2 dispersion correction has a long tail): Please check the input option 'cell_dofree'.<br>
<br>
  'cell_factor = 15' does not make much sense; your threhould variables 'etot_conv_thr' and 'forc_conv_thr' look very loose to me. Are you sure about "disk_io = 'high'"?<br>
<br>
    Greetings,<br>
<br>
       apsi<br>
<br>
-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=<wbr>*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=<wbr>-=*=-=*=-=*=-=*=-<br>
  Ari Paavo Seitsonen / <a href="mailto:Ari.P.Seitsonen@iki.fi" target="_blank">Ari.P.Seitsonen@iki.fi</a> / <a href="http://www.iki.fi/~apsi/" rel="noreferrer" target="_blank">http://www.iki.fi/~apsi/</a><br>
    Ecole Normale Supérieure (ENS), Département de Chimie, Paris<br>
    Mobile (F) : +33 789 37 24 25    (CH) : +41 79 71 90 935<div class="HOEnZb"><div class="h5"><br>
<br>
<br>
On Mon, 6 Mar 2017, Anindya Bose wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear Sir,<br>
I have used this code but I didn't get the direct band gap structure of monolayer WSe2.<br>
&CONTROL<br>
  calculation='vc-relax',<br>
  outdir='monolayer WSe2',<br>
  prefix='calc',<br>
  pseudo_dir='/home/anindya/Desk<wbr>top/pseudopotential',<br>
  verbosity='high',<br>
  disk_io='high',<br>
  wf_collect=.true.,<br>
  nstep = 50,<br>
  etot_conv_thr=1d-02,<br>
  forc_conv_thr=1d-02,<br>
/<br>
<br>
&SYSTEM<br>
  ibrav=4,<br>
  celldm(1)=6.2020811408d0, celldm(3)=4.3566118221d0,<br>
  nat=3,<br>
  ntyp=2,<br>
  ecutwfc=40.0d0,<br>
  ecutrho =130.0d0,<br>
  nbnd=200,<br>
  vdw_corr='Grimme-D2',<br>
  force_symmorphic=.true.,<br>
  input_dft='PBE',<br>
  occupations='fixed',<br>
/<br>
<br>
&ELECTRONS<br>
  diagonalization='david',<br>
  conv_thr=1d-08,<br>
  mixing_mode='plain',<br>
  mixing_beta=0.700d0,<br>
/<br>
<br>
&ions<br>
  ion_dynamics ='bfgs',<br>
/<br>
 <br>
&cell<br>
  cell_dynamics ='bfgs',<br>
  cell_factor=15,<br>
/<br>
<br>
ATOMIC_SPECIES<br>
  Se 78.960000d0 Se_ONCV_PBE-1.0_r.oncvpsp.upf<br>
  W 183.840000d0 W_ONCV_PBE-1.1_r.oncvpsp.upf<br>
<br>
ATOMIC_POSITIONS {alat}<br>
  W    0.5000000000d0   0.2886751346d0   2.1782995738d0<br>
  Se   0.5000000000d0  -0.2886751346d0   1.6692977456d0<br>
  Se   0.5000000000d0  -0.2886751346d0   2.6873014019d0<br>
<br>
K_POINTS {automatic}<br>
  10 10 1 0 0 0<br>
<br>
I have enclosed my output band structure with this mail.Can you please help me in this regard.How can I get the correct band diagram of monolayer<br>
WSe2(direct bandgap). I have used all the codes for spin calculation and I got the spin orbit coupling but didn't find perfect band diagram.I will<br>
be waiting for your response.<br>
</blockquote>
</div></div><br>______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br></blockquote></div><br></div>