Mason–Weaver equation explained

The Mason–Weaver equation (named after Max Mason and Warren Weaver) describes the sedimentation and diffusion of solutes under a uniform force, usually a gravitational field.^[1] Assuming that the gravitational field is aligned in the z direction (Fig. 1), the Mason–Weaver equation may be written

(\partial c)/(\partial t) = D (\partial ^{2} c)/(\partial z ^{2}) + sg (\partial c)/(\partial z)

where t is the time, c is the solute concentration (moles per unit length in the z-direction), and the parameters D, s, and g represent the solute diffusion constant, sedimentation coefficient and the (presumed constant) acceleration of gravity, respectively.

D (\partial c)/(\partial z) + s g c = 0 at the top and bottom of the cell, denoted as z _{a} and z _{b} , respectively (Fig. 1). These boundary conditions correspond to the physical requirement that no solute pass through the top and bottom of the cell, i.e., that the flux there be zero. The cell is assumed to be rectangular and aligned with the Cartesian axes (Fig. 1), so that the net flux through the side walls is likewise zero. Hence, the total amount of solute in the cellN _{tot} = \int ^{z _{a}} _{z _{b}} dz c(z, t) is conserved, i.e., dN _{tot/dt} = 0 .

Derivation of the Mason–Weaver equation

f v , the force of gravity m g and the buoyant force \rho V g , where g is the acceleration of gravity, V is the solute particle volume and \rho is the solvent density. At equilibrium (typically reached in roughly 10 ns for molecular solutes), the particle attains a terminal velocity v _{term} where the three forces are balanced. Since V equals the particle mass m times its partial specific volume \bar{\nu} , the equilibrium condition may be written asf v _{term} = m ( 1 - \bar{\nu} \rho) g \stackrel{ def }\ m_b g

where

m _{b} is the buoyant mass.s \stackrel{ def }\ m_b / f = v_\text/g. Since the drag coefficient f is related to the diffusion constant D by the Einstein relationD = (k _{B} T)/(f) ,

the ratio of s and D equals

(s)/(D) = (m _{b})/(k _{B} T)

where

k _{B} is the Boltzmann constant and T is the temperature in kelvins.

The flux J at any point is given by

J = - D (\partial c)/(\partial z) - v _{term} c = - D (\partial c)/(\partial z) - s g c.

The first term describes the flux due to diffusion down a concentration gradient, whereas the second term describes the convective flux due to the average velocity

v _{term} of the particles. A positive net flux out of a small volume produces a negative change in the local concentration within that volume(\partial c)/(\partial t) =- (\partial J)/(\partial z) .

Substituting the equation for the flux J produces the Mason–Weaver equation

(\partial c)/(\partial t) = D (\partial ^{2} c)/(\partial z ^{2}) + sg (\partial c)/(\partial z) .

The dimensionless Mason–Weaver equation

The parameters D, s and g determine a length scale

z _{0} z _{0} \stackrel{ def }\ \frac

and a time scale

t _{0} t _{0} \stackrel{ def }\ \frac

Defining the dimensionless variables

\zeta \stackrel{ def }\ z/z_0 and \tau \stackrel{ def }\ t/t_0, the Mason–Weaver equation becomes(\partial c)/(= \partial \tau) (\partial ^{2} c)/(\partial \zeta ^{2}) + (\partial c)/(\partial \zeta)

subject to the boundary conditions

(\partial c)/(\partial \zeta) + c = 0 at the top and bottom of the cell, \zeta _{a} and \zeta _{b} , respectively.

Solution of the Mason–Weaver equation

This partial differential equation may be solved by separation of variables. Defining

c(\zeta,\tau) \stackrel{ def }\ e^ T(\tau) P(\zeta), we obtain two ordinary differential equations coupled by a constant \beta (dT)/(d \tau) + \beta T = 0 (d ^{2} P)/(d \zeta ^{2}) + \left[ \beta - (1)/(4) \right] P = 0

where acceptable values of

\beta are defined by the boundary conditions(dP)/(d\zeta) + (1)/(2) P = 0

at the upper and lower boundaries,

\zeta _{a} and \zeta _{b} , respectively. Since the T equation has the solution T(\tau) = T _{0} e ^{-\beta} , where T _{0} is a constant, the Mason–Weaver equation is reduced to solving for the function P(\zeta) .

The ordinary differential equation for P and its boundary conditions satisfy the criteriafor a Sturm–Liouville problem, from which several conclusions follow. First, there is a discrete set of orthonormal eigenfunctions

