<div dir="ltr">Hello everyone,<br><div>Thanks to the developers for providing such a versatile open-source platform for us to use.</div><div>My question today is about using a finite electric field on a 2D buckled structure; I want to see the SOC + electric field in version 6.4.1. I did relax, with the following input without SOC. </div><div>1. Without SOC, when I try to do scf, nscf, and bands, do I need the lelfiled flags still included? or just a normal <span aria-invalid="grammar" class="Lm ng" style="border-style:none;background:none">bands</span> calculation? Because in the nscf step, if I have the flags, I can't use the specific band path, including the high symmetry points (XGMX). I just have to use a dense mesh in the automatic option. </div><div>2. When I do SOC + lelfiled, do I relax with lelfiled and SOC again?</div><div>3. In SOC + lelfield, for the band structure calculation, do I keep the <span aria-invalid="spelling" class="LI ng" style="border-style:none;background:none">lelfield</span> + <span aria-invalid="grammar" class="Lm ng" style="border-style:none;background:none">SOC</span> flags or just the <span aria-invalid="grammar" class="Lm ng" style="border-style:none;background:none">SOC</span> flags?</div><div><br></div><div><span style="background-color:transparent">&CONTROL</span></div><div>    calculation = "relax"<br>    outdir      = "./t/"<br>    prefix      = "aelec"<br>    pseudo_dir  = "./"<br>    lelfield    = .true.,<br>    verbosity   = "high"<br>    nstep       = 100<br>/<br><br>&SYSTEM<br>    celldm(1)   =  12, <br>    degauss     =  0.01,<br>    ecutwfc     =  90,<br>    ecutrho     =  360,  <br>    ibrav       =  0,<br>    nat         =  8,<br>    ntyp        =  2,<br>/<br><br>&ELECTRONS<br>   conv_thr         =  1.00000e-08<br>   electron_maxstep =  600<br>   mixing_beta      =  1.00000e-01<br>   efield_cart(1)   =  0.d0,<br>   efield_cart(2)   =  0.d0,<br>   efield_cart(3)   =  0.001<br>/<br><br>&IONS<br>    ion_dynamics  = 'bfgs'<br>/<br><br>CELL_PARAMETERS {alat= 12.93834100}<br>   1.003011396   0.000005821   0.000000000<br>   0.000005821   1.003011396   0.000000000<br>   0.000000000   0.000000000   2.722061412<br><br>ATOMIC_SPECIES<br>Ga  69.723000   Ga.pbe-dnl-kjpaw_psl.1.0.0.UPF<br>P   30.973800   P.pbe-n-kjpaw_psl.1.0.0.UPF<br><br>ATOMIC_POSITIONS (crystal)<br><br>K_POINTS {automatic}<br> 6 6 1 0 0 0<br></div><div><br></div><div><div>I have been stuck on this for a while, searched the forum as well, <span class="gmail-Lm gmail-ng" style="background:none;border:0px;outline:none"><span aria-invalid="grammar" class="Lm ng" style="border-style:none;background:none">couldn't</span></span> find a proper workflow. Kindly help and guide. Kindly point out if I am doing anything else wrong.</div><div><br></div></div><div>Barsha Pal</div><div>Research Scholar</div><div>Department of Physics</div><div>Assam University Silchar</div><div>Silchar, India.</div></div>