- Posts: 8
polarizable water - vapor causes box to expand in semiisotropic simulation
- mboy@seas.upenn.edu
-
Topic Author
- Offline
- Fresh Boarder
Less
More
1 year 10 months ago - 1 year 10 months ago #9287
by mboy@seas.upenn.edu
polarizable water - vapor causes box to expand in semiisotropic simulation was created by mboy@seas.upenn.edu
Hi all,
I am attempting to use MARTINI to simulate a polyelectrolyte polymer brush using polarizable MARTINI water. It involves using an explicit substrate with position restraints on the substrate and the ends of the polymer brush, and periodic boundary conditions in xyz. I run the following steps:
1) initialize the box using packmol
2) steep energy minimization
3) NVT annealing at 400K to 300K
4) NPT at 300K
Thermostat and barostat settings are shown below:
dt = 0.010
tcoupl = v-rescale
tc-grps = INIT_PMCG PCL_W SUBS
tau_t = 1.0 1.0 1.0
ref_t = 300 300 300
Pcoupl = berendsen
Pcoupltype = semiisotropic ; allows for different pressure in x-y and z directions
tau_p = 6.00 ;
compressibility = 0 3e-4 ; x-y in compressible, compressible in z direction
ref_p = 1.0 1.0 ;
However, when run the NPT step, vapor bubbles of low density seem to form in the polarizable water, and the z dimension of the box keeps expanding until the simulation crashes.
Has anyone seen this issue in membrane or other semi-isotropic pressure simulations? Do you have any suggestions to prevent vapor bubble formation during equilibration?
Any help would be much appreciated!
Thanks,
Mike
I am attempting to use MARTINI to simulate a polyelectrolyte polymer brush using polarizable MARTINI water. It involves using an explicit substrate with position restraints on the substrate and the ends of the polymer brush, and periodic boundary conditions in xyz. I run the following steps:
1) initialize the box using packmol
2) steep energy minimization
3) NVT annealing at 400K to 300K
4) NPT at 300K
Thermostat and barostat settings are shown below:
dt = 0.010
tcoupl = v-rescale
tc-grps = INIT_PMCG PCL_W SUBS
tau_t = 1.0 1.0 1.0
ref_t = 300 300 300
Pcoupl = berendsen
Pcoupltype = semiisotropic ; allows for different pressure in x-y and z directions
tau_p = 6.00 ;
compressibility = 0 3e-4 ; x-y in compressible, compressible in z direction
ref_p = 1.0 1.0 ;
However, when run the NPT step, vapor bubbles of low density seem to form in the polarizable water, and the z dimension of the box keeps expanding until the simulation crashes.
Has anyone seen this issue in membrane or other semi-isotropic pressure simulations? Do you have any suggestions to prevent vapor bubble formation during equilibration?
Any help would be much appreciated!
Thanks,
Mike
Last edit: 1 year 10 months ago by mboy@seas.upenn.edu.
Please Log in or Create an account to join the conversation.
- mboy@seas.upenn.edu
-
Topic Author
- Offline
- Fresh Boarder
Less
More
- Posts: 8
1 year 10 months ago #9299
by mboy@seas.upenn.edu
Replied by mboy@seas.upenn.edu on topic polarizable water - vapor causes box to expand in semiisotropic simulation
Running the first npt simulation at 1000 bar reference pressure, then allowing the pressure to rescale to 1 bar reference pressure eliminated the vapor bubble.
Cheers,
Mike
Cheers,
Mike
Please Log in or Create an account to join the conversation.
- riccardo
-
- Offline
- Platinum Boarder
1 year 10 months ago #9300
by riccardo
Replied by riccardo on topic polarizable water - vapor causes box to expand in semiisotropic simulation
Thanks for following up in your own thread with a solution, Mike!
Please Log in or Create an account to join the conversation.
Time to create page: 0.092 seconds