Pith. sign in

REVIEW 3 major objections 5 minor 84 references

Inverse Thermodynamics: Designing Interactions for Targeted Phase Behavior

T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

Pith's one-line read The paper establishes design rules that place azeotropes at prescribed compositions in binary patchy mixtures, including an 'ideal' design that is azeotropic at every concentration.

desk verdict Extends prior azeotrope work with a clean 'ideal azeotrope' rule and an energy-tuning protocol, but the bond-energy validation compares theory at T=0.08 to simulation at T=0.098. read the letter →

arxiv 2506.11856 v1 pith:2VPFNSTQ submitted 2025-06-13 cond-mat.soft

classification cond-mat.soft
keywords inversethermodynamicsazeotropypatchyparticlesWertheimperturbationtheoryGibbsensemblesimulationphasebehaviordesignbinarymixturesbond-energytuning
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper aims to turn azeotropy—the condition where coexisting liquid and vapor have identical composition—into a designable property of mixtures of patchy colloids, spherical particles decorated with a small number of directional binding sites. Using Wertheim first-order perturbation theory, a closed-form theory of fluids with directional bonds, as the design engine, it derives two control knobs: bond topology, which sets which patches bind, and bond energies, which set how strongly they bind. The central result is that if every patch has exactly one complementary bonding partner on each species, the law of mass action becomes independent of composition, so the mixture shows an azeotrope at every concentration $x_A \in [0,1]$—a state the authors call ideal azeotropy. They also show that changing a single bond energy from $\epsilon$ to $\epsilon'=1.35$ in the N2c8 mixture—two species with eight distinct patch types—moves its azeotropic point from $x=0.5$ to $x=0.6$, and Gibbs-ensemble simulations confirm both predictions. If correct, this provides a route to prescribe not just the structure a mixture assembles into, but the thermodynamic window in which it does so.

What carries the argument

The central object is a graph representation of a patchy-particle mixture: species are nodes, each connected to its patch nodes, and edges between patches mark complementary bonds with edge weights equal to bond energies. The thermodynamic engine is Wertheim first-order perturbation theory, whose mass-action equation (7) gives the probability $X_\alpha^{(i)}$ that patch $\alpha$ on species $i$ is unbonded. The design condition for an ideal azeotrope is that every patch has one complementary partner on each species, so the mass-action equation becomes independent of molar fraction and the mixture's free energy reduces to a single-component form; the energy-tuning variant iterates the bond energy while tracking the chemical-potential signature of coexistence to move the azeotropic point to a target concentration.

What would settle it

Run Gibbs-ensemble simulations of the energy-shifted N2c8 mixture at $T=0.098$ and $\rho=0.1$ with starting concentrations $x=0.1$, $0.3$, $0.5$, $0.7$, and $0.9$: if any box pair other than $x=0.6$ keeps its initial concentration in both phases, the predicted azeotropic shift fails. For the ideal-azeotrope design, repeat the simulation at $T=0.138$ with initial densities away from $\rho=0.2$, say $0.15$ and $0.25$; if coexisting concentrations drift from their starting values, the concentration independence of Eq. (8) is not realized.

Watch

Extended reading notes

Core claim

On the paper's own terms, the discovery is that azeotropic demixing can be programmed from the microscopic interaction graph rather than left as a property that emerges. The argument runs through the law of mass action for patch bonding: the azeotropic solution $X_\alpha^{(i)} = X$ makes the mixture's Wertheim free energy reduce to that of a single-component fluid. When each patch has $N_s$ complementary partners, one per species, the mass-action equation collapses to $X + \phi X^2 \Delta - 1 = 0$, which contains no molar fraction, so the system is azeotropic at every composition. The paper further shows that bond-energy asymmetry can shift the azeotropic concentration: in the N2c8 system, increasing one self-complementary bond energy to $\epsilon'=1.35$ relocates the azeotrope from $x=0.5$ to $x=0.6$, a prediction verified by Gibbs-ensemble Monte Carlo simulations, and chemical-potential measurements via the $S_0$ method confirm the near-ideal, single-component-like behavior of the always-azeotropic design.

Load-bearing premise

The load-bearing premise is that the mass-action equation (7) from first-order thermodynamic perturbation theory gives quantitatively accurate bonding statistics and free energies for these directionally bonded patchy colloids; if that approximate theory misestimates bonding at the simulated or experimental state points, the designed azeotropic compositions would drift or vanish.

Editorial extensions

