<div>Dear Pwscf users,</div>
<div>        I did a variable cell calculations for a base centered monoclinic crystal.  The cell displayed in xcrysden is correct, and PWscf can find the symmetry.  However, in the optimizing process, the cell parameter changed unexpectedly, which may be wrong.  In the user's guide of PWscf, the base centered monoclinic cell is defined as follows:</div>

<div>base centered monoclinic<br>=============================<br>   v1 = (  a/2,         0,                -c/2),<br>   v2 = (b*cos(gamma), b*sin(gamma), 0),<br>   v3 = (  a/2,         0,                  c/2),<br>where gamma is the angle between axis a and b<br>
=============================</div>
<div>In the first step of optimize process, the cell is changed as following:</div>
<div>CELL_PARAMETERS (alat)<br>   0.483414893   0.003687739  -0.271981740<br>  -1.013180906   0.556415406   0.000000000<br>   0.483414893   0.003687739   0.271981740</div>
<div>My question is : why the quantities of second row of v1, v3 vector do not retain zero!</div>
<div> </div>
<div>Any suggestion and comment is appriciated. Thank you in advance!</div>
<div> </div>
<div>The input file is attached!</div>
<div> &CONTROL<br>                       title = 'c2',<br>                      prefix = 'li',<br>                 calculation = 'vc-relax' ,<br>                restart_mode = 'from_scratch' ,<br>
                      outdir = './' ,<br>                      wfcdir = './' ,<br>                  pseudo_dir = './' ,<br>                     tstress = .true. ,<br>                     tprnfor = .true. ,<br>
               etot_conv_thr = 1.0D-5,<br>               forc_conv_thr = 1.0D-4,<br> /<br> &SYSTEM<br>                       ibrav = 13,<br>                   celldm(1) = 14.7756,<br>                   celldm(2) = 1.20528,<br>
                   celldm(3) = 0.56348,<br>                   celldm(4) = -0.873177,<br>                         nat = 12,<br>                        ntyp = 1,<br>                     ecutwfc = 200,<br>                       nosym = .false. ,<br>
 /<br> &ELECTRONS<br>            electron_maxstep = 100,<br>                    conv_thr = 1.D-10,<br> /<br> &IONS<br> ion_dynamics='bfgs',<br>/<br> &CELL<br> cell_dynamics='bfgs',<br> press=900,<br>
/<br>ATOMIC_SPECIES<br>   Li   6.914  Li.pz-s-mt.UPF<br>ATOMIC_POSITIONS (crystal)<br>Li   -0.46080   0.21740   1.49640<br>Li   -0.02770   0.02590   1.28350<br>Li   -0.80760   0.63690   0.91840<br>Li    0.23010   0.43900   0.29270<br>
Li    0.06870   0.87600   0.40410<br>Li   -0.04760   0.70260   1.56000<br>Li   -1.49640  -0.21740   0.46080<br>Li   -1.28350  -0.02590   0.02770<br>Li   -0.91840  -0.63690   0.80760<br>Li   -0.29270  -0.43900  -0.23010<br>
Li   -0.40410  -0.87600  -0.06870<br>Li   -1.56000  -0.70260   0.04760<br>K_POINTS automatic<br>  8 8 8   1 1 1<br></div>
<div>-- <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><br></div>