Brownian dynamics (BD) is an accurate and computationally efficient mesoscale simulation technique used to study the dynamics and material properties of synthetic and biological polymeric solutions with different architectures, under both equilibrium and nonequilibrium (with the presence of an external field) conditions. A common model used to describe the physics of a polymer is the bead-spring model, where beads resemble the centers of hydrodynamic resistance and they are connected by a network of springs, where the spring tension is generally related nonlinearly to the spring extension. In comparison to a fully resolved molecular dynamics model, this bead-spring model abstracts away the fine (unnecessary) details of real macromoleculaes.
The stochastic differential equation governing the positions of beads in a bead-spring chain is [Öttinger (1989)],
where dt is the time step, Pe is the Peclet number, representing the strength of the flow relative to diffusive motion, is the transposed, gradient of the imposed velocity field, and
is an increment of a 3
dimensional Wiener process. The net non-hydrodynamic, non-Brownian force acting on each bead is
where is all conservative forces, i.e., the net entropic spring force, the excluded volume force preventing beads from intersecting, and bending forces.
is any external force that may be applied. The hydrodynamic interactions (HI) between beads are captured by
which relates forces acting on a bead at position
to a velocity perturbation on a bead at location
, i.e.,
The form of is available analytically for polymers in simple domains such as a polymer in free space or a polymer near a plane wall, and for more complex domains this may be determined numerically. As a result of the fluctuation-dissipation theorem, we require that
The diffusivity matrix is a function of the current bead coordinates
, and thus must be recalculated at every time step of the simulation. As a result, the matrix
must also be recalculated through a decomposition of
as given by the equation given above.