Dear users and developers of CP2K,<div><br></div><div>I tried starting the simulation of the crystal structure using cp2k. To do that, I used an input file:</div><div><br></div><div>&GLOBAL<br>  PROJECT crys_pbc_opt<br>  RUN_TYPE GEO_OPT <br>  PRINT_LEVEL MEDIUM<br>  EXTENDED_FFT_LENGTHS .TRUE.<br>&END GLOBAL<br><br>&MOTION <br>&GEO_OPT<br>  <br>  ! KEEP_SPACE_GROUP .TRUE.<br>  OPTIMIZER LBFGS    <br><br>&END GEO_OPT<br>&END MOTION <br><br>&FORCE_EVAL<br>  METHOD Quickstep              ! Electronic structure method (DFT,...)<br>  &DFT<br>    BASIS_SET_FILE_NAME  BASIS_MOLOPT<br>    POTENTIAL_FILE_NAME  POTENTIAL<br><br>    &POISSON                    ! Solver requested for non periodic calculations<br>      PERIODIC XYZ<br>    &END POISSON<br>    &SCF                        ! Parameters controlling the convergence of the scf. This section should not be changed. <br>      SCF_GUESS ATOMIC<br>      EPS_SCF 1.0E-6<br>      MAX_SCF 300<br>    &END SCF<br>    &XC                        ! Parameters needed to compute the electronic exchange potential <br>      &XC_FUNCTIONAL PBE<br>      &END XC_FUNCTIONAL<br>    &END XC<br>  &END DFT<br><br>  &SUBSYS<br>    &CELL<br>      CELL_FILE_FORMAT CIF<br>      ABC 9.2568 140.3208 18.9674<br>      PERIODIC XYZ<br>    &END CELL<br>    &TOPOLOGY                    ! Section used to center the atomic coordinates in the given box. Useful for big molecules<br>      COORD_FILE_FORMAT cif<br>      COORD_FILE_NAME  ./cell.cif<br>      ! &CENTER_COORDINATES<br>        ! CENTER_POINT 5. 5. 5.<br>      ! &END<br>    &END<br>    &KIND H<br>      ELEMENT H<br>      BASIS_SET DZVP-MOLOPT-GTH<br>      POTENTIAL GTH-PBE-q1<br>    &END KIND<br>    &KIND O<br>      ELEMENT O<br>      BASIS_SET DZVP-MOLOPT-GTH<br>      POTENTIAL GTH-PBE-q6<br>    &END KIND<br>    &KIND C<br>      ELEMENT C<br>      BASIS_SET DZVP-MOLOPT-GTH<br>      POTENTIAL GTH-PBE-q4<br>    &END KIND<br>    &KIND N<br>      ELEMENT N<br>      BASIS_SET DZVP-MOLOPT-GTH<br>      POTENTIAL GTH-PBE-q5<br>    &END KIND<br>  &END SUBSYS<br>&END FORCE_EVAL<br></div><div><br></div><div> ##############################################<br><br></div><div>The simulation starts, but only one step of SCF wavefunction optimization is passed, and then the calculation stops without errors. </div><div><br></div><div> Re-scaling the density matrix to get the right number of electrons<br>                  # Electrons              Trace(P)               Scaling factor<br>                         8576              8575.456                        1.000<br><br><br> SCF WAVEFUNCTION OPTIMIZATION<br><br>  Step     Update method      Time    Convergence         Total energy    Change<br>  ------------------------------------------------------------------------------<br>     1 P_Mix/Diag. 0.40E+00 ******     1.37713667    -13502.7304100109 -1.35E+04<br></div><div><br></div><div>How do you think, what happens? How could I change the input file to solve this problem?</div><div><br></div><div>Best regards, </div><div>Victor</div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">суббота, 20 августа 2022 г. в 17:41:05 UTC+3, Victor Nazarychev: <br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><p>Dear users and developers of CP2K,<br><br>I am a new user of CP2K. I want to use CP2K to simulate the DFT of a polymer crystal, which contains ~2500 atoms (C, N, O, H). This initial crystal conformation was suggested by the experimenters, and I need to optimize this structure to find the minimal energy state.<br><br>I'm considering using CP2K due to its ability to do periodic DFT calculations. My question is related to what DFT functional and basis sets can be used to optimize a periodic crystal structure if such a massive system is considered? Maybe some there exists some examples of an input file for the problem of optimizing crystal systems using periodic boundary conditions?<br><br>For these calculations, I have a workstation with an AMD Ryzen 5950x and 64 GB of memory. Also, is it possible to speed up the calculation with the GPU?<br><br>Best wishes,<br><br>Victor<br></p></blockquote></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/a2a11a0e-43d2-4110-986d-4598dc1b0f3dn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/a2a11a0e-43d2-4110-986d-4598dc1b0f3dn%40googlegroups.com</a>.<br />