<div dir="ltr"><div>A quick glance at the code seems to indicate that cell_dofree='ibrav+something'  is implemented only for BFGS minimization, not for variable-cell molecular dynamics</div><div><br></div><div>Paolo<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Mar 7, 2022 at 5:14 PM Tianhua Wang via users <<a href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div>Dear Quantum ESPRESSO developers,</div><div>  I am doing the cell optimization for cassiterite (SnO<font size="1">2</font>) using the variable cell shape molecular dynamics (i.e., <i>cell_dynamics='damp-w'</i>) in QE 7.0. I had set <i>ibrav=6</i> for the crystal of cassiterite, and wanted to constrain <i>cell_dofree='ibrav+shape'</i> during the optimization. However, I found the none-zero values of v1_y&z, v2_x&z, v3_x&y in the output file, e.g.,</div><div><br></div><div><div><i>CELL_PARAMETERS (alat=  8.95238858)</i></div><div><i>   0.999715701   <u>0.000000089</u>  <u>-0.000000001</u></i></div><div><i>   <u>0.000000089</u>   0.999715701  <u>-0.000000001</u></i></div><div><i>  <u>-0.000000001</u>  <u>-0.000000001</u>   0.672987836</i></div></div><div><br></div><div>which was not absolutely consistent with the initial ibrav choice (i.e., <i>ibrav=6</i>). Then I tested the performance of <i>cell_dofree='ibrav'</i>, and found the same phenomenon. It seems that the <i>cell_dofree</i> command does not strictly constrain the initial ibrav choice. I also tested <i>cell_dofree='xyz+shape'</i> although I had guessed that it did not work.</div><div>  I am wondering if there is any way to strictly achieve <i>cell_dofree='ibrav+shape'</i> (which is equal to <i>cell_dofree='xyz+shape'</i> in this case) during the cell optimization. Perhaps another solution is to fix the angles during the cell optimization, but how to achieve? My input file is attached for reference.</div><div>  Thank you for your time and I am looking forward to your reply.</div><div>Sincerely,</div><div>Tianhua Wang</div></div><div><br></div><div><hr style="margin:0px 0px 10px;border-color:currentcolor currentcolor rgb(228,229,230);border-style:none none solid;border-width:0px 0px 1px;height:0px;line-height:0;font-size:0px;padding:20px 0px 0px;width:50px" align="left"><div style="font-size:14px;font-family:Verdana;color:rgb(0,0,0)"><div>Tianhua Wang</div><div>State Key Laboratory for Mineral Deposits Research</div><div>School of Earth Sciences and Engineering</div><div>Nanjing University</div><div>163 Xianlin Avenue, Nanjing</div><div>Jiangsu 210023, China</div><div><br></div><div><br></div></div></div>_______________________________________________<br>
The Quantum ESPRESSO community stands by the Ukrainian<br>
people and expresses its concerns about the devastating<br>
effects that the Russian military offensive has on their<br>
country and on the free and peaceful scientific, cultural,<br>
and economic cooperation amongst peoples<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 206, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>