[lammps-users] Domino effect with bonds

Hi everyone,

I’m trying to model a system which starts with a bond between the right brown bead and the right blue bead. Then,

  • when the left brown bead gets closer to the blue one, a bond forms between the left brown and the left blue bead with a probability of 30%
    -when there is a bond between the left blue bead and the left brown bead, then we break the bond between the right brown bead and the right blue bead.

Is this possible to do? If so, suggestions are more than welcome.

Thanks!

Untitled.png

This is probably possible using fix bond/react (https://lammps.sandia.gov/doc/fix_bond_react.html).

Untitled.png