Simulation code for: Diffusiophoretic Assembly of Hard Spheres via Reaction-Diffusion Instabilities

This repository contains Fortran source code for fast, particle-resolved simulations of polydisperse hard-sphere colloids in a 2D domain. The biochemical signal is generated by solving reaction–diffusion equations for Turing patterns, which serve as a blueprint for particle motion. Here, we used Brusselator model for the reaction. Using this blueprint, the code directly resolves Lagrangian particles driven by diffusiophoresis and Brownian dynamics, while also undergoing hard-sphere interactions. In addition, the code includes a finite-difference solver for the continuum description of colloidal concentration (advection–diffusion).