<div dir="ltr"><div><div><div><div><div><div>Dear QE users and developers<br></div>I do vc-relax with BFGS algorithm as follows:<br><br>&CONTROL<br>  calculation = "vc-relax",<br>  pseudo_dir  = '/home/uriargaman/espresso-5.0.3/pseudo/',<br>  outdir      = '/home/uriargaman/tmp/tmp28898/',<br>  tstress=.true.<br>  verbosity='high',<br> /<br>&system<br>    ibrav=  4, celldm(1) =5.548454954360, celldm(3) =1.581979631777413 , nat= 2, ntyp= 1,<br>    ecutwfc =80, ecutrho=1200, <br>    occupations='smearing', smearing="methfessel-paxton", degauss=0.01,<br> /<br>&ELECTRONS<br>  conv_thr    = 1.D-9,<br>  electron_maxstep=1000<br>  mixing_beta = 0.7D0,<br>/<br>&IONS<br>   ion_dynamics="bfgs",<br>/<br>&cell<br>  cell_dynamics="bfgs",<br>  press=0.6,<br>  press_conv_thr=0.01,<br>/<br>ATOMIC_SPECIES<br>  Ti   47.9  Ti.pbe-sp-van_ak.UPF<br>ATOMIC_POSITIONS (crystal)<br>   Ti      0.666666666    0.333333333    0.250000000    0  0  0<br>   Ti      0.333333333    0.666666666    0.750000000<br>K_POINTS {automatic}<br> 24 24 24 0 0 0<br><br></div>In the last iteration (after begin final coordinates) I got stress:<br>          total   stress  (Ry/bohr**3)                   (kbar)     P=    0.71<br>   0.00000485  -0.00000000  -0.00000000          0.71     -0.00     -0.00<br>  -0.00000000   0.00000485   0.00000000         -0.00      0.71      0.00<br>  -0.00000000   0.00000000   0.00000481         -0.00      0.00      0.71<br><br>and the stress in the iteration before (last iteration before  begin final coordinates):<br>  total   stress  (Ry/bohr**3)                   (kbar)     P=    0.73<br>   0.00000499  -0.00000000   0.00000000          0.73     -0.00      0.00<br>  -0.00000000   0.00000499   0.00000000         -0.00      0.73      0.00<br>   0.00000000   0.00000000   0.00000495          0.00      0.00      0.73<br><br>I did'nt reach the pressure threshold although I have very good convergence with respect to all numerical parameters. The threshold is achieved only for hydrostatic condition but not for the value of the stress (P=0.73 kbar instead of 0.6 kbar). Is the parameter "press_conv_thr" is just for the hydrostatic condition and not for the absolute value of the stress? How can I get more accurate stress?<br></div>Thank you very much<br></div>Uri Argaman<br></div>Ben-Gurion University<br></div>Israel<br></div>