Shuangxing Dai wrote:
Thank you for your help. I checked all the force field, all of them are
with default non-bond params in LJ. So I still need to try new solution.
I add my non bond parameters in the .top file generated. Force field
files were not changed. Then I use grompp, the error was:
That wasn't the question. Do you *need* the rest of the force field? If
you're only doing ZnO then you don't need the rest. IIRC you've never
told us what you're trying to simulate, which tends to frustrate any
effort to help you.
'[ nonbond_params ]'
Invalid order for directive nonbond_params
Yep, these can't appear anywhere after a [molecule] directive. See
chapter 5 of the manual.
Rather confused by setting up Buckingham potential parameters...
Make sure you read the advice you're being given:
If you need other atoms as well, you will need to find
a consistent force field, probably or completely LJ
or completely Buckingham.
and
> No, because it's not a reasonable goal. Your new atomtypes are
talking a
> different vdW "language". That doesn't work.
and
> > And what
> > force field to choose if not ffgmx?
>
> That's your homework I'm afraid. You'll have to read the literature.
Cheers,
Mark
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/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/mailing_lists/users.php