Tuesday, July 7, 2015

diffusion - Expression for "rotational diffusivity"; orientation random-walk of thin rod-like particles?



From this answer and from the Stokes-Einstein equation the diffusivity of a particle of radius R in a fluid of viscosity η is


D=kBT6πηR


where ξ=6πηR is a coefficient of friction Stokes' law such that for velocity v the viscous drag force is


FD=ξv.


The diffusion in 1 dimension would then be given as


x2=2Dt.


As suggested in this answer to that question, the diffusivity is usually so strongly limited by the drag force that it does not depend much on the density of the particle.


Question: Is there an orientation analog to the positional diffusion? For example, if the particle were a long thin rod, would the direction of its axis move in a random-walk type process? If so, is there an analogous orientation diffusivity, perhaps something like Drot where θ2=2Drott based on the length of the thin rod and likewise independent of particle density?


update: As pointed out by @KyleKanos there is indeed the concept of "rotational diffusivity". This question mentions that for a spherical particle:


DrotkBTζfkBT(8πη)(r)3



and @KyleKanos's answer begins to describe how it can be used.


What I'm looking for here instead is Drot for a long thin rod.



Answer



Rotational diffusion was solved by Edwardes / Perrin / Langevin at their times. It is normally made for particles of Ellipsoid shape, since it is easier to solve analitycally. A long and thin ellipsoid is your rod. Try Google "ellipsoid brownian motion". The model is for macroscopic particles immersed in a viscous fluid, but should work also for long molecules in a gas.


Solution, thin ellipsoid of axis a (roughly equivalent to 2a rod) and (small) radius b:


D=kTCr


Where Cr is the friction for the ellipsoid rotating in the fluid; for our thin rod,


Cr16πη3Pr


Where again, the Pr is an elliptic integral that represent the particle size, at the moment I don't manage to calculate:


Pr=0dx(a2+x)32(b2+x)



All this discussion work for orientation the mayor axis of a thin rod (rotation around the two perpendicular axes); if the particle has a stranger shape, I guess we must consider all components and tensors.


Here some paper where the model si discussed in english: Paper Wegener et al, A general ellipsoid cannot always .. Koening, Brownian motion of an ellipsoid And the book The Langevin equation by Coffey has a chapter. Wikipedia Perrin friction factors


No comments:

Post a Comment

classical mechanics - Moment of a force about a given axis (Torque) - Scalar or vectorial?

I am studying Statics and saw that: The moment of a force about a given axis (or Torque) is defined by the equation: $M_X = (\vec r \times \...