pith. sign in

arxiv: 1906.10510 · v1 · pith:DAC34XCOnew · submitted 2019-06-25 · ❄️ cond-mat.soft · physics.chem-ph

Energy Landscapes for Digital Alchemy

John W. R. Morgan , Sharon C. Glotzer This is my paper

Pith reviewed 2026-05-25 16:12 UTC · model grok-4.3

classification ❄️ cond-mat.soft physics.chem-ph
keywords energy landscapesdigital alchemyalchemical landscapepotential parametersgeometrical optimisationlandscape frustrationtransition statessmall clusters
0
0 comments X

The pith

Treating potential parameters as degrees of freedom creates an alchemical energy landscape on which the lowest-energy minimum is easy to locate and minima concentrate in particular regions.

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

The paper extends energy landscape methods to digital alchemy by treating the parameters of the interparticle potential as additional degrees of freedom alongside particle positions. For small clusters, geometrical optimisation is used to locate minima and transition states on this combined landscape. The results establish that parameters yielding the global lowest energy minimum are straightforward to identify, that the minima concentrate in specific areas of parameter space, and that the landscape is more frustrated than the fixed-parameter case due to high barriers separating competing low-energy structures. Transition states are classified according to whether they become minima or remain saddles once parameters are fixed, with the former showing a significant alchemical component.

Core claim

In the alchemical energy landscape, where potential parameters vary as degrees of freedom, it is straightforward to find the parameters that give the lowest energy minimum for small clusters, the distribution of minima is concentrated in particular areas, and the landscape is more frustrated in terms of competition between low-energy structures separated by high barriers. Transition states on this landscape are classified by whether they become minima or transition states when the potential parameters are fixed; those that become minima have a significant alchemical component while those that remain transition states are characterised mainly by atomic displacements.

What carries the argument

The alchemical landscape defined by treating the parameters of the potential as degrees of freedom in addition to atomic coordinates.

If this is right

  • The parameters giving the lowest energy minimum are easy to locate by geometrical optimisation on the alchemical landscape.
  • Minima on the alchemical landscape concentrate in particular areas of parameter space.
  • The alchemical landscape is more frustrated, with greater competition between low-energy structures separated by high barriers.
  • Transition states on the alchemical landscape divide into those with a significant alchemical component that become minima when parameters are fixed and those that remain transition states characterised by atomic displacements.

Where Pith is reading between the lines

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

  • The concentration of minima suggests that inverse design of interactions could be simplified by searching the alchemical landscape rather than fixing parameters in advance.
  • The increased frustration may imply that pathways between different stable structures require larger parameter changes than expected from fixed-potential landscapes.
  • Extending the classification of transition states to larger clusters could reveal whether the alchemical component remains a useful diagnostic for design problems.

Load-bearing premise

Geometrical optimisation methods locate the relevant minima and transition states throughout the combined configuration-plus-parameter space without missing important low-energy regions or misclassifying saddles.

What would settle it

An independent global search over parameter space that locates a lower-energy minimum outside the regions identified by the geometrical optimisation, or that uncovers additional low-barrier transition states connecting to unvisited minima.

Figures

Figures reproduced from arXiv: 1906.10510 by John W. R. Morgan, Sharon C. Glotzer.

