
<div dir="ltr">Hi,<div>After doing a vc-relax calculation, I used the relaxed cell parameters and coordinates to run a scf calculation with the same convergence criteria and cutoffs. However, the system's stress and pressure are shown to be very different. A similar issue I encountered on the QE forum, was fixed by tweaking the convergence against the Ecut value. But before I even started the computations, I verified convergence against Ecut and Kpoints.Â  I have used the NC PseudoDojo PP. Any assistance in this respect would be greatly appreciated.<br></div><div><br></div><div><b>Here is the output stress from the last step of vc-relax:</b></div><div>Â  Â  Â  Â  Â  total Â  stress Â (Ry/bohr**3) Â  Â  Â  Â  Â  Â  Â  Â  Â  (kbar) Â  Â  P= Â  Â  Â  Â 0.01<br>Â  Â 0.00000007 Â -0.00000001 Â  0.00000000 Â  Â  Â  Â  Â  Â 0.01 Â  Â  Â  -0.00 Â  Â  Â  Â 0.00<br>Â  -0.00000001 Â  0.00000007 Â -0.00000001 Â  Â  Â  Â  Â  -0.00 Â  Â  Â  Â 0.01 Â  Â  Â  -0.00<br>Â  Â 0.00000000 Â -0.00000001 Â  0.00000008 Â  Â  Â  Â  Â  Â 0.00 Â  Â  Â  -0.00 Â  Â  Â  Â 0.01<br><br></div><div><b>Here is the output stress from scf:</b><br></div><div>Â  Â  Â  Â  Â total Â  stress Â (Ry/bohr**3) Â  Â  Â  Â  Â  Â  Â  Â  Â  (kbar) Â  Â  P= Â  Â  Â  -0.81<br>Â  Â 0.00000173 Â  0.00001357 Â -0.00001496 Â  Â  Â  Â  Â  Â 0.26 Â  Â  Â  Â 2.00 Â  Â  Â  -2.20<br>Â  Â 0.00001357 Â  0.00000039 Â  0.00000247 Â  Â  Â  Â  Â  Â 2.00 Â  Â  Â  Â 0.06 Â  Â  Â  Â 0.36<br>Â  -0.00001496 Â  0.00000247 Â -0.00001869 Â  Â  Â  Â  Â  -2.20 Â  Â  Â  Â 0.36 Â  Â  Â  -2.75<br></div><div><br></div><div><b>vc-relaxÂ input:</b></div><div>&CONTROL<br>Â  calculation = 'vc-relax'<br>Â  etot_conv_thr = Â  1.0000000000d-06<br>Â  forc_conv_thr = Â  1.0000000000d-05<br>Â  outdir = './out/'<br>Â  restart_mode = "restart"<br>Â  pseudo_dir = '/home/ikharagpur/Desktop/Aritri/HAP_QE_NEW_FINAL/HAP/convergance/ecut'<br>Â  tprnfor = .true.<br>Â  tstress = .true.<br>Â  verbosity = 'high'<br>Â  nstep = 200<br>/<br>&SYSTEM<br>Â  degauss = Â  7.3498618000d-03<br>Â  ecutrho = Â  6.40000000000d+02<br>Â  ecutwfc = Â  7.00000000000d+01<br>Â  ibrav = 0<br>Â  nat = 41<br>Â  nosym = .false.<br>Â  ntyp = 6<br>Â  occupations = 'smearing'<br>Â  smearing = 'cold'<br>/<br>&ELECTRONS<br>Â  conv_thr = Â  1.000000000d-06<br>Â  electron_maxstep = 200<br>Â  mixing_beta = Â  4.0000000000d-01<br>/<br><br>&IONS<br>Â  Â  ion_dynamics = "bfgs"<br>/<br><br>&CELL<br>press_conv_thr = 2.0000000000d-01<br>/<br>ATOMIC_SPECIES<br>Ca Â  Â  40.078 Ca.upf<br>Na Â  Â  22.98977 Â Na.upf<br>C Â  Â  Â 12.01070 Â C.upf<br>O Â  Â  Â 15.99940 Â O.upf<br>P Â  Â  Â 30.973761 P.upf<br>F Â  Â  Â 18.99840 Â F.upf<br>ATOMIC_POSITIONS (angstrom)<br>Ca Â  Â  Â  Â  Â  Â 4.9287293961 Â  Â  Â  Â 2.5312304672 Â  Â  Â  -0.0636442378<br>Ca Â  Â  Â  Â  Â  Â 4.9349185508 Â  Â  Â  Â 2.9027082111 Â  Â  Â  Â 3.5370417241<br>Ca Â  Â  Â  Â  Â  -0.1577580449 Â  Â  Â  Â 5.2894025247 Â  Â  Â  Â 3.4702936429<br>Ca Â  Â  Â  Â  Â  -0.0739768540 Â  Â  Â  Â 5.6370186571 Â  Â  Â  Â 6.8292427360<br>Ca Â  Â  Â  Â  Â  Â 8.2864929497 Â  Â  Â  Â 1.8261397245 Â  Â  Â  Â 1.8328431870<br>Ca Â  Â  Â  Â  Â  Â 2.2609502006 Â  Â  Â  -0.0600744649 Â  Â  Â  Â 1.7143079190<br>Ca Â  Â  Â  Â  Â  Â 2.5803500800 Â  Â  Â  Â 8.0301042324 Â  Â  Â  Â 5.1805666090<br>Ca Â  Â  Â  Â  Â  Â 1.2884866870 Â  Â  Â  Â 2.1466631854 Â  Â  Â  Â 5.2813513288<br>Ca Â  Â  Â  Â  Â  Â 3.5158951682 Â  Â  Â  Â 6.0372944022 Â  Â  Â  Â 1.5643802295<br>Na Â  Â  Â  Â  Â  -3.4620658996 Â  Â  Â  Â 6.2795759096 Â  Â  Â  Â 5.1320586858<br>P Â  Â  Â  Â  Â  Â  3.0867436650 Â  Â  Â  Â 4.9588322631 Â  Â  Â  Â 5.0692930838<br>P Â  Â  Â  Â  Â  Â -2.6763991933 Â  Â  Â  Â 5.0837191424 Â  Â  Â  Â 1.6232868296<br>P Â  Â  Â  Â  Â  Â  7.4815069785 Â  Â  Â  Â 2.9665452844 Â  Â  Â  Â 5.2783593876<br>P Â  Â  Â  Â  Â  Â  3.6993356186 Â  Â  Â  Â 0.3106309483 Â  Â  Â  Â 5.0653678851<br>P Â  Â  Â  Â  Â  Â  1.0698739046 Â  Â  Â  Â 7.9066715577 Â  Â  Â  Â 1.7873359347<br>C Â  Â  Â  Â  Â  Â  1.2319164303 Â  Â  Â  Â 3.3343087228 Â  Â  Â  Â 1.7769185095<br>O Â  Â  Â  Â  Â  Â  5.4234243858 Â  Â  Â  Â 4.2067912520 Â  Â  Â  Â 1.5295885888<br>O Â  Â  Â  Â  Â  Â -0.6403864365 Â  Â  Â  Â 3.8380742071 Â  Â  Â  Â 5.3624856051<br>O Â  Â  Â  Â  Â  Â -0.8563492736 Â  Â  Â  Â 6.9192183026 Â  Â  Â  Â 5.0299634841<br>O Â  Â  Â  Â  Â  Â  5.4654395719 Â  Â  Â  Â 1.2958222863 Â  Â  Â  Â 1.8055312311<br>O Â  Â  Â  Â  Â  Â  1.5979572302 Â  Â  Â  Â 5.3945572754 Â  Â  Â  Â 5.0209743070<br>O Â  Â  Â  Â  Â  Â -1.3937275362 Â  Â  Â  Â 4.1855912819 Â  Â  Â  Â 1.6316785468<br>O Â  Â  Â  Â  Â  Â  6.1807498417 Â  Â  Â  Â 3.8078762598 Â  Â  Â  Â 5.2946494123<br>O Â  Â  Â  Â  Â  Â  5.0962113435 Â  Â  Â  Â 1.0053557461 Â  Â  Â  Â 5.0068747439<br>O Â  Â  Â  Â  Â  Â -0.2479127385 Â  Â  Â  Â 7.1008839512 Â  Â  Â  Â 1.9479371612<br>O Â  Â  Â  Â  Â  Â  3.2273831465 Â  Â  Â  Â 3.4634648722 Â  Â  Â  Â 5.5060881250<br>O Â  Â  Â  Â  Â  Â -2.6816919427 Â  Â  Â  Â 5.9171673155 Â  Â  Â  Â 2.9218077268<br>O Â  Â  Â  Â  Â  Â  7.4667745856 Â  Â  Â  Â 2.1504763747 Â  Â  Â  Â 3.9587589893<br>O Â  Â  Â  Â  Â  Â  3.0395685571 Â  Â  Â  Â 0.7662604690 Â  Â  Â  Â 6.4089715877<br>O Â  Â  Â  Â  Â  Â  1.6807850560 Â  Â  Â  Â 7.4400441153 Â  Â  Â  Â 0.4356710059<br>O Â  Â  Â  Â  Â  Â  3.7496780255 Â  Â  Â  Â 5.0285541608 Â  Â  Â  Â 3.6678802859<br>O Â  Â  Â  Â  Â  Â  7.4304843715 Â  Â  Â  Â 1.9745826851 Â  Â  Â  Â 6.4795274695<br>O Â  Â  Â  Â  Â  Â -2.5490400608 Â  Â  Â  Â 6.0618753239 Â  Â  Â  Â 0.4335569310<br>O Â  Â  Â  Â  Â  Â  2.1013955811 Â  Â  Â  Â 7.5530328240 Â  Â  Â  Â 2.8790462527<br>O Â  Â  Â  Â  Â  Â  2.8601969330 Â  Â  Â  Â 0.8778458001 Â  Â  Â  Â 3.8781941654<br>O Â  Â  Â  Â  Â  Â  3.7949586106 Â  Â  Â  Â 5.8956756067 Â  Â  Â  Â 6.0712536830<br>O Â  Â  Â  Â  Â  Â  0.8670346278 Â  Â  Â  Â 3.0547681296 Â  Â  Â  Â 2.9882678476<br>O Â  Â  Â  Â  Â  Â  1.1759918381 Â  Â  Â  Â 2.3994316994 Â  Â  Â  Â 0.8866049917<br>O Â  Â  Â  Â  Â  Â  1.5424153852 Â  Â  Â  Â 4.5568348312 Â  Â  Â  Â 1.4666970156<br>F Â  Â  Â  Â  Â  Â -0.0049782532 Â  Â  Â  -0.1503358719 Â  Â  Â  Â 1.6561706198<br>F Â  Â  Â  Â  Â  Â  0.1059848801 Â  Â  Â  Â 0.2892097568 Â  Â  Â  Â 5.0465194963<br><br><br>K_POINTS automatic<br>2 2 3 0 0 0<br>CELL_PARAMETERS {angstrom}<br>Â  9.397300 Â  0.000000 Â  0.000000<br>Â -4.698650 Â  8.138301 Â  0.000000<br>Â  0.000000 Â  0.000000 Â  6.878200<br></div><div><br></div><div><b>SCF input:</b></div><div>&CONTROL<br>Â  calculation = 'scf'<br>Â  etot_conv_thr = Â  1.0000000000d-06<br>Â  forc_conv_thr = Â  1.0000000000d-05<br>Â  outdir = './out/'<br>Â  restart_mode = "restart"<br>Â  pseudo_dir = '/home/ikharagpur/Desktop/Aritri/HAP_QE_NEW_FINAL/HAP/convergance/ecut'<br>Â  tprnfor = .true.<br>Â  tstress = .true.<br>Â  verbosity = 'high'<br>Â  nstep = 200<br>/<br>&SYSTEM<br>Â  degauss = Â  7.3498618000d-03<br>Â  ecutrho = Â  6.40000000000d+02<br>Â  ecutwfc = Â  7.00000000000d+01<br>Â  ibrav = 0<br>Â  nat = 41<br>Â  nosym = .false.<br>Â  ntyp = 6<br>Â  occupations = 'smearing'<br>Â  smearing = 'cold'<br>/<br>&ELECTRONS<br>Â  conv_thr = Â  1.000000000d-06<br>Â  electron_maxstep = 200<br>Â  mixing_beta = Â  4.0000000000d-01<br>/<br><br>&IONS<br>Â  Â  ion_dynamics = "bfgs"<br>/<br><br>&CELL<br>press_conv_thr = 2.0000000000d-01<br>/<br>ATOMIC_SPECIES<br>Ca Â  Â  40.078 Ca.upf<br>Na Â  Â  22.98977 Â Na.upf<br>C Â  Â  Â 12.01070 Â C.upf<br>O Â  Â  Â 15.99940 Â O.upf<br>P Â  Â  Â 30.973761 P.upf<br>F Â  Â  Â 18.99840 Â F.upf<br>CELL_PARAMETERS (angstrom)<br>Â  Â 9.403522322 Â  0.089346951 Â -0.122107790<br>Â  -4.624574075 Â  8.051846973 Â  0.082501283<br>Â  -0.089342394 Â  0.017896376 Â  6.858790895<br><br>ATOMIC_POSITIONS (angstrom)<br>Ca Â  Â  Â  Â  Â  Â 4.9770967179 Â  Â  Â  Â 2.5653460075 Â  Â  Â  -0.1337751051<br>Ca Â  Â  Â  Â  Â  Â 4.9372354934 Â  Â  Â  Â 2.9499612895 Â  Â  Â  Â 3.4521186383<br>Ca Â  Â  Â  Â  Â  -0.2034442691 Â  Â  Â  Â 5.2530858364 Â  Â  Â  Â 3.4846530148<br>Ca Â  Â  Â  Â  Â  -0.1137416011 Â  Â  Â  Â 5.6308719994 Â  Â  Â  Â 6.8378219159<br>Ca Â  Â  Â  Â  Â  Â 8.2977611076 Â  Â  Â  Â 1.8800297962 Â  Â  Â  Â 1.7736446849<br>Ca Â  Â  Â  Â  Â  Â 2.2413533700 Â  Â  Â  -0.0004311837 Â  Â  Â  Â 1.6492932113<br>Ca Â  Â  Â  Â  Â  Â 2.5998562512 Â  Â  Â  Â 8.0245043365 Â  Â  Â  Â 5.1748616941<br>Ca Â  Â  Â  Â  Â  Â 1.2506199411 Â  Â  Â  Â 2.1667262606 Â  Â  Â  Â 5.2671066291<br>Ca Â  Â  Â  Â  Â  Â 3.5541793955 Â  Â  Â  Â 6.0330065573 Â  Â  Â  Â 1.5227272245<br>Na Â  Â  Â  Â  Â  -3.4956016665 Â  Â  Â  Â 6.2398706047 Â  Â  Â  Â 5.1457184840<br>P Â  Â  Â  Â  Â  Â  3.0566889679 Â  Â  Â  Â 4.9793535438 Â  Â  Â  Â 5.0167609882<br>P Â  Â  Â  Â  Â  Â -2.6527279052 Â  Â  Â  Â 5.0483379029 Â  Â  Â  Â 1.6427735806<br>P Â  Â  Â  Â  Â  Â  7.4441823231 Â  Â  Â  Â 3.0411198798 Â  Â  Â  Â 5.2001935411<br>P Â  Â  Â  Â  Â  Â  3.6424844026 Â  Â  Â  Â 0.3624968438 Â  Â  Â  Â 4.9862395101<br>P Â  Â  Â  Â  Â  Â  1.1282750360 Â  Â  Â  Â 7.8820131992 Â  Â  Â  Â 1.8006951844<br>C Â  Â  Â  Â  Â  Â  1.2510149622 Â  Â  Â  Â 3.3000132805 Â  Â  Â  Â 1.7883657663<br>O Â  Â  Â  Â  Â  Â  5.4429703417 Â  Â  Â  Â 4.2737740659 Â  Â  Â  Â 1.4446893386<br>O Â  Â  Â  Â  Â  Â -0.6738195196 Â  Â  Â  Â 3.8099927194 Â  Â  Â  Â 5.4003635803<br>O Â  Â  Â  Â  Â  Â -0.8628277153 Â  Â  Â  Â 6.8847306005 Â  Â  Â  Â 5.0381510390<br>O Â  Â  Â  Â  Â  Â  5.4357139734 Â  Â  Â  Â 1.3538452729 Â  Â  Â  Â 1.7347050153<br>O Â  Â  Â  Â  Â  Â  1.5682645853 Â  Â  Â  Â 5.4071648285 Â  Â  Â  Â 4.9732803722<br>O Â  Â  Â  Â  Â  Â -1.3803282481 Â  Â  Â  Â 4.1392385455 Â  Â  Â  Â 1.6342972619<br>O Â  Â  Â  Â  Â  Â  6.1476486660 Â  Â  Â  Â 3.8856145701 Â  Â  Â  Â 5.2283043788<br>O Â  Â  Â  Â  Â  Â  5.0436443940 Â  Â  Â  Â 1.0499318227 Â  Â  Â  Â 4.9116948440<br>O Â  Â  Â  Â  Â  Â -0.1813703769 Â  Â  Â  Â 7.0691733500 Â  Â  Â  Â 1.9897343088<br>O Â  Â  Â  Â  Â  Â  3.1868368926 Â  Â  Â  Â 3.4940352445 Â  Â  Â  Â 5.4742190961<br>O Â  Â  Â  Â  Â  Â -2.6384743048 Â  Â  Â  Â 5.8674555732 Â  Â  Â  Â 2.9521034467<br>O Â  Â  Â  Â  Â  Â  7.4112431425 Â  Â  Â  Â 2.2431045948 Â  Â  Â  Â 3.8687306171<br>O Â  Â  Â  Â  Â  Â  2.9998260819 Â  Â  Â  Â 0.8176325894 Â  Â  Â  Â 6.3372934848<br>O Â  Â  Â  Â  Â  Â  1.7154095157 Â  Â  Â  Â 7.4080787828 Â  Â  Â  Â 0.4394134783<br>O Â  Â  Â  Â  Â  Â  3.7349447320 Â  Â  Â  Â 5.0315628428 Â  Â  Â  Â 3.6206036999<br>O Â  Â  Â  Â  Â  Â  7.3859083710 Â  Â  Â  Â 2.0358197649 Â  Â  Â  Â 6.3901175611<br>O Â  Â  Â  Â  Â  Â -2.5132765701 Â  Â  Â  Â 6.0348969987 Â  Â  Â  Â 0.4629313985<br>O Â  Â  Â  Â  Â  Â  2.1946381316 Â  Â  Â  Â 7.5655991791 Â  Â  Â  Â 2.8668076838<br>O Â  Â  Â  Â  Â  Â  2.7966145386 Â  Â  Â  Â 0.9308171326 Â  Â  Â  Â 3.8087744773<br>O Â  Â  Â  Â  Â  Â  3.7629853133 Â  Â  Â  Â 5.9368419567 Â  Â  Â  Â 6.0029497429<br>O Â  Â  Â  Â  Â  Â  0.8283603096 Â  Â  Â  Â 3.0530829701 Â  Â  Â  Â 2.9898800237<br>O Â  Â  Â  Â  Â  Â  1.2402701808 Â  Â  Â  Â 2.3363280977 Â  Â  Â  Â 0.9260179599<br>O Â  Â  Â  Â  Â  Â  1.5764439696 Â  Â  Â  Â 4.5094877888 Â  Â  Â  Â 1.4558604255<br>F Â  Â  Â  Â  Â  Â -0.0267498884 Â  Â  Â  -0.1708338921 Â  Â  Â  Â 1.6329080094<br>F Â  Â  Â  Â  Â  Â  0.0483918544 Â  Â  Â  Â 0.3169071048 Â  Â  Â  Â 5.0200898238<br><br>K_POINTS automatic<br>2 2 3 0 0 0<br></div><div>Regards,</div><div>Aritri RoyÂ </div><div>Junior Research Fellow,</div><div>Indian Institute of Technology Kharagpur</div><div>West Bengal, India</div><div><br></div></div>