pith. sign in

arxiv: 2606.13586 · v1 · pith:6FGTS6CDnew · submitted 2026-06-11 · ❄️ cond-mat.soft · cond-mat.stat-mech· physics.chem-ph

Exploratory digital alchemy for colloidal crystal discovery

Shih-Kuang (Alex) Lee , Sun-Ting Tsai , Sharon C. Glotzer This is my paper

Pith reviewed 2026-06-27 05:20 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.stat-mechphysics.chem-ph
keywords colloidal self-assemblydigital alchemymetadynamicsFrank-Kasper phasesforward designparticle designfree-energy landscapes
0
0 comments X

The pith

Exploratory Digital Alchemy discovers stable colloidal crystal phases by removing any pre-specified target structure.

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

The paper presents Exploratory Digital Alchemy as a forward design method that starts from the existing Digital Alchemy framework but drops the requirement to name a target crystal in advance. An exploration bias drawn from metadynamics is added to drive sampling across possible structures. Demonstrations on two-dimensional Lennard-Jones-Gauss particles map free-energy landscapes at varying temperatures, while three-dimensional oscillating pair potentials yield parameter sets that stabilize a range of Frank-Kasper phases. A reader would care because conventional inverse design can only optimize toward a guessed structure and therefore cannot guarantee that the result is the most stable or only stable outcome.

Core claim

Exploratory Digital Alchemy releases the target-crystal constraint of Digital Alchemy and augments the scheme with a metadynamics-style bias, thereby enabling parameter exploration that identifies stable and metastable colloidal crystals without prior specification of the assembled structure.

What carries the argument

Exploratory Digital Alchemy (EDA), which fuses the digital-alchemy ensemble with a metadynamics bias to sample configuration space across particle-design parameters.

If this is right

  • A wide range of oscillating pair potential parameters can be scanned to locate those that stabilize metastable Frank-Kasper phases.
  • Free-energy landscapes for Lennard-Jones-Gauss particles can be constructed at multiple temperatures without fixing the target lattice.
  • The combined framework can be applied to study alchemical reactions inside a generalized statistical ensemble.
  • Particle designs can be optimized for self-assembly outcomes whose structures are unknown at the start of the search.

Where Pith is reading between the lines

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

  • The same bias-augmented scheme could be tested on other pair potentials or on systems with three-body interactions to check whether additional crystal families become discoverable.
  • If the method scales to larger particle numbers, it may reduce reliance on exhaustive enumeration when screening colloidal building blocks for new materials.
  • Integration with other enhanced-sampling techniques might further accelerate identification of low-free-energy assemblies.

Load-bearing premise

The metadynamics-style bias produces sampling of configuration space that is sufficiently unbiased to reveal truly stable structures rather than artifacts of the added potential.

What would settle it

A direct free-energy comparison in which a crystal phase found by EDA is shown to have higher free energy than a different phase that the method never reported under the same potential parameters.

read the original abstract

Digital Alchemy (DA), introduced by Van Anders et al., is a statistical mechanics-based generalized thermodynamic ensemble method that employs computer simulations to optimize colloidal particle design. This approach applies the principles of statistical mechanics to predict and tailor particle attributes that lead to desired self-assembled structures or material properties. However, as an inverse design method, its main limitation is that the target structure must be known \textit{a priori}. Therefore, the optimal design from DA does not guarantee the targeted structure is the most or the only stable one. This highlights the importance of forward design with an exploratory scheme for optimizing novel colloid designs, which becomes more suitable in such cases. In this paper, we introduce Exploratory Digital Alchemy (EDA), an enhanced forward design scheme that begins by releasing the constraint of the target crystal from DA, followed by an exploration-oriented bias that has been extensively used in enhanced sampling methods such as metadynamics (MetaD). We demonstrate the utility of EDA through examples involving particles interacting via a two-dimensional Lennard-Jones Gauss potential (LJGP) and a three-dimensional oscillating pair potential (OPP). We applied EDA to study the free energy landscapes given different potential parameters of LJGP at different temperatures. With the exploratory scheme, we've also successfully identified a wide range of OPP potential parameters that stabilize metastable Frank-Kasper phases. Our approach fuses the standard DA framework with metadynamics, which could potentially be useful for studying alchemical reactions in a generalized ensemble.

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

3 major / 2 minor

Summary. The manuscript introduces Exploratory Digital Alchemy (EDA) as an extension of Digital Alchemy (DA) that removes the a priori target-structure constraint and augments the ensemble with a metadynamics-style bias to enable forward exploration of colloidal self-assembly. It applies the method to two-dimensional Lennard-Jones-Gauss potentials (LJGP) at varying temperatures and to three-dimensional oscillating pair potentials (OPP), claiming that EDA identifies parameter sets stabilizing a range of Frank-Kasper phases.

