<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body>
<p>Dear Geng,</p>
<p>You should decrease eps_scf to 1e-7 or 1e-8. The default 1e-6 is
good for the initial stage of the optimization, when the
forces/pressure are still large. Later on this should be adjusted,
especially if you request such tight convergence criteria in
&CELL_opt. Decreasing this value reduces the noise in the
forces and improves the convergence of the cell optimization.</p>
<p>Cheers,</p>
<p>Fabian<br>
</p>
<br>
<div class="moz-cite-prefix">On 03.01.2022 21:40, GENG YUAN wrote:<br>
</div>
<blockquote type="cite"
cite="mid:d824ebcb-7f45-41ec-9255-25fe50047760n@googlegroups.com">
<meta http-equiv="content-type" content="text/html; charset=UTF-8">
Here are the screenshots.
<div><img alt="beforecellopt_view1.png" data-iml="2190296"
src="cid:c479680e-3000-4027-8361-deb4c4072de3"
moz-do-not-send="true" width="422px" height="282.166px"><img
alt="aftercellopt_view1.png" data-iml="2233329.8000001907"
src="cid:bdf8e830-997f-4c90-a48d-efcd24bf5831"
moz-do-not-send="true" width="417px" height="283.208px"><br>
</div>
<div>Top view before cell optimization
Top view after
cell optimization</div>
<div><br>
</div>
<div><img alt="beforecellopt_view2.png"
data-iml="2288950.6000003815"
src="cid:c4653aca-771f-46f2-abd8-5aaadd07226d"
moz-do-not-send="true" width="405px" height="256.076px"><img
alt="aftercellopt_view2.png" data-iml="2324640.200000286"
src="cid:8249d18c-3e36-4210-821b-1ea849f3b910"
moz-do-not-send="true" width="406px" height="259.467px"><br>
</div>
<div>Another view before cell optimization
Another view after
cell optimization</div>
<div><br>
</div>
<div>Sincerely,</div>
<div>Geng</div>
<div><br>
</div>
<div class="gmail_quote">
<div dir="auto" class="gmail_attr">在2022年1月3日星期一 UTC-5
15:34:56<GENG YUAN> 写道：<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><br>
</div>
<div>Hope you are going great and happy new year.</div>
<div><br>
</div>
<div>I am running the cell optimization for a single bulk Pt
(134 atoms) with fcc 111 surface. By checking my output
file, I notice that the cell optimization takes 446 steps to
converge (which is almost close to the maximum number of
optimization steps:500). Besides, a certain degree of
structure deformation is observed during optimization
(please see the attached files and screenshots).</div>
<div><br>
</div>
<div>I am wondering if there are any other ways to let the
cell optimization converge faster (in other words, converge
in fewer steps) and whether the above-mentioned deformation
looks reasonable?</div>
<div>Below is the input: (attached are input file, output
file, original coordinate, trajectory, and the screenshots
of the trajectory)</div>
<div>------------------------------------------------------------</div>
<div><font size="2">&GLOBAL<br>
PROJECT Pt_cellopt<br>
RUN_TYPE CELL_OPT<br>
PRINT_LEVEL LOW<br>
&END GLOBAL<br>
&FORCE_EVAL<br>
METHOD QS<br>
STRESS_TENSOR ANALYTICAL<br>
&SUBSYS<br>
&CELL<br>
ABC 13.7 13.7 13.7<br>
SYMMETRY CUBIC<br>
&CELL_REF<br>
ABC 13.7*1.5 13.7*1.5 13.7*1.5<br>
&END <br>
&END CELL<br>
&TOPOLOGY<br>
COORD_FILE_NAME ./Pt134.xyz<br>
COORDINATE xyz<br>
&END<br>
&KIND Pt<br>
BASIS_SET DZVP-MOLOPT-SR-GTH<br>
POTENTIAL GTH-PBE-q18<br>
&END KIND<br>
&END SUBSYS<br>
&DFT<br>
BASIS_SET_FILE_NAME ./BASIS_MOLOPT<br>
POTENTIAL_FILE_NAME ./GTH_POTENTIALS<br>
&QS<br>
EPS_DEFAULT 1.0E-12<br>
&END QS<br>
&MGRID<br>
CUTOFF 400<br>
NGRIDS 5<br>
REL_CUTOFF 40<br>
&END MGRID<br>
&SCF<br>
SCF_GUESS ATOMIC<br>
EPS_SCF 1.0E-06<br>
MAX_SCF 200<br>
ADDED_MOS 100<br>
&OUTER_SCF<br>
MAX_SCF 50<br>
EPS_SCF 1.0E-6<br>
&END<br>
&DIAGONALIZATION T<br>
ALGORITHM STANDARD<br>
&END DIAGONALIZATION<br>
&MIXING T<br>
METHOD BROYDEN_MIXING<br>
ALPHA 0.1<br>
NBROYDEN 8<br>
&END MIXING<br>
&SMEAR ON<br>
METHOD FERMI_DIRAC<br>
ELECTRONIC_TEMPERATURE [K] 500<br>
&END SMEAR<br>
&PRINT<br>
&RESTART ON<br>
&END RESTART<br>
&END PRINT<br>
&END SCF<br>
&XC<br>
&XC_FUNCTIONAL<br>
&LIBXC<br>
FUNCTIONAL XC_GGA_X_PBE_R<br>
&END<br>
&LIBXC<br>
FUNCTIONAL XC_GGA_C_PBE<br>
&END<br>
&END XC_FUNCTIONAL<br>
&VDW_POTENTIAL<br>
POTENTIAL_TYPE PAIR_POTENTIAL<br>
&PAIR_POTENTIAL<br>
LONG_RANGE_CORRECTION<br>
PARAMETER_FILE_NAME dftd3.dat<br>
TYPE DFTD3<br>
REFERENCE_FUNCTIONAL revPBE<br>
&END PAIR_POTENTIAL<br>
&END VDW_POTENTIAL<br>
&XC_GRID<br>
XC_DERIV NN50_SMOOTH<br>
&END XC_GRID<br>
&END XC<br>
&END DFT<br>
&END FORCE_EVAL<br>
&MOTION<br>
&CELL_OPT<br>
EXTERNAL_PRESSURE 0<br>
TYPE DIRECT_CELL_OPT<br>
KEEP_SYMMETRY TRUE<br>
MAX_DR 1.0E-05<br>
MAX_FORCE 1.0E-05<br>
RMS_DR 1.0E-05 <br>
RMS_FORCE 1.0E-05<br>
MAX_ITER 500<br>
OPTIMIZER BFGS<br>
&END CELL_OPT<br>
&END MOTION </font><br>
</div>
<div><br>
</div>
<div>Many thanks in advance,</div>
<div>Sincerely,</div>
<div>Geng</div>
<div><br>
</div>
</blockquote>
</div>
-- <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"
moz-do-not-send="true" class="moz-txt-link-freetext">cp2k+unsubscribe@googlegroups.com</a>.<br>
To view this discussion on the web visit <a
href="https://groups.google.com/d/msgid/cp2k/d824ebcb-7f45-41ec-9255-25fe50047760n%40googlegroups.com?utm_medium=email&utm_source=footer"
moz-do-not-send="true">https://groups.google.com/d/msgid/cp2k/d824ebcb-7f45-41ec-9255-25fe50047760n%40googlegroups.com</a>.<br>
</blockquote>
</body>
</html>
<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/715e86cc-3f4e-6881-ecdc-a8bd6b12777d%40gmail.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/715e86cc-3f4e-6881-ecdc-a8bd6b12777d%40gmail.com</a>.<br />