Hi Margherita, can you attach your xyz file and your full input file as well?<br /><br />Thanks,<br /><br />Michela<br /><br /><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, November 11, 2024 at 5:44:37 AM UTC-5 Margherita wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">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></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 visit <a href="https://groups.google.com/d/msgid/cp2k/784ae1bd-4bb3-4874-8f87-3b0148d8fb1en%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/784ae1bd-4bb3-4874-8f87-3b0148d8fb1en%40googlegroups.com</a>.<br />