[gmx-users] Liquid-liquid interface and electric field

[Gschwend Grégoire]

Tank you for the comment on the parameters. The system is neutral.  

I do no think that the problem is the energy barrier because in that case the ions would be blocked at the interface but they would at least move in the direction imposed by the field (and not randomly move in their respective phases). Do you know if their is a way to check that the field is really applied? 
________________________________________
De : gromacs.org_gmx-users-bounces at maillist.sys.kth.se <gromacs.org_gmx-users-bounces at maillist.sys.kth.se> de la part de David van der Spoel <spoel at xray.bmc.uu.se>
Envoyé : mercredi 16 novembre 2016 17:02
À : gmx-users at gromacs.org
Objet : Re: [gmx-users] Liquid-liquid interface and electric field

On 16/11/16 11:54, Gschwend Grégoire wrote:
> Thank you for yor answer.
>
> So, with and without electric field the ions move, sometimes form pairs or separate, etc. Everything looks "normal". The temperature of the box is stable at 300K, and the pressure around 1bar. The size of the box is stable as well. The simulation lasts 2 ns. I think that this simulation time should be enough to see an ion transfer. For instance, the simulation length was 200 ps with a field of 0.5 V/nm in this publication:  J. Phys. Chem. A., 103, 10274-10279 (1999),  (http://pubs.acs.org/doi/abs/10.1021/jp991967w).
>
> Please find below the parameters that I use for the simulation, maybe this could help...
>
For your parameters, please consider setting
fourierspacing = 0.12
and use just one group for t coupling. I hope your simulation system is
neutral in order to prevent artifacts.

It could be that the energy barrier is too high which you could find out
using a PMF calculation.

> constraint_algorithm    = lincs
> constraints               = all-bonds
> lincs_iter                = 1
> lincs_order               = 4
>
> integrator    = md
> nsteps                = 2000000
> dt                = 0.001
>
> cutoff-scheme   = Verlet
> ns_type                   = grid
> nstlist                   = 10
> rcoulomb          = 1.0
> rvdw              = 1.0
> coulombtype       = PME
> pme_order         = 4
> fourierspacing        = 0.16
>
> tcoupl                = V-rescale
> tc-grps               = Aqueous Organic
> tau_t         = 0.5 0.5
> ref_t         = 300 300
> pcoupl                        = Berendsen
> pcoupltype            = semiisotropic
> tau_p                 = 5.0
> ref_p                 = 1.0 1.0
> compressibility     = 4.5e-5 4.5e-5
>
> E-1;E-1;E-0
>
> refcoord_scaling    = com
> pbc           = xyz
> DispCorr      = EnerPres
> gen_vel               = no
>
> ________________________________________
> De : gromacs.org_gmx-users-bounces at maillist.sys.kth.se <gromacs.org_gmx-users-bounces at maillist.sys.kth.se> de la part de David van der Spoel <spoel at xray.bmc.uu.se>
> Envoyé : mercredi 16 novembre 2016 11:09
> À : gmx-users at gromacs.org
> Objet : Re: [gmx-users] Liquid-liquid interface and electric field
>
> On 16/11/16 09:51, Gschwend Grégoire wrote:
>> Hello everyone,
>>
>>
>> I am trying to simulate an ion transfer at an interface between two immiscible electrolytes solutions. I have prepared a 5x5x5 nm box with the solvents and electrolytes. The system is stable and well equilibrated. The problem now is that when I apply an electric field in the z direction (normal to the interface)  it seems that nothing happens. By nothing I mean nothing different that with a simulation without electric field. I have tried many values for the magnitude of the field, even unrealistic values as high as 5 V/nm (while the ions should transfer around 0.5 V/nm max).
>>
>>
>> In the simulation I use periodic boudary conditions in all three directions. Thus, I was wondering if this was compatible with an electric field in the z direction. Indeed, the PBC implies that the potential is the same at z=0 and z=box_size, which implies that there is no field. Of couse the field could be periodic, but in the manual it is said that the frequency of the cosine is zero.
>>
>>
>> Should I use walls and PBC in x and y directions only?
> No, that should not be necessary.
> Please be more specific than "nothing happens". Does the ion not move at
> all? Then the temperature is zero. How long is the simulation?
>>
>>
>> Thank you for your help and suggestions,
>>
>>
>> Grégoire
>>
>
>
> --
> David van der Spoel, Ph.D., Professor of Biology
> Dept. of Cell & Molec. Biol., Uppsala University.
> Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205.
> spoel at xray.bmc.uu.se    http://folding.bmc.uu.se
> --
> Gromacs Users mailing list
>
> * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-request at gromacs.org.
>


--
David van der Spoel, Ph.D., Professor of Biology
Dept. of Cell & Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205.
spoel at xray.bmc.uu.se    http://folding.bmc.uu.se
--
Gromacs Users mailing list

* Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-request at gromacs.org.