- Posts: 98
Inserting position restraints
- Sithara Perera
- Topic Author
- Visitor
6 years 11 months ago #7225
by Sithara Perera
Inserting position restraints was created by Sithara Perera
Hi!,
I' making a 512 DPPC bi-layer out of 128 bilayer available in www.cgmartini.nl . The force field used is martini.
But once the energy minimization is done, the bi-layer crashes and accumulates to globules.
So I made an position.itp file using the command
gmx genrestr -f dppc221.gro -n test.ndx -o lipid_posre.itp -fc 1000 1000 1000
And included in the topology file.
Topology file also include
include "martini_v2.0.itp"
#include "martini_v2.0_lipids.itp" from martini as well.
But then when the command
gmx grompp -f minim.mdp -c dppc221.gro -p dppc_bilayer.top -o em2.tpr
is given , following error occurs.
Fatal error:
[ file posre.itp, line 29 ]:
Atom index (25) in position_restraints out of bounds (1-24).
This probably means that you have inserted topology section "position_restraints"
in a part belonging to a different molecule than you intended to.
In that case move the "position_restraints" section to the right molecule.
How can I fix this?
I' making a 512 DPPC bi-layer out of 128 bilayer available in www.cgmartini.nl . The force field used is martini.
But once the energy minimization is done, the bi-layer crashes and accumulates to globules.
So I made an position.itp file using the command
gmx genrestr -f dppc221.gro -n test.ndx -o lipid_posre.itp -fc 1000 1000 1000
And included in the topology file.
Topology file also include
include "martini_v2.0.itp"
#include "martini_v2.0_lipids.itp" from martini as well.
But then when the command
gmx grompp -f minim.mdp -c dppc221.gro -p dppc_bilayer.top -o em2.tpr
is given , following error occurs.
Fatal error:
[ file posre.itp, line 29 ]:
Atom index (25) in position_restraints out of bounds (1-24).
This probably means that you have inserted topology section "position_restraints"
in a part belonging to a different molecule than you intended to.
In that case move the "position_restraints" section to the right molecule.
How can I fix this?
Please Log in or Create an account to join the conversation.
- bart
- Offline
- Admin
Less
More
6 years 11 months ago #7226
by bart
Replied by bart on topic Inserting position restraints
Heya,
I am not sure, maybe someone else would know, but to make a 512 lipids bilayer from scratch would not take more than a few minutes. You can use 'insane', a program developed by our lab which you can find under downloads (its part of the complex bilayer tutorial), or you could even do random assembly and be done within an hour. Of course this doesn't answer you question, but it might still help you continue you research. For sure I would not walk down the path of using constraints anyway. I could advise on doing the lipids tutorials anyway if you intend to use martini for that purpose. We just released the GMX 5 version.
Cheers
I am not sure, maybe someone else would know, but to make a 512 lipids bilayer from scratch would not take more than a few minutes. You can use 'insane', a program developed by our lab which you can find under downloads (its part of the complex bilayer tutorial), or you could even do random assembly and be done within an hour. Of course this doesn't answer you question, but it might still help you continue you research. For sure I would not walk down the path of using constraints anyway. I could advise on doing the lipids tutorials anyway if you intend to use martini for that purpose. We just released the GMX 5 version.
Cheers
Please Log in or Create an account to join the conversation.
Time to create page: 0.099 seconds