P _{k(\zeta)} that satisfy the ordinary differential equation and boundary conditions. Second, the corresponding eigenvalues \beta _{k} are real, bounded below by a lowest eigenvalue \beta _{0} and grow asymptotically like k ^{2} where the nonnegative integer k is the rank of the eigenvalue. (In our case, the lowest eigenvalue is zero, corresponding to the equilibrium solution.) Third, the eigenfunctions form a complete set; any solution for c(\zeta, \tau) can be expressed as a weighted sum of the eigenfunctionsc(\zeta, \tau) = \sum ^{inf ty} _{k =0} c _{k} P _{k(\zeta)} e ^{- \beta _{k\tau}}

where

c _{k} are constant coefficients determined from the initial distribution c(\zeta, \tau =0 ) c _{k} = \int ^{\zeta _{b}} _{\zeta _{a}} d\zeta c(\zeta, \tau =0 ) e ^{\zeta/2} P _{k(\zeta)}

At equilibrium,

\beta =0 (by definition) and the equilibrium concentration distribution ise ^{-\zeta/2} P _{0(\zeta)} = B e ^{-\zeta} = B e ^{- m _{b} gz/k _{B} T}

which agrees with the Boltzmann distribution. The

P _{0(\zeta)} function satisfies the ordinary differential equation and boundary conditions at all values of \zeta (as may be verified by substitution), and the constant B may be determined from the total amount of soluteB = N _{tot} \left( (sg)/(D) \right) \left( (1)/(e ^{- \zeta _{b}} - e ^{- \zeta _{a}}) \right)

To find the non-equilibrium values of the eigenvalues

\beta _{k} , we proceed as follows. The P equation has the form of a simple harmonic oscillator with solutions P(\zeta) = e ^{i\omega _{k\zeta}} where\omega _{k} = \pm \sqrt{\beta _{k} - (1)/(4 })

Depending on the value of

\beta _{k} , \omega _{k} is either purely real (\beta _{k\geq (1)/(4)} ) or purely imaginary (\beta _{k} < (1)/(4) ). Only one purely imaginary solution can satisfy the boundary conditions, namely, the equilibrium solution. Hence, the non-equilibrium eigenfunctions can be written asP(\zeta) = A \cos{\omega _{k} \zeta} + B \sin{\omega _{k} \zeta}

where A and B are constants and

\omega is real and strictly positive.\rho and phase \varphi as new variables,u \stackrel{ def }\ \rho \sin(\varphi) \ \stackrel\ Pv \stackrel{ def }\ \rho \cos(\varphi) \ \stackrel\ - \frac 1 \omega\left(\frac \right)\rho \stackrel{ def }\ u^2 + v^2\tan(\varphi) \stackrel{ def }\ v / u

the second-order equation for P is factored into two simple first-order equations

(d\rho)/(d\zeta) = 0 (d\varphi)/(d\zeta) = \omega

Remarkably, the transformed boundary conditions are independent of

\rho and the endpoints \zeta _{a} and \zeta _{b} \tan(\varphi _{a)} = \tan(\varphi _{b)} = (1)/(2 \omega _{k})

Therefore, we obtain an equation

\varphi _{a} - \varphi _{b} + k\pi = k\pi = \int ^{\zeta _{a}} _{\zeta _{b}} d\zeta (d\varphi)/(d\zeta) = \omega _{k} (\zeta _{a} - \zeta _{b)}

giving an exact solution for the frequencies

\omega _{k} \omega _{k} = (k\pi)/(\zeta _{a} - \zeta _{b})

The eigenfrequencies

\omega _{k} are positive as required, since \zeta _{a} > \zeta _{b} , and comprise the set of harmonics of the fundamental frequency \omega _{1} \stackrel{ def }\ \pi/(\zeta_a - \zeta_b). Finally, the eigenvalues \beta _{k} can be derived from \omega _{k} \beta _{k} = \omega ^{2} _{k} + (1)/(4)

Taken together, the non-equilibrium components of the solution correspond to a Fourier series decomposition of the initial concentration distribution

c(\zeta, \tau =0 ) multiplied by the weighting function e ^{\zeta/2} . Each Fourier component decays independently as e ^{- \beta _{k\tau}} , where \beta _{k} is given above in terms of the Fourier series frequencies \omega _{k} .

See also

• Lamm equation
• The Archibald approach, and a simpler presentation of the basic physics of the Mason–Weaver equation than the original.^[2]

Notes and References

1. Mason . M . Weaver W . 1924 . The Settling of Small Particles in a Fluid . . 23 . 3 . 412–426 . 10.1103/PhysRev.23.412 . 1924PhRv...23..412M.
2. Archibald . William J. . The Process of Diffusion in a Centrifugal Field of Force . Physical Review . American Physical Society (APS) . 53 . 9 . 1938-05-01 . 0031-899X . 10.1103/physrev.53.746 . 746–752. 1938PhRv...53..746A .