Dear all,<div><br /></div><div>I have been trying to optimize a unit cell of the 2D semiconductor CrSBr using GPW method implemented in CP2K. The initial cell parameters were a = 3.6 A, b = 4.8 A and c = 25.0 A (obtained from previous calculations with Quantum Espresso and in good agreement with experimental data). After 128 optimization cycles the geometry is optimized and the obtained cell parameters were a = 3.64 A and b = 6.60 A, which leads to a meaningless structure since there is no longer a 2D structure. Also the geometry of the unit cell is really distorted. Any suggestions to improve the accuracy of my calculation? Thank you in advanced.</div><div><br /></div><div>I used this input:</div><div><br /></div><div>&GLOBAL<br />  PROJECT test2<br />  RUN_TYPE CELL_OPT<br />  PRINT_LEVEL MEDIUM<br />&END GLOBAL<br /><br />&FORCE_EVAL<br />  METHOD Quickstep              ! Electronic structure method (DFT,...)<br />  STRESS_TENSOR ANALYTICAL<br />  &DFT<br />    BASIS_SET_FILE_NAME  /BASIS_MOLOPT<br />    POTENTIAL_FILE_NAME  /GTH_POTENTIALS<br />    MULTIPLICITY 7<br />    &POISSON                    ! Solver requested for non periodic calculations<br />      PERIODIC XYZ<br />    &END POISSON<br />     &XC<br />      &XC_FUNCTIONAL PBE<br />      &END XC_FUNCTIONAL<br />     &END XC<br />    &SCF<br />      SCF_GUESS SPARSE<br />      MAX_SCF 100<br />      &OT<br />      PRECONDITIONER FULL_KINETIC<br />      MINIMIZER CG<br />      &END OT<br />      &OUTER_SCF<br />       MAX_SCF 50<br />      &END<br />    &END SCF<br />    &MGRID<br />      NGRIDS 4<br />      CUTOFF 550<br />      REL_CUTOFF 90<br />    &END MGRID<br />    &QS<br />       METHOD GPW<br />    &END<br />  &END DFT<br />  &SUBSYS<br />    &CELL<br />      ABC 3.601367151 4.818258107 25.0<br />      ALPHA_BETA_GAMMA 90 90 90<br />      PERIODIC XYZ<br />      SYMMETRY ORTHORHOMBIC<br />    &END CELL<br />    &TOPOLOGY                    ! Section used to center the atomic coordinates in the given box. Useful for big molecules<br />      COORD_FILE_FORMAT xyz<br />      COORD_FILE_NAME  ./crsbr.xyz<br />    &END<br />    &KIND Cr<br />      ELEMENT Cr<br />      BASIS_SET DZVP-MOLOPT-SR-GTH<br />      POTENTIAL GTH-PBE-q14<br />      MAGNETIZATION 3.0<br />    &END KIND<br />    &KIND S<br />      ELEMENT S<br />      BASIS_SET TZV2P-MOLOPT-GTH<br />      POTENTIAL GTH-PBE-q6<br />    &END KIND<br />    &KIND Br<br />      ELEMENT Br<br />      BASIS_SET DZVP-MOLOPT-SR-GTH<br />      POTENTIAL GTH-PBE-q7<br />    &END KIND<br />  &END SUBSYS<br />&END FORCE_EVAL<br />&MOTION<br />  &CELL_OPT<br />    TYPE DIRECT_CELL_OPT<br />    MAX_DR    1.0E-03<br />    MAX_FORCE 1.0E-03<br />    RMS_DR    1.0E-03<br />    RMS_FORCE 1.0E-03<br />    MAX_ITER 200<br />    OPTIMIZER BFGS<br />    KEEP_SYMMETRY<br />  &END CELL_OPT<br /> &GEO_OPT<br />    MAX_DR    1.0E-03<br />    MAX_FORCE 1.0E-03<br />    RMS_DR    1.0E-03<br />    RMS_FORCE 1.0E-03<br />    OPTIMIZER BFGS<br />  &END GEO_OPT<br />&END MOTION<br /></div><div><br /></div><div><br /></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/ffb8824d-4222-4ba7-8a85-f0a99c49ee46n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/ffb8824d-4222-4ba7-8a85-f0a99c49ee46n%40googlegroups.com</a>.<br />