Simulations of concentrated suspensions of rigid fibers: Relationship between short-time diffusivities and the long-time rotational diffusion

Abstract

Brownian dynamics simulations of the behavior of suspensions of fibers demonstrate that the scaling of the rotational diffusivity with respect to the number density (nL3) is a sensitive function of the thickness and the parameter L2D(R0)/D(T0), where D(R0) is the rotational diffusivity at infinite dilution, D(T0) is the average center-of-mass diffusivity at infinite dilution, and L is the fiber length. Existing theories for the long-time rotational diffusivities of rigid fibers in the semidilute and concentrated regimes fail to accurately account for the relationship with the dilute values of the rotational and translational diffusivities of the various physical models used to simulate the fibers. The concentration regime studied in this work ranges from a number density of nL3 approximately 0-150, which is below the transition from an isotropic to nematic state. The effect of the fiber thickness was studied by performing simulations of rods with aspect ratios (fiber length over diameter) of 25, 50, and 500, as well as performing projections for infinitely thin fibers. The excluded volume of the rods was enforced through the use of short-range potentials. For a rod with an aspect ratio of 50 with a parameter of L2D(R0)/D(T0)=9, which corresponds to a slender-body model of the individual fibers, the rotational diffusivity (D(R)) scales as D(R)/D(R0) approximately (nL3)(-1.9) in the concentration regime of 70 < or = nL3 < or = 150. Similarly with a parameter of L2D(R0)/D(T0)=4, corresponding to a rigid-dumbbell model, the rotational diffusivity scales as D(R)/D(R0) approximately (nL3)(-1.1) over the same range of concentrations. For rods with aspect ratios of 25, it is observed that a difference in the scaling is seen for L2D(R0)/D(T0) approximately < 8, with higher values of this ratio exhibiting essentially the same scaling. Additional values of the ratio L2D(R0)/D(T0) were investigated to determine the overall behavior of the suspension dynamics with respect to this parameter. These findings resolve discrepancies between simulation results for rotational diffusivities reported by previous investigators and provide new insights for the development of an accurate theory for the diffusivity of rigid rods suspended in solution.