<div dir="ltr">Thanks Matt,<br><br>Its working correctly now.<br><br>On Tuesday, April 17, 2018 at 4:34:13 PM UTC-7, Matt W wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hi,<div><br></div><div>you should either specify </div><div><br></div><div>ABC and alpha, beta, gamma</div><div><br></div><div>or</div><div><br></div><div>give the whole set of vectors, not both.</div><div><br></div><div>If the set of A, B, C vectors are given any angles are ignored.</div><div><br></div><div>Also, I suspect  your supercell likely isn't large enough. By default cp2k only uses gamma point to sample the Brillouin zone. You should probably use a 2x2x1 supercell to get reasonably accurate results.</div><div><br></div><div>Matt</div><div><br>On Wednesday, April 18, 2018 at 12:03:43 AM UTC+1, Aniruddha Dive wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi All,<br><br>I am trying to do cell optimization of a SiS2 crystal structure which has a triclinic symmetry. Below is the section of input file where I specify the cell vectors and the angles:<br><br> &SUBSYS<br>    &CELL<br>      <b>SYMMETRY TRICLINIC<br>      A 5.58515   0.00000   0.00000<br>      B 0.00000   6.21427   0.00000<br>      C 0.00000   0.00000  10.92676<br>      ALPHA_BETA_GAMMA 116.86 113.96 90.00<br>      PERIODIC XYZ</b><br>    &END<br>    ! unit cells that are orthorhombic are more efficient with CP2K<br>    ! atom coordinates can be in the &COORD section,<br>    ! or provided as an external file.<br>    &TOPOLOGY<br>       COORD_FILE_FORMAT PDB<br>       COORD_FILE_NAME SiS2_computed.pdb<br>       CONNECTIVITY OFF<br>       &CENTER_COORDINATES<br>       &END CENTER_COORDINATES<br>    &END<br><br>When I try to run the calculations, in the output file I can see the correct vectors a,b and c but the angles are still 90.0 90.0 and 90.0. Below is the portion from output file:<br><br> <b>CELL_TOP| Volume [angstrom^3]:                 <wbr>                         379.242<br> CELL_TOP| Vector a [angstrom     5.585     0.000     0.000    |a| =       5.585<br> CELL_TOP| Vector b [angstrom     0.000     6.214     0.000    |b| =       6.214<br> CELL_TOP| Vector c [angstrom     0.000     0.000    10.927    |c| =      10.927<br> CELL_TOP| Angle (b,c), alpha [degree]:                     <wbr>              90.000<br> CELL_TOP| Angle (a,c), beta  [degree]:                     <wbr>              90.000<br> CELL_TOP| Angle (a,b), gamma [degree]:                     <wbr>              90.000<br> CELL_TOP| Requested initial symmetry:                     <wbr>            TRICLINIC<br> CELL_TOP| Numerically orthorhombic:                 <wbr>                        YES<br><br> CELL| Volume [angstrom^3]:                 <wbr>                             379.242<br> CELL| Vector a [angstrom]:       5.585     0.000     0.000    |a| =       5.585<br> CELL| Vector b [angstrom]:       0.000     6.214     0.000    |b| =       6.214<br> CELL| Vector c [angstrom]:       0.000     0.000    10.927    |c| =      10.927<br> CELL| Angle (b,c), alpha [degree]:                     <wbr>                  90.000<br> CELL| Angle (a,c), beta  [degree]:                     <wbr>                  90.000<br> CELL| Angle (a,b), gamma [degree]:                     <wbr>                  90.000<br> CELL| Requested initial symmetry:                     <wbr>                TRICLINIC<br> CELL| Numerically orthorhombic:                 <wbr>                            YES<br><br> CELL_REF| Volume [angstrom^3]:                 <wbr>                         379.242<br> CELL_REF| Vector a [angstrom     5.585     0.000     0.000    |a| =       5.585<br> CELL_REF| Vector b [angstrom     0.000     6.214     0.000    |b| =       6.214<br> CELL_REF| Vector c [angstrom     0.000     0.000    10.927    |c| =      10.927<br> CELL_REF| Angle (b,c), alpha [degree]:                     <wbr>              90.000<br> CELL_REF| Angle (a,c), beta  [degree]:                     <wbr>              90.000<br> CELL_REF| Angle (a,b), gamma [degree]:                     <wbr>              90.000<br> CELL_REF| Requested initial symmetry:                     <wbr>            TRICLINIC<br> CELL_REF| Numerically orthorhombic:                 <wbr>                        YES</b><br><br>Because of this the complete geometry of the cell has changed during the optimization giving me weird results. Please find attached input files as well as the coordinate file.<br><br>Kindly let me know how to fix the issue.<br><br>Regards,<br>Aniruddha M Dive<br>PhD Candidate, Washington State University<br><br><br></div></blockquote></div></div></blockquote></div>