We derive an expression for the kinetic coefficient of growth of crystals from solution as diffusion over an energy barrier U, followed by unimpeded incorporation (73,74). The barrier U may be of electrostatic origin (74); however, for the ferritin/apoferritin pair, it probably accounts for the need to expel the water molecules attached to the incoming molecules and to the growth site (82). Repulsive potentials owing to water structuring at hydrophobic and hydrophilic surface patches can have significant strength and range (83,84).

In Fig. 12, the resulting potential relief is schematically depicted. The potential reaches its maximum value Umax at the crossing of the increasing branch, owing to the repulsion between the incoming solute molecules and the crystal surface at medium separations, and the receding branch, which corresponds to the short-range attraction required if the molecules should enter the growth site. We position the beginning of the coordinate axis x = 0 at the location of this maximum. We assign a finite curvature of U(x) about this maximum and link it with the expulsion of the last few solvent molecules as an incoming solute molecule joins the crystal. The finite curvature assumption follows previous solutions to similar problems (85,86). The distance d, used as the upper integration limit, is bound from below by the range of interaction of the solute molecules with the surface, which can be a few solute molecular sizes, and from above by the distances between the solute molecules in the solution bulk, n-1/3 ^ 0.2 mm. Thus, d can be chosen significantly longer than the molecular sizes and kink lengths, and using a one-dimensional model is justified. Since the rate of diffusion over a sharp barrier only depends on the curvature around the maximum (85,86), the choice of d does not affect the result

To calculate the flux J of molecules with concentration n that, driven by a concentration gradient, overcome a barrier to reach the surface, we orient the coordinate x perpendicular to a growing surface and denote the potential relief close to this surface as U(x). From the generalized Fick law, J = (nD/kBT)dyddx, with |(T,x) = |0(T) + kBT ln[yn(x)] + U(x) and 1 (10), J is linked to U(x), n(x) and the gradient of n as (73,86)

dx dx I

with D being the Stokes diffusion coefficient of the molecules. In search of a steady J = const, we integrate Eq. 17 with two sets of boundary conditions:

(1) that at a certain distance from the surface S, x >8, U = 0, and n = n8; and

Dividing by D and multiplying both sides by exp[U(x)/kBT], we get

J exp[U(x)/kT] = — n(x) exp[U(x)/k„T] D B dx\ B

Integrating from x = 0 to x = 8, using the boundary conditions at x = 0 and x >8, we get

an analog to equation 9.51 in ref. 74 and the Fuchs expression for coagulation of particles interacting through U(x).

If U(x) has a sharp maximum at x = 0, we can represent it with a symmetric function around the point of the maximum. As shown below, in many cases \d2U/dx2\ < a, and this justifies the assumption of a sharp maximum. We use only the first two members of the Taylor series: U(x) = Umax - 1/2\d2U/dx2\x2. The minus sign stems from d2U/dx2 < 0 at the maximum. Then, the integral

" kBT

" kBT

2 |
d2( U/kBT) | |||||

n |
dx2 |
The approximate equality above is based on 5 >> [1/2\d2U/dx2\] 1/2, the half-width of the Gaussian function in Eq. 20 Finally, " kBT Note that only half of the flux J from Eq. 21 contributes to growth: on top of the barrier, the force driving the molecules into the crystal is zero, and a molecule has equal chances of getting incorporated, or going back to the solution (87). With this, and introducing the parameter L as the radius of curvature of U(x)/kBT at its maximum
Was this article helpful? |

## Post a comment