<meta http-equiv="Content-Type" content="text/html; charset=GB18030"><div><div>Dear Quantum ESPRESSO developers,</div><div>  Below is my E-mail to the QE users forum. It can be seen that this letter is not a SPAM but a sincere consultation of a problem. It is very confused and depressed for me to see that this E-mail is marked as a SPAM in the Mail Archive of QE. Could you please check the E-mail again and remove this mark?</div></div><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><div><br></div><div><br></div><div style="font-size: 12px;font-family: Arial Narrow;padding:2px 0 2px 0;">------------------ Original ------------------</div><div style="font-size: 12px;background:#efefef;padding:8px;"><div><b>From:</b>                                                                                                                        "Tianhua Wang"                                                                                    <njuwth@qq.com>;</div><div><b>Date:</b> Tue, Mar 8, 2022 00:13 AM</div><div><b>To:</b> "users"<users@lists.quantum-espresso.org>;<wbr></div><div></div><div><b>Subject:</b> A problem about cell_dofree='ibrav'</div></div><div><br></div><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 align="left" style="margin: 0 0 10px 0;border: 0;border-bottom:1px solid #E4E5E6;height:0;line-height:0;font-size:0;padding: 20px 0 0 0;width: 50px;"><div style="font-size:14px;font-family:Verdana;color:#000;"><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></div>