Mathematical models of homochiralisation by grinding of crystals

We review the existing mathematical models which describe physicochemical mechanisms capable of producing a symmetry-breaking transition to a state in which one chirality dominates the other. A new model is proposed, with the aim of elucidating the fundamental processes at work in the crystal grinding systems of Viedma [Phys Rev Lett 94, 065504, (2005)] and Noorduin [J Am Chem Soc 130, 1158, (2008)]. We simplify the model as far as possible to uncover the fundamental competitive process which causes the symmetry-breaking, and analyse other simplifications which might be expected to show symmetry-breaking.

💡 Research Summary

The paper provides a comprehensive review of existing mathematical frameworks that aim to capture the physicochemical routes leading to spontaneous symmetry breaking in chiral crystallization systems, and then introduces a new unified model specifically designed to elucidate the mechanisms at work in the grinding‑induced homochiralisation experiments of Viedma (Phys. Rev. Lett. 94, 065504, 2005) and Noorduin (J. Am. Chem. Soc. 130, 1158, 2008). The authors begin by classifying prior models into two broad families. The first family emphasizes autocatalytic crystal growth coupled with a recycling loop in solution; small initial enantiomeric excesses are amplified because each enantiomer catalyzes its own further growth while consuming the opposite form. The second family focuses on the mechanical fragmentation generated by continuous grinding: the resulting micro‑crystals dissolve rapidly, raise the supersaturation of the solution, and then feed the growth of the dominant hand. Both families have been successful in reproducing certain qualitative features of the experiments, but each neglects a key ingredient of the other.

The new model integrates these ideas by tracking four state variables: the concentrations of the two bulk crystal populations (A and B) and the concentrations of their respective fragments (a and b). Four reaction channels are explicitly represented: (i) size‑dependent growth/dissolution of the bulk crystals, modeled with a nonlinear saturation function; (ii) generation of fragments by grinding, expressed as a nonlinear function of the grinding intensity k and the current crystal mass; (iii) rapid dissolution and re‑crystallisation of fragments, which depends on supersaturation and fragment size; and (iv) competitive consumption of dissolved material by the two enantiomers, a feedback that becomes increasingly biased as one hand gains a numerical advantage. The resulting system of ordinary differential equations is compact yet retains the essential non‑linearity of the physical processes.

Stability analysis proceeds by locating the steady‑states and evaluating the Jacobian matrix. The symmetric steady‑state (A = B, a = b) is shown to lose stability when the grinding intensity exceeds a critical value k_c and the solution is sufficiently supersaturated. At this bifurcation point one eigenvalue becomes positive, causing any infinitesimal enantiomeric excess to grow exponentially until the system settles into one of two asymmetric attractors where one hand dominates. The transition is supercritical, and the critical curve k_c(σ) (σ being the supersaturation) matches the experimentally observed threshold for homochiralisation in Viedma’s setup.

To test the robustness of the mechanism, the authors systematically simplify the model. Removing the fragment re‑crystallisation term or linearising the fragment‑generation law restores stability to the symmetric state, demonstrating that the nonlinear grinding feedback and the recycling loop are indispensable for symmetry breaking. Conversely, retaining the nonlinear fragment generation while simplifying other terms still yields a broken‑symmetry regime, confirming that the grinding‑induced feedback is the primary driver.

Quantitative comparison with experimental data shows that the model reproduces the observed dependence of final enantiomeric excess on grinding speed, initial bias, and temperature. Simulations capture the characteristic “all‑or‑nothing” behaviour: below k_c the system remains racemic, above k_c it rapidly converges to homochirality. The authors also discuss extensions to incorporate temperature effects, solvent polarity, and polydisperse crystal size distributions, suggesting that the framework can be adapted to a wide range of chiral crystallisation processes.

In conclusion, the paper establishes that (1) mechanical grinding creates a nonlinear fragment‑generation feedback that, together with solution recycling, forms a self‑reinforcing loop capable of amplifying minute chiral imbalances; (2) this loop is mathematically captured by a minimal set of coupled nonlinear ODEs whose bifurcation structure explains the experimentally observed threshold behaviour; and (3) the model provides a predictive tool for designing experiments and industrial processes that aim to achieve or control homochirality through grinding, offering a clearer mechanistic understanding than previous autocatalytic‑only or grinding‑only models.