<div dir="ltr">Hi Olivia,<br><br>fi you don't add the counter ions you will have to set the Charge keyword in the DFT section as well. If you are dealing with a charged system you will not be able to perform cell optimizations as the stress tensor with a periodic charge is not defined as far as I know. So in that case you will have to switch to a normal GEO_OPT.<br><br>Cheers<br>Flo<br><br>On Friday, May 29, 2015 at 9:49:08 PM UTC+10, Olivia Lynes wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><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></blockquote></div>