Journal of Chemical Physics, Vol.114, No.18, 7989-7992, 2001
Needlelike motion of prolate ellipsoids in the sea of spheres
Molecular dynamics simulations of translational motion of isolated prolate ellipsoids in the sea of spheres have been carried out for several different values of the aspect ratio (kappa), obtained by changing either the length or the diameter of the ellipsoids, at several different solvent densities. The interaction among the spheres is given by the Lennard-Jones pair potential while that between spheres and ellipsoids is given by a modified Gay-Berne potential. Both the mean-square displacements of the center of mass of the ellipsoids and their orientational time correlation function have been calculated. It is found that at short to intermediate times, the motion of ellipsoids is anisotropic and primarily needlelike-the molecules prefer to move parallel to their long axis. The ratio of these two diffusion constants (D-parallel to and D-perpendicular to) approaches kappa, suggesting a decoupling of D-parallel to from the length of the ellipsoid. The diffusion becomes isotropic in the long time with the total diffusion coefficient given by D-parallel to+2D(perpendicular to). The crossover from the anisotropic to the isotropic diffusion is surprisingly sharp and clear in most cases.