Dear pwscf users:<br>&nbsp;&nbsp;&nbsp;&nbsp; I checked VCSexample in V4.0, and found a disagreement in the *.in and *.out. For example:<br>In the <a href="http://As.vcs00.in">As.vcs00.in</a> file,the system is defined as follows:<br>++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++<br>
+&nbsp;&nbsp; ibrav = 0 ,<br>+&nbsp;&nbsp; A = 3.70971016 ,<br>+&nbsp;&nbsp; B = 3.70971016 ,<br>+&nbsp;&nbsp; C = 3.70971016 ,<br>+&nbsp;&nbsp; cosAB = 0.49517470 ,<br>+&nbsp;&nbsp; cosAC = 0.49517470 ,<br>+&nbsp;&nbsp; cosBC = 0.49517470 ,<br>+&nbsp;&nbsp; <br>+&nbsp;&nbsp; CELL_PARAMETERS cubic <br>+&nbsp;&nbsp;&nbsp; 0.58012956&nbsp; 0.00000000&nbsp; 0.81452422<br>
+&nbsp;&nbsp; -0.29006459&nbsp; 0.50240689&nbsp; 0.81452422<br>+&nbsp;&nbsp; -0.29006459 -0.50240689&nbsp; 0.81452422<br>++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++&nbsp;&nbsp; <br>In the beginning of As.vcs00.out file,the system is changed.<br>
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++<br>+&nbsp;&nbsp;&nbsp;&nbsp; bravais-lattice index&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0<br>+&nbsp;&nbsp;&nbsp;&nbsp; lattice parameter (a_0)&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 7.0103&nbsp; a.u.<br>+&nbsp;&nbsp;&nbsp;&nbsp; unit-cell volume&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp; 245.3705 (a.u.)^3<br>
+&nbsp;&nbsp;&nbsp;&nbsp; <br>+&nbsp;&nbsp;&nbsp;&nbsp; celldm(1)=&nbsp;&nbsp; 7.010336&nbsp; celldm(2)=&nbsp;&nbsp; 1.000000&nbsp; celldm(3)=&nbsp;&nbsp; 1.000000<br>+&nbsp;&nbsp;&nbsp;&nbsp; celldm(4)=&nbsp;&nbsp; 0.495175&nbsp; celldm(5)=&nbsp;&nbsp; 0.000000&nbsp; celldm(6)=&nbsp;&nbsp; 0.000000<br>+<br>+&nbsp;&nbsp;&nbsp;&nbsp; crystal axes: (cart. coord. in units of a_0)<br>
+&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; a(1) = (&nbsp; 0.580130&nbsp; 0.000000&nbsp; 0.814524 )&nbsp; <br>+&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; a(2) = ( -0.290065&nbsp; 0.502407&nbsp; 0.814524 )&nbsp; <br>+&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; a(3) = ( -0.290065 -0.502407&nbsp; 0.814524 )<br>+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
a=3.70971016A -&gt; celldm(1)=7.010336a.u. and celldm(2)=celldm(3)=1<br>celldm(4)=cos(bc)=0.495175<br>but why celldm(5),celldm(6)are equal to 0?<br clear="all"><br>-- <br>Y. C. Cheng<br>Department of Phyics<br>Nanjing University<br>
Nanjing 210093<br>P. R. China<br>Tel: 86-25-83592907<br>Email: <a href="mailto:yccheng.nju@gmail.com">yccheng.nju@gmail.com</a>