Static Structure and the Bulk Modulus

In addition to gauging long-ranged order, the static structure factor

where is the mean number density of particles at temperature

Errors in *S*(*q*) due to the limited extent of the field of view
were minimized by applying the Hann windowing function
to the structural data[35] before calculating *S*(*q*).
We estimated the resolution limit for the measured
structure factors by calculating
*S*(*q*) for simulated distributions of points using
the same lattice constants and vacancy concentrations
as the actual crystals and including normally distributed
random displacements large enough to account
for the crystals' Lindemann factors.
These simulated defects constitute small departures from an otherwise
perfect and thus ideally rigid lattice.
The crystals' long-wavelength
structure factors should exceed the simulated results
because of order-disrupting elastic and plastic deformations.
We also expect the experimental structure factors to
verge below estimates based on the bulk moduli calculated
in Section 5.
The anticipated ranges of values for
appear as shaded regions in Fig. 7.

As can be seen in Fig. 7(a), data for the dilute, weakly interacting crystal fall quite close to the anticipated range. The small discrepancy might be accounted for by the vacancies in the field of view. Given the spheres in the field of view, these would tend to depress the bulk modulus by a bit more than two percent both because of the reduction in density and because of the loss of supporting interactions. With these caveats, the linear interaction theory appears to be consistent with both the dynamics and the elasticity of the dilute crystal.

The long-wavelength behavior of *S*(*q*) for the dense,
comparatively rigid crystal, on the other
hand, falls significantly
outside the anticipated range.
In fact, the bulk modulus corresponding to
is only
,
less than a third
of the expected value.
In other words,
the crystal is much softer than expected.
At least part of this difference can be attributed to
structural defects such as vacancies.
This cannot be the whole story, however.
The rigid crystal has only
vacancies
out of
spheres.
Disorder's effect on the dense crystal's bulk modulus should be proportionately
smaller than it was in the dilute crystal.
And yet the discrepancy is much larger.

The discrepancy for the dense crystal also cannot be ascribed
to inadequacies in the charge renormalization formula in
Eq. (25).
If this were the case, then alternative values of *Z* and could be found which would be consistent with both the
dynamics in the potential of mean force and the elasticity of the crystal.
No such values can be found.
Combining Eq. (24) for the potential
of mean force with Eq. (30) for the bulk modulus
yields physically meaningful values for
only if

(33) |

Given the measured values of , the smallest possible bulk modulus for the rigid crystal is

Finally, we cannot simply dismiss the DLVO theory and seek consistency in other models for the pair potential. The DLVO theory successfully describes the electrostatic interactions between isolated pairs of the spheres used in the present study [12]. It very nearly accounts for the elastic properties of the more dilute crystal. When viewed in this light, its inability to describe the dense crystal's properties indicates not just a failure of the DLVO theory, but rather a failure of linear superposition. This is more serious since the linear superposition approximation underlies all linearized theories for colloidal electrostatic interactions. Not surprisingly, it seems to fail when used to describe strongly interacting many-body systems. Under these conditions, the nonlinearity ignored in the derivation of Eq. (21) apparently results in non-additive electrostatic interactions which significantly affect the crystals' elastic properties. This interpretation is consistent with recent numerical [36] and analytical [29] studies of nonlinear coupling in charge-stabilized suspensions. Such many-body effects should be less important in more weakly coupled systems. Thus, the dilute crystal's properties are adequately described by the linear theory while the dense crystal's are not. Perhaps the most troubling point is that there seems to be no simple way to anticipate when the linearized description will suffice and when it will not.