If this is right

  • A mixture designed with the ideal-azeotrope topology behaves as a single-component fluid in coexistence: bubble and dew curves lie parallel to the concentration axis and tie-lines stay vertical, regardless of the starting concentration.
  • By placing the azeotrope at the stoichiometric composition of a target crystal, the coexisting liquid can be made to match the crystal's composition, which the authors argue optimizes two-step nucleation pathways.
  • The bond-energy strategy works without changing the bonding topology, so azeotrope position can be tuned in designs where connectivity is fixed by other constraints.
  • Validation requires only a few Gibbs-ensemble runs that check whether concentration stays constant in coexisting boxes, rather than a full phase-diagram computation.
  • Chemical potentials measured from structure factors are nearly ideal at all concentrations for the always-azeotropic mixture, so composition changes barely alter the excess chemical potential.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Beyond the paper: the ideal-azeotropy condition should generalize to more than two species simply by requiring every patch to have one complementary partner on each species, suggesting a family of mixtures that are azeotropic across all mixing ratios.
  • Beyond the paper: the energy-tuning iteration amounts to an inverse-design protocol in which the azeotrope location is used as a measurable target to calibrate interaction parameters, a route that could transfer to DNA-coated colloids where bond strengths are adjustable.
  • Beyond the paper: an always-azeotropic mixture preserves composition during condensation or evaporation, which could make it useful for processes requiring compositional uniformity, such as drying of multicomponent dispersions.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The manuscript presents an 'inverse thermodynamics' framework for binary mixtures of patchy particles, in which interaction topologies and bond energies are designed to realize targeted azeotropic behavior. Using Wertheim first-order perturbation theory and Gibbs-ensemble simulations, the authors identify two design strategies: (i) a 'bond-exclusivity' topology yielding an ideal azeotrope, where the azeotropic condition holds at every concentration (Eq. 8), and (ii) a bond-energy modulation algorithm that is claimed to shift the azeotropic composition of the N2c8 mixture from x=0.5 to x=0.6 by setting a self-complementary bond energy to epsilon'=1.35. The ideal-azeotrope prediction is validated by Gibbs-ensemble simulations showing concentration-independent coexistence densities and temperatures, and by S0-method chemical-potential calculations. The energy-shift prediction is tested by Gibbs-ensemble simulations at T=0.098, while the theoretical phase diagrams of Fig. 8 are computed at T=0.08.

Significance. If the results hold, the paper extends inverse design from structural targets to thermodynamic conditions, which is a conceptually valuable step for programmable self-assembly. The ideal-azeotrope construction is elegant and is supported by independent simulations, not by post-hoc fitting. The use of Wertheim theory and Gibbs-ensemble validation is methodologically sound, and the connection to practical self-assembly (matching azeotrope composition to crystal stoichiometry) gives the work practical relevance. However, the evidence for the energy-control route is limited: only one target composition is demonstrated, and the theory/simulation comparison is performed at different temperatures. The claim of 'any prescribed composition' is therefore broader than what is currently established.

major comments (3)
  1. [Azeotropy from bond-energy control; Figs. 8 and 9] The theoretical prediction that setting epsilon' = 1.35 shifts the N2c8 azeotrope from x=0.5 to x=0.6 is presented in Fig. 8 at T=0.08, but the validating Gibbs-ensemble simulations in Fig. 9 are run at T=0.098. Since azeotropic composition generally varies with temperature, observing concentration stability at x=0.6 at T=0.098 does not quantitatively confirm the T=0.08 prediction. To support the central claim of bond-energy control, provide either a theoretical phase diagram at T=0.098 or a Gibbs-ensemble simulation at T=0.08, or otherwise discuss why the azeotropic composition is expected to be temperature-independent in this range.
  2. [Abstract and Conclusions; 'Azeotropy from bond-energy control'] The manuscript repeatedly states that the method can place an azeotrope at 'any prescribed composition,' but the energy-control route is demonstrated for a single target (x=0.6). The iterative algorithm described in the text and Fig. 7 is a trial-and-error procedure without an explicit proof of convergence or reachability for arbitrary x_goal. Demonstrating at least one additional target composition (e.g., x=0.4 or x=0.7) or providing a formal argument for the generality of the approach would be necessary to substantiate the 'any' claim.
  3. [Ideal Azeotropy, Eq. (8)] The derivation of the concentration-independent mass-action equation (8) is stated rather than shown. The text says that when each patch has Ns complementary patches, one per species, the mass-balance equation becomes independent of molar fractions, but the intermediate steps (summation over complementary patches, cancellation of x_i terms) are omitted. A short derivation would clarify the conditions under which the ideal-azeotrope solution X_alpha^(i)=X is self-consistent, especially because the argument is the basis for the 'always azeotropic' claim.
