Thanks for your reply,  Marcella. <div>I have found a solution. Once I change coordination input file into .xyz rather than .pdb, the initalzation time decreases a lot. I didn't know why, but it works.<br></div><div>Output file is attached, it may provide more information.</div><div><br></div><div><br></div><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">在2020年9月5日星期六 UTC+8 下午9:37:18<Marcella Iannuzzi> 写道：<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div><br></div>The direct comparing  of systems with different electronic properties is not straightforward.<div>There is no strong reason to expect the same behaviour. </div><div>CeO2 is known to be a difficult case,  see the  the literature. You might need to adapt the model.</div><div><br></div><div>Without output and timings it is anyway hard to say something and give advise. </div><div>Do you restart the MD reading the wavefunction from the restart file as initial guess?</div><div><br></div><div>Kind regards</div><div>Marcella</div><div><br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Saturday, September 5, 2020 at 5:02:20 AM UTC+2 <a href data-email-masked rel="nofollow">dm...@gmail.com</a> 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">Hi,<div>I started an AIMD calculation of CeO2(110) slab, including 144 atoms and 2016 electrons. With current configuration, it costs more than 300s for one MD step. Compared to another system, Pd20 cluster on TiO2(110), including 451 atoms and 3810 electrons, with the same DFT configuration, it costs about 100s for one step. I'm using cp2k 7.1, OT method and ASPC as extrapolation method, and i found that it only cost less than 30s for SCF converged, yet 200s needed to start next MD step.</div><div>Is there any problems with CeO2 systems or my input files?</div><div>Thanks!</div></blockquote></div></blockquote></div>