Figure 1
Figure 1. Figure 1: The density of minima as a function of the alchemica [PITH_FULL_IMAGE:figures/full_fig_p017_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: The OPP potential for the three sets of fixed paramet [PITH_FULL_IMAGE:figures/full_fig_p018_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Disconnectivity graphs for the alchemical landsc [PITH_FULL_IMAGE:figures/full_fig_p019_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: The variation of the frustration indices with temp [PITH_FULL_IMAGE:figures/full_fig_p020_4.png] view at source ↗
read the original abstract

We apply energy landscape methods to digital alchemy, defining a system in which the parameters of the potential are treated as degrees of freedom. Using geometrical optimisation, we locate minima and transition states on the landscape for small clusters. We show that it is easy to find the parameters that give the lowest energy minimum, and that the distribution of minima on the alchemical landscape is concentrated in particular areas. We also conclude that the alchemical landscape is more frustrated, in terms of competition between low energy structures separated by high barriers. Transition states on the alchemical landscape are classified by whether they become minima or transition states when the potential parameters are fixed. Those that become minima have a significant alchemical component, while those that remain as transition states can be characterised mainly in terms of atomic displacements.

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 paper applies energy landscape methods to 'digital alchemy' by augmenting the configuration space of small atomic clusters with continuous potential parameters as additional degrees of freedom, thereby defining an alchemical landscape. Geometrical optimization (basin-hopping style) is used to locate minima and transition states in this joint space. The central claims are that parameters yielding the lowest-energy minimum are easy to locate, that the distribution of minima is concentrated in particular regions of parameter space, and that the alchemical landscape is more frustrated (low-energy structures separated by high barriers) than the conventional configurational landscape. Transition states are further classified according to whether they remain transition states or become minima when parameters are fixed, with the former dominated by atomic displacements and the latter having a significant alchemical component.

Significance. If the optimization results are complete, the work supplies a concrete demonstration that parameter space can be explored on equal footing with configuration space for cluster systems, together with a useful classification of saddles that distinguishes alchemical from configurational character. This could inform inverse-design strategies in soft-matter systems where interaction parameters are tunable. The approach itself is a straightforward extension of existing landscape methods, so its significance hinges on whether the reported concentration and frustration metrics survive more exhaustive sampling.

major comments (2)
  1. [Methods] Methods (optimization procedure): No benchmarks are reported—multiple random restarts, exhaustive enumeration for the smallest cluster sizes (N≤5), or recovery of known global minima for standard potentials—to establish that the geometrical optimization has not missed lower-energy basins or misclassified saddles in the augmented configuration-plus-parameter space. This directly undermines the claims of easy location of lowest-energy parameters, concentration of minima, and increased frustration, all of which presuppose exhaustive sampling.
  2. [Results] Results (frustration and concentration claims): The assertion that the alchemical landscape is 'more frustrated' is presented without a quantitative metric (e.g., barrier-height distributions or disconnectivity-graph statistics) or a side-by-side comparison against the fixed-parameter configurational landscape for the same clusters and potentials, rendering the comparative conclusion difficult to assess.
minor comments (2)
  1. [Abstract] Abstract: Specific numerical details—cluster sizes examined, functional form of the potential, number of independent optimization runs, and convergence criteria—are omitted, preventing the reader from gauging the scope and robustness of the reported findings.
  2. [Figures and text] Figure captions and text: Notation for the alchemical coordinates and the precise definition of 'frustration' should be clarified to avoid ambiguity when comparing landscapes.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and their constructive comments. We address the major comments below and will incorporate revisions to improve the presentation of our results.

read point-by-point responses

  1. Referee: [Methods] Methods (optimization procedure): No benchmarks are reported—multiple random restarts, exhaustive enumeration for the smallest cluster sizes (N≤5), or recovery of known global minima for standard potentials—to establish that the geometrical optimization has not missed lower-energy basins or misclassified saddles in the augmented configuration-plus-parameter space. This directly undermines the claims of easy location of lowest-energy parameters, concentration of minima, and increased frustration, all of which presuppose exhaustive sampling.

    Authors: We agree that benchmarks are important to validate the optimization procedure. In the revised manuscript, we will add results from multiple random restarts for the cluster sizes considered, as well as exhaustive enumeration for N≤5 where feasible, and demonstrate recovery of known global minima for standard potentials. These additions will confirm that lower-energy basins have not been missed and support the reported findings on the location of lowest-energy minima and the distribution of minima. revision: yes

  2. Referee: [Results] Results (frustration and concentration claims): The assertion that the alchemical landscape is 'more frustrated' is presented without a quantitative metric (e.g., barrier-height distributions or disconnectivity-graph statistics) or a side-by-side comparison against the fixed-parameter configurational landscape for the same clusters and potentials, rendering the comparative conclusion difficult to assess.

    Authors: We acknowledge that the manuscript lacks a quantitative comparison for the frustration claim. We will revise the results section to include quantitative metrics such as barrier height distributions and disconnectivity graph statistics for both the alchemical landscape and the corresponding fixed-parameter configurational landscapes. This will provide a direct side-by-side comparison and allow readers to assess the increased frustration in the alchemical case. revision: yes

Circularity Check

0 steps flagged

No significant circularity; claims are direct outputs of geometrical optimization

full rationale

The paper defines an alchemical landscape by treating potential parameters as continuous degrees of freedom and applies standard geometrical optimization (basin-hopping style) to locate minima and transition states. Reported results—ease of locating lowest-energy parameters, concentration of minima, and increased frustration—are direct numerical outputs of this search procedure. No equations reduce by construction to fitted inputs, no predictions are statistically forced from the same data, and no load-bearing self-citations or uniqueness theorems are invoked. The derivation chain is self-contained against external benchmarks and receives a low circularity score.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

The central claims rest on the assumption that standard geometrical optimisation reliably samples the joint configuration-parameter space; no free parameters are introduced in the abstract, and the only invented framing is the alchemical landscape itself.

axioms (1)
  • domain assumption Geometrical optimisation locates the relevant minima and transition states in the extended alchemical space.
    Invoked to justify the reported locations of minima and the classification of transition states.
invented entities (1)
  • alchemical landscape no independent evidence
    purpose: To treat potential parameters as explicit degrees of freedom alongside particle coordinates.
    New framing introduced to combine parameter search with energy landscape exploration; no independent evidence outside the computational results is provided.

pith-pipeline@v0.9.0 · 5654 in / 1323 out tokens · 39707 ms · 2026-05-25T16:12:19.347590+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

45 extracted references · 45 canonical work pages

  1. [1]

    Zhang, Z.; Glotzer, S. C. Self-assembly of patchy particles. 2004, 4, 1407--1413

    work page 2004

  2. [2]

    C.; Solomon, M

    Glotzer, S. C.; Solomon, M. J. Anisotropy of building blocks and their assembly into complex structures. 2007, 6, 557--562

    work page 2007

  3. [3]

    B.; Kretzschmar, I

    Pawar, A. B.; Kretzschmar, I. Fabrication, assembly, and application of patchy particles. 2007, 31, 150--168

    work page 2007

  4. [4]

    Sacanna, S.; Pine, D. J. Shape-anisotropic colloids: Building blocks for complex assemblies. 2011, 16, 96--105

    work page 2011

  5. [5]

    Wales, D. J. Energy landscapes; Cambridge University Press, 2003

    work page 2003

  6. [6]

    Morgan, J. W. R.; Chakrabarti, D.; Dorsaz, N.; Wales, D. J. Designing a B ernal spiral from patchy colloids. 2013, 7, 1246--1256

    work page 2013

  7. [7]

    Supracolloidal reconfigurable polyhedra via hierarchical self-assembly

    Morphew, D.; Chakrabarti, D. Supracolloidal reconfigurable polyhedra via hierarchical self-assembly. 2016, 12, 9633--9640

    work page 2016

  8. [8]

    N.; Wales, D

    Fejer, S. N.; Wales, D. J. Design of a kagome lattice from soft anisotropic particles. 2015, 11, 6663--6668

    work page 2015

  9. [9]

    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. 2015, 9, 9542--9553

    work page 2015

  10. [10]

    X.; van Anders, G.; Newman, R

    Du, C. X.; van Anders, G.; Newman, R. S.; Glotzer, S. C. Shape-driven solid-solid transitions in colloids. 2017, 114, E3892--E3899

    work page 2017

  11. [11]

    K.; van Anders, G.; Dodd, P

    Cersonsky, R. K.; van Anders, G.; Dodd, P. M.; Glotzer, S. C. Relevance of packing to colloidal self-assembly. 2018,

    work page 2018

  12. [12]

    Mihalkovi c , M.; Henley, C. L. Empirical oscillating potentials for alloys from ab initio fits and the prediction of quasicrystal-related structures in the A l- C u- S c system. 2012, 85, 092102

    work page 2012

  13. [13]

    F.; Phillips, C

    Engel, M.; Damasceno, P. F.; Phillips, C. L.; Glotzer, S. C. Computational self-assembly of a one-component icosahedral quasicrystal. 2015, 14, 109--116

    work page 2015

  14. [14]

    D.; Ford, J

    Stoddard, S. D.; Ford, J. Numerical experiments on the stochastic behavior of a L ennard- J ones gas system. 1973, 8, 1504--1512

    work page 1973

  15. [15]

    A.; Lorenze, C

    Anderson, J. A.; Lorenze, C. D.; Travesset, A. General purpose molecular dynamics simulations fully implemented on graphics processing units. 2008, 227, 5342--5359

    work page 2008

  16. [16]

    D.; Anderson, J

    Glaser, J.; Nguyen, T. D.; Anderson, J. A.; Liu, P.; Spiga, F.; Millan, J. A.; Morse, D. C.; Glotzer, S. C. Strong scaling of general-purpose molecular dynamics simulations on GPUs. 2015, 192, 97--107

    work page 2015

  17. [17]

    Li, Z.; Scheraga, H. A. Monte C arlo-minimization approach to the multiple-minima problem in protein folding. 1987, 84, 6611--6615

    work page 1987

  18. [18]

    Li, Z.; Scheraga, H. A. Structure and free energy of complex thermodynamic systems. 1988, 179, 333

    work page 1988

  19. [19]

    J.; Doye, J

    Wales, D. J.; Doye, J. P. K. Global optimization by basin-hopping and the lowest energy structures of L ennard- J ones clusters containing up to 110 atoms. 1997, 101, 5111--5116

    work page 1997

  20. [20]

    Wales, D. J. GMIN : A program for basin-hopping global optimisation. http://www-wales.ch.cam.ac.uk/software.html

    work page

  21. [21]

    A.; Wales, D

    Trygubenko, S. A.; Wales, D. J. A doubly nudged elastic band method for finding transition states. 2004, 120, 2082--2094

    work page 2004

  22. [22]

    A.; Wales, D

    Trygubenko, S. A.; Wales, D. J. Erratum: A doubly nudged elastic band method for finding transition states [J. Chem. Phys. 120, 2082 (2004)]. 2004, 120, 7820--7820

    work page 2082

  23. [23]

    Optimization methods for finding minimum energy paths

    Sheppard, D.; Terrell, R.; Henkelman, G. Optimization methods for finding minimum energy paths. 2008, 128, 134106

    work page 2008

  24. [24]

    Quantum and thermal effects in H_2 dissociative adsorption: Evaluation of free energy barriers in multidimensional quantum systems

    Mills, G.; J\'onsson, H. Quantum and thermal effects in H_2 dissociative adsorption: Evaluation of free energy barriers in multidimensional quantum systems. 1994, 72, 1124--1127

    work page 1994

  25. [25]

    Mills, G.; J\'onsson, H.; Schenter, G. K. Reversible work transition state theory: application to dissociative adsorption of hydrogen. 1995, 324, 305--337

    work page 1995

  26. [26]

    J\'onsson, H.; Mills, G.; Jacobsen, K. W. Classical and quantum dynamics in condensed phase simulations; World Scientific: Singapore, 1998

    work page 1998

  27. [27]

    P.; J\'onsson, H

    Henkelman, G.; Uberuaga, B. P.; J\'onsson, H. A climbing image nudged elastic band method for finding saddle points and minimum energy paths. 2000, 113, 9901--9904

    work page 2000

  28. [28]

    Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points

    Henkelman, G.; J\'onsson, H. Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points. The Journal of Chemical Physics 2000, 113, 9978--9985

    work page 2000

  29. [29]

    J.; Wales, D

    Munro, L. J.; Wales, D. J. Defect migration in crystalline silicon. 1999, 59, 3969--3980

    work page 1999

  30. [30]

    Wales, D. J. OPTIM : A Program for optimizing geometries and calculating reaction pathways. http://www-wales.ch.cam.ac.uk/software.html

    work page

  31. [31]

    Wales, D. J. PATHSAMPLE : A program for generating connected stationary point databases and extracting global kinetics. http://www-wales.ch.cam.ac.uk/software.html

    work page

  32. [32]

    J.; Doye, J

    Wales, D. J.; Doye, J. P. K.; Miller, M. A.; Mortenson, P. N.; Walsh, T. R. Energy landscapes: From clusters to biomolecules. 2000, 115, 1--111

    work page 2000

  33. [33]

    J.; Munro, L

    Kumeda, Y.; Wales, D. J.; Munro, L. J. Transition states and rearrangement mechanisms from hybrid eigenvector-following and density functional theory: Application to C_ 10 H_ 10 and defect migration in crystalline silicon. 2001, 341, 185--194

    work page 2001

  34. [34]

    Wales, D. J. Exploring energy landscapes. 2018, 69, 401--425

    work page 2018

  35. [35]

    D.; Onuchic, J

    Bryangelson, J. D.; Onuchic, J. N.; Socci, N. D.; Wolynes, P. G. Funnels, pathways, and the energy landscape of protein folding: A synthesis. 1995, 21, 167--195

    work page 1995

  36. [36]

    N.; Luthey-Schulten, Z.; Wolynes, P

    Onuchic, J. N.; Luthey-Schulten, Z.; Wolynes, P. G. Theory of protein folding: The energy landscape perspective. 1997, 48, 545--600

    work page 1997

  37. [37]

    K.; Stevenson, J

    de Souza, V. K.; Stevenson, J. D.; Niblett, S. P.; Farrell, J. D.; Wales, D. J. Defining and quantifying frustration in the energy landscape: Applications to atomic and molecular clusters, biomolecules, jammed and glassy systems. 2017, 146, 124103

    work page 2017

  38. [38]

    H.; Weber, T

    Stillinger, F. H.; Weber, T. A. Packing structures and transitions in liquids and solids. 1984, 225, 983--989

    work page 1984

  39. [39]

    J.; Doye, J

    Wales, D. J.; Doye, J. P. K. Stationary points and dynamics in high-dimensional systems. 2003, 119, 12409--12416

    work page 2003

  40. [40]

    Morgan, J. W. R.; Wales, D. J. Energy landscapes of planar colloidal clusters. 2014, 6, 10717--10726

    work page 2014

  41. [41]

    M.; Karplus, M

    Becker, O. M.; Karplus, M. The topology of multidimensional potential energy surfaces: theory and application to peptide structure and kinetics. 1997, 106, 1495--1517

    work page 1997

  42. [42]

    J.; Miller, M

    Wales, D. J.; Miller, M. A.; Walsh, T. A. Archetypal energy landscapes. 1998, 394, 758--760

    work page 1998

  43. [43]

    Wales, D. J. Coexistence in small inert gas clusters. 1993, 78, 151--171

    work page 1993

  44. [44]

    Criterion for the action of applied stress in the martensitic transformation

    Patel, J.; Cohen, M. Criterion for the action of applied stress in the martensitic transformation. 1953, 1, 531--538

    work page 1953

  45. [45]

    ,kFdu^up&E@'d@Kr)2XVh]O<=Q2`SI4aMG'(0oS6AP`M#^ r,&S_s G@n0Q2Dj0,>66s+5-N.?AiI>iR=F(7c>sI>5nF/b&-P914+rN*+m+]f>/b@o[i

    Murrell, J. N.; Laidler, K. J. Symmetries of activated complexes. 1968, 64, 371--377 mcitethebibliography opp.tex0000664000000000000000000006646213504416614011120 0ustar rootroot [journal=jctcce,manuscript=article] achemso [version=3] mhchem [T1] fontenc amsmath subfigure booktabs psfrag xspace placeins acs articletitle = true table John W.\ R.\ Morgan jw...

    work page doi:10.5281/zenodo.3256234 1968