[gmx-users] Re: Tabulated potentials and normal nonbonded interaction at the same time

Thomas Schlesier schlesi at uni-mainz.de
Thu Sep 16 12:57:35 CEST 2010


>> Hi,
>> first, thank you for you're answers.
>> I looked more into the stuff.
>>
>> I think i would use the following setup:
>>
>> Since the CG-solvent is uncharged, i can use the normal coulomb-
>> interaction (PME) for the system.
>>
>> For the vdw-part, i would use 'vdw_type=user' and:
>> energy_grps = CG Mol
>> energygrp_tables = CG CG (CG self-interaction) CG Mol (CG-
>> Molecule interaction)
>> (CG = CG-solvent; Mol=Molecule)
>> So i would need two tables (for each interaction one).
>> With this i should get all the vdw-interaction where the CG-
>> solvent is involved.
>>
>> But since i need normal vdw-interactions for the Molecule self-
>> interaction, i would need in the .mdp-file
>> 'vdw_type = user' AND 'vdw_type = cut-off'
>>
>> Or am i wrong (hope so)? Because the only way i can think of to 
>> solve the problem, would be that i make for every vdw-self-
>> interaction of the Molecule an own table, which wouldn't be good 
>> (since i read somewhere that tables are slower then to 
>> calculated the interaction normally).
>>
>> So how could i solve the problem, or this there no problem 
>> because i didn't understand the whole stuff correctly (seems so 
>> because of Marks reply)?
> 
> I haven't ever actually done it, but I'm expecting "vdw_type= cut-off" and "energygrp_tables = CG CG CG Mol" to work if you supply the two tables. What I can see in the manual neither confirms nor denies this, but it is cost-free to try :-)
> 
> Mark

Ok, i will try it and report back if it works.
Thomas




More information about the gromacs.org_gmx-users mailing list