<div>Dear all,</div>  <div>After many vc-relax calculations(first ibrav=1,then using ibrav=0, new cell_parameters,  new positions), I obtained <BR>the following results. It is obvious the cell angle parameters have been changed and the structure isn't the cubic structure.</div>  <div><BR> 1. If I still set ibrav=1 in the next scf calculation, how to read the new lattice constant? <BR>I have attempted to use celldm(1)*(a11+a22+a33)/3 as the new celldm(1),<BR>(axx is the diagonal element of cell_parameters) ,the output atomic_positions(alat) as the new atomic_positions,and neglect all of the nondiagonal elements.<BR>After the first self-consistent,I gain stresses with very large value. How can I adapt the input parameters to avoid this condition and get smaller elements?</div>  <div><BR>2.Must the k-points are same in the vc-relax and scf?</div>  <div> </div>  <div>any help will be appreciated!</div>  <div><BR>Best!</div>  <div> </div>  <div>Niu
 Li<BR>Harbin Institute of Technology<BR>China</div>  <div>==================================================<BR>!    total energy              =  -729.30084995 Ry<BR>     Harris-Foulkes estimate   =  -729.30084995 Ry<BR>     estimated scf accuracy    <        8.1E-10 Ry</div>  <div>     convergence has been achieved</div>  <div>     Forces acting on atoms (Ry/au):</div>  <div>     atom   1 type  1   force =     0.00000028   -0.00000032   -0.00000017<BR>          ...<BR>     atom  62 type  1   force =    -0.00000135    0.00000206   
 0.00000190<BR>     atom  63 type  1   force =    -0.00000012   -0.00000115    0.00000144<BR>     atom  64 type  1   force =    -0.00000077   -0.00000041   -0.00000140</div>  <div>     Total force =     0.000024     Total SCF correction =     0.000091</div>  <div><BR>     entering subroutine stress ...</div>  <div>          total   stress  (Ry/bohr**3)                   (kbar)     P=    0.06<BR>   0.00000025   0.00000101   0.00000074         
 0.04      0.15      0.11<BR>   0.00000101   0.00000338  -0.00000180          0.15      0.50     -0.27<BR>   0.00000074  -0.00000180  -0.00000237          0.11     -0.27     -0.35</div>  <div><BR> Wentzcovitch Damped Dynamics: convergence achieved, Efinal=  -729.30084995</div>  <div>------------------------------------------------------------------------</div>  <div>     Final estimate of lattice vectors (input alat units)<BR>   0.988103212   0.021498496  -0.002655708<BR>   0.021493430   1.013103454  -0.000556700<BR>  -0.002655491  -0.000556401   0.983829644<BR>  final unit-cell volume =  
 2649.6086 (a.u.)^3<BR>  input alat =      13.9103 (a.u.)</div>  <div>CELL_PARAMETERS (alat)<BR>   0.988103212   0.021498496  -0.002655708<BR>   0.021493430   1.013103454  -0.000556700<BR>  -0.002655491  -0.000556401   0.983829644</div>  <div>ATOMIC_POSITIONS (angstrom)<BR>C        0.292559767  -0.288781798  -0.891844186<BR>C        1.355532741  -0.743068572   3.628414295<BR>...<BR>C        3.388679954   3.539286798   4.255326154</div>  <div>=========================================</div><p>
                <hr size=1><a href="http://cn.mail.yahoo.com" target=blank>ÇÀ×¢ÑÅ»¢Ãâ·ÑÓÊÏä-3.5GÈÝÁ¿£¬20M¸½¼þ£¡</a>