<div>Dear all,</div><div><br /></div>I wanted to do a simulation of water molecules inducing catalyst surface remodeling, so I wanted to implement the addition of a limiting potential to bring water molecules close to the catalyst surface, leading to catalyst surface remodeling.<br /><br />This was done by referring to the SI section of this literature (doi.org/10.1021/acscatal.8b01232) <br /><div><br /></div><div>I have set the limiting potential in the inp file as follows. Where 1-141 is specified for all atoms in the water molecule, Z0 is taken as the z-value of the catalyst surface, and B is approximated by the calculation of the reference formula.</div><div><i>  &EXTERNAL_POTENTIAL<br />    ATOMS_LIST 1..141<br />    FUNCTION -1 * (1 - ((Z - Z0) / B)^12) / (1 - ((Z - Z0) / B)^24)<br />    VALUES [eV] 0.20 [angstrom] 5.0<br />    PARAMETERS B Z0<br />  &END EXTERNAL_POTENTIAL<br />&END FORCE_EVAL</i><br /></div><div>Running 2.7ps with NVT resulted in some of the water molecules going up and some going down, and no dissociation of the water molecules was found to occur. So, my question is : 1.Does this mean that the effect of the catalyst (Cu2O) surface on the adsorption of water molecules is stronger than the effect of the confining potential (i.e., running NVT without confining potential also gives this result, and adding confining potential is not very useful?) 2. What should I do if I want all water molecules above the vacuum layer to be close to the catalyst surface, or even if water dissociation occurs?<i><br /></i></div><div><br /></div><div><div>Your help would be greatly appreciated!</div><div><br /></div><div>All the best,</div></div><div>chorin.</div><div><br /></div><div><i><br /></i></div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/7bf8fff8-936f-4d3f-86d9-f6b65b7cde49n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/7bf8fff8-936f-4d3f-86d9-f6b65b7cde49n%40googlegroups.com</a>.<br />