<div>Dear all,</div><div>I'm trying to study interaction of  a Sulfur atom with a Gold surface under different electric fields and calculate its binding energy. I'm using saw-like potentials for adding electric field. I have two questions:</div>

<div>1- Potential parameters : As I read the manual and pw-forum, saw like potentials should not change slope where atoms exist but with this condition, required parameters would not be unique and one can assign different values to emaxpos and eopreg. What would be the difference if I choose ,for example, <i><b>emaxpos=0.45 & eopreg=0.1</b></i> or <i><b>emaxpos=0.49 & eopreg=0.02</b></i>? What  happens if I have a compact unit cell and the place where slope changes is near atoms? </div>


<div><br></div><div>2- Larger electric fields: When I do my calculations for larger electric fields, it needs more scf steps to converge as for some cases I need to change "electron_maxstep" to a higher value than 100. Is there a trick to speed up convergency?</div>

<div><br></div><div>Any help would be greatly appreciated.</div><div><br></div><div><br></div><div>below is a sample of the input files which I use:</div><div><br></div><div><div> &CONTROL</div><div>                 calculation = 'scf' ,</div>

<div>             restart_mode = 'from_scratch' ,</div><div>                         outdir = '/home/farzad/tmp/' ,</div><div>                pseudo_dir = '/home/farzad/Desktop/espresso-5.0.1-GPU/pseudo/' ,</div>

<div>                       tstress = .false. ,</div><div>                       tprnfor = .true. ,</div><div>                        tefield = .true. ,</div><div> /</div><div> &SYSTEM</div><div>                                  ibrav = 14,</div>

<div>                           celldm(1) = 15.6688,</div><div>                           celldm(2) = 1,</div><div>                           celldm(3) = 1.5,</div><div>                           celldm(4) = 0,</div><div>

                           celldm(5) = 0,</div><div>                           celldm(6) = 0,</div><div>                                     nat = 9,</div><div>                                   ntyp = 2,</div><div>                              ecutwfc = 40 ,</div>

<div>                              ecutrho = 320 ,</div><div>                       occupations = 'smearing' ,</div><div>                             degauss = 0.01 ,</div><div>                            smearing = 'marzari-vanderbilt' ,</div>

<div>                                  nspin = 2 ,</div><div>   starting_magnetization(1) = 0.02,</div><div>   starting_magnetization(2) = 0.01,</div><div>                 exxdiv_treatment = 'gygi-baldereschi' ,</div>

<div>                                      edir = 3 ,</div><div>                              emaxpos = 0.55 ,</div><div>                                 eopreg = 0.1 ,</div><div>                                   eamp = 0.04 ,</div>

<div> /</div><div> &ELECTRONS</div><div>                electron_maxstep = 500 ,</div><div>                    diagonalization = 'cg' ,</div><div> /</div><div>ATOMIC_SPECIES</div><div>   Au  196.96657  Au.pw91-van_ak.UPF </div>

<div>    S   32.06500  S.pw91-van_ak.UPF </div><div>ATOMIC_POSITIONS bohr </div><div>   Au      0.000000000    0.000000000    0.000000000    </div><div>   Au      3.917200000    3.917200000    0.000000000    </div><div>   Au      7.834400000    0.000000000    0.000000000    </div>

<div>   Au     11.751600000    3.917200000    0.000000000    </div><div>   Au      0.000000000    7.834400000    0.000000000    </div><div>   Au      7.834400000    7.834400000    0.000000000    </div><div>   Au      3.917200000   11.751600000    0.000000000    </div>

<div>   Au     11.751600000   11.751600000    0.000000000    </div><div>    S      7.834400000    3.917200000    2.050000000    </div><div>K_POINTS automatic </div><div>  4 4 4   0 0 0 </div><div><br></div><div><br></div>

<div>Best Regards,</div><div>----------------------------</div><div>Farzad Shirazian</div><div>Sharif University of Technology</div></div>