pith. sign in

arxiv: 2605.12870 · v2 · pith:ES64IQ25new · submitted 2026-05-13 · ❄️ cond-mat.soft · cond-mat.stat-mech· physics.chem-ph

PACSim: A Flexible Simulation Framework for Polymer-Attenuated Coulombic Self-Assembly

Pith reviewed 2026-05-20 22:08 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.stat-mechphysics.chem-ph
keywords polymer-attenuated Coulombic self-assemblymolecular dynamics simulationcolloidal crystalsOpenMM frameworkcharged colloidspolymer brushself-assembly simulation
0
0 comments X

The pith

PACSim is an OpenMM-based framework for molecular dynamics simulations of polymer-attenuated Coulombic self-assembly in charged colloids.

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

The paper presents PACSim as an open-source tool that lets researchers run molecular dynamics simulations of how charged spherical particles coated with polymer brushes assemble into crystals. In PACS, crystallization depends on particle concentration, charge, size, and the salt level in the surrounding solution. The framework supports a range of these experimentally relevant conditions while allowing different interaction potentials and connections to enhanced sampling or machine-learning tools. A sympathetic reader would value the ability to predict assembly outcomes and inspect particle-level mechanisms before or alongside physical experiments.

Core claim

PACSim enables MD simulation studies of assembly by PACS across a range of experimentally relevant scenarios. It is built on top of OpenMM to support the implementation of different interaction potentials as well as integration with other tools such as enhanced-sampling and machine-learning frameworks, and it provides particle-level insight into the assembly processes.

What carries the argument

PACSim, an open-source MD simulation framework built on OpenMM that implements interaction potentials for polymer-coated charged colloids and enables flexible exploration of assembly conditions.

If this is right

  • Researchers can test how changes in colloid concentration, charge, size, or salt concentration affect whether crystals form and which structures appear.
  • The same framework supports studies that combine standard MD with enhanced sampling methods to reach longer assembly timescales.
  • Integration with machine-learning tools allows training on simulation trajectories to identify key assembly mechanisms or guide further simulations.
  • Methodological improvements, such as new ways to represent polymer-brush repulsion or electrostatic screening, can be developed and tested inside the PACSim environment.

Where Pith is reading between the lines

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

  • Virtual parameter sweeps in PACSim could reduce the number of trial-and-error synthesis runs needed to discover new colloidal crystal lattices.
  • The modular design suggests straightforward extension to related self-assembly problems that also combine electrostatic attraction with steric repulsion.
  • Coupling PACSim outputs to automated experimental feedback loops would create closed-loop design of colloidal materials.

Load-bearing premise

The coarse-grained interaction potentials in PACSim capture enough of the physics of polymer-coated charged colloids that simulated assembly results remain predictive of laboratory experiments.

What would settle it

Direct comparison of crystal structures or assembly pathways produced by PACSim for a chosen set of colloid concentration, charge, size, and salt concentration against the structures observed in the corresponding physical PACS experiment.

Figures

Figures reproduced from arXiv: 2605.12870 by Glen M. Hocky, John P. Marquardt, Michael S. Chen, Nicole Smina, Philipp H\"ollmer, Stefano Sacanna, Steven van Kesteren.

