From: Athreya, Nagendra Bala Murali (nathreya_at_illinois.edu)
Date: Tue Aug 11 2020 - 19:50:48 CDT
I am trying to simulate an ion channel under electric bias with asymmetric concentration gradient conditions (where X,Y-directions can be periodic, however Z-direction needs to be non-periodic for ions). The ions on one side of the channel is 10mM and the other side is 1M.
When I tried to simulate with PBC (using wrapAll = ON), the ions move to the other side and the whole system neutralizes over time.
When I turn wrapAll OFF, and use a tclBC script to push ions from the z-boundaries (as suggested here: https://www.ks.uiuc.edu/Research/namd/mailing_list/namd-l.2015-2016/1610.html), the ions just expand in all directions but everything else remains in the initial defined cuboid.
I am not sure what is going on. I have attached the minimization, equilibration, and script when voltage is applied. Am I missing something?
I would appreciate any help I can get.
This archive was generated by hypermail 2.1.6 : Fri Dec 31 2021 - 23:17:09 CST