Brownian dynamics simulation of diffusion-influenced bimolecular reactions

A method is developed and tested for extracting diffusion‐controlled rate constants for condensed phase bimolecular reactions from Brownian dynamics trajectory simulations. This method will be useful when highly detailed model systems are employed, such as those required to explore the complicated range of interactions between enzymes and their substrates. The method is verified by comparing with exact analytical results for simple cases of spheres with uniform reactivity subject to various centrosymmetric Coulombic and Oseen slip hydrodynamic interactions. The utility of the method is illustrated for more complicated cases involving anisotropic reactivity and rotational diffusion.