Figure 1
Figure 1. Figure 1: FIG. 1. Tuning parameters of the PACS potential [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. (A) The [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Specification of a cluster template in a LAMMPS [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Configuration file for [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Specification of a crystal structure in a CIF file that [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Configuration file for [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Configuration file for [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Example Python implementation of the electrostatic [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Configuration file excerpt for [PITH_FULL_IMAGE:figures/full_fig_p010_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Cluster sizes during an MD simulation seeded with [PITH_FULL_IMAGE:figures/full_fig_p010_12.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Stability of crystalline seeds during [PITH_FULL_IMAGE:figures/full_fig_p010_10.png] view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14. Late-time snapshots from MD simulations with [PITH_FULL_IMAGE:figures/full_fig_p011_14.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13. Configuration file excerpts for adding a substrate ex [PITH_FULL_IMAGE:figures/full_fig_p011_13.png] view at source ↗
Figure 16
Figure 16. Figure 16: FIG. 16. Debye length [PITH_FULL_IMAGE:figures/full_fig_p012_16.png] view at source ↗
Figure 15
Figure 15. Figure 15: FIG. 15. Configuration file excerpt for [PITH_FULL_IMAGE:figures/full_fig_p012_15.png] view at source ↗
Figure 18
Figure 18. Figure 18: FIG. 18. Average coordination number during a well [PITH_FULL_IMAGE:figures/full_fig_p013_18.png] view at source ↗
Figure 17
Figure 17. Figure 17: FIG. 17. MD step time as a function of number of particles [PITH_FULL_IMAGE:figures/full_fig_p012_17.png] view at source ↗
read the original abstract

Polymer-Attenuated Coulombic Self-Assembly (PACS) is a flexible experimental approach for generating crystals from simple colloidal building blocks. The central components are charged spherical particles coated with a polymer brush that prevents irreversible aggregation. Whether oppositely charged colloids crystallize, and which structures they form, depends on several factors, including colloid concentration, charge, and size, as well as the salt concentration of the solution. Molecular dynamics (MD) simulations are a powerful tool for predicting the outcomes of PACS assembly experiments and also provide particle-level insight into the assembly processes. Here, we present an open-source simulation framework, PACSim, that enables MD simulation studies of assembly by PACS across a range of experimentally relevant scenarios. PACSim is built on top of OpenMM, a flexible MD simulation framework that readily supports the implementation of different interaction potentials, as well as integration with other tools such as enhanced-sampling and machine-learning frameworks. We describe the motivation for PACSim, outline its features, report methodological advancements inspired by this framework, and provide examples of its use.

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 / 2 minor

Summary. The manuscript introduces PACSim, an open-source MD simulation framework built on OpenMM for studying Polymer-Attenuated Coulombic Self-Assembly (PACS) of charged colloids coated with polymer brushes. It describes the motivation from experimental PACS, outlines framework features including support for different interaction potentials, reports methodological advancements, and provides usage examples for simulating assembly across experimentally relevant regimes of concentration, charge, size, and salt concentration.

Significance. If the implemented effective potentials prove accurate, PACSim would offer a useful extensible platform for the soft-matter community, enabling particle-resolved simulations of PACS and straightforward coupling to enhanced sampling or machine-learning methods. The OpenMM foundation is a clear strength for flexibility and reproducibility. However, the absence of any quantitative validation against experimental PACS data limits the immediate significance of the contribution.

major comments (2)
  1. [Abstract and §3] Abstract and §3 (Examples of Use): the central claim that PACSim 'enables MD simulation studies of assembly by PACS across a range of experimentally relevant scenarios' and supports predictive work is not accompanied by any direct, quantitative benchmarks (e.g., simulated vs. measured critical salt concentrations or observed lattice symmetries for a known oppositely-charged pair). Without such checks the predictive utility remains untested.
  2. [§2] §2 (Implementation): the polymer-brush and screened-Coulomb potentials are presented without sensitivity analysis or comparison to more detailed models, leaving open the possibility that many-body or conformation-dependent effects controlling real crystallization are missed.
minor comments (2)
  1. [§2] Notation for the effective interaction parameters (e.g., brush thickness, effective charge) is introduced without a consolidated table, making it difficult to reproduce the example runs.
  2. [Figures] Figure captions should explicitly state the OpenMM integrator settings and cutoff distances used in the reported trajectories.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive comments on our manuscript introducing the PACSim framework. We address each major comment below and indicate the revisions made to strengthen the presentation of the work.

read point-by-point responses
  1. Referee: [Abstract and §3] Abstract and §3 (Examples of Use): the central claim that PACSim 'enables MD simulation studies of assembly by PACS across a range of experimentally relevant scenarios' and supports predictive work is not accompanied by any direct, quantitative benchmarks (e.g., simulated vs. measured critical salt concentrations or observed lattice symmetries for a known oppositely-charged pair). Without such checks the predictive utility remains untested.

    Authors: We agree that direct quantitative benchmarks against specific experimental datasets (such as critical salt concentrations or lattice symmetries) are not included. The manuscript presents PACSim primarily as an extensible simulation framework built on OpenMM to enable studies across experimentally relevant regimes, with the examples in §3 illustrating qualitative assembly behaviors consistent with PACS phenomenology. We have revised the abstract and §3 to clarify the scope of the claims, emphasizing the framework's role in supporting future predictive and validation work rather than claiming immediate predictive utility. A new paragraph has been added discussing strategies for quantitative comparison with experiment. revision: partial

  2. Referee: [§2] §2 (Implementation): the polymer-brush and screened-Coulomb potentials are presented without sensitivity analysis or comparison to more detailed models, leaving open the possibility that many-body or conformation-dependent effects controlling real crystallization are missed.

    Authors: The effective potentials were selected to capture the essential physics of polymer-attenuated Coulombic interactions while enabling efficient large-scale simulations. We acknowledge the value of sensitivity analysis. In the revised manuscript we have added a dedicated subsection to §2 that performs sensitivity analysis on key parameters including brush thickness, grafting density, and Debye length, showing their influence on assembly outcomes. We also discuss the approximations inherent to the effective-potential approach and note that PACSim's modular design permits future extensions to more detailed or explicit-polymer models. revision: yes

Circularity Check

0 steps flagged

No circularity: software framework paper with no derivation chain or fitted predictions

full rationale

The paper introduces PACSim, an OpenMM-based MD simulation framework for studying polymer-attenuated Coulombic self-assembly. It outlines motivation, features, methodological advancements, and usage examples but presents no mathematical derivations, parameter fittings, predictions, or uniqueness theorems. No load-bearing steps reduce to self-defined inputs, self-citations, or ansatzes. The central claim concerns enabling simulations across scenarios rather than deriving results from first principles within the work. This is a standard tool-description paper whose assumptions about potential validity are external to any internal loop.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The contribution is a software framework rather than new physical theory; it inherits standard molecular-dynamics assumptions and OpenMM capabilities without introducing new fitted parameters or postulated entities in the abstract.

axioms (1)
  • domain assumption Molecular dynamics with suitably chosen interaction potentials can reproduce the self-assembly behavior of polymer-brush-coated charged colloids.
    Central premise required for the framework to be useful for predicting experimental outcomes.

pith-pipeline@v0.9.0 · 5757 in / 1219 out tokens · 39484 ms · 2026-05-20T22:08:03.104415+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

  • IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    The pairwise PACS interaction can be modeled by combining screened electrostatics with polymer-brush repulsion... VP(h) ... derived within the Alexander–de Gennes model... VE(h) ... analogous to the screened Coulomb potential appearing in DLVO theory

  • IndisputableMonolith/Foundation/BranchSelection.lean branch_selection unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    PACSim is built on top of OpenMM... custom force framework... analytical expressions... polymer-brush repulsion VP and electrostatic interaction VE are always present

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

55 extracted references · 55 canonical work pages

  1. [1]

    Gasser, E

    U. Gasser, E. R. Weeks, A. Schofield, P. N. Pusey, and D. A. Weitz, Real-space imaging of nucleation and growth in colloidal crystallization, Science 292, 258 (2001)

  2. [2]

    Hueckel, G

    T. Hueckel, G. M. Hocky, and S. Sacanna, Total synthesis of colloidal matter, Nat. Rev. Mater. 6, 1053 (2021)

  3. [3]

    S. Zang, A. W. Hauser, S. Paul, G. M. Hocky, and S. Sacanna, Enabling three-dimensional real-space analy- sis of ionic colloidal crystallization, Nat. Mater. 23, 1131 (2024)

  4. [4]

    Ducrot, M

    ´E. Ducrot, M. He, G.-R. Yi, and D. J. Pine, Colloidal al- loys with preassembled clusters and spheres, Nat. Mater. 16, 652 (2017)

  5. [5]

    Y. Wang, I. C. Jenkins, J. T. McGinley, T. Sinno, and J. C. Crocker, Colloidal crystals with diamond symmetry at optical lengthscales, Nat. Commun. 8, 14173 (2017)

  6. [6]

    Zhang, A. S. Keys, T. Chen, and S. C. Glotzer, Self- assembly of patchy particles into diamond structures through molecular mimicry, Langmuir 21, 11547 (2005)

  7. [7]

    P. F. Damasceno, M. Engel, and S. C. Glotzer, Predic- tive self-assembly of polyhedra into complex structures, Science 337, 453 (2012)

  8. [8]

    Dhakal, K

    S. Dhakal, K. L. Kohlstedt, G. C. Schatz, C. A. Mirkin, and M. Olvera de la Cruz, Growth dynamics for DNA- guided nanoparticle crystallization, ACS Nano 7, 10948 (2013)

  9. [9]

    H. Fang, M. F. Hagan, and W. B. Rogers, Two-step crystallization and solid–solid transitions in binary col- loidal mixtures, Proc. Natl. Acad. Sci. U.S.A. 117, 27927 (2020)

  10. [10]

    Hueckel, G

    T. Hueckel, G. M. Hocky, J. Palacci, and S. Sacanna, Ionic solids from common colloids, Nature 580, 487 (2020)

  11. [11]

    Dijkstra and E

    M. Dijkstra and E. Luijten, From predictive modelling to machine learning and reverse engineering of colloidal self-assembly, Nat. Mater. 20, 762 (2021)

  12. [12]

    S. Zang, S. Paul, C. W. Leung, M. S. Chen, T. Hueckel, G. M. Hocky, and S. Sacanna, Direct observation and control of non-classical crystallization pathways in binary colloidal systems, Nat. Commun. 16, 3645 (2025)

  13. [13]

    R. J. Hunter, Foundations of Colloid Science , 2nd ed. (Oxford University Press, 2001)

  14. [14]

    J. N. Israelachvili, Intermolecular and Surface Forces, 3rd ed. (Academic Press, 2011)

  15. [15]

    C. A. Mirkin, R. L. Letsinger, R. C. Mucic, and J. J. Storhoff, A DNA-based method for rationally assembling nanoparticles into macroscopic materials, Nature 382, 607 (1996)

  16. [16]

    A. P. Alivisatos, K. P. Johnsson, X. Peng, T. E. Wilson, C. J. Loweth, M. P. Bruchez, and P. G. Schultz, Orga- nization of ’nanocrystal molecules’ using DNA, Nature 382, 609 (1996)

  17. [17]

    W. B. Rogers, W. M. Shih, and V. N. Manoharan, Using DNA to program the self-assembly of colloidal nanoparti- cles and microparticles, Nat. Rev. Mater.1, 16008 (2016)

  18. [18]

    M. E. Leunissen, C. G. Christova, A.-P. Hynninen, C. P. Royall, A. I. Campbell, A. Imhof, M. Dijkstra, R. van Roij, and A. van Blaaderen, Ionic colloidal crystals of oppositely charged particles, Nature 437, 235 (2005)

  19. [19]

    van Kesteren, N

    S. van Kesteren, N. Smina, S. Zang, C. W. Leung, G. M. Hocky, and S. Sacanna, Light-controlled colloidal crys- tallization, Chem 12, 102917 (2026)

  20. [20]

    J. A. Anderson, J. Glaser, and S. C. Glotzer, HOOMD- blue: A Python package for high-performance molecu- lar dynamics and hard particle Monte Carlo simulations, Comput. Mater. Sci. 173, 109363 (2020)

  21. [21]

    A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolin- tineanu, W. M. Brown, P. S. Crozier, P. J. in ’t Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton, LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales, Comput. Phys. Commun....

  22. [22]

    Eastman, R

    P. Eastman, R. Galvelis, R. P. Pel´ aez, C. R. A. Abreu, S. E. Farr, E. Gallicchio, A. Gorenko, M. M. Henry, F. Hu, J. Huang, A. Kr¨ amer, J. Michel, J. A. Mitchell, V. S. Pande, J. P. Rodrigues, J. Rodriguez-Guerra, A. C. Simmonett, S. Singh, J. Swails, P. Turner, Y. Wang, I. Zhang, J. D. Chodera, G. De Fabritiis, and T. E. Mark- land, OpenMM 8: Molecula...

  23. [23]

    J. P. Gales, M. J. Kim, G. M. Hocky, and D. J. Pine, Crystallization of non-convex colloids: the roles of parti- cle shape and entropy, Soft Matter 21, 7021 (2025)

  24. [24]

    Mitra, C

    G. Mitra, C. Chang, A. McMullen, D. Puchall, J. Brujic, and G. M. Hocky, A coarse-grained simulation model for colloidal self-assembly via explicit mobile binders, Soft 15 Matter 19, 4223 (2023)

  25. [25]

    T. Shu, G. Mitra, J. Alberts, M. P. Viana, E. D. Levy, G. M. Hocky, and L. J. Holt, Mesoscale molecular assem- bly is favored by the active, crowded cytoplasm, PRX Life 2, 033001 (2024)

  26. [26]

    T. R. Holoman, B. Prajwal, G. M. Hocky, and T. M. Truskett, Simulating dynamic bonding in soft materials, Curr. Opin. Colloid Interface Sci. 83, 102019 (2026)

  27. [27]

    Bonomi, G

    M. Bonomi, G. Bussi, C. Camilloni, G. A. Tribello, P. Ba- nas, A. Barducci, M. Bernetti, P. G. Bolhuis, S. Bot- taro, D. Branduardi, R. Capelli, P. Carloni, M. Ceriotti, A. Cesari, H. Chen, W. Chen, F. Colizzi, S. De, M. D. L. Pierre, D. Donadio, V. Drobot, B. Ensing, A. L. Fer- guson, M. Filizola, J. S. Fraser, H. Fu, P. Gasparotto, F. L. Gervasio, F. Gi...

  28. [28]

    G. A. Tribello, M. Bonomi, G. Bussi, C. Camilloni, B. I. Armstrong, A. Arsiccio, S. Aureli, F. Ballabio, M. Bernetti, L. Bonati, S. G. Brookes, Z. F. Brotza- kis, R. Capelli, M. Ceriotti, K.-T. Chan, P. Cossio, S. Dasetty, D. Donadio, B. Ensing, A. L. Ferguson, G. Fraux, J. D. Gale, F. L. Gervasio, T. Giorgino, N. S. Herringer, G. M. Hocky, S. E. Hoff, M....

  29. [29]

    Paszke, S

    A. Paszke, S. Gross, F. Massa, A. Lerer, J. Bradbury, G. Chanan, T. Killeen, Z. Lin, N. Gimelshein, L. Antiga, A. Desmaison, A. K¨ opf, E. Yang, Z. DeVito, M. Raison, A. Tejani, S. Chilamkurthy, B. Steiner, L. Fang, J. Bai, and S. Chintala, PyTorch: an imperative style, high- performance deep learning library, NeurIPS 33, 8026 (2019)

  30. [30]

    P. F. Zubieta Rico, L. Schneider, G. R. P´ erez-Lemus, R. Alessandri, S. Dasetty, T. D. Nguyen, C. A. Men´ endez, Y. Wu, Y. Jin, Y. Xu, S. Varner, J. A. Parker, A. L. Fer- guson, J. K. Whitmer, and J. J. de Pablo, PySAGES: flexible, advanced sampling methods accelerated with GPUs, npj Comput. Mater. 10, 35 (2024)

  31. [31]

    Alexander, Adsorption of chain molecules with a po- lar head a scaling description, J

    S. Alexander, Adsorption of chain molecules with a po- lar head a scaling description, J. Phys. France 38, 983 (1977)

  32. [32]

    de Gennes, Polymers at an interface; a simplified view, Adv

    P. de Gennes, Polymers at an interface; a simplified view, Adv. Colloid Interface Sci. 27, 189 (1987)

  33. [33]

    Derjaguin and L

    B. Derjaguin and L. Landau, Theory of the stability of strongly charged lyophobic sols and of the adhesion of strongly charged particles in solutions of electrolytes, Prog. Surf. Sci. 43, 30 (1993), reprint of the 1941 Acta Physicochim. U.R.S.S. article

  34. [34]

    E. J. W. Verwey and J. T. G. Overbeek, Theory of the Stability of Lyophobic Colloids (Elsevier, 1948)

  35. [35]

    Asakura and F

    S. Asakura and F. Oosawa, Interaction between particles suspended in solutions of macromolecules, J. Polym. Sci. 33, 183 (1958)

  36. [36]

    Vrij, Polymers at interfaces and the interactions in colloidal dispersions, Pure Appl

    A. Vrij, Polymers at interfaces and the interactions in colloidal dispersions, Pure Appl. Chem. 48, 471 (1976)

  37. [37]

    Dijkstra, J

    M. Dijkstra, J. M. Brader, and R. Evans, Phase be- haviour and structure of model colloid-polymer mixtures, J. Phys.: Condens. Matter 11, 10079 (1999)

  38. [38]

    Dijkstra, Computer simulations of charge and steric stabilised colloidal suspensions, Curr

    M. Dijkstra, Computer simulations of charge and steric stabilised colloidal suspensions, Curr. Opin. Colloid In- terface Sci. 6, 372 (2001)

  39. [39]

    Sacanna, W

    S. Sacanna, W. T. M. Irvine, P. M. Chaikin, and D. J. Pine, Lock and key colloids, Nature 464, 575 (2010)

  40. [40]

    C. N. Likos, K. A. Vaynberg, H. L¨ owen, and N. J. Wag- ner, Colloidal stabilization by adsorbed gelatin, Lang- muir 16, 4100 (2000)

  41. [41]

    S. J. O’Shea, M. E. Welland, and T. Rayment, An atomic force microscope study of grafted polymers on mica, Langmuir 9, 1826 (1993)

  42. [42]

    R. Hogg, T. W. Healy, and D. W. Fuerstenau, Mutual coagulation of colloidal dispersions, Trans. Faraday Soc. 62, 1638 (1966)

  43. [43]

    H. N. W. Lekkerkerker and R. Tuinier, Colloids and the Depletion Interaction (Springer, 2011)

  44. [44]

    A. V. Petukhov, R. Tuinier, and G. J. Vroege, Entropic patchiness: Effects of colloid shape and depletion, Curr. Opin. Colloid Interface Sci. 30, 54 (2017)

  45. [45]

    Stukowski, Visualization and analysis of atomistic simulation data with OVITO—the Open Visualization Tool, Model

    A. Stukowski, Visualization and analysis of atomistic simulation data with OVITO—the Open Visualization Tool, Model. Simul. Mater. Sci. Eng. 18, 015012 (2009)

  46. [46]

    Humphrey, A

    W. Humphrey, A. Dalke, and K. Schulten, VMD: Visual molecular dynamics, J. Mol. Graph. 14, 33 (1996)

  47. [47]

    van Kesteren, S

    S. van Kesteren, S. Zang, G. M. Hocky, and S. Sacanna, Structure selection in ionic colloidal crystals via indepen- dent charge tuning, In review (2026)

  48. [48]

    Hjorth Larsen, J

    A. Hjorth Larsen, J. Jørgen Mortensen, J. Blomqvist, I. E. Castelli, R. Christensen, M. Du lak, J. Friis, M. N. Groves, B. Hammer, C. Hargus, E. D. Hermes, P. C. Jennings, P. Bjerre Jensen, J. Kermode, J. R. Kitchin, E. Leonhard Kolsbjerg, J. Kubal, K. Kaasb- jerg, S. Lysgaard, J. Bergmann Maronsson, T. Max- son, T. Olsen, L. Pastewka, A. Peterson, C. Ros...

  49. [49]

    A. Jain, S. P. Ong, G. Hautier, W. Chen, W. D. Richards, S. Dacek, S. Cholia, D. Gunter, D. Skinner, G. Ceder, and K. A. Persson, Commentary: The Materials Project: A materials genome approach to accelerating materials innovation, APL Mater. 1, 011002 (2013)

  50. [50]

    Zhang, X

    Z. Zhang, X. Liu, K. Yan, M. E. Tuckerman, and J. Liu, Unified efficient thermostat scheme for the canonical ensemble with holonomic or isokinetic constraints via molecular dynamics, J. Phys. Chem. A 123, 6056 (2019)

  51. [51]

    Eastman and V

    P. Eastman and V. S. Pande, Constant constraint ma- 16 trix approximation: A robust, parallelizable constraint method for molecular simulations, J. Chem. Theory Comput. 6, 434 (2010)

  52. [52]

    Eastman and V

    P. Eastman and V. Pande, Accelerating development and execution speed with just-in-time GPU code generation, in GPU Computing Gems Jade Edition , edited by W. W. Hwu (Morgan Kaufmann, Boston, 2012) pp. 399–407

  53. [53]

    B. Peng, F. Smallenburg, A. Imhof, M. Dijkstra, and A. van Blaaderen, Colloidal clusters by using emulsions and dumbbell-shaped particles: Experiments and simu- lations, Angew. Chem. Int. Ed. 52, 6709 (2013)

  54. [54]

    M. He, J. P. Gales, ´E. Ducrot, Z. Gong, G.-R. Yi, S. Sacanna, and D. J. Pine, Colloidal diamond, Nature 585, 524 (2020)

  55. [55]

    Bussi and A

    G. Bussi and A. Laio, Using metadynamics to explore complex free-energy landscapes, Nat. Rev. Phys. 2, 200 (2020)