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>

<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/1fd00a66-b8f5-418d-ac5b-4c160bcde24fn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/1fd00a66-b8f5-418d-ac5b-4c160bcde24fn%40googlegroups.com</a>.<br />