On 5/19/19 11:32 PM, Nidhin Thomas wrote:
Hello GROMACS users, I tried to apply a static electric field across a water box. I used anisotropic pressure coupling. Details of the pressure coupling is given below. pcoupl = Parrinello-Rahman pcoupltype = anisotropic tau_p = 5.0 compressibility = 4.5e-5 4.5e-5 4.5e-5 0 0 0 ref_p = 1.0 1.0 1.0 0 0 0 ; constraints = h-bonds constraint_algorithm = LINCS continuation = yes ; electric-field-y = 0.1 0 0 0 I see that water box size reduced in the direction of application of electric field (Y) and increased in X and Z direction. Original Size = (3.90 x 3.90 x 3.90) nm Final Size = (4.39 x 2.66 x 5.16) nm Could someone please help me understand why the water box is shrinking in the direction of electric field ? can someone please let me know how external electric field is actually implemented in GROMACS ?
Anisotropic pressure coupling intrinsically leads to the distortion of the box shape. Run a control simulation of water with no electric field and you will see the same thing.
-Justin -- ================================================== Justin A. Lemkul, Ph.D. Assistant Professor Office: 301 Fralin Hall Lab: 303 Engel Hall Virginia Tech Department of Biochemistry 340 West Campus Dr. Blacksburg, VA 24061 jalem...@vt.edu | (540) 231-3129 http://www.thelemkullab.com ================================================== -- 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-requ...@gromacs.org.