<div dir="ltr"><div><br></div>Dear Anna, <div><br></div><div><br></div><div>which type of calculations did you carry out in preparation of the MD run?</div><div>Is your DFT model reproducing the correct lattice structure?<br></div><div>Did you optimise the slab and check that the electronic structure of the bulk and at the surface is what expected from reference data?</div><div>I do not know this specific system, but the parameters in the DFT section seem to be not very well tuned. Is the SCF converging? Do you get reasonable forces?<br></div><div><br></div><div>Did you model the adsorption of individual water molecules to evaluate the binding energy and compare to references? </div><div>Should the layer of molecules be ice-like or do you want to simulate liquid water?</div><div>If the coordinates are obtained from some structural data or geometry optimisation, did you anneal slowly from 0K to the final temperature?</div><div>Why are you running NPE? Did you already thermalised at constant volume?</div><div><br></div><div>Kind regards</div><div>Marcella</div><div><br></div><div><br></div><div><br><br>On Thursday, March 21, 2019 at 3:55:19 PM UTC+1, Anna Anic wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hey everybody,<div><br></div><div>I am trying to do Ab-initio MD for a metal-water interface. I am using a symmetric slab and positioned a bilayer of water at both surface sites. </div><div>Now, if I run the simulation my water layer just decomposes into the vacuum region. That is why I would like to implement an vacuum potential barrier at some point in the z-direction to prevent my water from just "flying away".</div><div>I tried using the &EXTERNAL_POTENTIAL keyword and defined a harmonic function as a barrier.</div><div>But when I run the simulation the hydrogen atoms from the H2O-bilayer  are "blown away".</div><div>The question is, what did I do wrong? Is the idea legit? Does anyone have a proposal how to fix the problem or how to do it differentey.</div><div><br></div><div>Greetings,</div><div><br></div><div>Anna</div><div><br></div><div>(I attached my input file,geometry and trajectory )</div></div></blockquote></div></div>