- Posts: 1
CG simulations of a membrane protein with an applied electric field?
- mprevost
- Topic Author
- Offline
- Fresh Boarder
Less
More
5 years 7 months ago #7791
by mprevost
CG simulations of a membrane protein with an applied electric field? was created by mprevost
Hello,
I am trying to simulate a CG system containing a membrane protein subjected to an electric field.
I constructed the system with Charmm-gui and ran the first minimisation and equilibrations with the Charm-gui grimaces input file without problems.
A following equilibration with an imposed electric field of Ez=10 V/nm crashes. Reducing the electric field to 0.5 solves the problem. However the molecular system shrinks along the Z direction (direction perpendicular to the membrane) to an almost 2D system!
Could you tell me what's wrong? Can one simulate a CG system with an applied electric field?
Thanks on advance for your help.
I am trying to simulate a CG system containing a membrane protein subjected to an electric field.
I constructed the system with Charmm-gui and ran the first minimisation and equilibrations with the Charm-gui grimaces input file without problems.
A following equilibration with an imposed electric field of Ez=10 V/nm crashes. Reducing the electric field to 0.5 solves the problem. However the molecular system shrinks along the Z direction (direction perpendicular to the membrane) to an almost 2D system!
Could you tell me what's wrong? Can one simulate a CG system with an applied electric field?
Thanks on advance for your help.
Please Log in or Create an account to join the conversation.
- bart
- Offline
- Admin
Less
More
- Posts: 98
5 years 7 months ago #7803
by bart
Replied by bart on topic CG simulations of a membrane protein with an applied electric field?
Could it be that you box is set to semiisotropic and that the imposed electric field introduces pressure in z? Then the box would adept to that, squeezing it smaller and smaller. A solution could be to to do nvt or to set the box to isotropic, or maybe the best, to meassure the imposed pressure and set the z pressure to the new imposed level for the barostat and keep it at npt semiisotropic.
Please Log in or Create an account to join the conversation.
Time to create page: 0.086 seconds