Dear Onofrio, <div><br></div><div>To me you system looks quite well behaved. </div><div>The temperature is not low, maybe it helps reducing the time-step of activating a variable time-step, to avoid too large moves between two integration steps. </div><div>Best</div><div>Marcella</div><div><br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, August 30, 2022 at 2:14:03 PM UTC+2 onofr...@gmail.com 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 community,<p>I have done a
second-generation Car-Parrinello (SGCP) MD simulation for a system of
6 toluene molecules on a 4-layer Cu slab (6x6) at 900K, as follows:</p>
<p>1. the initial
system (including 3 adsorbed toluene + 3 toluene in gas phase, see
attached image) was first geometrically optimized</p>
<p>2. to equilibrate
the geometrically-optimized system I used NVE + velocity rescaling
looking at the conservation of temperature per region (1°region: Cu
slab; 2°region: toluene molecules) by using the &THERMAL_REGION
section with DO_LANGEVIN_DEFAULT  F and &DEFINE_REGION. See
attached “2_step___Equilibration” plot (4 ps of BOMD run)</p>
<p>3. 4 ps of sgcpMD
run. I tuned GAMMA and NOISY_GAMMA per region monitoring again the
conservation of temperature per region (see attached plot
“3_step___sgcpMD” and text files).</p>
<p><br>

</p>
<p>However, after about the
first 400 sgcpMD steps the SCF runs started to converge in more than 10 steps
(between 11 and 16) from 3 steps, as can be seen in the attached .out file. I have
tried different settings in the &OT section and I haven’t found
any solution to keep the number of steps per SCF run small. 
</p>
<p>Note: 
</p>
<p>- FULL_KINETIC as
preconditioner seems to work better than FULL_ALL.</p><p>- every time I
restart the same sgcpMD run the first SCF runs converge in few steps
(typically 3), but then the number increase to more than 10
steps.</p>
<p>How to keep the
number of steps per SCF run small? This behavior could be attributed
to a bad initial geometry or exclusively to the &OT section
settings?<br></p>
<p>Any help would be
greatly appreciated.<br></p>
<p>Thanks.<br>
<br>
Onofrio</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/2e9059b2-9149-479f-b1c2-0acc0f9bd047n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/2e9059b2-9149-479f-b1c2-0acc0f9bd047n%40googlegroups.com</a>.<br />