<div dir="ltr">Dear researchers<div><br></div><div>I have been trying to observe the effects of electric field on the bandstrcuture of multilayered TMDC materials. I have been able to calculate the effect of electric field on bandstructure of Monolayer MoS2. However I am kind of stuck at the multilayer case. The scf calculation would not converge. I am attaching the code here: </div><div><br></div><div><div>&CONTROL</div><div>      calculation = 'scf',</div><div>     restart_mode = 'from_scratch',</div><div>           outdir = 'D:\Quantum ESPRESSO 64-bit 5.1.2-mpich2\MoSSeAAA187EF\outdir',</div><div>           prefix = 'MoSSeAAA187EF',</div><div>       <span class="" style="white-space:pre">       </span>   lelfield=.true.,</div><div>               nberrycyc= 1</div><div>/</div><div>&SYSTEM</div><div>            ibrav = 0,<br></div><div>                a = 1.87,<br></div><div>              nat = 9,<br></div><div>           ntyp = 3,<br></div><div>      occupations = 'fixed',</div><div>        !smearing = 'methfessel-paxton',<br></div><div>      degauss = 0.001,<br></div><div>          !nspin=2,</div><div>      ecutwfc = 70,<br></div><div>     !ecutrho = 160,<br></div><div>      nbnd = 50,<br></div><div>           vdw_corr='grimme-d2',</div><div>        !london      = .true.</div><div>        london_s6   = 0.75    ! this is the default</div><div>        london_rcut = 200.00  ! this is the default</div><div>       !starting_magnetization=0.6, </div><div>       !input_dft='vdw-DF'</div><div>/</div><div><br></div><div>&ELECTRONS</div><div>    diagonalization='david',<br></div><div>    conv_thr =  1.0d-8,</div><div>    mixing_beta = 0.5,</div><div>    startingwfc='random',</div><div>    efield_cart(1)=0.d0,efield_cart(2)=0.d0,efield_cart(3)=0.01d0</div><div>/<br></div><div><br></div><div>CELL_PARAMETERS </div><div>1.5 0.8660254038 0.00000</div><div>1.5 -0.8660254038 0.00000</div><div>0.0000 0.00000 25.00000  </div><div>ATOMIC_SPECIES</div><div><br></div><div>S 32.065 S.pbe-hgh.UPF</div><div>Mo 95.94 Mo.pbe-sp-hgh.UPF</div><div>Se 78.96 Se.pbe-hgh.UPF</div><div><br></div><div><br></div><div>ATOMIC_POSITIONS (angstrom)</div><div>S        0.000000000   0.000000000   2.038611513<br></div><div>S        0.000000000   0.000000000   5.121719265</div><div>Se       0.000000000  -0.000000000   8.874475491</div><div>Se       0.000000000   0.000000000  12.265594331</div><div>S        0.000000000  -0.000000000  16.018257986</div><div>S        0.000000000   0.000000000  19.101377884</div><div>Mo       1.870000000   0.000000000   3.580650275</div><div>Mo       1.870000000   0.000000000  10.570030193</div><div>Mo       1.870000000   0.000000000  17.559283062</div><div><br></div><div><br></div><div>K_POINTS {automatic}</div><div>8 8 1   0 0 0</div></div><div><br></div><div>Am I missing something here? I am quite new to quantum espresso. Any help would be much appreciated. Thanks in advance.</div><div><br></div><div>Sincerely yours</div><div>Kanak Datta</div><div>EEE, BUET</div></div>