Elastic fluctuations as observed in a confocal slice

Recent confocal experiments on colloidal solids motivate a fuller study of the projection of three-dimensional fluctuations onto a two-dimensional confocal slice. We show that the effective theory of a projected crystal displays several exceptional features, such as non-standard exponents in the dispersion relations. We provide analytic expressions for the effective two-dimensional elastic properties which allow one to work back from sliced experimental observations to three-dimensional elastic constants.

💡 Research Summary

The paper addresses a practical problem that arises in modern confocal microscopy of colloidal crystals: experimentalists often record particle motions only within a single two‑dimensional (2D) focal plane, typically the (111) plane of a three‑dimensional (3D) crystal, because fast scanning of the full volume is technically demanding. The authors develop a rigorous theoretical framework that connects the elastic fluctuations observed in such a 2D slice to the underlying 3D elastic constants of the bulk material.

Starting from the standard linear elasticity description of a cubic crystal, the elastic energy is expressed in terms of the displacement field u_i and its gradient tensor u_{ij}. The elastic matrix D̂_{ik}(k) (Eq. 1) contains the Lamé coefficients λ, the shear modulus μ, and an anisotropy parameter ν. Its inverse defines the static elastic Green function Ĝ_{ij}(k) (Eq. 2), which governs the thermal correlation ⟨u_i(k)u_j(−k)⟩ = β^{−1}Ĝ_{ij}(k) (Eq. 4). For each wave vector k the matrix D̂ has three eigenvalues d_i(k); the auxiliary variables ω_i^2(k)=d_i(k) play the role of squared “frequencies” in the harmonic theory.

To obtain the effective theory for the observed slice, the authors project the 3D Green function onto the plane orthogonal to N = (1,1,1)/√3. In real space the projected Green function Q_{αβ}(x) (Eq. 8) retains an almost isotropic form with scalar prefactors C_0(x) and C_{00}(x) that are determined by contracting the full 3D Green function with the projection operators δ_{ij}−N_iN_j and \bar{x}i\bar{x}j/x^2 (Eqs. 9a‑b). Fourier transforming Q{αβ}(x) yields a 2D elastic matrix \hat D{αβ}(q) whose inverse \hat Q_{αβ}(q) has the structure (Eq. 10) with coefficients that scale as 1/q. This scaling is the origin of the anomalous dispersion: unlike standard 3D elasticity where ω∝q, the projected system exhibits ω∝√q.

Because the exact anisotropic Green function is analytically intractable except for high‑symmetry directions, the authors adopt FedoroV’s isotropic approximation. They define effective isotropic Lamé parameters λ̃ = λ+ν/5 and μ̃ = μ+ν/5 (Eq. 11) and replace the full Green function by the isotropic form (Eq. 12). This yields simple expressions for the scalar coefficients A_0, A_{00} and, after projection, for the 2D coefficients E_0 and E_{00}. Consequently, the effective dispersion relations for the slice become (Eq. 14):

• Transverse mode: ω_⊥^2 = 2 μ̃ q
• Longitudinal mode: ω_∥^2 = 4 μ̃(λ̃+2μ̃)/(λ̃+3μ̃) q

Both branches scale linearly with the magnitude of the 2D wave vector q, i.e., ω∝√q, a non‑Debye behavior that is a direct consequence of dimensional reduction.

The theoretical predictions are tested against large‑scale molecular dynamics simulations of a face‑centered cubic (FCC) crystal containing 4 147 200 particles at volume fraction φ = 0.57. The full 3D dispersion curves (ω∝k) are used to extract the bulk elastic constants λ, μ, ν = 42.8, 51.8, −53.8 (in units where particle diameter, mass, and k_BT are unity). The same simulation data are then sliced along the (111) plane, and the 2D displacement correlations are analyzed. The measured 2D dispersion follows the √q law, and the numerical prefactors agree with the analytic formulas (Eqs. 11‑14) within 3 %. Additional simulations of a 2D hexagonal lattice projected onto a line, where the full isotropic 2D elasticity can be applied without approximation, confirm the same ω^2∝q scaling.

The paper concludes that the projection of 3D elastic fluctuations onto a confocal slice leads inevitably to anomalous dispersion relations, and that the derived analytic expressions provide a practical route for experimentalists to infer bulk elastic constants from purely 2D measurements. This framework is broadly applicable to other systems where a lower‑dimensional observation window is imposed, such as Stokesian hydrodynamics or thin‑film mechanics, and it clarifies that non‑standard exponents observed in sliced data do not necessarily signal exotic or glassy behavior but are a generic consequence of dimensional reduction.