<p dir="ltr">You cannot simulate the system as periodic. You have to use a large super cell with vacuum spaces along the three spatial directions. Moreover, you can simulate the solution only via importing a large number of atomic positions of the order of 25000 at least, which does not seem to be feasible by QE even in case of using random generators unless you have a supercomputer facility. </p>
<div class="gmail_extra"><br><div class="gmail_quote">On May 6, 2017 9:09 AM, "D J Anick" <<a href="mailto:david.anick@rcn.com">david.anick@rcn.com</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello plane wavers,<br>
<br>
I am interested in using QE for a molecular dynamics simulation of an aqueous solution containing a solute, modeling it as a 3-D periodic cell.  Principal questions would be about solvation shell geometries, distribution of configurations adopted by the solute, and H-bond duration / stability.<br>
<br>
Does anyone know of examples, benchmarks, or recommendations?  Would the X.blyp-van.ak pseudopotentials be appropriate?  Any guidance or advice about parameter settings for this application would be much appreciated.<br>
<br>
Thank you in advance,<br>
David Anick<br>
david.anick###<a href="http://rcn.com" rel="noreferrer" target="_blank">rcn.com</a><br>
______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br>
</blockquote></div></div>