Dear CP2K users, <div>I am trying to run CELL_OPT calculations. </div><div>The calculations start running correctly, but after one or a few steps it stops without any warning or error message. I stops right before the next SCF cycle. </div><div><br /></div><div>This is the part of my input:</div><div><br /></div><div><font face="Courier New">&GLOBAL<br />   PRINT_LEVEL MEDIUM<br />   PROJECT_NAME interface<br />   RUN_TYPE CELL_OPT<br />&END GLOBAL<br /><br />&MOTION<br />  &CONSTRAINT<br />    &FIXED_ATOMS<br />      COMPONENTS_TO_FIX xyz<br />      LIST 256..325<br />    &END FIXED_ATOMS<br />  &END CONSTRAINT<br /><br />  &CELL_OPT<br />    OPTIMIZER BFGS<br />    TYPE DIRECT_CELL_OPT<br />    KEEP_ANGLES   .TRUE. <br />    KEEP_SYMMETRY .TRUE.          <br />  &END CELL_OPT<br /><br />  &GEO_OPT<br />    TYPE MINIMIZATION<br />    OPTIMIZER BFGS<br />  &END GEO_OPT<br />&END MOTION<br /><br />&FORCE_EVAL<br />  METHOD QS<br />  STRESS_TENSOR ANALYTICAL<br /><br />  &DFT<br />    BASIS_SET_FILE_NAME ./BASIS_SETS<br />    POTENTIAL_FILE_NAME ./GTH_POTENTIALS<br />    <br />    &MGRID<br />      NGRIDS 4<br />      CUTOFF 300<br />      REL_CUTOFF 60<br />    &END MGRID<br /><br />    &QS<br />      EPS_DEFAULT 1.0E-14<br />    &END QS<br /><br />    &SCF<br />      SCF_GUESS RESTART<br />      EPS_SCF 1.0E-07<br />      MAX_SCF 500<br />      ADDED_MOS 1000<br />      CHOLESKY INVERSE<br /><br />      &SMEAR ON<br /><span style="white-space: pre;">      </span>     METHOD FERMI_DIRAC<br /><span style="white-space: pre;">     </span>     ELECTRONIC_TEMPERATURE 300<br />      &END SMEAR <br /><br />      &DIAGONALIZATION<br /><span style="white-space: pre;">    </span>     ALGORITHM STANDARD <br />      &END DIAGONALIZATION <br /><br />      &MIXING <br />        METHOD BROYDEN_MIXING<br />        ALPHA 0.1<br />        BETA 0.5<br />        NBROYDEN 8<br />      &END MIXING<br />    &END SCF<br /><br />    &XC<br />      &XC_FUNCTIONAL PBE<br />      &END XC_FUNCTIONAL<br />    &END XC</font></div><div><font face="Courier New"><br />  &END DFT </font></div><div><font face="Courier New"> .</font></div><div><font face="Courier New"> .</font></div><div><font face="Courier New"> .</font></div><div><font face="Courier New">&END FORCE_EVAL</font></div><div><br /></div><div>The output file just stops like so:</div><div><br /></div><div><font face="Courier New"> Number of electrons:                                                       2980<br /> Number of occupied orbitals:                                               1535<br /> Number of molecular orbitals:                                              2490<br /><br /> Number of orbital functions:                                               6500<br /> Number of independent orbital functions:                                   6500<br /><br /> Parameters for the always stable predictor-corrector (ASPC) method:<br /><br />  ASPC order: 0<br /><br />  B(1) =   1.000000<br /><br /> Extrapolation method: ASPC<br /><br /><br /> SCF WAVEFUNCTION OPTIMIZATION<br /><br />  Step     Update method      Time    Convergence         Total energy    Change<br />  ------------------------------------------------------------------------------</font></div><div><br /></div><div>The system is an interface between two solids and the cell contains approximately 400 atoms. I am using cp2k-2023.2. Does anyone have any ideas on why this might be happening and how to resolve the issue? </div><div>Thank you in advance, </div><div>Margherita </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 visit <a href="https://groups.google.com/d/msgid/cp2k/db028335-c083-4fbc-b11c-ec406058559dn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/db028335-c083-4fbc-b11c-ec406058559dn%40googlegroups.com</a>.<br />