Abstract
We study the pair dynamics of two self-propelled sphere dimers in the chemically active medium in which a cubic autocatalytic chemical reaction takes place. Concentration gradient around the dimer, created by reactions occurring on the catalytic sphere surface and responsible for the self-propulsion, is greatly influenced by the chemical activities of the environment. Consequently, the pair dynamics of two dimers mediated by the concentration field are affected. In the particle-based mesoscopic simulation, we combine molecular dynamics (MD) for potential interactions and reactive multiparticle collision dynamics (RMPC) for solvent flow and bulk reactions. Our results indicate three different configurations between a pair of dimers after the collision, i.e., two possible scenarios of bound dimer pairs and one unbound dimer pair. A phase diagram is sketched as a function of the rate coefficients of the environment reactions. Since the pair interactions are the basic elements of larger scale systems, we believe the results may shed light on the understanding of the collective dynamics.