Phase separation kinetics and cluster dynamics in two-dimensional active dumbbell systems

Molecular dynamics simulations were employed to investigate the phase separation process of a two-dimensional active Brownian dumbbell model. We evaluated the time dependence of the typical size of the dense component using the scaling properties of the structure factor, along with the averaged number of clusters and their radii of gyration. The growth observed is faster than in active disk models, and this effect is further enhanced under stronger activity. Next, we focused on studying the hexatic order of the clusters. The length associated with the orientational order increases algebraically with time and faster than for spherical active Brownian particles. Under weak active forces, most clusters exhibit a uniform internal orientational order. However, under strong forces, large clusters consist of domains with different orientational orders. We demonstrated that the latter configurations are not stable, and given sufficient time to evolve, they eventually achieve homogeneous configurations as well. No gas bubbles are formed within the clusters, even when there are patches of different hexatic order. Finally, attention was directed towards the geometry and motion of the clusters themselves. By employing a tracking algorithm, we showed that clusters smaller than the typical size at the observation time exhibit regular shapes, while larger ones display fractal characteristics. In between collisions or break-ups, the clusters behave as solid bodies. Their centers of mass undergo circular motion, with radii increasing with the cluster size. The angular velocity of the center of mass equals that of the constituents with respect to their center of mass. These observations were rationalised with a simple mechanical model.