pith. sign in

arxiv: 1906.10680 · v1 · pith:LA2R2L23new · submitted 2019-06-25 · ❄️ cond-mat.soft

Softness suppresses fivefold symmetry and enhances crystallization of binary Laves phases in nearly hard spheres

Tonnishtha Dasgupta , Gabriele M. Coli , Marjolein Dijkstra This is my paper

Pith reviewed 2026-05-25 15:47 UTC · model grok-4.3

classification ❄️ cond-mat.soft
keywords Laves phasesfivefold symmetrysoft potentialsbinary mixturescolloidal crystallizationhard spheresphotonic crystals
0
0 comments X

The pith

Softness in interparticle potentials suppresses fivefold symmetry and enables Laves phase crystallization in binary nearly hard spheres.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes that introducing softness into the repulsive interactions between colloidal particles reduces the amount of fivefold symmetric local ordering in the fluid state of a binary mixture. This reduction allows the system to crystallize more readily into Laves phases, which are otherwise kinetically trapped in hard-sphere systems. A reader would care because these Laves phases serve as precursors for open crystal structures like diamond and pyrochlore that are valuable for photonic materials but difficult to assemble directly.

Core claim

Using computer simulations, the authors show that softness in the interparticle potential suppresses the degree of fivefold symmetry in the binary fluid phase and enhances crystallization of Laves phases in nearly hard spheres, structures that had never been observed to form spontaneously in hard-sphere fluid mixtures due to slow dynamics.

What carries the argument

Softness in the repulsive part of the pair potential, which reduces fivefold symmetry in the fluid and thereby lowers kinetic barriers to Laves ordering.

If this is right

  • Laves phases become accessible via spontaneous crystallization in simulations of binary mixtures with soft potentials.
  • Selective removal of one sublattice from the Laves phase yields diamond and pyrochlore structures.
  • The enhancement holds for nearly hard spheres, indicating the effect applies close to the hard-sphere limit.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same softness-induced symmetry reduction may aid assembly in other frustrated colloidal systems beyond Laves phases.
  • Experiments with microgel particles or other mildly soft colloids could test whether the simulated effect appears in the lab.
  • If finite-size effects are ruled out in larger runs, the mechanism suggests a general design rule for tuning potentials to bypass fivefold traps.

Load-bearing premise

The chosen soft potentials and simulation protocols produce dynamics and local structure representative of real nearly hard-sphere colloids, with the observed reduction in fivefold symmetry as the operative mechanism.

What would settle it

An observation that fivefold symmetry levels remain identical between hard and soft potentials yet crystallization rates still increase would falsify the claimed causal link.

Figures

Figures reproduced from arXiv: 1906.10680 by Gabriele M. Coli, Marjolein Dijkstra, Tonnishtha Dasgupta.

