<div dir="ltr">Dear Kondrin<div>Thank you for your suggestion. I actually followed the S Ponce lecture notes from the following link to calculate el-ph coupling constant. </div><div><a href="https://indico.ictp.it/event/8301/session/96/contribution/534/material/0/0.pdf">https://indico.ictp.it/event/8301/session/96/contribution/534/material/0/0.pdf</a>. But in your message here <a href="https://lists.quantum-espresso.org/pipermail/users/2022-September/049463.html">https://lists.quantum-espresso.org/pipermail/users/2022-September/049463.html</a> you suggest to include la2f = .true for q2r, matdyn.x files as well. After that we should post process the files to get lambda value without using lambda.x. Am I right in understanding your message clearly. ?? <br></div><div>Regards</div><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><div><br></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>