REVIEW 2 major objections 5 minor 38 references
The underappreciated role of nonspecific interactions in the crystallization of DNA-coated colloids
T0 review · 2 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read This paper shows that nonspecific steric and van der Waals interactions, not only designed DNA binding, determine the equilibrium crystal structure of binary DNA-coated colloids, producing a tunable family of body-centered tetragonal…
desk verdict New continuous BCT polymorph family tuned by nonspecific interactions, with strong experiment-simulation agreement; the theory has a calibrated parameter that needs sensitivity testing. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the dimensionless BCT parameter C, which interpolates between the cubic CsCl (C=0) and CuAu (C=1) unit cells. The argument is carried by the microscopic pair-potential model of DNA-coated colloids (reference [15]), which self-consistently computes the free energy of DNA hybridization, polymer-brush steric repulsion, and van der Waals attraction as a function of particle separation; these potentials feed constant-pressure molecular dynamics that predict equilibrium C. A thermodynamic perturbation theory then decomposes the crystal free energy into a harmonic reference crystal with A-B bonds, an anharmonic correction from same-type A-A/B-B repulsion, and a perturbative van der Waals term, showing how the balance of these small driving forces selects C. The anharmonic term and the vdW term are the decisive machinery: one favors low C, the other high C, and their near-cancellation explains both the non-cubic equilibrium and the sensitivity to design parameters.
What would settle it
Measure the polymer grafting density (for example, by quantifying DNA and PEO coverage per particle) for each particle diameter and molecular weight, feed the measured densities into the microscopic pair-potential model, and check whether the predicted equilibrium C values still match the experimentally observed BCT lattices; if the C-MW trend depends strongly on the assumed constant density, the central claim is falsified.
Extended reading notes
Core claim
For a binary suspension of same-sized, micron-scale DNA-coated colloids, the equilibrium crystal is generically a non-cubic body-centered tetragonal (BCT) lattice, continuously parameterized by C in [0,1], rather than one of the two cubic end members. The paper's central discovery is that this C value is determined by a competition between specific DNA-mediated attraction and nonspecific forces of similar magnitude: steric repulsion between polymer brushes pushes C down, van der Waals attraction between same-type particle pairs pushes C up, and the harmonic phonon entropy of the reference crystal also favors higher C. Tuning polymer molecular weight, particle diameter, or the mixing fraction of complementary DNA strands shifts the balance and continuously moves the lattice between CsCl-like (C~0) and CuAu-like (C~1) structures. Simulations using pair potentials from a microscopic model reproduce the experimental C values nearly quantitatively, and the perturbation theory shows the same trend. The paper concludes that nonspecific interactions are not a nuisance to be ignored but a control knob for programmable self-assembly.
Load-bearing premise
The quantitative match between simulation and experiment assumes that the polymer brush density is identical across all particle sizes and polymer molecular weights, even though each batch is synthesized separately; if grafting density varies with molecular weight or particle size, the predicted pair potentials and C values shift.
Editorial extensions
If this is right
- Same-sized binary DNA-coated colloids can spontaneously form a one-parameter family of BCT crystals, including continuous transitions between CsCl and CuAu, not just the two cubic extremes.
- Changing the polymer brush molecular weight or the particle diameter shifts the equilibrium C parameter by altering how strongly van der Waals attraction acts at the brush-contact distance.
- Van der Waals attraction is not required for CuAu-like BCT crystals; harmonic phonon entropy alone favors high C, so BCT stability survives even when vdW is removed.
- Introducing a DNA-mediated attraction between same-type particles also tunes C continuously from 0.1 to 0.8 as the strand-mixing fraction increases, confirming that any same-type attraction of comparable range can play the vdW role.
- The coexistence theory predicts melting temperatures with absolute errors under 2 degrees Celsius, so the same pair potentials plus perturbation theory can be used to predict both structure and thermal stability.
Reading between the lines
- A general design rule suggested by the paper, though not stated as such: any weak attraction active when same-type particles approach (vdW, DNA hybridization, depletion) should push the equilibrium toward CuAu-like structures, whereas thicker or denser brushes push toward CsCl-like structures.
- The flat free-energy landscape near the midpoint of the continuous transition may explain the batch-to-batch and crystal-to-crystal scatter seen in prior DNA-coated-colloid experiments, and why certain predicted cubic lattices have been hard to realize.
- One testable extension is to replace vdW attraction with a tunable depletion interaction; the model predicts the same continuous C shift, which could be checked in the same 600 nm and 67 kDa system.
- For nanometer-scale DNA-coated particles, where vdW is much weaker, the theory implies a baseline C set by phonon entropy and brush repulsion; dedicated simulations could predict whether BCT intermediates should appear there too.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports a combined experimental, simulation, and theoretical study of binary DNA-coated colloids. Experimentally, 600-nm and 430-nm particles with PEO brushes of four molecular weights form body-centered tetragonal (BCT) crystals whose continuously varying aspect-ratio parameter C is characterized by confocal imaging and RDF matching. Constant-pressure simulations using independent pair potentials from Ref. [15] reproduce the experimental C values (R = 0.95) without fitted parameters. A perturbation-theory model, using a harmonic reference binary crystal and perturbative corrections for van der Waals attraction and anharmonic A-A/B-B repulsion, reproduces the qualitative trends (R = 0.88) and attributes structure selection to a balance among weak harmonic, anharmonic, and vdW contributions. A separate experiment with mixed DNA strands shows a continuous, tunable transition from CsCl-like to CuAu-like structures. Simulation, analysis, and theory scripts are provided on GitHub.
Significance. Strengths: the experimental survey is systematic, the simulation-experiment agreement is quantitative and achieved with independently computed pair potentials, melting temperatures are predicted with R = 0.98, and the paper makes falsifiable predictions (a continuous BCT family and the result that vdW attraction is not required for high C values). The theoretical model is a useful interpretative tool, but its quantitative C predictions are less secure because they rely on a calibrated Gaussian approximation for g_AA,ref; this weakness does not affect the experimental or simulation conclusions. If the theoretical parametrization is made robust, the work will be an important step toward including nonspecific interactions in the design of programmable colloidal crystals.
major comments (2)
- [Appendix B, Eq. (B6), and Fig. 3B] The anharmonic entropy term that stabilizes low-C BCT crystals uses g_AA,ref approximated as a Gaussian whose width is set by a single parameter alpha, chosen from zero-pressure simulations at C = 0. The paper states that this term is essential for stabilizing low-C BCT crystals, yet no sensitivity analysis is reported. Because the predicted equilibrium C is the minimum of mu_coex(C) set by the competition between this anharmonic term and the vdW attraction, a different alpha or an alpha that varies with C, molecular weight, or particle diameter would shift the theoretical C values and the reported R = 0.88 correlation. Please provide a sensitivity analysis over a physically plausible range of alpha, or determine alpha independently across the phase diagram, and specify the units or dimensionality of alpha.
- [Supplementary Information, Sec. II (Solving the crystal structures)] The experimental structure classification uses a look-up table containing only BCC/FCC/BCT facets, so BCT is the only non-cubic structure family the classifier can return. Because the central experimental claim is that the assembled crystals are BCT, this identification is circular unless alternative binary lattices (e.g., simple tetragonal, orthorhombic, or wurtzite-type arrangements) are explicitly ruled out by the RDF data. Please include such alternative structures in the comparison or state clearly that the BCT assignment is an assumption whose credibility rests on the independent simulation agreement shown in Fig. 3.
minor comments (5)
- [Sec. II.E and Appendix B] The symbol alpha is used for the strand-mixing fraction in Sec. II.E and for the Gaussian width scaling factor in Appendix B; this notation conflict should be resolved.
- [Appendix B, Eq. (B6)] The meanings of x and d_AA in Eq. (B6) are not defined: please state the coordinate used in the integrals and define d_AA as the minimum A-A/B-B pair distance, and clarify why the upper limit of both integrals is d_AA.
- [Fig. 2] The caption says the pair potentials are shown as black curves, but the figure uses distinct styling for the two molecular weights; please label the curves directly or in the legend.
- [Fig. 3B and Fig. S3] The reported Pearson correlation coefficients are based on a small number of systems; please state the number of data points for each R value and, for R = 0.88, provide an uncertainty estimate such as a bootstrap confidence interval.
- [Sec. II.C] The sentence stating that C decreases from roughly 0.6 to 0.2 as MW increases applies only to the 600-nm particles; for the 430-nm particles the C values are already low and nearly independent of MW, so please qualify the statement.
Circularity Check
No significant circularity: the central experimental and simulation claims rest on an independent pair-potential model, and the theory's one calibrated width parameter is anchored to a C=0 reference simulation rather than to the experimental C values being predicted.
full rationale
The paper's main claim is supported by a chain that is not self-referential. The pair potentials are taken from an independent microscopic model (Ref. [15], Cui et al., not authored by the present authors), and the constant-pressure simulations are parameter-free predictions from those potentials: they initialize multiple BCT structures and let the box dimensions relax under the predicted pair potentials, then compare the resulting C values with experimental measurements. The experimental C values are determined geometrically from confocal imaging and a BCT lookup table, not from the simulation or theory. The theoretical model does contain one calibrated quantity: the width parameter α in the approximate A-A/B-B radial distribution function, g_AA,ref(x) ∝ exp[−(α k_AB/2)(x−d_AA/2)^2], chosen as α≈1×10^-4 from zero-pressure simulations of the reference crystal with C=0 (Appendix B). This anharmonic term is stated to be essential for stabilizing low-C BCT crystals, so the quantitative C predictions of the theory do depend on this choice; a different α could shift the predicted minima of µ_coex(C). However, this is a model-approximation caveat, not circularity: α is not fitted to the experimental C values being predicted, and the main experimental/simulation comparison does not use α at all. No load-bearing self-citation, imported uniqueness claim, or renaming of a known result is present. The theory's α calibration is a legitimate correctness/robustness concern but does not make the derivation equivalent to its inputs.
Assumptions & free parameters
free parameters (2)
- alpha (Gaussian width scaling factor for g_AA,ref) =
approximately 1 x 10^-4
- Crystal-structure lookup step =
0.025 in C
assumptions (6)
- domain assumption The microscopic pair-potential model of DNA-coated colloids (Cui et al., Ref. [15]) accurately describes the experimental systems, including vdW, steric, and DNA hybridization free energies.
- domain assumption Polymer graft density is constant across all particle sizes and PEO molecular weights.
- domain assumption The observed crystals represent the equilibrium structure at room temperature.
- standard math The A-B pair potential can be approximated as a harmonic spring for the reference crystal free energy.
- ad hoc to paper The A-A/B-B radial distribution function in the reference crystal is Gaussian with width controlled by alpha.
- standard math The fluid phase is described by the Carnahan-Starling hard-sphere equation of state with a mean-field attraction term.
Cite this review
Pith. "Pith review of The underappreciated role of nonspecific interactions in the crystallization of DNA-coated colloids." pith.science (2026). https://pith.science/paper/E5IRMQRL
@misc{pith2026250102220,
author = {Pith},
title = {Pith review of: The underappreciated role of nonspecific interactions in the crystallization of DNA-coated colloids},
year = {2026},
howpublished = {\url{https://pith.science/paper/E5IRMQRL}},
note = {Machine review of arXiv:2501.02220}
}
read the original abstract
Over the last decade, the field of programmable self-assembly has seen an explosion in the diversity of crystal lattices that can be synthesized from DNA-coated colloidal nanometer- and micrometer-scale particles. The prevailing wisdom has been that a particular crystal structure can be targeted by designing the DNA-mediated interactions, to enforce binding between specific particle pairs, and the particle diameters, to control the packing of the various species. In this article, we show that other ubiquitous nonspecific interactions can play equally important roles in determining the relative stability of different crystal polymorphs and therefore what crystal structure is most likely to form in an experiment. For a binary mixture of same-sized DNA-coated colloidal micrometer-scale particles, we show how changing the magnitudes of nonspecific steric and van der Waals interactions gives rise to a family of binary body-centered tetragonal crystals, including both cesium-chloride and copper-gold crystals. Simulations using pair potentials that account for these interactions reproduce our experimental observations quantitatively, and a theoretical model reveals how a subtle balance between specific and nonspecific forces determines the equilibrium crystal structure. These results highlight the importance of accounting for nonspecific interactions in the crystal-engineering design process.
Figures
Reference graph
Works this paper leans on
-
[15]
F. Cui, S. Marbach, J. A. Zheng, M. Holmes-Cerfon, and D. J. Pine, Nat. Commun. 13, 1 (2022)
work page 2022
-
[1]
G. M. Whitesides and B. Grzybowski, Science 295, 2418 (2002)
work page 2002
-
[2]
C. A. Mirkin, R. L. Letsinger, R. C. Mucic, and J. J. Storhoff, Nature 382, 607 (1996)
work page 1996
-
[3]
A. P. Alivisatos, K. P. Johnsson, X. Peng, T. E. Wilson, C. J. Loweth, M. P. Bruchez Jr, and P. G. Schultz, Nature 382, 609 (1996)
work page 1996
-
[4]
W. M. Jacobs and W. B. Rogers, arXiv preprint arXiv:2409.08988 (2024)
arXiv 2024
-
[5]
S. Y. Park, A. K. R. Lytton-Jean, B. Lee, S. Weigand, G. C. Schatz, and C. A. Mirkin, Nature 451, 553 (2008)
work page 2008
-
[6]
D. Nykypanchuk, M. M. Maye, D. van der Lelie, and O. Gang, Nature 451, 549 (2008)
work page 2008
-
[7]
R. J. Macfarlane, B. Lee, M. R. Jones, N. Harris, G. C. Schatz, and C. A. Mirkin, Science 334, 204 (2011)
work page 2011
Show all 38 references
-
[8]
C. R. Laramy, M. N. O’Brien, and C. A. Mirkin, inSpher- ical Nucleic Acids (Jenny Stanford Publishing, 2020) pp. 227–289
2020
-
[9]
P. L. Biancaniello, A. J. Kim, and J. C. Crocker, Phys. Rev. Lett. 94, 058302 (2005)
2005
-
[10]
Y. Wang, Y. Wang, X. Zheng, ´E. Ducrot, J. S. Yodh, M. Weck, and D. J. Pine, Nat. Commun. 6, 7253 (2015)
2015
-
[11]
M. He, J. P. Gales, ´E. Ducrot, Z. Gong, G.-R. Yi, S. Sacanna, and D. J. Pine, Nature 585, 524 (2020)
2020
-
[12]
Hensley, W
A. Hensley, W. M. Jacobs, and W. B. Rogers, Proc. Natl. Acad. Sci. U.S.A. 119, e2114050118 (2022)
2022
-
[13]
Hensley, T
A. Hensley, T. E. Videbæk, H. Seyforth, W. M. Jacobs, and W. B. Rogers, Nat. Commun. 14, 4237 (2023)
2023
-
[14]
W. B. Rogers and J. C. Crocker, Proc. Natl. Acad. Sci. U.S.A. 108, 15687 (2011)
2011
-
[16]
Liang, N
Y. Liang, N. Hilal, P. Langston, and V. Starov, Adv. Coll. Int. Sci. 134, 151 (2007)
2007
-
[17]
J. N. Israelachvili, Intermolecular and surface forces (Academic press, 2011)
2011
-
[18]
R. Mao, B. Minevich, D. McKeen, Q. Chen, F. Lu, O. Gang, and J. Mittal, Proc. Natl. Acad. Sci. U.S.A. 120, e2302037120 (2023)
2023
-
[19]
Y. Wang, Y. Wang, X. Zheng, ´E. Ducrot, M.-G. Lee, G.- R. Yi, M. Weck, and D. J. Pine, J. Am. Chem. Soc. 137, 10760 (2015)
2015
-
[20]
A. V. Tkachenko, Proc. Natl. Acad. Sci. U.S.A. 113, 10269 (2016)
2016
-
[21]
J. S. Oh, Y. Wang, D. J. Pine, and G.-R. Yi, Chemistry of Materials 27, 8337 (2015)
2015
-
[22]
J. C. Crocker and D. G. Grier, J. Coll. Int. Sci. 179, 298 (1996)
1996
-
[23]
M. T. Casey, R. T. Scarlett, W. Benjamin Rogers, I. Jenkins, T. Sinno, and J. C. Crocker, Nat. Commun. 3, 1209 (2012)
2012
-
[24]
H. Hu, P. S. Ruiz, and R. Ni, Phys. Rev. Lett. 120, 048003 (2018)
2018
-
[25]
SantaLucia Jr and D
J. SantaLucia Jr and D. Hicks, Annu. Rev. Biophys. Biomol. Struct. 33, 415 (2004)
2004
-
[26]
H. Fang, M. F. Hagan, and W. B. Rogers, Proc. Natl. Acad. Sci. U.S.A. 117, 27927 (2020)
2020
-
[27]
Zhang, S
Y. Zhang, S. Pal, B. Srinivasan, T. Vo, S. Kumar, and O. Gang, Nature materials 14, 840 (2015)
2015
-
[28]
R. V. Thaner, Y. Kim, T. I. Li, R. J. Macfarlane, S. T. Nguyen, M. Olvera de la Cruz, and C. A. Mirkin, Nano letters 15, 5545 (2015)
2015
-
[29]
W. B. Rogers, W. M. Shih, and V. N. Manoharan, Nat. Rev. Mater. 1, 16008 (2016)
2016
-
[30]
Elser, Phys
V. Elser, Phys. Rev. E 89, 052404 (2014)
2014
-
[31]
THE UNDERAPPRECIATED ROLE OF NONSPECIFIC INTERACTIONS IN THE CRYSTALLIZATION OF DNA-COATED COLLOIDS
J.-P. Hansen and I. R. McDonald, Theory of simple liq- uids: with applications to soft matter (Academic press, 2013). 1 SUPPLEMENTARY INFORMATION FOR “THE UNDERAPPRECIATED ROLE OF NONSPECIFIC INTERACTIONS IN THE CRYSTALLIZATION OF DNA-COATED COLLOIDS” I. MATERIALS AND METHODS ...
2013 arXiv
-
[32]
J. S. Oh, Y. Wang, D. J. Pine, and G.-R. Yi, High-density PEO-b-DNA brushes on polymer particles for colloidal superstructures, Chem. Mater. 27, 8337 (2015)
2015
-
[33]
Hensley, T
A. Hensley, T. E. Videbæk, H. Seyforth, W. M. Jacobs, and W. B. Rogers, Macroscopic photonic single crystals via seeded growth of DNA-coated colloids, Nat. Commun. 14, 4237 (2023)
2023
-
[34]
F. Cui, S. Marbach, J. A. Zheng, M. Holmes-Cerfon, and D. J. Pine, Comprehensive view of microscopic interactions between DNA-coated colloids, Nat. Commun. 13, 1 (2022)
2022
-
[35]
Hensley, W
A. Hensley, W. M. Jacobs, and W. B. Rogers, Self-assembly of photonic crystals by controlling the nucleation and growth of dna-coated colloids, Proc. Natl. Acad. Sci. U.S.A. 119, e2114050118 (2022)
2022
-
[36]
J. C. Crocker and D. G. Grier, Methods of digital video microscopy for colloidal studies, J. Coll. Int. Sci. 179, 298 (1996)
1996
-
[37]
Chandler, Introduction to Modern Statistical Mechanics (Oxford University Press, 1987)
D. Chandler, Introduction to Modern Statistical Mechanics (Oxford University Press, 1987)
1987
-
[38]
Hansen and I
J.-P. Hansen and I. R. McDonald, Theory of simple liquids: with applications to soft matter (Academic press, 2013). 4 IV. ADDITIONAL FIGURES NaCl concentration (mM) 0 100 200 300 400 500 600 0.0 0.2 0.4 0.6 0.8 1.0 600nm 6.5kDa PEO C FIG. S1. Effect of NaCl concentration on th...
2013
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.