Combining 'fix bond/react' and 'if'

Not sure if that is of any help, but I remember that I had a small modification of of bond/react to accept lammps variables as input to the reaction probability. But this was a per reaction and not a per atom value if I remember correctly. I can look that up if necessary.

Regards
Wolfgang

Hi Wolfgang,

I’m not sure that capability helps in the current case, however, it sounds useful. Would you be willing to submit a pull request on Github, or send to me?

Thanks,
Jake

Hi Jake

Good opportunity to learn a little bit on how to use git. I’ll try myself first.

Regards

Wolfgang