I am using the latest version of LAMMPS, and I am using these statements in LAMMPS input script to create bonds between atom type 1 and atome type 2 and also between these two types of atoms themselves:
However, an error is generated, saying that:
Fix bond/create induced too many angles/dihedrals/impropers per atom
Any opinion about what might cause this kind of error?
Thanks!
Best,
Jim Miao
Physics Dept.
Case Western Reserve Univ.
"latest" is a very inaccurate description. it can change within an hour
what is the latest version. currently the latest is 6 June 2014. is that
the version you are using?
and I am using these statements in LAMMPS input script to create bonds
between atom type 1 and atome type 2 and also between these two types of
atoms themselves:
However, an error is generated, saying that:
Fix bond/create induced too many angles/dihedrals/impropers per atom
Any opinion about what might cause this kind of error?
probably you didn't reserve sufficient space in LAMMPS internal data
structures to accommodate all those bonds/angles/dihedrals.
Also, if I set the “Rmin” in the “fix bond/create” to be 3.5 rather than 4.0, it works.
I checked the created bonds, dihedrals, and anlges, they are correct.
This error happens when Rmin is a little large.
If you believe lammps is not working correctly, please provide a simple and complete test case, and we can look into it. Until then, we have to assume that the error is in your input deck. It is near impossible to debug things like this based on a few fragments of information.
Please check the documentation very carefully. Perhaps there is something that you have overlooked.
The exact error massage is :
ERROR: Fix bond/create induced too many angles/dihedrals/impropers per atom (…/fix_bond_create.cpp:679)
I find this error happens when we issue these three fixes together:
There are several ways for this error to be trigger, so its hard to know without seeing your input. The first check is to make sure there is really enough extra bonds, angles, dihedrals for all of the bonding possibilities. Also, since you are using 3 fixes to make the same type of bond, a bond can be formed by one fix and not seen by the other two fixes. Each fix counts the bonds per atom at the beginning of the run, then increments as new bonds are added. So one fix could be trying to form a bond with an atom that was already ‘maxed out’ by a different fix.