<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Hi,<div><br></div><div>I am a chemical crystallographer attempting to better understand and test usage of "supposedly straightforward" QE vc-relax in optimizing cell constants.  I am using form I of paracetamol as a test case. The starting structure is from the 20K neutron data.   For this test calculation I am actually following guidelines which were outlined in an article on the subject (<span style="color:rgb(117,117,117);line-height:1.5">J. Chem. Theory Comput. 2014, 10, 3423−3437</span><span style="line-height:1.5">), the dispersion correction is vdW-DF2.</span></div><div><span style="line-height:1.5"><br></span></div><div><span style="line-height:1.5">the experimental lattice parameters are as follows :</span></div><div><span style="line-height:1.5"><br></span></div><div><div>_cell_length_a                   12.667(4)</div><div>_cell_length_b                   9.166(3)</div><div>_cell_length_c                   7.073(3)</div><div>_cell_angle_alpha                90</div><div>_cell_angle_beta                 115.51(2)</div><div>_cell_angle_gamma                90</div><div>_cell_volume                     741.156</div></div><div><br></div><div>the output form the optimization is</div><div><br></div><div><div>_cell_length_a           12.994686</div><div>_cell_length_b           9.115636</div><div>_cell_length_c           6.461433</div><div>_cell_angle_alpha        90.000000</div><div>_cell_angle_beta         115.510024</div><div>_cell_angle_gamma        90.000000</div></div><div><br></div><div>which is strange because the angle beta is unchanged and the a and c lengths have become inaccurate. I believe the deviation is too high. </div><div><span style="line-height:1.5"><br></span></div><div><span style="line-height:1.5">attached is the .in file</span></div><div><span style="line-height:1.5">I have  missed something, and cannot figure out what that is? </span></div><div><span style="line-height:1.5">can anyone help me master the art of predicting reliable cell constants using  QE </span><span style="line-height:1.5">vc-relax?</span></div><div><br></div><div>any help is greatly appreciated. </div><div><br></div><div>-Eric</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div>







</div></div></div></div></div>