Parallelization of Gillespie algorithm based on binary words

We present an improvement of the Gillespie Exact Stochastic Simulation Algorithm, which leverages a bitwise representation of variables to perform independent simulations in parallel. We show that the subsequent gain in computational yield is significant, and it may allow to perform simulations of non-well mixed chemical systems. We illustrate this idea with simulations of Frank model, originally introduced to explain the emergence of homochirality in prebiotic systems.