Hi there everyone,<div><br></div><div>First time posting here, thank you in advance! </div><div>I am trying to do a structural optimization of FAU zeolite, which has a rhombohedral unit cell (a = b = c = 17.9 A, alpha = beta = gamma = 60 degrees). </div><div><br></div><div>I tried running the input file attached, but unfortunately, the crystallinity of the structure collapsed. </div><div><br></div><div><img alt="structure_collapse.png" width="423px" height="194.073px" data-atf="true" src="cid:31987916-0f0a-4849-99b1-edd75818c3ef"><br></div><div><br></div><div>I believe that the reason is related to the line in my input file "SYMMETRY RHOMBOHEDRAL", since the structure is preserved if this line is excluded. However, without SYMMETRY RHOMBOHEDRAL, the a, b and c parameters of the cell deviate from 17.9 A. I looked through the CP2K manual, but was unfortunately unable to find any information which might help solve the problem.<br></div><div><br></div><div>A part of my input file which is attached:</div><div><div>&CELL</div><div>    &CELL_REF</div><div>      A     12.6592017  12.6592017   0.0000000 </div><div>      B     12.6592017   0.0000000  12.6592017 </div><div>      C     0.0000000  12.6592017  12.6592017</div><div>      MULTIPLE_UNIT_CELL  1 1 1</div><div>      ALPHA_BETA_GAMMA  60 60 60</div><div>      SYMMETRY  RHOMBOHEDRAL</div><div>    &END CELL_REF</div><div>    A     12.6592017  12.6592017   0.0000000 </div><div>    B     12.6592017   0.0000000  12.6592017 </div><div>    C     0.0000000  12.6592017  12.6592017</div><div>    MULTIPLE_UNIT_CELL  1 1 1</div><div>    ALPHA_BETA_GAMMA  60 60 60</div><div>    SYMMETRY  RHOMBOHEDRAL</div><div>     &END CELL</div></div><div><br></div><div>Would anyone happen to know of a way to do a structural optimisation for a rhombohedral system while preserving the a, b and c parameters and without the structure falling apart?</div><div><br></div><div>Thank you in advance,</div><div>Tomoki</div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/ce488bba-4350-4be6-9849-3d52ff083bcbn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/ce488bba-4350-4be6-9849-3d52ff083bcbn%40googlegroups.com</a>.<br />