Molecular Dynamics Simulation of the Motion of Colloidal Nanoparticles in a Solute Concentration Gradient and a Comparison to the Continuum Limit

Chemical-mechanical transduction mechanisms which can actuate the movement of colloids through liquids are highly sought after as engines to propel miniaturized micro-and nanobots. One mechanism involves harnessing the long-range van der Waals attractive forces between the colloid and solute molecules dissolved in the liquid around the particle. If a concentration gradient of this solute is applied across the particle, then the imbalance in the van der Waals attraction drives the particle towards the higher concentration of solute. We present a molecular dynamics simulation using Lennard-Jones interactions between molecules of the solvent, solute, and colloid cluster which include short-range repulsive and long-range attractive potentials. The simulations demonstrate that a solute gradient can propel nanosized colloids, and that the velocity decreases with the colloid size. The solute-colloid short-range repulsive interactions are observed to be restricted to a region of specifically adsorbed solutes on the particle surface which are symmetrically adsorbed and do not contribute to the motion. The size of this region provides a cutoff for a continuum level description of the motion, and with this cutoff, continuum calculations are in excellent agreement with the molecular dynamics simulation results, completing a description of the propulsion from the nano-to the microscale.