minor comments (5)
  1. [Methods, Eq. (6) and the following paragraph] The text says 'X_alpha^(i) is defined by the law of mass action in Eq. (6)', but Eq. (6) is the free-energy expression; the mass-action equation is later introduced as Eq. (7). Please correct the cross-reference.
  2. [Methods, second paragraph] There is a typo: 'simualations' should be 'simulations'.
  3. [Fig. 8 caption] The caption reads 'Phase diagrams of N2c8 mixtures for different bond-energy .' with a missing word; it should be 'for different bond energies'.
  4. [Fig. 9 caption and text] The text states that initial concentrations x=0.2, 0.4, 0.6, 0.8 were simulated, but Fig. 9 appears to show results only for x=0.2, 0.4, and 0.6. Please clarify whether x=0.8 was simulated and, if so, include it in the figure or remove it from the text.
  5. [Fig. 4 caption] The caption says 'temperature density phase diagrams' but the axes are likely temperature and density; please make the axis labels explicit in the caption.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the azeotrope predictions are derived from Wertheim theory and independently checked by Gibbs-ensemble simulations; the only self-citation supplies a starting point and is not load-bearing, though the energy-shift validation compares T=0.08 theory with T=0.098 simulation.

full rationale

The central claims are not circular. The ideal-azeotropy prediction follows from the mass-action equation (Eq. 7): under the stated graph condition that each patch has one complementary partner on each species, and with the symmetric bonding-probability condition X_alpha^(i) = X, the equation reduces to Eq. 8, X + phi X^2 Delta - 1 = 0, which contains no molar fraction x_A. The authors then predict concentration-independent coexistence, and this is tested by independent Gibbs-ensemble simulations (Figs. 3 and 4) and by S0-method chemical-potential measurements (Fig. 6), not by refitting the theory. The energy-control result is likewise not a fitted input renamed as a prediction: epsilon' = 1.35 is chosen as a designed input, the azeotrope location x = 0.6 is a theoretical output, and the Gibbs-ensemble simulation is a separate check. The paper's reliance on the authors' prior work [45] for the known equimolar azeotrope of N2c8 is a minor, non-load-bearing self-citation: it provides the starting point of the iterative algorithm, but the shifted azeotrope is verified independently. One validation gap exists and is a correctness concern rather than circularity: the Wertheim phase diagrams in Fig. 8 are computed at T = 0.08, while the confirming Gibbs-ensemble simulations in Fig. 9 are run at T = 0.098, with no temperature-matched theory curve or simulation at the prediction temperature; since azeotrope composition can depend on temperature, this weakens quantitative confirmation but does not make the derivation equivalent to its inputs. No circular step of the enumerated kinds is present.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

No new physical entities are postulated. The graph designs are new interaction topologies, not new particles or forces; they are specified by patch-patch binding rules and energies within the known Kern-Frenkel model.

free parameters (3)
  • patch binding energy scale epsilon' for N2c8 self-complementary pairs = 1.35 times baseline epsilon
    Chosen by the design algorithm to target the azeotropic composition x=0.6; it is a design variable rather than an empirical fit, but the quantitative prediction depends on this value.
  • Kern-Frenkel square-well width delta = not reported in text
    The pair potential in Eq. 1 depends on delta; the numerical value is not given, presumably inherited from prior work, but needed for reproduction.
  • patch angular width theta_max = not reported in text
    The orientation function F in Eq. 2 depends on theta_max; the value is not stated in the paper.
assumptions (5)
  • domain assumption Wertheim first-order perturbation theory provides an accurate free energy for patchy particle mixtures.
    Used in Eq. 6 and throughout the theoretical predictions; it is an approximate statistical-mechanical theory, not exact.
  • domain assumption The Kern-Frenkel potential faithfully represents the interactions of patchy colloids, including DNA-functionalized particles.
    Adopted in Methods Eqs. 1-2 as the simulation model; experimental correspondence is asserted but not demonstrated here.
  • domain assumption The law of mass action (Eq. 7) correctly describes bonding probabilities with at most one bond per patch.
    Central to deriving the ideal azeotrope and the phase diagrams; this is the Wertheim single-bond-per-site approximation.
  • domain assumption Isochoric thermodynamics (refs 70,71) can be used to draw binodals from the Wertheim free energy.
    Invoked in Methods to construct pressure-concentration and density-concentration diagrams.
  • domain assumption The S0 method relates zero-wavenumber structure factors to excess chemical potentials.
    Used in Methods and Fig. 6 to measure chemical potentials; it relies on Kirkwood-Buff theory.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Inverse Thermodynamics: Designing Interactions for Targeted Phase Behavior." pith.science (2026). https://pith.science/paper/2VPFNSTQ

