Recognition: unknown
PACSim: A Flexible Simulation Framework for Polymer-Attenuated Coulombic Self-Assembly
Pith reviewed 2026-05-14 18:50 UTC · model grok-4.3
The pith
PACSim is an open-source molecular dynamics framework for simulating polymer-attenuated Coulombic self-assembly of colloids.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors present PACSim as a flexible open-source MD simulation framework built on OpenMM that enables studies of assembly by PACS across a range of experimentally relevant scenarios through implementation of interaction potentials capturing polymer-brush attenuation and Coulombic forces, with support for integration with enhanced-sampling and machine-learning tools.
What carries the argument
PACSim, the simulation framework that implements customizable interaction potentials for polymer-brush attenuation and Coulombic forces on top of OpenMM.
If this is right
- Simulations can predict whether oppositely charged colloids crystallize and which structures form for given concentrations, charges, sizes, and salt levels.
- Particle-level trajectories provide insight into assembly mechanisms that complement experimental data.
- The framework supports modeling across multiple experimentally relevant parameter ranges without requiring new code for each case.
- Integration with enhanced sampling and machine learning methods becomes straightforward for studying rare assembly events.
Where Pith is reading between the lines
- Direct side-by-side comparison of PACSim trajectories against confocal microscopy movies from the same parameters could confirm or refute the model's fidelity.
- The open-source code could be extended by users to test hypothetical particle sizes or brush lengths before synthesis.
- Screening wide ranges of salt concentration in simulation first could guide efficient experimental trials to locate new crystal phases.
- Coupling PACSim output to structure-prediction algorithms might accelerate identification of stable assemblies from simple building blocks.
Load-bearing premise
The interaction potentials implemented in PACSim accurately capture the polymer-brush attenuation and Coulombic forces present in real PACS experiments without post-hoc adjustments.
What would settle it
Running PACSim with parameters matching a published PACS experiment and obtaining crystal structures or assembly kinetics that differ from the experimental observations.
Figures
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.
Referee Report
Summary. The manuscript presents PACSim, an open-source molecular dynamics simulation framework built on OpenMM for investigating Polymer-Attenuated Coulombic Self-Assembly (PACS) of charged colloidal particles coated with polymer brushes. It outlines the motivation for the tool, describes its features including flexible implementation of polymer-brush repulsion and screened Coulomb potentials, reports methodological advancements, and provides examples of simulations across varying colloid concentrations, charges, sizes, and salt conditions.
Significance. If validated, PACSim would offer a reusable, extensible platform for MD studies of PACS that integrates with enhanced sampling and machine-learning tools, potentially aiding prediction of assembly outcomes and particle-level insights in colloidal crystallization. The open-source nature and OpenMM foundation are strengths for reproducibility and extensibility in soft-matter simulations.
major comments (2)
- [Abstract] Abstract: The central claim that PACSim 'enables MD simulation studies of assembly by PACS across a range of experimentally relevant scenarios' is load-bearing for the paper's contribution but is not supported by any direct validation; no comparisons of simulated lattice constants, phase behavior, or critical salt concentrations to published PACS experiments are reported.
- [Examples of its use] Examples section: The provided simulation examples demonstrate code usage and trajectory generation but supply no quantitative benchmarks (e.g., radial distribution functions, assembly yields, or error bars) against experimental data, leaving the fidelity of the implemented polymer-brush attenuation and Coulomb potentials untested.
minor comments (2)
- [Framework features] The description of custom force implementations in OpenMM could include explicit code snippets or pseudocode for the polymer-brush repulsion term to improve reproducibility for users.
- [Methodological advancements] No discussion of computational performance (e.g., scaling with particle number or system size) is provided, which would help assess practicality for experimentally relevant system sizes.
Simulated Author's Rebuttal
We thank the referee for their constructive review of our manuscript on PACSim. We address the major comments point by point below, clarifying the scope of the work as a framework presentation while outlining planned revisions.
read point-by-point responses
-
Referee: [Abstract] Abstract: The central claim that PACSim 'enables MD simulation studies of assembly by PACS across a range of experimentally relevant scenarios' is load-bearing for the paper's contribution but is not supported by any direct validation; no comparisons of simulated lattice constants, phase behavior, or critical salt concentrations to published PACS experiments are reported.
Authors: We agree that the manuscript does not provide direct experimental validation such as comparisons of lattice constants, phase behavior, or critical salt concentrations. The primary aim is to introduce the open-source PACSim framework built on OpenMM, including its flexible implementation of polymer-brush and screened-Coulomb potentials, along with examples demonstrating its use across relevant parameter ranges (colloid concentration, charge, size, and salt). The abstract claim refers to this enabling capability rather than completed validation studies. To address the concern, we will revise the abstract to explicitly note that PACSim supplies the simulation infrastructure for such investigations, with the reported examples serving as demonstrations of functionality; comprehensive experimental benchmarking is planned as follow-on work. revision: partial
-
Referee: [Examples of its use] Examples section: The provided simulation examples demonstrate code usage and trajectory generation but supply no quantitative benchmarks (e.g., radial distribution functions, assembly yields, or error bars) against experimental data, leaving the fidelity of the implemented polymer-brush attenuation and Coulomb potentials untested.
Authors: The examples section is designed to illustrate practical code usage and the generation of trajectories under varied conditions. We acknowledge the absence of quantitative benchmarks against experiments. In the revised manuscript we will add radial distribution functions computed from the example trajectories, report assembly metrics where applicable, and include error estimates from replicate runs. We will also discuss how the implemented potentials reproduce expected theoretical behavior for polymer-brush repulsion and Debye-screened Coulomb interactions, thereby providing an internal consistency check even if full experimental matching remains beyond the current scope. revision: yes
Circularity Check
No circularity: PACSim is a tool-release paper with independent implementation description
full rationale
The manuscript presents an open-source MD framework (PACSim) built on OpenMM for simulating polymer-attenuated Coulombic self-assembly. Its central claim is the availability of the code and its features for exploring experimentally relevant scenarios. No derivation chain, fitted parameters renamed as predictions, or load-bearing self-citations appear in the abstract or described structure. The interaction potentials are implemented directly from standard forms (polymer brush repulsion plus screened Coulomb), with no reduction of outputs to inputs by construction. The reader's assessment of score 2 is consistent with minor self-citation norms but does not indicate circularity here.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption OpenMM supports flexible implementation of interaction potentials needed for PACS
Reference graph
Works this paper leans on
- [1]
-
[2]
T. Hueckel, G. M. Hocky, and S. Sacanna, Total synthesis of colloidal matter, Nat. Rev. Mater. 6, 1053 (2021)
work page 2021
-
[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)
work page 2024
- [4]
-
[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)
work page 2017
-
[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)
work page 2005
-
[7]
P. F. Damasceno, M. Engel, and S. C. Glotzer, Predic- tive self-assembly of polyhedra into complex structures, Science 337, 453 (2012)
work page 2012
- [8]
-
[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)
work page 2020
-
[10]
T. Hueckel, G. M. Hocky, J. Palacci, and S. Sacanna, Ionic solids from common colloids, Nature 580, 487 (2020)
work page 2020
-
[11]
M. Dijkstra and E. Luijten, From predictive modelling to machine learning and reverse engineering of colloidal self-assembly, Nat. Mater. 20, 762 (2021)
work page 2021
-
[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)
work page 2025
-
[13]
R. J. Hunter, Foundations of Colloid Science , 2nd ed. (Oxford University Press, 2001)
work page 2001
-
[14]
J. N. Israelachvili, Intermolecular and Surface Forces, 3rd ed. (Academic Press, 2011)
work page 2011
-
[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)
work page 1996
-
[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)
work page 1996
-
[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)
work page 2016
-
[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)
work page 2005
-
[19]
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)
work page 2026
-
[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)
work page 2020
-
[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....
work page 2022
-
[22]
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...
work page 2024
-
[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)
work page 2025
- [24]
-
[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)
work page 2024
-
[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)
work page 2026
-
[27]
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...
work page 2019
-
[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....
work page 2025
-
[29]
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)
work page 2019
-
[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)
work page 2024
-
[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)
work page 1977
-
[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)
work page 1987
-
[33]
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
work page 1993
-
[34]
E. J. W. Verwey and J. T. G. Overbeek, Theory of the Stability of Lyophobic Colloids (Elsevier, 1948)
work page 1948
-
[35]
S. Asakura and F. Oosawa, Interaction between particles suspended in solutions of macromolecules, J. Polym. Sci. 33, 183 (1958)
work page 1958
-
[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)
work page 1976
-
[37]
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)
work page 1999
-
[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)
work page 2001
-
[39]
S. Sacanna, W. T. M. Irvine, P. M. Chaikin, and D. J. Pine, Lock and key colloids, Nature 464, 575 (2010)
work page 2010
-
[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)
work page 2000
-
[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)
work page 1993
-
[42]
R. Hogg, T. W. Healy, and D. W. Fuerstenau, Mutual coagulation of colloidal dispersions, Trans. Faraday Soc. 62, 1638 (1966)
work page 1966
-
[43]
H. N. W. Lekkerkerker and R. Tuinier, Colloids and the Depletion Interaction (Springer, 2011)
work page 2011
-
[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)
work page 2017
-
[45]
A. Stukowski, Visualization and analysis of atomistic simulation data with OVITO—the Open Visualization Tool, Model. Simul. Mater. Sci. Eng. 18, 015012 (2009)
work page 2009
-
[46]
W. Humphrey, A. Dalke, and K. Schulten, VMD: Visual molecular dynamics, J. Mol. Graph. 14, 33 (1996)
work page 1996
-
[47]
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)
work page 2026
-
[48]
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...
work page 2017
-
[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)
work page 2013
- [50]
-
[51]
P. Eastman and V. S. Pande, Constant constraint ma- trix approximation: A robust, parallelizable constraint method for molecular simulations, J. Chem. Theory Comput. 6, 434 (2010)
work page 2010
-
[52]
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
work page 2012
-
[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)
work page 2013
-
[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)
work page 2020
-
[55]
G. Bussi and A. Laio, Using metadynamics to explore complex free-energy landscapes, Nat. Rev. Phys. 2, 200 (2020)
work page 2020
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.