Chiral symmetry breaking via crystallization of the glycine and alpha-amino acid system: a mathematical model

We introduce and numerically solve a mathematical model of the experimentally established mechanisms responsible for the symmetry breaking transition observed in the chiral crystallization experiments reported by I. Weissbuch, L. Addadi, L. Leiserowitz and M. Lahav, J. Am. Chem. Soc. 110 (1988), 561-567. The mathematical model is based on five basic processes: (1) The formation of achiral glycine clusters in solution, (2) The nucleation of oriented glycine crystals at the air/water interface in the presence of hydrophobic amino acids, (3) A kinetic orienting effect which inhibits crystal growth, (4) The enantioselective occlusion of the amino acids from solution, and (5) The growth of oriented host glycine crystals at the interface. We translate these processes into differential rate equations. We first study the model with the orienting process (2) without (3) and then combine both allowing us to make detailed comparisons of both orienting effects which actually act in unison in the experiment. Numerical results indicate that the model can yield a high percentage orientation of the mixed crystals at the interface and the consequent resolution of the initially racemic mixture of amino acids in solution. The model thus leads to separation of enantiomeric territories, the generation and amplification of optical activity by enantioselective occlusion of chiral additives through chiral surfaces of glycine crystals.

💡 Research Summary

The paper presents a quantitative kinetic model that reproduces the spontaneous mirror‑symmetry breaking observed in the glycine–α‑amino‑acid crystallization experiments of Weissbuch, Addadi, Leiserowitz and Lahav (JACS 1988). The authors identify five elementary processes that are experimentally documented: (1) formation of achiral glycine clusters in solution, (2) nucleation of oriented glycine crystals at the air–water interface in the presence of hydrophobic amino‑acid “templates”, (3) a kinetic inhibition effect exerted by soluble (hydrophilic) amino acids on the growth of embryonic nuclei, (4) enantioselective occlusion of the amino‑acid enantiomers into the growing crystal, and (5) subsequent growth of the oriented host crystal by uptake of glycine monomers.

These processes are translated into a set of ordinary differential equations. Glycine monomers (A₁) and clusters (A_r) are assumed to combine irreversibly, with an average seed size M that represents the typical nucleus. The nucleation step (process 2) is modeled as a first‑order reaction in which a single hydrophobic L‑ or D‑amino‑acid molecule at the interface triggers the formation of a crystal seed with a specific orientation: X (exposing the (010) face) or Y (exposing the (0 ¯10) face). Although real nucleation likely involves many molecules, the authors adopt a minimal “one‑template” assumption to keep the model tractable.

Enantioselective occlusion (process 4) is encoded by two coupled reactions: the X‑oriented crystal incorporates only D‑amino‑acid monomers, while the Y‑oriented crystal incorporates only L‑amino‑acid monomers. Each occlusion event adds a single amino‑acid unit to the crystal lattice, reflecting crystallographic observations that monomers, not clusters, are incorporated one at a time. The growth step (process 5) adds glycine monomers to the crystal face exposed to the solution, preserving the orientation already set by the template.

The kinetic inhibition (process 3) is introduced as a competing first‑order term that selectively slows the growth of nuclei whose exposed face is opposite to the predominant hydrophobic template. This term captures the experimentally observed “kinetic effect” whereby soluble amino acids inhibit the expansion of unfavored nuclei, thereby reinforcing the orientation bias introduced by the hydrophobic effect.

Two numerical experiments are performed. In the first, only the hydrophobic orientation effect (process 2) is active; the model shows that the racemic state is linearly unstable, and any infinitesimal excess of one enantiomer leads to the dominance of the corresponding crystal orientation. In the second, both the hydrophobic and kinetic inhibition effects are combined; the simulations reveal a stronger amplification of the favored orientation, with the fraction of crystals adopting the dominant face reaching 80–90 % under realistic parameter values.

Parameter sweeps demonstrate that the model’s behavior is robust to variations in template concentration, inhibition strength, and average seed size. Even a modest initial imbalance (e.g., L : D = 10 : 9) suffices to trigger a cascade in which the oriented crystals continuously occlude the same enantiomer, progressively depleting it from solution and enriching the opposite enantiomer. The resulting solution becomes highly enantiomerically enriched (up to >95 % of one handedness), in agreement with the experimental observations of optical activity generation without mechanical stirring or grinding.

The authors acknowledge simplifications: no fragmentation, no multi‑cluster aggregation, and a single‑type template assumption. Nevertheless, the model captures the essential feedback loop—template‑induced orientation, enantioselective occlusion, and kinetic inhibition—that drives symmetry breaking. By providing a minimal yet quantitative framework, the work bridges the gap between phenomenological chemistry experiments and dynamical systems theory, offering insight into how small stochastic fluctuations could be amplified into macroscopic homochirality in prebiotic environments. The study thus contributes a valuable theoretical tool for researchers investigating the origins of biological homochirality and the role of interfacial crystallization in early Earth chemistry.