@misc{pith2026250611856,
  author       = {Pith},
  title        = {Pith review of: Inverse Thermodynamics: Designing Interactions for Targeted Phase Behavior},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/2VPFNSTQ}},
  note         = {Machine review of arXiv:2506.11856}
}
read the original abstract

The traditional goal of inverse self-assembly is to design interactions that drive particles toward a desired target structure. However, achieving successful self-assembly also requires tuning the thermodynamic conditions under which the structure is stable. In this work, we extend the inverse design paradigm to explicitly address this challenge by developing a framework for inverse thermodynamics, i.e. the design of interaction potentials that realize specific thermodynamic behavior. As a step in this direction, using patchy particle mixtures as a model system, we demonstrate how precise control over both bonding topology and bond energetics enables the programming of targeted phase behavior. In particular, we establish design principles for azeotropic demixing and show how to create mixtures that exhibit azeotropy at any prescribed composition. Our predictions are validated through Gibbs-ensemble simulations [A.Z. Panagiotopoulos,Molecular Physics 61, 813-826 (1987)]. These results highlight the necessity of coupling structural design with thermodynamic engineering, and provide a blueprint for controlling complex phase behavior in multi-component systems.

Figures

Figures reproduced from arXiv: 2506.11856 by the authors.

Figure 1
Figure 1. FIG. 1: N2c8 design [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: Ideal azeotropic binary mixture. [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3: Gibbs ensemble simulations of the always [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figures from the paper (7 more)
Figure 4
Figure 4. Figure 4: FIG. 4: Density concentration and temperature [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]
Figure 6
Figure 6. Figure 6: To determine the standard chemical potential for [PITH_FULL_IMAGE:figures/full_fig_p005_6.png]
Figure 5
Figure 5. Figure 5: FIG. 5: Nucleation of the always azeotropic binary mixture. [PITH_FULL_IMAGE:figures/full_fig_p006_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6: S0 method results [PITH_FULL_IMAGE:figures/full_fig_p006_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7: Chemical potential behaviors near and [PITH_FULL_IMAGE:figures/full_fig_p007_7.png]
Figure 9
Figure 9. Figure 9: FIG. 9: Gibbs ensemble simulations for the new [PITH_FULL_IMAGE:figures/full_fig_p008_9.png]
Figure 8
Figure 8. Figure 8: FIG. 8: Phase diagrams of N2c8 mixtures for dif [PITH_FULL_IMAGE:figures/full_fig_p008_8.png]

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

84 extracted references · 76 canonical work pages

  1. [1]

    Preisler, T

    Z. Preisler, T. Vissers, F. Smallenburg, and F. Sciortino, Crystals of janus colloids at various interaction ranges, The Journal of Chemical Physics 145 (2016)

  2. [2]

    T. A. Nguyen, A. Newton, D. J. Kraft, P. G. Bolhuis, and P. Schall, Tuning patchy bonds induced by critical casimir forces, Materials 10, 1265 (2017)

  3. [3]

    Kumar, G

    R. Kumar, G. M. Coli, M. Dijkstra, and S. Sastry, Inverse design of charged colloidal particle interactions for self assembly into specified crystal structures, The Journal of chemical physics 151 (2019)

  4. [4]

    Dasgupta, G

    T. Dasgupta, G. M. Coli, and M. Dijkstra, Tuning the glass transition: enhanced crystallization of the laves phases in nearly hard spheres, ACS nano14, 3957 (2020)

  5. [5]

    Dwyer, T

    T. Dwyer, T. C. Moore, J. A. Anderson, and S. C. Glotzer, Tunable assembly of host–guest colloidal crys- tals, Soft Matter 19, 7011 (2023)

  6. [6]

    Van Saarloos, V

    W. Van Saarloos, V. Vitelli, and Z. Zeravcic, Soft Mat- ter: Concepts, Phenomena, and Applications (Princeton University Press, 2024)

  7. [7]

    Sacanna and D

    S. Sacanna and D. J. Pine, Shape-anisotropic colloids: Building blocks for complex assemblies, Current opinion in colloid & interface science 16, 96 (2011)

  8. [8]

    Chakraborty, V

    I. Chakraborty, V. Meester, C. van der Wel, and D. J. Kraft, Colloidal joints with designed motion range and tunable joint flexibility, Nanoscale 9, 7814 (2017)

Show all 84 references
  1. [9]

    Chintha, S

    D. Chintha, S. K. Veesam, E. Boattini, L. Filion, and S. N. Punnathanam, Modeling of effective interactions between ligand coated nanoparticles through symmetry functions, The Journal of Chemical Physics 155 (2021)

  2. [10]

    Chakraborty, D

    I. Chakraborty, D. J. Pearce, R. W. Verweij, S. C. Matysik, L. Giomi, and D. J. Kraft, Self-assembly dy- namics of reconfigurable colloidal molecules, ACS nano 16, 2471 (2022)

  3. [11]

    Zhang, G

    Y. Zhang, G. Giunta, H. Liang, and M. Dijkstra, Shape- induced crystallization of binary dna-functionalized nanocubes, The Journal of Chemical Physics 158 (2023)

  4. [12]

    J. A. Diaz A, J. S. Oh, G.-R. Yi, and D. J. Pine, Photo- printing of faceted DNA patchy particles, Proceedings of the National Academy of Sciences 117, 10645 (2020)

  5. [13]

    McMullen, S

    A. McMullen, S. Hilgenfeldt, and J. Brujic, DNA self- organization controls valence in programmable colloid de- sign, Proceedings of the National Academy of Sciences 118, e2112604118 (2021)

  6. [14]

    H. Liu, M. Matthies, J. Russo, L. Rovigatti, R. P. Narayanan, T. Diep, D. McKeen, O. Gang, N. Stephanopoulos, F. Sciortino, et al., Inverse design of a pyrochlore lattice of dna origami through model-driven experiments, Science 384, 776 (2024)

  7. [15]

    W. Zhou, Y. Lim, H. Lin, S. Lee, Y. Li, Z. Huang, J. S. Du, B. Lee, S. Wang, A. S´ anchez-Iglesias,et al., Colloidal quasicrystals engineered with dna, Nature materials 23, 424 (2024)

  8. [16]

    D. J. Kraft, J. Hilhorst, M. A. Heinen, M. J. Hoogen- raad, B. Luigjes, and W. K. Kegel, Patchy polymer col- loids with tunable anisotropy dimensions, The Journal of Physical Chemistry B 115, 7175 (2011)

  9. [17]

    D. J. Kraft, R. Ni, F. Smallenburg, M. Hermes, K. Yoon, D. A. Weitz, A. van Blaaderen, J. Groenewold, M. Dijk- stra, and W. K. Kegel, Surface roughness directed self- assembly of patchy particles into colloidal micelles, Pro- ceedings of the National Academy of Sciences 109, 1...

  10. [18]

    W. Zhou, Y. Li, K. Je, T. Vo, H. Lin, B. E. Partridge, Z. Huang, S. C. Glotzer, and C. A. Mirkin, Space-tiled colloidal crystals from DNA-forced shape-complementary polyhedra pairing, Science 383, 312 (2024)

  11. [19]

    Russo, P

    J. Russo, P. Tartaglia, and F. Sciortino, Reversible gels of patchy particles: role of the valence, The Journal of chemical physics 131 (2009)

  12. [20]

    Swinkels, R

    P. Swinkels, R. Sinaasappel, Z. Gong, S. Sacanna, W. Meyer, F. Sciortino, and P. Schall, Networks of limited-valency patchy particles, Physical Review Letters 132, 078203 (2024)

  13. [21]

    Sciortino, A

    F. Sciortino, A. Giacometti, and G. Pastore, Phase di- agram of janus particles, Physical review letters 103, 237801 (2009)

  14. [22]

    Russo, J

    J. Russo, J. Tavares, P. Teixeira, M. Telo da Gama, and F. Sciortino, Reentrant phase diagram of network fluids, Physical review letters 106, 085703 (2011)

  15. [23]

    Rold´ an-Vargas, F

    S. Rold´ an-Vargas, F. Smallenburg, W. Kob, and F. Sciortino, Phase diagram of a reentrant gel of patchy particles, The Journal of chemical physics 139 (2013)

  16. [24]

    Bianchi, J

    E. Bianchi, J. Largo, P. Tartaglia, E. Zaccarelli, and F. Sciortino, Phase diagram of patchy colloids: Towards empty liquids, Physical review letters 97, 168301 (2006)

  17. [25]

    M. F. Hagan and G. M. Grason, Equilibrium mechanisms of self-limiting assembly, Reviews of modern physics 93, 025008 (2021)

  18. [26]

    Dijkstra and E

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

  19. [27]

    Semwal, C

    S. Semwal, C. Clowe-Coish, I. Saika-Voivod, and A. Yethiraj, Tunable colloids with dipolar and depletion interactions: Toward field-switchable crystals and gels, Physical Review X 12, 041021 (2022)

  20. [28]

    C. L. Kennedy, D. Sayasilpi, P. Schall, and J.-M. Meijer, Self-assembly of colloidal cube superstructures with crit- ical casimir attractions, Journal of Physics: Condensed Matter 34, 214005 (2022)

  21. [29]

    Marino, R

    E. Marino, R. A. LaCour, T. C. Moore, S. W. van Dongen, A. W. Keller, D. An, S. Yang, D. J. Rosen, G. Gouget, E. H. Tsai, et al. , Crystallization of bi- nary nanocrystal superlattices and the relevance of short- range attraction, Nature Synthesis 3, 111 (2024)

  22. [30]

    T. E. Videbæk, D. Hayakawa, G. M. Grason, M. F. Ha- gan, S. Fraden, and W. B. Rogers, Economical routes to size-specific assembly of self-closing structures, Science advances 10, eado5979 (2024)

  23. [31]

    Shelke, D

    Y. Shelke, D. J. Pearce, and D. J. Kraft, Optimal self- assembly pathways towards colloidal lattices with tun- able flexibility, arXiv preprint arXiv:2406.19260 (2024)

  24. [32]

    F. Lu, Y. Zhang, T. Dwyer, A. Michelson, T. C. Moore, H. Yan, K. Kisslinger, H. Zhang, X. Chen, S. C. Glotzer, et al., Octo-diamond crystal of nanoscale tetrahedra with interchanging chiral motifs, Nature materials , 1 (2025). 10

  25. [33]

    De Yoreo, S

    J. De Yoreo, S. Y. Schmid, B. Helfrecht, A. Stegmann, B. Legg, H. Pyles, J. Chen, J. R. Edison, M. Ziatdinov, Z. Preisler, et al. , Impact of colloidal forces and broken symmetry on the design of a protein-material interface, (2025)

  26. [34]

    L. B. Frechette, N. Sundararajan, F. Caballero, A. Tru- biano, and M. F. Hagan, Computer simulations show that liquid-liquid phase separation enhances self-assembly, arXiv preprint arXiv:2505.10400 (2025)

  27. [35]

    H. N. Lekkerkerker, W.-K. Poon, P. N. Pusey, A. Stroobants, and P. . Warren, Phase behaviour of col- loid+ polymer mixtures, Europhysics Letters 20, 559 (1992)

  28. [36]

    E. J. Meijer and D. Frenkel, Colloids dispersed in polymer solutions. a computer simulation study, The Journal of chemical physics 100, 6873 (1994)

  29. [37]

    Hemmer, E

    P. Hemmer, E. Velasco, L. Mederos, G. Navascu´ es, and G. Stell, Solid-solid transitions induced by repulsive in- teractions, The Journal of Chemical Physics 114, 2268 (2001)

  30. [38]

    Rold´ an-Vargas, F

    S. Rold´ an-Vargas, F. Smallenburg, W. Kob, and F. Sciortino, Gelling by heating, Scientific reports3, 2451 (2013)

  31. [39]

    W. B. Rogers and V. N. Manoharan, Programming col- loidal phase transitions with DNA strand displacement, Science 347, 639 (2015)

  32. [40]

    Bomboi, F

    F. Bomboi, F. Romano, M. Leo, J. Fernandez-Castanon, R. Cerbino, T. Bellini, F. Bordi, P. Filetici, and F. Sciortino, Re-entrant dna gels, Nature communica- tions 7, 13191 (2016)

  33. [41]

    Hvozd, Y

    T. Hvozd, Y. V. Kalyuzhnyi, V. Vlachy, and P. T. Cum- mings, Empty liquid state and re-entrant phase behav- ior of the patchy colloids confined in porous media, The Journal of Chemical Physics 156 (2022)

  34. [42]

    Russo, F

    J. Russo, F. Leoni, F. Martelli, and F. Sciortino, The physics of empty liquids: From patchy particles to water, Reports on Progress in Physics 85, 016601 (2022)

  35. [43]

    Sciortino and E

    F. Sciortino and E. Zaccarelli, Equilibrium gels of limited valence colloids, Current opinion in colloid & interface science 30, 90 (2017)

  36. [44]

    Beneduce, D

    C. Beneduce, D. EP Pinto, P. ˇSulc, F. Sciortino, and J. Russo, Two-step nucleation in a binary mixture of patchy particles, The Journal of Chemical Physics 158 (2023)

  37. [45]

    Beneduce, F

    C. Beneduce, F. Sciortino, P. ˇSulc, and J. Russo, Engi- neering azeotropy to optimize the self-assembly of col- loidal mixtures, ACS nano 17, 24841 (2023)

  38. [46]

    Russo and M

    G. Russo and M. Lattuada, Synthesis of patchy colloids with different chemical functionalities, Available at SSRN 4641669 (2024)

  39. [47]

    N. C. Seeman and H. F. Sleiman, Dna nanotechnology, Nature Reviews Materials 3, 1 (2017)

  40. [48]

    Xiong, S

    Y. Xiong, S. Yang, Y. Tian, A. Michelson, S. Xiang, H. Xin, and O. Gang, Three-dimensional patterning of nanoparticles by molecular stamping, ACS nano14, 6823 (2020)

  41. [49]

    Y. Li, W. Zhou, Y. Zhou, H. F. Cheng, B. Lee, X. Hu, E. W. Roth, V. P. Dravid, S. C. Glotzer, and C. A. Mirkin, Cocrystals combining order and correlated dis- order via colloidal crystal engineering with dna, Science Advances 11, eadu4919 (2025)

  42. [50]

    Zhang, C

    T. Zhang, C. Hartl, K. Frank, A. Heuer-Jungemann, S. Fischer, P. C. Nickels, B. Nickel, and T. Liedl, 3d dna origami crystals, Adv. Mater. 30, 1800273 (2018)

  43. [51]

    S. H. Park, H. Park, K. Hur, and S. Lee, Design of dna origami diamond photonic crystals, ACS Applied Bio Materials 3, 747 (2019)

  44. [52]

    Hayakawa, T

    D. Hayakawa, T. E. Videbaek, D. M. Hall, H. Fang, C. Sigl, E. Feigl, H. Dietz, S. Fraden, M. F. Hagan, G. M. Grason, et al., Geometrically programmed self-limited as- sembly of tubules using dna origami colloids, Proceedings of the National Academy of Sciences 119, e2207902119 (2022)

  45. [53]

    J. S. Kahn, D. C. Redeker, A. Michelson, A. Tkachenko, S. Hong, B. Minevich, and O. Gang, Arbitrary design of DNA-programmable 3d crystals through symmetry map- ping, ACS nano (2025)

  46. [54]

    R. Saha, D. Hayakawa, T. E. Videbaek, M. Price, W.-S. Wei, J. Pombo, D. Duke, G. Arya, G. M. Grason, W. B. Rogers, et al. , Modular programming of interaction and geometric specificity enables assembly of complex dna origami nanostructures, arXiv preprint arXiv:2502.05388 (2025)

  47. [55]

    Neophytou, F

    A. Neophytou, F. Sciortino, and J. Russo, Designing the self-assembly of disordered materials via color frustra- tion, Advanced Materials (2025 in pressing)

  48. [56]

    D. E. Pinto, P. Sulc, F. Sciortino, and J. Russo, Automat- ing blueprints for the assembly of colloidal quasicrystal clusters, ACS nano 19, 512 (2024)

  49. [57]

    D. E. Pinto, N. A. Ara´ ujo, P.ˇSulc, and J. Russo, Inverse design of self-folding 3d shells, Physical Review Letters 132, 118201 (2024)

  50. [58]

    C. Wang, A. P´ erez de Alba Ortz ´ ı, and M. Dijkstra, In- verse design method with enhanced sampling for complex open crystals: Application to novel zeolite self-assembly, ACS nano (2025)

  51. [59]

    C. Vega, E. Sanz, J. Abascal, and E. Noya, Determina- tion of phase diagrams via computer simulation: method- ology and applicationsto water, electrolytes and proteins, Journal of Physics: Condensed Matter 20, 153101 (2008)

  52. [60]

    de Las Heras, J

    D. de Las Heras, J. M. Tavares, and M. M. T. da Gama, Phase diagrams of binary mixtures of patchy colloids with distinct numbers of patches: the network fluid regime, Soft Matter 7, 5615 (2011)

  53. [61]

    Kern and D

    N. Kern and D. Frenkel, Fluid–fluid coexistence in col- loidal systems with short-ranged strongly directional at- traction, The Journal of chemical physics 118, 9882 (2003)

  54. [62]

    Bol, Monte Carlo simulations of fluid systems of wa- terlike molecules, Molecular Physics 45, 605 (1982)

    W. Bol, Monte Carlo simulations of fluid systems of wa- terlike molecules, Molecular Physics 45, 605 (1982)

  55. [63]

    Chen and J

    B. Chen and J. I. Siepmann, A novel Monte Carlo algo- rithm for simulating strongly associating fluids: Appli- cations to water, hydrogen fluoride, and acetic acid, The Journal of Physical Chemistry B 104, 8725 (2000)

  56. [64]

    Rovigatti, J

    L. Rovigatti, J. Russo, and F. Romano, How to simulate patchy particles, The European Physical Journal E 41, 59 (2018)

  57. [65]

    Cheng, Computing chemical potentials of solutions from structure factors, The Journal of Chemical Physics 157 (2022)

    B. Cheng, Computing chemical potentials of solutions from structure factors, The Journal of Chemical Physics 157 (2022)

  58. [66]

    Barrat and J.-P

    J.-L. Barrat and J.-P. Hansen, Basic concepts for simple and complex liquids (Cambridge University Press, 2003)

  59. [67]

    Wertheim, Fluids with highly directional attractive forces

    M. Wertheim, Fluids with highly directional attractive forces. i. statistical thermodynamics, J. Stat. Phys. 35, 19 (1984)

  60. [68]

    W. G. Chapman, G. Jackson, and K. E. Gubbins, Phase equilibria of associating fluids: Spherical molecules with multiple bonding sites, Mol. Phys. 65, 1 (1988). 11

  61. [69]

    Jonas, P

    H. Jonas, P. Schall, and P. Bolhuis, Extended wertheim theory predicts the anomalous chain length distributions of divalent patchy particles under extreme confinement, The Journal of Chemical Physics 157 (2022)

  62. [70]

    U. K. Deiters, Differential equations for the calculation of isopleths of multicomponent fluid mixtures, Fluid Phase Equilib. 447, 72 (2017)

  63. [71]

    I. H. Bell and U. K. Deiters, On the construction of bi- nary mixture p-x and t-x diagrams from isochoric ther- modynamics, AIChE J. 64, 2745 (2018)

  64. [72]

    Rovigatti, J

    L. Rovigatti, J. Russo, F. Romano, M. Matthies, L. Kroc, and P. ˇSulc, A simple solution to the problem of self- assembling cubic diamond crystals, Nanoscale 14, 14268 (2022)

  65. [73]

    A. Z. Panagiotopoulos, Direct determination of phase co- existence properties of fluids by Monte Carlo simulation in a new ensemble, Mol. Phys. 61, 813 (1987)

  66. [74]

    A. Z. Panagiotopoulos, N. Quirke, M. Stapleton, and D. Tildesley, Phase equilibria by simulation in the gibbs ensemble: Alternative derivation, generalization and ap- plication to mixture and membrane equilibria, Mol. Phys. 63, 527 (1988)

  67. [75]

    J. M. Smith, H. C. Van Ness, M. M. Abbott, and M. T. Swihart, Introduction to chemical engineering thermody- namics (McGraw-Hill Singapore, 1949)

  68. [76]

    W. J. Moore, Physical Chemistry (Prentice-Hall, 1962)

  69. [77]

    Whitelam and R

    S. Whitelam and R. L. Jack, The statistical mechanics of dynamic pathways to self-assembly, Annual review of physical chemistry 66, 143 (2015)

  70. [78]

    Bupathy, D

    A. Bupathy, D. Frenkel, and S. Sastry, Temperature pro- tocols to guide selective self-assembly of competing struc- tures, Proceedings of the National Academy of Sciences 119, e2119315119 (2022)

  71. [79]

    C. L. Bassani, G. Van Anders, U. Banin, D. Baranov, Q. Chen, M. Dijkstra, M. S. Dimitriyev, E. Efrati, J. Fa- raudo, O. Gang, et al. , Nanocrystal assemblies: Cur- rent advances and open problems, ACS nano 18, 14791 (2024)

  72. [80]

    Hayakawa, T

    D. Hayakawa, T. E. Videbæk, G. M. Grason, and W. B. Rogers, Symmetry-guided inverse design of self- assembling multiscale dna origami tilings, ACS nano 18, 19169 (2024)

  73. [81]

    Michelson, L

    A. Michelson, L. Shani, J. S. Kahn, D. C. Redeker, W.-I. Lee, K. R. DeOlivares, K. Kisslinger, N. Tiwale, H. Yan, A. Pattammattel, et al. , Scalable fabrication of chip-integrated 3d-nanostructured electronic devices via DNA-programmable assembly, Science Advances 11, eadt5620 (2025)

  74. [82]

    Romano, J

    F. Romano, J. Russo, L. Kroc, and P. ˇSulc, Designing patchy interactions to self-assemble arbitrary structures, Physical Review Letters 125, 118003 (2020)

  75. [83]

    Russo, F

    J. Russo, F. Romano, L. Kroc, F. Sciortino, L. Rovigatti, and P. ˇSulc, Sat-assembly: A new approach for designing self-assembling systems, Journal of Physics: Condensed Matter 34, 354002 (2022)

  76. [84]

    D. E. Pinto, P. ˇSulc, F. Sciortino, and J. Russo, De- sign strategies for the self-assembly of polyhedral shells, Proceedings of the National Academy of Sciences 120, e2219458120 (2023)

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.