Hi, LAMMPS users:
I found a bug maybe exists in the fix bond/create command. I have test fix bond/create and fix bond/break in a minimal model system, which is described as follows:
Two parallel polymer chains are set in the test model, each chain include 10 atoms, and they apart from 15 distance units. One of the chain include an atom site which can dynamically bonding and unbonding a new type bond with other atoms in another chain if distance constrains are satisfied.
I used extra/bond/per/atom 1 and extra/special/per/atom 25 to control the bonding number at the atom site, if the atom site has formed a bond with any atoms in another chain, it cannot form any bond furthermore. Compute property/local and dump/local were used to dump out the bond information. I found that:
(1) Extra/bond/per/atom 1 cannot constrain the numbers of new type bond no more than one at the atom site which formed with the atoms in another chain.
(2) If I use iparam and jparam parameters to constrain the bond number in the atom site, I found the information of bond/break’s cannot tell the command bond/create that the bonds have broken, and the atom site can form new bonds again. So, the maxbond in the iparam or jparam, it can only used once, if the bonds have broken, the maxbond can not update information.
The version of LAMMPS used in this test is 5May 2020.
The follows are input file, data file and output Bond list file.
Input file: