<div dir="ltr">Hi Flo,<div><br></div><div>Thanks for the clarification, I started with optimisations to get used to CP2K input/output and timings for larger systems. I'm moving MD when I add in the ions, I have a second question if you don't mind- If I'm adding in a +2 ion is it just a case of adding the CORE_CORRECTION tag into my KIND section for that atom?<br><br><br>Thank you!</div><div><br></div><div>Olivia </div><div><br>On Friday, May 29, 2015 at 7:52:41 AM UTC+1, Florian Schiffmann wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hi,<br><br>it might be that I misunderstand what you intend to do but for water boxes, i.e. a liquid, it is not meaningful to do optimizations. The energetic minimum should be the crystal or ice in your case. As liquids only exist at finite temperature you might have to sample using MD. <br>If you are working with ice structures then normal mode analysis can be used just the way everbody else described<br><br>Flo<br><br></div></blockquote></div></div>