
To assign 13 repulsion term (atoms separated by two bonds) of force field, what command should I use?

I’d like to include independent dihedral terms for cis (ijkl) and trans (ijkl’) cases where k is connected to j, l, and l’ (ijkl and ijkl’ are almost planar and the equil. dihedral angle is 0 for ijkl while it is 180 for ijkl’).
How can I define i,j,k,l,l’ for it?
Is it enough for it by employing two commands line assigning different dihedral types (e.g., 1,2) and using phi of 0, 180 respectively?
Thanks for your help in advance.
Best regards,
Luke
2011/7/26 JhonY. I. <[email protected]…8…>
 To assign 13 repulsion term (atoms separated by two bonds) of force field, what command should I use?
How are the atoms repelled? Is it via standard force field form? If so, then you can use the special_bonds command to tell LAMMPS to allow for something different than the default. If it’s not something that LAMMPS has builtin, you’ll need to modify the code to do it.
 I’d like to include independent dihedral terms for cis (ijkl) and trans (ijkl’) cases where k is connected to j, l, and l’ (ijkl and ijkl’ are almost planar and the equil. dihedral angle is 0 for ijkl while it is 180 for ijkl’).
How can I define i,j,k,l,l’ for it?
Is it enough for it by employing two commands line assigning different dihedral types (e.g., 1,2) and using phi of 0, 180 respectively?
Yes, you should be able to do it by defining two different types with different values for the corresponding coefficients in the potential.
–AEI