On 6/18/13 12:20 PM, Bastien Loubet wrote:
Dear GROMACS user,
I have been trying to define parameters for virtual sites in CHARMM lipid,
and I would like to know if there is something simillar to the [ bondtypes ]
for virtual sites.
Specifically the way I specify the parameters now is by having lines in the
topology such as (for a 3fout virtual sites):
[ virtual_sites3 ]
; aiajakal functc0c1c2
40393751 4 -0.26909 -0.26909 4.5519
But then I have to do that for every virtual sites in the molecule, even if
the parameters are the same.
What I would like to be able to do is:
[ ]
; aiajakal functc0c1c2
H C CC 4 -0.26909 -0.26909 4.5519
and then
[ virtual_sites3 ]
; aiajakal functc0c1c2
40393751 4
Where atom 40 is of type H and atom 39, 37 and 51 of type C and should
be are directive similar to [ bondtypes ] but for virtual site.
Does this directive exist ? I had no luck checking the manual.
Not likely, because virtual sites aren't generally built based on types, but
based on specific atoms. The workaround that came to mind is to make use of
#define macros, like how the Gromos96 force fields define bonded terms. You can
simply:
#define param -0.26909 -0.26909 4.5519
Then, for your virtual sites:
[ virtual_sites3 ]
; aiajakal functc0c1c2
40393751 4 param
Not quite as efficient as defining types, but since that's not an option, it
does save some typing.
-Justin
--
Justin A. Lemkul, Ph.D.
Research Scientist
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