Significance. A validated EDA framework would provide a practical route for discovering unanticipated stable or metastable colloidal crystals without presupposing their structure, extending the utility of alchemical ensembles beyond inverse design.

major comments (3)
  1. [Abstract; Results (OPP section)] Abstract and Results (OPP examples): the claim that EDA 'successfully identified a wide range of OPP potential parameters that stabilize metastable Frank-Kasper phases' is supported only by trajectories generated under the metadynamics bias; no unbiased long-time MD runs, free-energy calculations, or stability checks after bias removal are reported to confirm that the discovered phases are minima of the unbiased landscape.
  2. [Method (EDA construction)] Method (EDA bias implementation): the description of the metadynamics-style bias added to the DA ensemble does not include any diagnostic (e.g., bias-height convergence, reweighting tests, or comparison to unbiased sampling) demonstrating that the bias does not create or stabilize spurious configurations that would disappear once the bias is removed.
  3. [Results (LJGP subsection)] Results (LJGP free-energy landscapes): the statement that EDA was 'applied successfully' to LJGP landscapes supplies no quantitative metrics, error estimates, or comparison against standard unbiased sampling or known phase boundaries, leaving the utility claim without numerical support.
minor comments (2)
  1. [Method] Notation for the combined DA+MetaD ensemble is introduced without an explicit equation defining the modified partition function or the bias potential form; adding this would clarify the method.
  2. [Figures (OPP results)] Figure captions for the OPP parameter scans do not state the simulation length, bias parameters, or number of independent runs, making reproducibility difficult.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their thorough review and valuable feedback on our manuscript introducing Exploratory Digital Alchemy (EDA). We address each of the major comments below and outline the revisions we plan to make.

read point-by-point responses

  1. Referee: [Abstract; Results (OPP section)] Abstract and Results (OPP examples): the claim that EDA 'successfully identified a wide range of OPP potential parameters that stabilize metastable Frank-Kasper phases' is supported only by trajectories generated under the metadynamics bias; no unbiased long-time MD runs, free-energy calculations, or stability checks after bias removal are reported to confirm that the discovered phases are minima of the unbiased landscape.

    Authors: We acknowledge this limitation in the current presentation. The EDA method uses the metadynamics bias to explore the alchemical and configurational space, and the identified parameters are those for which the Frank-Kasper phases are visited and stabilized under the bias. However, to strengthen the claim of metastability in the unbiased ensemble, we will perform additional unbiased molecular dynamics simulations starting from the discovered configurations for selected parameter sets and report their stability over long times. We will also include free-energy estimates where feasible. revision: yes

  2. Referee: [Method (EDA construction)] Method (EDA bias implementation): the description of the metadynamics-style bias added to the DA ensemble does not include any diagnostic (e.g., bias-height convergence, reweighting tests, or comparison to unbiased sampling) demonstrating that the bias does not create or stabilize spurious configurations that would disappear once the bias is removed.

    Authors: The referee is correct that the method section lacks these diagnostics. In the revised manuscript, we will expand the Methods section to include convergence diagnostics for the bias potential, such as the time evolution of the bias height, and perform reweighting analyses or direct comparisons with unbiased sampling for key cases to demonstrate that the observed phases persist without the bias. revision: yes

  3. Referee: [Results (LJGP subsection)] Results (LJGP free-energy landscapes): the statement that EDA was 'applied successfully' to LJGP landscapes supplies no quantitative metrics, error estimates, or comparison against standard unbiased sampling or known phase boundaries, leaving the utility claim without numerical support.

    Authors: We agree that the LJGP results would benefit from more quantitative support. We will revise this section to include quantitative metrics, such as computed free energies with error estimates from block averaging or multiple independent runs, and where possible, compare the observed phase boundaries to those from standard unbiased simulations or literature values. revision: yes

Circularity Check

0 steps flagged

No circularity: method extension and simulation results are independent of inputs

full rationale

The paper presents EDA as a direct combination of existing DA (cited to Van Anders et al.) and standard metadynamics biasing, then reports simulation outcomes on LJGP and OPP potentials. No equations, parameters, or success criteria are defined in terms of the target outputs; the identification of Frank-Kasper phases is a reported simulation result rather than a fitted or renamed input. No self-citations are load-bearing for the central claim, and the derivation chain does not reduce any prediction to its own construction. The method is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the applicability of statistical mechanics ensembles to colloidal design and on the effectiveness of metadynamics bias for exploration; no specific numerical free parameters or new physical entities are introduced in the abstract.