Figure 1
Figure 1. Figure 1: Structure of the three types of Laves phases, showing the different stacking sequences of the large-sphere dimers, marked as ”aa”, ”bb” and ”cc”, when viewed along specific projection planes. The stacking of the large-sphere dimers is (i) “...aa-bb-cc...” for MgCu2, (ii) “...aa-bb...” for MgZn2, and (iii) “...aa-bb-cc-bb...” for MgNi2. packed MgCu2 crystal structure from a binary colloidal dispersion. By s… view at source ↗
Figure 2
Figure 2. Figure 2: Fluid-Laves Phase (LP) coexistence as denoted by the blue region of a binary mixture of WCA spheres with a diameter ratio q = 0.78 at a fixed composition xL = NL/(NL + NS) = 1/3 in the reduced temperature kBT / - reduced density ρσ3 L plane. In the limit of kBT / → 0, the system reduces to a binary mixture of hard spheres. using the topological cluster classification (TCC)19 for varying softness of the i… view at source ↗
Figure 3
Figure 3. Figure 3: Number fraction of particles Nc/N belonging to three different fivefold symmetry clusters as a function of the supersaturation β∆µ of the fluid phase of a binary mixture of WCA spheres at varying temperatures as labeled corresponding to different degrees of particle softness, and of a binary hard-sphere mixture. The three data sets correspond to pentagonal bipyramids (diamonds), defective icosahedra (bulle… view at source ↗
Figure 4
Figure 4. Figure 4: a) Probability distribution to observe a specific cluster topology P(N ) for the five-membered rings for the three ideal LPs, a supersaturated binary hard-sphere (BHS) fluid (β∆µ ' 0.533), and a binomial distribution, all at composition xL = 1/3. The 8 distinct cluster topologies are shown in b) with their index label N . is equilibrated further. Subsequently, the equilibrated configuration is simulated fo… view at source ↗
Figure 5
Figure 5. Figure 5: (a) The largest cluster size NCL with LP symmetry as a function of time t/τMD using the seeding approach in MD simulations of a binary mixture of WCA spheres in the NPT ensemble at temperature T ∗ = 0.2, composition xL = 1/3 and a diameter ratio q = 0.78 for varying pressures βP σ3 L with corresponding supersaturations β∆µ between brackets in order to estimate the critical pressure βPcσ 3 L . The initial s… view at source ↗
Figure 6
Figure 6. Figure 6: a) The height of the Gibbs free-energy barrier β∆Gc and b) the nucleation rate Jσ5 L /DL as a function of the chemical potential difference β∆µ between the fluid and the LP for a binary WCA mixture for the three different LPs and for temperatures kBT / = 0.025, 0.1, and 0.2. coexistence pressure βP σ3 L = 21.32. The results yield some interesting observations. At pressure βP σ3 L = 29, we find absence of … view at source ↗
Figure 7
Figure 7. Figure 7: a) Size of the largest crystalline cluster NCL for a binary mixture of WCA spheres with a diameter ratio q = 0.78 and temperature T ∗ = 0.2 as a function of time t/τMD for varying pressures βP σ3 L with corresponding supersaturations β∆µ between brackets using MD simulations in the NPT ensemble. Profiles are averaged over three independent simulations. b) Configuration of the MgCu2 phase, spontaneously for… view at source ↗
Figure 8
Figure 8. Figure 8: a) The reduced pressure βP σ∗3 L and b) the supersaturation β∆µ versus effective packing fraction η ∗ for a binary hard-sphere mixture and a binary WCA mixture for varying temperatures with a diameter ratio q = 0.78 and composition xL = 1/3. In c) we show the phase diagram of this binary WCA mixture in the reduced temperature kBT / - η ∗ plane. The yellow circles connected by a vertical dashed line denote… view at source ↗
Figure 9
Figure 9. Figure 9: a) The self-intermediate scattering function Fs(q, t) for the large spheres as a function of time tD0/σ2 L for a binary WCA mixture with a diameter ratio q = 0.78 at T ∗ = 0.2 for varying effective packing fractions η ∗ as obtained from MC simulations. b) The structural relaxation time τα as a function of | η ∗ − η ∗ c | /η∗ c for a binary mixture of WCA spheres at T ∗ = 0.025, 0.1, and 0.2 and a binary ha… view at source ↗
read the original abstract

Colloidal crystals with a diamond and pyrochlore structure display wide photonic band gaps at low refractive index contrasts. However, these low-coordinated and open structures are notoriously difficult to self-assemble from colloids interacting with simple pair interactions. To circumvent these problems, one can self-assemble both structures in a closely packed MgCu2 Laves phase from a binary mixture of colloidal spheres and then selectively remove one of the sublattices. Although Laves phases have been proven to be stable in a binary hard-sphere system, they have never been observed to spontaneously crystallize in such a fluid mixture in simulations nor in experiments of micron-sized hard spheres due to slow dynamics. Here we demonstrate, using computer simulations, that softness in the interparticle potential suppresses the degree of fivefold symmetry in the binary fluid phase and enhances crystallization of Laves phases in nearly hard spheres.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 0 minor

Summary. The manuscript claims that softness in the interparticle potential suppresses the degree of fivefold symmetry in the binary fluid phase of nearly hard spheres and thereby enhances spontaneous crystallization into stable MgCu2 Laves phases, which remain kinetically inaccessible in the hard-sphere limit.

Significance. If substantiated, the result would identify a structural route to overcome frustration in binary colloidal crystallization and enable practical self-assembly of Laves phases as precursors to open photonic crystals. The distinction between structural and purely kinetic effects of softness would be a useful contribution to the literature on colloidal self-assembly.

