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

Reply via email to