<div dir="ltr"><div>Dear <span name="Giuseppe Mattioli" class="gmail-gD">Mattioli,<br><br></span></div><span name="Giuseppe Mattioli" class="gmail-gD">Water is not a crystal and therefore it cannot be treated as a periodic str</span>ucture. The aforementioned example provides only a naive picture of water as an aqueous compartment leading to a crude approximation of the given system as well. Moreover, QE is not best suited for DFMD simulations of large, nonperiodic, liquid systems. <br></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><span style="color:rgb(0,0,0)"><span style="background-color:rgb(255,255,255)"><span style="font-size:12.8px">Regards,<br></span></span></span></div><div><span style="color:rgb(0,0,0)"><span style="background-color:rgb(255,255,255)"><span style="font-size:12.8px">Ashkan</span><i style="font-size:12.8px"><br><br></i></span></span></div><div><font size="2"><b><span style="color:rgb(0,0,0)"><span style="background-color:rgb(255,255,255)">Ashkan Shekaari</span></span></b></font></div><div><span style="color:rgb(0,0,0)"><span style="background-color:rgb(255,255,255)">Plasma Physics Research Center</span></span></div><div><span style="color:rgb(0,0,0)"><span style="background-color:rgb(255,255,255)"><span style="font-size:12.8px"><span style="font-size:12.8px">Science and Research Branch</span></span></span></span></div><div><span style="color:rgb(0,0,0)"><span style="background-color:rgb(255,255,255)">I A U, <span style="font-size:12.8px">14778-93855 Tehran, Iran.</span></span></span></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
<br><div class="gmail_quote">On Sat, May 6, 2017 at 12:21 PM, Giuseppe Mattioli <span dir="ltr"><<a href="mailto:giuseppe.mattioli@ism.cnr.it" target="_blank">giuseppe.mattioli@ism.cnr.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
Dear David and Ashkan<br>
<span class=""><br>
> You cannot simulate the system as periodic. You have to use a large super<br>
> cell with vacuum spaces along the three spatial directions. Moreover, you<br>
> can simulate the solution only via importing a large number of atomic<br>
> positions of the order of 25000 at least, which does not seem to be<br>
> feasible by QE even in case of using random generators unless you have a<br>
> supercomputer facility.<br>
<br>
</span>This comments are misleading, if not wrong in some cases.<br>
First of all you can for sure perform ab initio molecular dynamics<br>
(aimd) simulations of a water solution in a periodic box, and you do<br>
not need a huge number of atoms. One of the cp.x examples provides<br>
starting points, see here<br>
/Your_Path_To_Espresso_6.1/<wbr>Examples/CPV/example04<br>
<span class=""><br>
>> Does anyone know of examples, benchmarks, or recommendations?  Would the<br>
>> X.blyp-van.ak pseudopotentials be appropriate?  Any guidance or advice<br>
>> about parameter settings for this application would be much appreciated<br>
<br>
</span>Yes, you can use ultrasoft pseudopotentials such as X.blyp-van.ak.<br>
They used to be old and well-tested vanderbilt ultrasoft<br>
pseudopotentials generated many years ago by Axel Kohlmeyer (ak). But<br>
if you are not familiar with this terminology (ultrasoft,<br>
norm-conserving, van, mt, rrkj, paw, ...) you should start with some<br>
tutorial on pseudopotentials, because you must be sure to use<br>
converged plane-wave and density cutoffs for a given set of<br>
pseudopotentials.<br>
<br>
There is a very large number of options that must be set in the cp.x<br>
input. I cannot review them here. But there are also tons of<br>
literature on aimd simulations of water solutions, and you will easily<br>
find something that will help to choice a lot of parameters (box<br>
dimensions, NVT, NVP, NVE dynamics, thermostats, ...)<br>
<span class=""><br>
> Indeed, QE is not best suited for MD simulations and I strongly recommend<br>
> the gromacs package.<br>
<br>
</span>Of course Ashkan might be right on a couple of points: if you have<br>
very big molecules weakly interacting with the solvent, and you are<br>
interested in the morphology of the solute only, then you have to<br>
perform very long md simulations and to use a large number of water<br>
molecules, and this might be costly and time-consuming if performed at<br>
an aimd level. If you want to study proton exchanges, reactions,<br>
dipole dynamics, ... you need aimd. You surely know whether your<br>
scientific task strictly requires aimd rather than model-potential md.<br>
<br>
HTH<br>
Giuseppe<br>
<span class=""><br>
Quoting ashkan shekaari <<a href="mailto:shekaari@gmail.com">shekaari@gmail.com</a>>:<br>
<br>
> Indeed, QE is not best suited for MD simulations and I strongly recommend<br>
> the gromacs package.<br>
<br>
</span><div><div class="h5">> You cannot simulate the system as periodic. You have to use a large super<br>
> cell with vacuum spaces along the three spatial directions. Moreover, you<br>
> can simulate the solution only via importing a large number of atomic<br>
> positions of the order of 25000 at least, which does not seem to be<br>
> feasible by QE even in case of using random generators unless you have a<br>
> supercomputer facility.<br>
><br>
> On May 6, 2017 9:09 AM, "D J Anick" <<a href="mailto:david.anick@rcn.com">david.anick@rcn.com</a>> wrote:<br>
><br>
>> Hello plane wavers,<br>
>><br>
>> I am interested in using QE for a molecular dynamics simulation of an<br>
>> aqueous solution containing a solute, modeling it as a 3-D periodic cell.<br>
>> Principal questions would be about solvation shell geometries, distribution<br>
>> of configurations adopted by the solute, and H-bond duration / stability.<br>
>><br>
>> Does anyone know of examples, benchmarks, or recommendations?  Would the<br>
>> X.blyp-van.ak pseudopotentials be appropriate?  Any guidance or advice<br>
>> 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>
>><br>
<br>
<br>
</div></div>--<br>
******************************<wbr>**************************<br>
- Article premier - Les hommes naissent et demeurent<br>
libres et égaux en droits. Les distinctions sociales<br>
ne peuvent être fondées que sur l'utilité commune<br>
- Article 2 - Le but de toute association politique<br>
est la conservation des droits naturels et<br>
imprescriptibles de l'homme. Ces droits sont la liberté,<br>
la propriété, la sûreté et la résistance à l'oppression.<br>
******************************<wbr>**************************<br>
<span class="HOEnZb"><font color="#888888"><br>
    Giuseppe Mattioli<br>
    CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<br>
    v. Salaria Km 29,300 - C.P. 10<br>
    I 00015 - Monterotondo Stazione (RM)<br>
    Tel + 39 06 90672836 - Fax +39 06 90672316<br>
    E-mail: <<a href="mailto:giuseppe.mattioli@ism.cnr.it">giuseppe.mattioli@ism.cnr.it</a>><br>
    <a href="http://www.ism.cnr.it/en/staff/giuseppe-mattioli/" rel="noreferrer" target="_blank">http://www.ism.cnr.it/en/<wbr>staff/giuseppe-mattioli/</a><br>
    ResearcherID: F-6308-2012<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
<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></div></div></blockquote></div><br></div>