<div dir="ltr">Dear Xavier<div><font face="georgia, serif">Thank you for the detailed information. Yes I have already used pseudo potential from the </font><span style="color:rgb(0,0,0)"><font face="georgia, serif">PseudoDojo library and observed exactly what you have described. After the VC-relax the verification the verification scf calculation showed a completely different pressure tensor. I did not increase the cut-off till 200 Ry as you suggested but stopped at 90 Ry. Will give it a try with 200 Ry. Thanks. Regards</font></span></div><div><span style="color:rgb(0,0,0)"><font face="georgia, serif">K Ramesh Kumar</font></span></div><div><span style="color:rgb(0,0,0)"><font face="georgia, serif">Assistant Professor</font></span></div><div><span style="color:rgb(0,0,0)"><font face="georgia, serif">GITAM, University</font></span></div><div><span style="color:rgb(0,0,0)"><font face="georgia, serif">Vizag, India</font></span></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Oct 27, 2022 at 1:53 PM Ramesh Kumar Kamadurai <<a href="mailto:rkamadur@gitam.edu">rkamadur@gitam.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Dear users<div>I am trying to calculate el-Ph coupling constant and other properties for a superconductor. The input file is attached. I am a bit stuck with VC-Relax calculation as the system is not able to stabilize at 0.0 kbar. I followed all the suggestions given in forum</div><div>1. Increased Cut-off to maximum</div><div>2. Tried USPP,NC, Mixed and currently running with PBE-Sol</div><div>3. Sequentially done 2 VC Relax calculation but i am not able to sort out the issue.</div><div>&CONTROL<br>    calculation   = "vc-relax"<br>    etot_conv_thr =  4.00000e-05<br>    forc_conv_thr =  1.00000e-04<br>    max_seconds   =  1.72800e+05<br>    nstep         = 100<br>    outdir        = "/home/sjacobs/tmp/"<br>    prefix        = "PDZ-VCRELAX"<br>    pseudo_dir    = "C:\Users\RAMESH\.burai\.pseudopot"<br>    tprnfor       = .TRUE.<br>    tstress       = .TRUE.<br>    verbosity     = "high"<br>/<br><br>&SYSTEM<br>    a           =  6.36512e+00<br>    degauss     =  7.34986e-03<br>    ecutrho     =  6.40000e+02<br>    ecutwfc     =  8.00000e+01<br>    ibrav       = 2<br>    nat         = 4<br>    nosym       = .FALSE.<br>    nspin       = 1<br>    ntyp        = 3<br>    occupations = "smearing"<br>    smearing    = "marzari-vanderbilt"<br>/<br><br>&ELECTRONS<br>    conv_thr         =  8.00000e-10<br>    electron_maxstep = 80<br>    mixing_beta      =  4.00000e-01<br>    startingpot      = "atomic"<br>    startingwfc      = "atomic+random"<br>/<br><br>&IONS<br>    ion_dynamics = "bfgs"<br>/<br><br>&CELL<br>    cell_dofree    = "all"<br>    cell_dynamics  = "bfgs"<br>    press          =  0.00000e+00<br>    press_conv_thr =  5.00000e-01<br>/<br><br>K_POINTS {automatic}<br>12 12 12  0 0 0<br><br>ATOMIC_SPECIES<br>Ga     69.72300  Ga.pbesol-dnl-rrkjus_psl.1.0.0.UPF<br>Pd    106.42000  Pd_ONCV_PBEsol-1.0.upf<br>Zr     91.22400  zr_pbesol_v1.uspp.F.UPF<br><br>ATOMIC_POSITIONS {crystal}<br>Zr      0.000000   0.000000   0.000000<br>Ga      0.493162   0.493162   0.493162<br>Pd      0.739743   0.739743   0.739743<br>Pd      0.246581   0.246581   0.246581<br><br>The output stress is given below</div><div>entering subroutine stress ...<br>          total   stress  (Ry/bohr**3)                   (kbar)     P=    2.67<br>     entering subroutine stress ...<br>          total   stress  (Ry/bohr**3)                   (kbar)     P=   -0.41<br>     entering subroutine stress ...<br>          total   stress  (Ry/bohr**3)                   (kbar)     P=   -0.82<br>     entering subroutine stress ...<br>          total   stress  (Ry/bohr**3)                   (kbar)     P=    0.56<br>     entering subroutine stress ...<br>          total   stress  (Ry/bohr**3)                   (kbar)     P=    0.13<br>     entering subroutine stress ...<br>          total   stress  (Ry/bohr**3)                   (kbar)     P=   -0.08<br>     entering subroutine stress ...<br>          total   stress  (Ry/bohr**3)                   (kbar)     P=   -0.23<br>     entering subroutine stress ...<br>          total   stress  (Ry/bohr**3)                   (kbar)     P=   -0.35<br>     entering subroutine stress ...<br>          total   stress  (Ry/bohr**3)                   (kbar)     P=   -0.50<br>     entering subroutine stress ...<br>          total   stress  (Ry/bohr**3)                   (kbar)     P=   -0.63<br>     entering subroutine stress ...<br>          total   stress  (Ry/bohr**3)                   (kbar)     P=   -0.77<br>     entering subroutine stress ...<br>          total   stress  (Ry/bohr**3)                   (kbar)     P=   -0.83<br></div><div>My question is </div><div>What is the optimal pressure value to avoid negative frequencies for electron phonon calculation ? </div><div>Is there any mistake in my input file ?</div><div>Thanks in advance</div><div>Regards</div><div>Ramesh</div><div>Asst. Professor</div><div>GITAM University</div><div>Vizag, India</div><div><br></div></div>
</blockquote></div>