major comments (2)
  1. [Abstract] Abstract: the central claim is presented as a demonstration via computer simulations, yet no methods, potential parameters, system sizes, order parameters for fivefold symmetry, or crystallization metrics are supplied, so the support for the stated mechanism cannot be evaluated.
  2. [Abstract] Abstract/Results: the paper does not describe control simulations that hold diffusivity fixed (e.g., hard-sphere runs with auxiliary forces or temperature rescaling) while varying only the structural measure of fivefold symmetry, leaving the operative role of symmetry suppression untested relative to the kinetic acceleration produced by softness.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for highlighting these points about the abstract and the interpretation of our results. We respond to each major comment below.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the central claim is presented as a demonstration via computer simulations, yet no methods, potential parameters, system sizes, order parameters for fivefold symmetry, or crystallization metrics are supplied, so the support for the stated mechanism cannot be evaluated.

    Authors: Abstracts are concise summaries and do not typically contain full methodological details. The simulation methods, the specific softness parameter in the pair potential, system sizes (several thousand particles), order parameters used to quantify fivefold symmetry (local bond-orientational invariants), and crystallization metrics (time-dependent fractions of particles in Laves environments) are all specified in the Methods section and Results. The abstract therefore summarizes findings whose support is documented in the body of the paper. revision: no

  2. Referee: [Abstract] Abstract/Results: the paper does not describe control simulations that hold diffusivity fixed (e.g., hard-sphere runs with auxiliary forces or temperature rescaling) while varying only the structural measure of fivefold symmetry, leaving the operative role of symmetry suppression untested relative to the kinetic acceleration produced by softness.

    Authors: We agree that an explicit decoupling of structural and kinetic effects would strengthen the mechanistic claim. Our existing data show that fivefold symmetry is already suppressed in the metastable fluid before nucleation events occur, and that this structural change correlates with the onset of Laves crystallization across the range of softness parameters studied. We have added a paragraph in the revised manuscript that discusses this correlation and acknowledges the absence of fixed-diffusivity controls as a limitation of the present study. revision: partial

Circularity Check

0 steps flagged

No circularity: observational simulation results with no derivation chain or fitted predictions

full rationale

The paper reports direct outcomes from molecular dynamics simulations of binary colloidal mixtures under potentials of varying softness. The central claims (suppression of fivefold symmetry in the fluid and enhanced Laves-phase crystallization) are presented as empirical observations from comparing hard-sphere-like and softer potentials. No equations, ansatzes, fitted parameters renamed as predictions, or self-citation chains appear in the abstract or described content. The work contains no mathematical derivation that reduces to its own inputs; results are generated by running the simulations rather than by algebraic or definitional equivalence. This matches the default expectation for a non-circular simulation study.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review supplies no explicit free parameters, axioms, or invented entities; the simulation claim implicitly assumes standard molecular-dynamics methods and a representative soft potential.

pith-pipeline@v0.9.0 · 5690 in / 1009 out tokens · 24196 ms · 2026-05-25T15:47:54.076704+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

