We carry out a Monte Carlo simulation of stochastic effects for two models of intercellular calcium wave propagation in rat hepatocytes. Both models involve gap junction diffusion by a second messenger. We find that in general the stochastic effects improve agreement with experiment, for a reasonable choice of model parameters. Both stochastic models exhibit baseline fluctuations and variations in the peak heights of Ca$^{2+}$. In addition, we find for one model that there is a distribution of latency times, rather than a single latency time, with a distribution width which is comparable to the experimental observation of spike widths. We also find for the other model with low gap junction diffusion that it is possible for cell multiplets to oscillate independently initially, but to subsequently become synchronized.
Postcript version