Long-Range Hydrodynamic Interactions Implemented into LAMMPS (Parallel MD)

We use an explicit solvent method to study the interaction between particles and a fluid. Similar to the Particle-Mesh-Ewald schemes for electrodynamics, the long range hydrodynamic interactions are treated by interpolating the particle density onto a mesh. This is then coupled to the fluid, which we model using a thermal lattice Boltzmann scheme. Mass and momentum conserving noise in the lattice Boltzmann fluid scheme provides a thermostat for both the fluid and the particles. This work has been fully parallelized and implemented into LAMMPS, an open-source molecular dynamics code. We demonstrate the scheme with some examples from colloidal physics and flow over rough surfaces.