42 extracted references · 42 canonical work pages

  1. [1]

    Kalsin, A. M. et al. Electrostatic self-assembly of binary nanoparticle crystals with a diamond-like lattice. Science 312, 420–424 (2006)

    work page 2006

  2. [2]

    J., Chevalier, N

    Bishop, K. J., Chevalier, N. R. & Grzybowski, B. A. When and why like-sized, oppositely charged particles assemble into diamond-like crystals. The Journal of Physical Chemistry Letters 4, 1507–1511 (2013)

    work page 2013

  3. [3]

    V., Talapin, D

    Shevchenko, E. V., Talapin, D. V., Kotov, N. A., O’brien, S. & Murray, C. B. Structural diversity in binary nanoparticle superlattices. Nature 439, 55 (2006)

    work page 2006

  4. [4]

    Evers, W. H. et al. Entropy-driven formation of binary semiconductor-nanocrystal superlattices. Nano Letters 10, 4235–4241 (2010)

    work page 2010

  5. [5]

    in Frontiers in Colloid Science In Memoriam Professor Dr

    Yoshimura, S & Hachisu, S. in Frontiers in Colloid Science In Memoriam Professor Dr. Bun-ichi Tamamushi 59–70 (Springer, 1983)

    work page 1983

  6. [6]

    & Oki, K

    Hasaka, M., Nakashima, H. & Oki, K. Structure of the Laves phase observed in polystyrene latexes. Transactions of the Japan Institute of Metals 25, 65–72 (1984)

    work page 1984

  7. [7]

    H., Fukutomi, T

    Ma, G. H., Fukutomi, T. & Morone, N. Preparation and analysis of ordered structure of binary mixtures composed of poly (4-vinylpyridine) and polystyrene microgels.Journal of Colloid and Interface Science 168, 393–401 (1994)

    work page 1994

  8. [8]

    & Lasnier, B

    Gauthier, J.-P., Fritsch, E., Aguilar-Reyes, B., Barreau, A. & Lasnier, B. Phase de Laves dans la premi` ere opale CT bidisperse. Comptes Rendus Geoscience 336, 187–196 (2004)

    work page 2004

  9. [9]

    Cabane, B. et al. Hiding in plain view: colloidal self-assembly from polydisperse populations. Physical Review Letters 116, 208001 (2016)

    work page 2016

  10. [10]

    Formation of Laves phases in buoyancy matched hard sphere suspensions

    Schaertl, N, Botin, D, Palberg, T & Bartsch, E. Formation of Laves phases in buoyancy matched hard sphere suspensions. Soft Matter 14, 5130–5139 (24 2018)

    work page 2018

  11. [11]

    H., Vermolen, E

    Hynninen, A.-P., Thijssen, J. H., Vermolen, E. C., Dijkstra, M. & Van Blaaderen, A. Self-assembly route for photonic crystals with a bandgap in the visible region. Nature Materials 6, 202 (2007)

    work page 2007

  12. [12]

    & Dellago, C

    Jungblut, S. & Dellago, C. Crystallization of a binary Lennard-Jones mixture. The Journal of Chemical Physics 134, 104501 (2011)

    work page 2011

  13. [13]

    Frank, F. C. Supercooling of liquids. Proc. R. Soc. Lond. A 215, 43–46 (1952)

    work page 1952

  14. [14]

    D., Chandler, D

    Weeks, J. D., Chandler, D. & Andersen, H. C. Role of repulsive forces in determining the equilibrium structure of simple liquids. The Journal of Chemical Physics 54, 5237–5247 (1971)

    work page 1971

  15. [15]

    & Tanaka, H

    Kawasaki, T. & Tanaka, H. Formation of a crystal nucleus from liquid. Proceedings of the National Academy of Sciences 107, 14036–14041 (2010)

    work page 2010

  16. [16]

    & Dijkstra, M

    Filion, L., Ni, R., Frenkel, D. & Dijkstra, M. Simulation of nucleation in almost hard-sphere colloids: The discrepancy between experiment and simulation persists. The Journal of Chemical Physics 134, 134901 (2011). 18

    work page 2011

  17. [17]

    & Speck, T

    Richard, D. & Speck, T. Crystallization of hard spheres revisited. I. Extracting kinetics and free energy landscape from forward flux sampling. The Journal of Chemical Physics 148, 124110 (2018)

    work page 2018

  18. [18]

    & Speck, T

    Richard, D. & Speck, T. Crystallization of hard spheres revisited. II. Thermodynamic modeling, nucleation work, and the surface of tension. The Journal of Chemical Physics 148, 224102 (2018)

    work page 2018

  19. [19]

    R., Eggers, J

    Malins, A., Williams, S. R., Eggers, J. & Royall, C. P. Identification of structure in condensed matter with the topological cluster classification. The Journal of Chemical Physics 139, 234506 (2013)

    work page 2013

  20. [20]

    & Royall, C

    Taffs, J. & Royall, C. P. The role of fivefold symmetry in suppressing crystallization. Nature Communications 7, 13225 (2016)

    work page 2016

  21. [21]

    & Royall, C

    Wood, N., Russo, J., Turci, F. & Royall, C. P. Coupling of sedimentation and liquid structure: Influence on hard sphere nucleation. The Journal of Chemical Physics 149, 204506 (2018)

    work page 2018

  22. [22]

    R., Vega, C., Valeriani, C

    Espinosa, J. R., Vega, C., Valeriani, C. & Sanz, E. Seeding approach to crystal nucleation. The Journal of Chemical Physics 144, 034501 (2016)

    work page 2016

  23. [23]

    R., Sampedro, P., Valeriani, C., Vega, C

    Espinosa, J. R., Sampedro, P., Valeriani, C., Vega, C. & Sanz, E. Lattice mold technique for the calculation of crystal nucleation rates. Faraday Discussions 195, 569–582 (2017)

    work page 2017

  24. [24]

    & McDonald, I

    Hansen, J.-P. & McDonald, I. R. Theory of simple liquids (Elsevier, 1990)

    work page 1990

  25. [25]

    & Ashcroft, N

    Rosenfeld, Y. & Ashcroft, N. Theory of simple classical fluids: Universality in the short-range structure. Physical Review A 20, 1208 (1979)

    work page 1979

  26. [26]

    Generalized Lindemann melting law

    Ross, M. Generalized Lindemann melting law. Physical Review 184, 233 (1969)

    work page 1969

  27. [27]

    R., Costigliola, L., Bailey, N

    Pedersen, U. R., Costigliola, L., Bailey, N. P., Schrøder, T. B. & Dyre, J. C. Thermodynamics of freezing and melting. Nature Communications 7, 12386 (2016)

    work page 2016

  28. [28]

    & Speck, T

    Richard, D. & Speck, T. The role of shear in crystallization kinetics: From suppression to enhancement. Scientific Reports 5, 14610 (2015)

    work page 2015

  29. [29]

    Ni, R., Stuart, M. A. C. & Dijkstra, M. Pushing the glass transition towards random close packing using self- propelled hard spheres. Nature Communications 4, 2704 (2013)

    work page 2013

  30. [30]

    Recent tests of the mode-coupling theory for glassy dynamics

    G¨ otze, W. Recent tests of the mode-coupling theory for glassy dynamics. Journal of Physics: Condensed Matter 11, A1 (1999)

    work page 1999

  31. [31]

    A., Jadrich, R

    Lindquist, B. A., Jadrich, R. B. & Truskett, T. M. Communication: From close-packed to topologically close- packed: Formation of Laves phases in moderately polydisperse hard-sphere mixtures. The Journal of Chemical Physics 148, 191101 (2018)

    work page 2018

  32. [32]

    K., Varela-Rosales, N

    Bommineni, P. K., Varela-Rosales, N. R., Klement, M. & Engel, M. Complex Crystals from Size-Disperse Spheres. Physical Review Letters 122, 128005 (2019)

    work page 2019

  33. [33]

    & Dijkstra, M

    Avvisati, G., Dasgupta, T. & Dijkstra, M. Fabrication of Colloidal Laves Phases via Hard Tetramers and Hard Spheres: Bulk Phase Diagram and Sedimentation Behavior. ACS Nano 11, 7702–7709 (2017). 19

    work page 2017

  34. [34]

    & Pine, D

    Ducrot, ´E., He, M., Yi, G.-R. & Pine, D. J. Colloidal alloys with preassembled clusters and spheres. Nature Materials 16, 652 (2017)

    work page 2017

  35. [35]

    A., Lorenz, C

    Anderson, J. A., Lorenz, C. D. & Travesset, A. General purpose molecular dynamics simulations fully implemented on graphics processing units. Journal of Computational Physics 227, 5342 –5359 (2008)

    work page 2008

  36. [36]

    Glaser, J. et al. Strong scaling of general-purpose molecular dynamics simulations on GPUs. Computer Physics Communications 192, 97 –107 (2015)

    work page 2015

  37. [37]

    J., Tobias, D

    Martyna, G. J., Tobias, D. J. & Klein, M. L. Constant pressure molecular dynamics algorithms. The Journal of Chemical Physics 101, 4177–4189 (1994)

    work page 1994

  38. [38]

    M., Trizac, E, Pronk, S & Frenkel, D

    Polson, J. M., Trizac, E, Pronk, S & Frenkel, D. Finite-size corrections to the free energies of crystalline solids. The Journal of Chemical Physics 112, 5339 (2000)

    work page 2000

  39. [39]

    & Noya, E

    Vega, C, Sanz, E, Abascal, J. & Noya, E. Determination of phase diagrams via computer simulation: methodology and applications to water, electrolytes and proteins. Journal of Physics: Condensed Matter 20, 153101 (2008)

    work page 2008

  40. [40]

    & Smit, B

    Frenkel, D. & Smit, B. Understanding molecular simulation: from algorithms to applications (Elsevier, 2001)

    work page 2001

  41. [41]

    J., Nelson, D

    Steinhardt, P. J., Nelson, D. R. & Ronchetti, M. Bond-orientational order in liquids and glasses. Physical Review B 28, 784 (1983)

    work page 1983

  42. [42]

    Computational Sciences for Energy Research

    Lechner, W. & Dellago, C. Accurate determination of crystal structures based on averaged local bond order parameters. The Journal of Chemical Physics 129, 114707 (2008). Acknowledgements T.D. and M.D. acknowledge financial support from the Industrial Partnership Programme, “Computational Sciences for Energy Research” (Grant no. 13CSER025), of the Netherlan...

    work page 2008