<DIV>&nbsp;Dear pwscf users</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; I run&nbsp;vc-relax to&nbsp;&nbsp;optimize the lattice constant of ZnO unit cell,my input file and out file as follows:</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </DIV>
<DIV>&nbsp;&amp;CONTROL<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; title = ZnOVC ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; calculation = 'vc-relax' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; restart_mode = 'from_scratch' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pseudo_dir = '/root/espresso-3.2.3/pseudo/' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; prefix = ZnO ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; etot_conv_thr = 1.D-5 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nstep = 100 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tstress = .true. ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tprnfor = .true. ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; dt = 100 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tefield = .false. ,<BR>&nbsp;/<BR>&nbsp;&amp;SYSTEM<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ibrav = 4,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; celldm(1) = 6.14161,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; celldm(3) = 1.6022,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nat = 4,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ntyp = 2,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutwfc = 40 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutrho = 300 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nosym = .true. ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nbnd = 25,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; occupations = 'smearing' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; degauss = 0.002 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; smearing = 'methfessel-paxton' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; lda_plus_u = .false. ,<BR>&nbsp;/<BR>&nbsp;&amp;ELECTRONS<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; conv_thr = 1.D-8 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; startingpot = 'atomic' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; startingwfc = 'atomic' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mixing_mode = 'plain' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mixing_beta = 0.6 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; diagonalization = 'david' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; diago_full_acc = .true. ,<BR>&nbsp;/<BR>&nbsp;&amp;IONS<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ion_dynamics = 'damp' ,<BR>&nbsp;/<BR>&nbsp;&amp;CELL<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; cell_dynamics = 'damp-w' ,<BR>&nbsp;/<BR>ATOMIC_SPECIES<BR>&nbsp;&nbsp; Zn&nbsp;&nbsp; 65.00000&nbsp; Zn.pz-van_ak.UPF <BR>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp; 16.00000&nbsp; O.pz-rrkjus.UPF <BR>ATOMIC_POSITIONS crystal <BR>&nbsp;&nbsp; Zn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.333333333&nbsp;&nbsp;&nbsp; 0.666666667&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1 <BR>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.333333333&nbsp;&nbsp;&nbsp; 0.666666667&nbsp;&nbsp;&nbsp; 0.381700000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1 <BR>&nbsp;&nbsp; Zn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.666666667&nbsp;&nbsp;&nbsp; 0.333333333&nbsp;&nbsp;&nbsp; 0.500000000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1 <BR>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.666666667&nbsp;&nbsp;&nbsp; 0.333333300&nbsp;&nbsp;&nbsp; 0.881700000&nbsp;&nbsp;&nbsp; 1&nbsp; 1&nbsp; 1 <BR>K_POINTS automatic <BR>&nbsp; 4 4 4&nbsp;&nbsp; 0 0 0 </DIV>
<DIV>&nbsp;</DIV>
<DIV>--------------------------------------------------------------------------------</DIV>
<DIV>&nbsp; Program PWSCF&nbsp;&nbsp;&nbsp;&nbsp; v.3.2.2&nbsp; starts ...<BR>&nbsp;&nbsp;&nbsp;&nbsp; Today is 18Aug2008 at 23: 7:36 </DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp; Ultrasoft (Vanderbilt) Pseudopotentials</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp; Current dimensions of program pwscf are:</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp; ntypx = 10&nbsp;&nbsp; npk = 40000&nbsp; lmax =&nbsp; 3<BR>&nbsp;&nbsp;&nbsp;&nbsp; nchix =&nbsp; 6&nbsp; ndmx =&nbsp; 2000&nbsp; nbrx = 14&nbsp; nqfx =&nbsp; 8</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp; Title: <BR>&nbsp;&nbsp;&nbsp;&nbsp; ZnOVC&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </DIV>
<DIV><BR>&nbsp;&nbsp;&nbsp;&nbsp; bravais-lattice index&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4<BR>&nbsp;&nbsp;&nbsp;&nbsp; lattice parameter (a_0)&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 6.1416&nbsp; a.u.<BR>&nbsp;&nbsp;&nbsp;&nbsp; unit-cell volume&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp; 321.4357 (a.u.)^3<BR>&nbsp;&nbsp;&nbsp;&nbsp; number of atoms/cell&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4<BR>&nbsp;&nbsp;&nbsp;&nbsp; number of atomic types&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2<BR>&nbsp;&nbsp;&nbsp;&nbsp; kinetic-energy cutoff&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 40.0000&nbsp; Ry<BR>&nbsp;&nbsp;&nbsp;&nbsp; charge density cutoff&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp; 300.0000&nbsp; Ry<BR>&nbsp;&nbsp;&nbsp;&nbsp; convergence threshold&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.0E-08<BR>&nbsp;&nbsp;&nbsp;&nbsp; beta&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.6000<BR>&nbsp;&nbsp;&nbsp;&nbsp; number of iterations used =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 8&nbsp; plain&nbsp;&nbsp;&nbsp;&nbsp; mixing<BR>&nbsp;&nbsp;&nbsp;&nbsp; Exchange-correlation&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; SLA&nbsp; PW&nbsp;&nbsp; PBE&nbsp; PBE (1434)<BR>&nbsp;&nbsp;&nbsp;&nbsp; nstep&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 100</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp; celldm(1)=&nbsp;&nbsp; 6.141610&nbsp; celldm(2)=&nbsp;&nbsp; 0.000000&nbsp; celldm(3)=&nbsp;&nbsp; 1.602200<BR>&nbsp;&nbsp;&nbsp;&nbsp; celldm(4)=&nbsp;&nbsp; 0.000000&nbsp; celldm(5)=&nbsp;&nbsp; 0.000000&nbsp; celldm(6)=&nbsp;&nbsp; 0.000000</DIV>
<DIV>&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; 1.000000&nbsp; 0.000000&nbsp; 0.000000 )&nbsp; <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; a(2) = ( -0.500000&nbsp; 0.866025&nbsp; 0.000000 )&nbsp; <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; a(3) = (&nbsp; 0.000000&nbsp; 0.000000&nbsp; 1.602200 )&nbsp; </DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp; reciprocal axes: (cart. coord. in units 2 pi/a_0)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; b(1) = (&nbsp; 1.000000&nbsp; 0.577350&nbsp; 0.000000 )&nbsp; <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; b(2) = (&nbsp; 0.000000&nbsp; 1.154701&nbsp; 0.000000 )&nbsp; <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; b(3) = (&nbsp; 0.000000&nbsp; 0.000000&nbsp; 0.624142 )&nbsp; </DIV>
<DIV><BR>.......</DIV>
<DIV>.......</DIV>
<DIV>outfile:</DIV>
<DIV><BR>&nbsp;&nbsp;&nbsp;&nbsp; Final estimate of lattice vectors (input alat units)<BR>&nbsp;&nbsp; 1.011331020&nbsp;&nbsp; 0.000000000&nbsp;&nbsp; 0.000000000<BR>&nbsp; -0.505665507&nbsp;&nbsp; 0.875838356&nbsp;&nbsp; 0.000000001<BR>&nbsp;&nbsp; 0.000000000&nbsp; -0.000000002&nbsp;&nbsp; 1.627761723<BR>&nbsp; final unit-cell volume =&nbsp;&nbsp;&nbsp; 334.0064 (a.u.)^3<BR>&nbsp; input alat =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 6.1416 (a.u.)</DIV>
<DIV>CELL_PARAMETERS (alat)<BR>&nbsp;&nbsp; 1.011331020&nbsp;&nbsp; 0.000000000&nbsp;&nbsp; 0.000000000<BR>&nbsp; -0.505665507&nbsp;&nbsp; 0.875838356&nbsp;&nbsp; 0.000000001<BR>&nbsp;&nbsp; 0.000000000&nbsp; -0.000000002&nbsp;&nbsp; 1.627761723</DIV>
<DIV>ATOMIC_POSITIONS (crystal)<BR>Zn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.333333347&nbsp;&nbsp; 0.666666677&nbsp;&nbsp; 0.000438189<BR>O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.333333356&nbsp;&nbsp; 0.666666677&nbsp;&nbsp; 0.380132292<BR>Zn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.666666690&nbsp;&nbsp; 0.333333345&nbsp;&nbsp; 0.500438960<BR>O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.666666695&nbsp;&nbsp; 0.333333331&nbsp;&nbsp; 0.880132483</DIV>
<DIV>&nbsp;</DIV>
<DIV>from the outfile I can see the c/a=1.627761723 and from the final unit-cell volume =334.0064 (a.u.)^3,I can compute the a=6.189678(Bohr)=3.2754(A)¡£From the value of caculating ,we can see that the c/a and a are both bigger than the experiment&nbsp;,&nbsp;I don't konw whether the result is right or wrong&nbsp;.</DIV>
<DIV>&nbsp;</DIV>
<DIV>Anybody who know about the question and give me the useful advice will be appreciated .</DIV>
<DIV>&nbsp;</DIV>
<DIV>Sincerely </DIV>
<DIV>&nbsp;</DIV>
<DIV><BR>&nbsp;</DIV>
<DIV>&nbsp;</DIV>