<div>Dear Marcella,</div><div>Thanks for your reply!</div><div>I have used VASP for doing AIMD simulations, the CP2K input parameter comes from the VASP input file.</div><div>Because VASP NPT simulation is not stable, I want to use CP2K NPT ensemble for doing AIMD simulations.</div><div>As you suggested, I change a small system (ZrO2)26(Y2O3)3, and the system has 93 atoms and about 10 angstroms per side.</div><div>In fact, the PAW energy cutoff for VASP is 400 eV, and I change the cutoff for CP2K calculation (in attachment).</div><div>In the input set of CP2K, the basis set and potential is DZVP-GTH-PADE, does the DZVP-GTH-PADE belong to LDA ?</div><div>I don't know how to set input for annealing temperature, annealing velocity, and thermostat.</div><div>In fact, I want to anneal from melting state (1200 K) to 900 K with about the annealing velocity 300K/1.5ps (2*10^14 K/s),</div><div>and then do NVT ensemble simulation at 900 K from the structure of the NPT (1200 ->900 K NPT simulations), until reach 300 K, and then calculate the</div><div>phonon density of states at 300 K.</div><div>That is NPT( 1200->900 K), and then NVT(900 K), and NPT(900->600 K), and then NVT(600 K), and NPT(600-300 K), and then NVT(300 K) simulations.</div><div>And how to set input file for annealing temperature, annealing velocity, and thermostat?</div><div>Thanks,</div><div>Zhijian</div><br>