<meta http-equiv="Content-Type" content="text/html; charset=GB18030"><div><span style="font-size: 14px;">Dear QE developer</span></div><div><span style="font-size: 14px;"><br></span></div><div><span style="font-size: 14px;">I did the structural relaxation calculation for Pb at pressure of 3kbar£¬ this is relax.in</span></div><div><span style="font-size: 14px;"><br></span></div><div><span style="font-size: 14px;">&CONTROL</span></div><div><span style="font-size: 14px;">  prefix='PbTaSe2',</span></div><div><span style="font-size: 14px;">  calculation='vc-relax', pseudo_dir='../../', outdir='./tmp', verbosity='high',</span></div><div><span style="font-size: 14px;">  tprnfor=.true., tstress=.true., forc_conv_thr=1.0d-4, nstep=100,</span></div><div><span style="font-size: 14px;">/</span></div><div><span style="font-size: 14px;">&SYSTEM</span></div><div><span style="font-size: 14px;">  ibrav= 0, nat= 4, ntyp= 3, </span></div><div><span style="font-size: 14px;">  occupations = 'smearing', smearing = 'gauss', degauss = 1.0d-2,</span></div><div><span style="font-size: 14px;">  ecutwfc = 50, ecutrho = 500,</span></div><div><span style="font-size: 14px;">/</span></div><div><span style="font-size: 14px;">&ELECTRONS</span></div><div><span style="font-size: 14px;">  conv_thr = 1.0d-8</span></div><div><span style="font-size: 14px;">  mixing_beta = 0.7d0</span></div><div><span style="font-size: 14px;">/</span></div><div><span style="font-size: 14px;">&IONS</span></div><div><span style="font-size: 14px;">/</span></div><div><span style="font-size: 14px;">&CELL</span></div><div><span style="font-size: 14px;">  cell_dynamics  = "bfgs"</span></div><div><span style="font-size: 14px;">  press          =  2.70000e+00</span></div><div><span style="font-size: 14px;">  press_conv_thr =0.1</span></div><div><span style="font-size: 14px;">/</span></div><div><span style="font-size: 14px;">ATOMIC_SPECIES</span></div><div><span style="font-size: 14px;">  Ta 180.9479 Ta.UPF</span></div><div><span style="font-size: 14px;">  Pb 207.2 Pb.UPF</span></div><div><span style="font-size: 14px;">  Se 78.96 Se.UPF</span></div><div><span style="font-size: 14px;">CELL_PARAMETERS (angstrom)</span></div><div><span style="font-size: 14px;">   3.4834700000  0.0000000000  0.0000000000</span></div><div><span style="font-size: 14px;">  -1.7417350000  3.0167735133  0.0000000000</span></div><div><span style="font-size: 14px;">   0.0000000000  0.0000000000  9.4755300000</span></div><div><span style="font-size: 14px;">ATOMIC_POSITIONS (crystal)</span></div><div><span style="font-size: 14px;">  Pb  0.0000000000  0.0000000000  0.0000000000</span></div><div><span style="font-size: 14px;">  Se  0.0000000000  0.0000000000  0.3246710000</span></div><div><span style="font-size: 14px;">  Ta  0.3333330000  0.6666670000  0.5000000000</span></div><div><span style="font-size: 14px;">  Se  0.0000000000  0.0000000000  0.6753290000</span></div><div><span style="font-size: 14px;">K_POINTS {automatic}</span></div><div><span style="font-size: 14px;">  15 15 5 0 0 0</span></div><div><span style="font-size: 14px;"><br></span></div><div><span style="font-size: 14px;">In relax.in, press=3.0kbar. But after pw.x is done, relax.out show that total stress  (Ry/bohr**3)  (kbar)     P=-0.83, and This is inconsistent with press=3. this is part of relax.out</span></div><div><span style="font-size: 14px;"><br></span></div><div><span style="font-size: 14px;">[rencong@mu01 phonon]$ grep 'Computing stress (Cartesian axis) and pressure' -A 15 relax.out</span></div><div><br></div><div>....</div><div><br></div><div><span style="font-size: 14px;">     Computing stress (Cartesian axis) and pressure</span></div><div><span style="font-size: 14px;"><br></span></div><div><span style="font-size: 14px;">          total   stress  (Ry/bohr**3)                   (kbar)     P=        2.93</span></div><div><span style="font-size: 14px;">   0.00002030  -0.00000000   0.00000000            2.99       -0.00        0.00</span></div><div><span style="font-size: 14px;">  -0.00000000   0.00002030  -0.00000000           -0.00        2.99       -0.00</span></div><div><span style="font-size: 14px;">  -0.00000000   0.00000000   0.00001921           -0.00        0.00        2.83</span></div><div><span style="font-size: 14px;"><br></span></div><div><span style="font-size: 14px;">     kinetic stress (kbar)  26877.23      0.00      0.00</span></div><div><span style="font-size: 14px;">                                0.00  26877.23     -0.00</span></div><div><span style="font-size: 14px;">                                0.00      0.00  26796.57</span></div><div><span style="font-size: 14px;"><br></span></div><div><span style="font-size: 14px;">     local   stress (kbar)-103340.53      0.03      0.00</span></div><div><span style="font-size: 14px;">                                0.03-103340.50      0.00</span></div><div><span style="font-size: 14px;">                                0.00      0.00  89394.48</span></div><div><span style="font-size: 14px;"><br></span></div><div><span style="font-size: 14px;">     nonloc. stress (kbar)   7797.15      0.00      0.00</span></div><div><span style="font-size: 14px;">--</span></div><div><span style="font-size: 14px;">     Computing stress (Cartesian axis) and pressure</span></div><div><span style="font-size: 14px;"><br></span></div><div><span style="font-size: 14px;">          total   stress  (Ry/bohr**3)                   (kbar)     P=        2.96</span></div><div><span style="font-size: 14px;">   0.00002009  -0.00000000   0.00000000            2.96       -0.00        0.00</span></div><div><span style="font-size: 14px;">  -0.00000000   0.00002009  -0.00000000           -0.00        2.96       -0.00</span></div><div><span style="font-size: 14px;">   0.00000000  -0.00000000   0.00002028            0.00       -0.00        2.98</span></div><div><span style="font-size: 14px;"><br></span></div><div><span style="font-size: 14px;">     kinetic stress (kbar)  26879.49      0.00      0.00</span></div><div><span style="font-size: 14px;">                                0.00  26879.49     -0.00</span></div><div><span style="font-size: 14px;">                                0.00     -0.00  26798.86</span></div><div><span style="font-size: 14px;"><br></span></div><div><span style="font-size: 14px;">     local   stress (kbar)-103299.75      0.03      0.00</span></div><div><span style="font-size: 14px;">                                0.03-103299.72      0.00</span></div><div><span style="font-size: 14px;">                                0.00      0.00  89348.84</span></div><div><span style="font-size: 14px;"><br></span></div><div><span style="font-size: 14px;">     nonloc. stress (kbar)   7797.79     -0.00      0.00</span></div><div><span style="font-size: 14px;">--</span></div><div><span style="font-size: 14px;">     Computing stress (Cartesian axis) and pressure</span></div><div><span style="font-size: 14px;"><br></span></div><div><span style="font-size: 14px;">          total   stress  (Ry/bohr**3)                   (kbar)     P=       -0.83</span></div><div><span style="font-size: 14px;">  -0.00000653  -0.00000000  -0.00000000           -0.96       -0.00       -0.00</span></div><div><span style="font-size: 14px;">   0.00000000  -0.00000653   0.00000000            0.00       -0.96        0.00</span></div><div><span style="font-size: 14px;">   0.00000000   0.00000000  -0.00000384            0.00        0.00       -0.56</span></div><div><span style="font-size: 14px;"><br></span></div><div><span style="font-size: 14px;">     kinetic stress (kbar)  26883.07      0.00     -0.00</span></div><div><span style="font-size: 14px;">                                0.00  26883.07     -0.00</span></div><div><span style="font-size: 14px;">                               -0.00     -0.00  26802.38</span></div><div><span style="font-size: 14px;"><br></span></div><div><span style="font-size: 14px;">     local   stress (kbar)-103303.36      0.03      0.00</span></div><div><span style="font-size: 14px;">                                0.03-103303.33      0.00</span></div><div><span style="font-size: 14px;">                                0.00      0.00  89345.22</span></div><div><span style="font-size: 14px;"><br></span></div><div><span style="font-size: 14px;">     nonloc. stress (kbar)   7795.83     -0.00      0.00</span></div>