<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
Hello Ramesh,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof ContentPasted0 ContentPasted1">
Norm-conserving (NC) pseudopotentials (PP) require higher energy cutoffs. You need NC PP for electron-phonon calculations. NC PP from PseudoDojo website are very good, at least for what I've tested on organic materials. Then you don't need ecutrho for NC PP,
which are made to work with default ecutrho = 4 x ecutwfc.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof ContentPasted0 ContentPasted1">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof ContentPasted0 ContentPasted1 ContentPasted2">
So, you should try first to increase your cutoff, up to 200 Ry should do the job to have the asymptotic behavior with PseudoDojo NC PP (unfortunately, other NC PP may require much higher cutoffs to get the asymptotic behavior). Then you can choose the cutoff
where the asymptotic error is just below 1 meV/atom.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof ContentPasted0 ContentPasted1 ContentPasted2">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof ContentPasted0 ContentPasted1 ContentPasted2">
After you do a vc-relax simulation, always check the verification step at the end of the output. If the pressure tensor is different from the pressure tensor of the last optimization step, then it means your cutoff is not large enough, or your system has converging
issues (initial geometry, k-points...).<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof ContentPasted0 ContentPasted1 ContentPasted2">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof ContentPasted0 ContentPasted1 ContentPasted2 ContentPasted3 ContentPasted4 ContentPasted5">
I think your forc_conv_thr of 1e-4 is good. You can tighten a little bit etot_conv_thr to 1e-5 and press_conv_thr to 0.1 kbar.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof ContentPasted0 ContentPasted1 ContentPasted2 ContentPasted3 ContentPasted4 ContentPasted5">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof ContentPasted0 ContentPasted1 ContentPasted2 ContentPasted3 ContentPasted4 ContentPasted5">
Regards,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof ContentPasted0 ContentPasted1 ContentPasted2 ContentPasted3 ContentPasted4 ContentPasted5">
Xavier<br>
</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> users <users-bounces@lists.quantum-espresso.org> on behalf of Ramesh Kumar Kamadurai via users <users@lists.quantum-espresso.org><br>
<b>Sent:</b> Thursday, October 27, 2022 3:23 AM<br>
<b>To:</b> users@lists.quantum-espresso.org <users@lists.quantum-espresso.org><br>
<b>Subject:</b> [QE-users] Electron Phonon Calculation: tolerance limit for P</font>
<div> </div>
</div>
<div>
<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>
</div>
</body>
</html>