I'm building an agent-based model for convection-reaction simulations. Basically my particles are moving at a certain speed in their environment and when they encounter receptors they can bind to them with a certain probability. I want this probability to be a function of their speed -the faster you go the less chance you have to bind-.
To have realistic numbers, I'm looking at binding rates, expressed in second-1 or molecule-1.second-1. Do you have any idea how I could turn these rates into individual probabilities, if it's even possible ? Or if I should look more into having the receptors pick some particles around them every time step, and the slower particles would have a greater cumulated probability of being picked, since they spend more time around the receptor ?
I feel like I might have to look at something like Gillespie algorithm, but not sure how, because I really want to keep the agent-based structure.