Hello to all,<br><br>I have tried for the vc-relax of Mo2C hexagonal 
system, job is not done properly it come out in between and the atomic 
positions <br><br>what I have given in the input file and what I am getting in the out file is different. <br>
(i) Why the crystal positions changes. what I think that the nature of the atomic postions must be retained.<br><br>In Input file<br>ATOMIC_POSITIONS (alat)<br>
C        0.00  0.00 0.00<br>
Mo       0.3333333   0.666666   0.25<br>
Mo       0.6666666   0.333333   0.75<br>
<br><br><br>In out file we get<br><br>CELL_PARAMETERS (alat=  5.67296009)<br>   0.946518633   0.001215132  -0.004855406<br>  -0.474470641   0.819090796  -0.000958606<br>   0.009278926   0.007480907   1.741942673<br><br>ATOMIC_POSITIONS (alat)<br>

C        0.156590059   0.041446237  -0.212516996<br>Mo       0.157761886   0.589915381   0.269109338<br>Mo       0.636664528   0.321113609   1.041609035<br><br><br><br><br>(ii) why cell_parameter changes in the out file . (I think that for hexagonal system it must be what given in the input file)<br>

<br><br>my input file<br><br>&amp;control<br>    calculation=&#39;vc-relax&#39;<br>    restart_mode=&#39;from_scratch&#39;,<br>    pseudo_dir = &#39;/home/ashish/softwares/espresso/pseudo/&#39;,<br>    outdir=&#39;/home/ashish/dft/espresso_mo2c/test1/tmp_espresso/&#39;<br>
    prefix=&#39;mo2c&#39;<br>    tprnfor = .true.,<br>    tstress = .true.<br> /<br> &amp;system<br>    ibrav=4, celldm(1) =5.672960087 , celldm(3)=1.573617588 , nat=3, ntyp=2,<br>    nspin = 1,  starting_magnetization(1)=0.7,<br>
    ecutwfc = 24.0, ecutrho = 288.0, <br>    occupations=&#39;smearing&#39;, smearing=&#39;methfessel-paxton&#39;, degauss=0.02<br> /<br> &amp;electrons<br>    diagonalization=&#39;cg&#39;<br>    conv_thr = 1.0e-8<br>    mixing_beta = 0.7<br>
/<br>&amp;ions<br>  ion_dynamics=&#39;bfgs&#39;<br> /<br>&amp;cell<br>cell_dynamics=&#39;bfgs&#39;<br>/<br>ATOMIC_SPECIES {crystal}<br> C  12.01 C.pw91-van_ak.UPF  <br> Mo 95.94 Mo.pw91-n-van.UPF <br>ATOMIC_POSITIONS<br> C 0.0 0.0 0.0<br>
 Mo 1/3 2/3 1/4<br> Mo 2/3 1/3 3/4<br><br>K_POINTS automatic<br> 4 4 4   0 0 0 <br><br><br> can any one help me. What is the mistake (in my input file or else)<br><br><br>with best wishes<br><font color="#888888">Ashish<br>
<br></font>