Abstract
The random motion of the diffusing reactants in continuum space was simulated directly based on the first-passage algorithm. As a result, accurate results for the reaction rate of a cluster of spherical sinks were obtained. The asymptotic expressions for the reaction rate varying with the number of sinks were confirmed for specified configurations.