axioms (1)
  • domain assumption Statistical mechanics principles can be used to optimize particle attributes via computer simulations in a generalized ensemble.
    This underpins both the original DA and the proposed EDA extension.

pith-pipeline@v0.9.1-grok · 5800 in / 1316 out tokens · 29256 ms · 2026-06-27T05:20:04.486104+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

34 extracted references · 1 canonical work pages · 1 internal anchor

  1. [1]

    A., Engel, M

    Boles, M. A., Engel, M. & Talapin, D. V . Self-assembly of colloidal nanocrystals: From intricate structures to functional materials.Chemical Reviews116, 11220–11289 (2016)

    2016

  2. [2]

    & Wei, G

    Wang, L., Gong, C., Yuan, X. & Wei, G. Controlling the self-assembly of biomolecules into functional nanomaterials through internal interactions and external stimulations: A review.Nanomaterials9, 285 (2019)

    2019

  3. [3]

    A.et al.Sustainable formulation polymers for home, beauty and personal care: challenges and opportunities.Chemical Science14, 12926–12940 (2023)

    Picken, C. A.et al.Sustainable formulation polymers for home, beauty and personal care: challenges and opportunities.Chemical Science14, 12926–12940 (2023)

    2023

  4. [4]

    S., Leonforte, F., Greaves, A., Rubio, R

    Luengo, G. S., Leonforte, F., Greaves, A., Rubio, R. G. & Guzman, E. Physico-chemical challenges on the self- assembly of natural and bio-based ingredients on hair surfaces: towards sustainable haircare formulations.Green Chemistry25, 7863–7882 (2023)

    2023

  5. [5]

    J.et al.Engineering precision nanoparticles for drug delivery.Nature Reviews Drug Discovery20, 101–124 (2021)

    Mitchell, M. J.et al.Engineering precision nanoparticles for drug delivery.Nature Reviews Drug Discovery20, 101–124 (2021)

    2021

  6. [6]

    S., Dodd, P

    Van Anders, G., Klotsa, D., Karas, A. S., Dodd, P. M. & Glotzer, S. C. Digital alchemy for materials design: Colloids and beyond.Acs Nano9, 9542–9553 (2015)

    2015

  7. [7]

    M., Dshemuchadse, J

    Geng, Y ., van Anders, G., Dodd, P. M., Dshemuchadse, J. & Glotzer, S. C. Engineering entropy for the inverse design of colloidal crystals from hard shapes.Science Advances5, eaaw0514 (2019)

    2019

  8. [8]

    Y ., Moore, T

    Rivera-Rivera, L. Y ., Moore, T. C. & Glotzer, S. C. Inverse design of triblock janus spheres for self-assembly of complex structures in the crystallization slot via digital alchemy.Soft Matter19, 2726–2736 (2023)

    2023

  9. [9]

    C., van Anders, G

    Zhou, P., Proctor, J. C., van Anders, G. & Glotzer, S. C. Alchemical molecular dynamics for inverse design. Molecular Physics117, 3968–3980 (2019)

    2019

  10. [10]

    & Glotzer, S

    Zhou, P. & Glotzer, S. C. Inverse design of isotropic pair potentials using digital alchemy with a generalized fourier potential.The European Physical Journal B94, 243 (2021)

    2021

  11. [11]

    & Glotzer, S

    Geng, Y ., van Anders, G. & Glotzer, S. C. Synthesizable nanoparticle eigenshapes for colloidal crystals.Nanoscale 13, 13301–13309 (2021)

    2021

  12. [12]

    Y ., Lee, S.-K

    Zhang, M. Y ., Lee, S.-K. A., Glotzer, S. C. & Lindsey, R. K. A generalized machine-learning framework for developing alchemical many-body interaction models for polymer-grafted nanoparticles.Journal of Chemical Theory and Computation21, 9853–9867 (2025). 10

    2025

  13. [13]

    & Parrinello, M

    Laio, A. & Parrinello, M. Escaping free-energy minima.Proceedings of the national academy of sciences99, 12562–12566 (2002)

    2002

  14. [14]

    T., Bussi, G

    Hsu, W.-T., Piomponi, V ., Merz, P. T., Bussi, G. & Shirts, M. R. Alchemical metadynamics: Adding alchemical variables to metadynamics to enhance sampling in free energy calculations.Journal of Chemical Theory and Computation19, 1805–1817 (2023)

    2023

  15. [15]

    & Trebin, H.-R

    Engel, M. & Trebin, H.-R. Self-assembly of monatomic complex crystals and quasicrystals with a double- well interaction potential.Physical Review Letters98, 225505 (2007). URL https://link.aps.org/doi/10.1103/ PhysRevLett.98.225505

    2007

  16. [16]

    F., Phillips, C

    Dshemuchadse, J., Damasceno, P. F., Phillips, C. L., Engel, M. & Glotzer, S. C. Moving beyond the constraints of chemistry via crystal structure discovery with isotropic multiwell pair potentials.Proceedings of the National Academy of Sciences118, e2024034118 (2021)

    2021

  17. [17]

    & Parrinello, M

    Barducci, A., Bussi, G. & Parrinello, M. Well-tempered metadynamics: a smoothly converging and tunable free- energy method.Physical Review Letters100, 020603 (2008)

    2008

  18. [18]

    & Parrinello, M

    Valsson, O., Tiwary, P. & Parrinello, M. Enhancing important fluctuations: Rare events and metadynamics from a conceptual viewpoint.Annual Review of Physical Chemistry67, 159–184 (2016)

    2016

  19. [19]

    & Laio, A

    Bussi, G. & Laio, A. Using metadynamics to explore complex free-energy landscapes.Nature Reviews Physics2, 200–212 (2020)

    2020

  20. [20]

    C., Stillinger, F

    Rechtsman, M. C., Stillinger, F. H. & Torquato, S. Optimized interactions for targeted self-assembly: application to a honeycomb lattice.Physical Review Letters95, 228301 (2005)

    2005

  21. [21]

    & Henley, C

    Mihalkovi ˇc, M. & Henley, C. Empirical oscillating potentials for alloys from ab initio fits and the prediction of quasicrystal-related structures in the al-cu-sc system.Physical Review B—Condensed Matter and Materials Physics85, 092102 (2012)

    2012

  22. [22]

    Elenius, M.et al.Structural model for octagonal quasicrystals derived from octagonal symmetry elements arising in β-mn crystallization of a simple monatomic liquid.Physical Review B—Condensed Matter and Materials Physics 79, 144201 (2009)

    2009

  23. [23]

    W.)Handbook of computational statistics: Concepts and methods549–569 (Springer, 2011)

    inMultivariate density estimation and visualization(ed.Scott, D. W.)Handbook of computational statistics: Concepts and methods549–569 (Springer, 2011)

    2011

  24. [24]

    W.Density estimation for statistics and data analysis(Routledge, 2018)

    Silverman, B. W.Density estimation for statistics and data analysis(Routledge, 2018)

    2018

  25. [25]

    W., Rosenbluth, M

    Metropolis, N., Rosenbluth, A. W., Rosenbluth, M. N., Teller, A. H. & Teller, E. Equation of state calculations by fast computing machines.The Journal of Chemical Physics21, 1087–1092 (1953)

    1953

  26. [26]

    J., Tobias, D

    Martyna, G. J., Tobias, D. J. & Klein, M. L. Constant pressure molecular dynamics algorithms.The Journal of Chemical Physics101, 10–1063 (1994)

    1994

  27. [27]

    & Parrinello, M

    Tiwary, P. & Parrinello, M. A time-independent free energy estimator for metadynamics.The Journal of Physical Chemistry B119, 736–742 (2015)

    2015

  28. [28]

    A., Glaser, J

    Anderson, J. A., Glaser, J. & Glotzer, S. C. Hoomd-blue: A python package for high-performance molecular dynamics and hard particle monte carlo simulations.Computational Materials Science173, 109363 (2020)

    2020

  29. [29]

    Ramasubramani, V .et al.freud: A software suite for high throughput analysis of particle simulation data.Computer Physics Communications254, 107275 (2020)

    2020

  30. [30]

    Visualization and analysis of atomistic simulation data with ovito–the open visualization tool

    Stukowski, A. Visualization and analysis of atomistic simulation data with ovito–the open visualization tool. Modelling and Simulation in Materials Science and Engineering18, 015012 (2009). 11

    2009

  31. [31]

    C., Schröder-Turk, G

    Mickel, W., Kapfer, S. C., Schröder-Turk, G. E. & Mecke, K. Shortcomings of the bond orientational order parameters for the analysis of disordered particulate matter.The Journal of Chemical Physics138(2013)

    2013

  32. [32]

    & Dellago, C

    Lechner, W. & Dellago, C. Accurate determination of crystal structures based on averaged local bond order parameters.The Journal of Chemical Physics129(2008)

    2008

  33. [33]

    Fijan, D., Rashidi, M. R. W. & Glotzer, S. C. Quantifying local point-group-symmetry order in complex particle systems.arXiv preprint arXiv:2509.12665(2025)

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  34. [34]

    S., Tsai, S.-T

    Zhao, C. S., Tsai, S.-T. & Glotzer, S. C. Hybrid monte carlo metadynamics (hybridmc-metad).The Journal of Chemical Physics164(2026). 12

    2026