<div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">I am just starting to use this software so i have undoubtedly made a simple mistake, however I cannot seem to get SCF convergence after 100 itterations. Here is the input file I am using:</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><div class="gmail_default"> &CONTROL</div><div class="gmail_default">                       title = MoS2_relax ,</div><div class="gmail_default">                 calculation = 'relax' ,</div><div class="gmail_default">                  wf_collect = .true. ,</div><div class="gmail_default">                      outdir = '/home/users/drubin/MoS2/temp/' ,</div><div class="gmail_default">                      wfcdir = '/home/users/drubin/MoS2/cpu_temp/' ,</div><div class="gmail_default">                  pseudo_dir = '/home/users/drubin/MoS2/pseudo/' ,</div><div class="gmail_default">                      prefix = 'MoS2_relax' ,</div><div class="gmail_default">                   verbosity = 'high' ,</div><div class="gmail_default">                     tprnfor = .true. ,</div><div class="gmail_default"> /</div><div class="gmail_default"> &SYSTEM</div><div class="gmail_default">                       ibrav = 4,</div><div class="gmail_default">                   celldm(1) = 5.972,</div><div class="gmail_default">                   celldm(3) = 24.472,</div><div class="gmail_default">                         nat = 6,</div><div class="gmail_default">                        ntyp = 2,</div><div class="gmail_default">                     ecutwfc = 48 ,</div><div class="gmail_default">                     ecutrho = 407 ,</div><div class="gmail_default">            exxdiv_treatment = 'gygi-baldereschi' ,</div><div class="gmail_default"> /</div><div class="gmail_default"> &ELECTRONS</div><div class="gmail_default">             diagonalization = 'david' ,</div><div class="gmail_default">            diago_david_ndim = 8,</div><div class="gmail_default"> /</div><div class="gmail_default"> &IONS</div><div class="gmail_default"> /</div><div class="gmail_default">ATOMIC_SPECIES</div><div class="gmail_default">   Mo   95.96000  Mo.pbe-spn-kjpaw_psl.0.2.UPF</div><div class="gmail_default">    S   32.06600  S.pbe-n-kjpaw_psl.0.1.UPF</div><div class="gmail_default">ATOMIC_POSITIONS crystal</div><div class="gmail_default">   Mo      0.333333300    0.666666700    0.250000000</div><div class="gmail_default">   Mo      0.666666700    0.333333333    0.750000000</div><div class="gmail_default">    S      0.666666700    0.333333333    0.121000000</div><div class="gmail_default">    S      0.333333333    0.666666700    0.878999999</div><div class="gmail_default">    S      0.333333333    0.666666700    0.620999990</div><div class="gmail_default">    S      0.666666700    0.333333333    0.379000000</div><div class="gmail_default">K_POINTS automatic</div><div class="gmail_default">  3 3 3   1 1 1</div><div><br></div><div><br></div><div>Any help is greatly appreciated. </div><div><br></div><div>Dan R.</div></div></div>