pith. sign in

arxiv: 1907.03148 · v2 · pith:FAPQCXLGnew · submitted 2019-07-06 · ❄️ cond-mat.soft · cond-mat.stat-mech

Morphology Formation in Binary Mixtures upon Gradual Destabilisation

Charley Schaefer , Stefan Paquay , Tom C. B. McLeish This is my paper

Pith reviewed 2026-05-25 01:33 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.stat-mech
keywords phase separationbinary mixturecritical pointquench ratedynamic exponentbrownian motionmonte carlomolecular dynamics
0
0 comments X

The pith

In 2D binary mixtures slowly destabilised near the critical point the characteristic length scale of the structure follows the 4/15 power of the quench rate rather than the mean-field 1/6.

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

The paper examines spontaneous liquid-liquid phase separation in binary mixtures that are slowly destabilised close to the critical point using kinetic Monte Carlo and molecular dynamics simulations of a binary surface fluid. It establishes that the characteristic length scale decreases with the quench rate to the 4/15 dynamic critical exponent in two dimensions. This is in contrast to the mean-field theory's 1/6 power. The result shows that the cluster formation dynamics from undriven Brownian motion is more sensitive to the destabilisation rate than mean-field approaches predict. The authors extend the discussion to three-dimensional systems involving liquid crystals and passive or active particles.

Core claim

Using simulations, the authors show that the characteristic length scale of the emerging structure in a two-dimensional binary mixture decreases with the 4/15 dynamic critical exponent of the quench rate rather than the mean-field 1/6th power when the mixture is gradually destabilised near the critical point.

What carries the argument

Dynamic critical exponent of 4/15 for the dependence of structure length scale on quench rate in Brownian motion governed cluster formation.

If this is right

  • The dynamics of cluster formation are much more sensitive on the rate of destabilisation than expected from mean-field theory.
  • Phase separation takes place close to the critical point in slowly destabilised mixtures.
  • The finding has implications for 3D systems with ordering liquid crystals as well as phase-separating passive or active particles.

Where Pith is reading between the lines

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

  • This sensitivity may affect models of phase separation in biological systems where mixtures are slowly destabilised.
  • In three dimensions the exponent might differ and require new simulations to determine.
  • Experiments with controlled quench rates in colloidal systems could test the predicted scaling.

Load-bearing premise

The kinetic Monte Carlo and molecular dynamics simulations accurately capture the cluster formation dynamics governed by thermodynamically undriven Brownian motion under gradual destabilisation near the critical point.

What would settle it

A direct measurement of the scaling exponent between structure length scale and quench rate in a two-dimensional binary fluid system that yields a value close to one sixth instead of four fifteenths would contradict the central claim.

Figures

Figures reproduced from arXiv: 1907.03148 by Charley Schaefer, Stefan Paquay, Tom C. B. McLeish.

Figure 4
Figure 4. Figure 4: This figure shows the time evolution of the characteris [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
read the original abstract

Spontaneous liquid-liquid phase separation is commonly understood in terms of phenomenological mean-field theories. These theories correctly predict the structural features of the fluid at sufficiently long time scales and wavelengths. However, these conditions are not met in various examples in biology and materials science where the mixture is slowly destabilised, and phase separation takes place close to the critical point. Using kinetic Monte Carlo and molecular dynamics simulations of a binary surface fluid under these conditions, we show that the characteristic length scale of the emerging structure decreases, in 2D, with the 4/15 dynamic critical exponent of the quench rate rather than the mean-field 1/6th power. Hence, the dynamics of cluster formation governed by thermodynamically undriven Brownian motion is much more sensitive on the rate of destabilisation than expected from mean-field theory. We discuss the expected implications of this finding to 3D systems with ordering liquid crystals, as well as phase-separating passive or active particles.

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 / 1 minor

Summary. The paper investigates morphology formation in binary mixtures under gradual destabilisation near the critical point using kinetic Monte Carlo and molecular dynamics simulations of a 2D binary surface fluid. It claims that the characteristic length scale of the emerging structure scales with the quench rate according to the dynamic critical exponent 4/15 rather than the mean-field prediction of 1/6, implying greater sensitivity of cluster formation dynamics (governed by undriven Brownian motion) to the destabilisation rate than expected from mean-field theory. Implications for 3D systems involving liquid crystals or passive/active particles are discussed.

Significance. If the central scaling result holds with proper verification, the work would challenge the applicability of mean-field theories to slowly quenched phase separation near criticality, with relevance to biological and materials systems. The use of simulations to identify a non-mean-field dynamic exponent is potentially valuable, but the absence of detailed methods and analysis in the provided information limits assessment of whether this constitutes a robust, falsifiable prediction.

major comments (3)
  1. [Methods] Methods section: The kinetic Monte Carlo and molecular dynamics simulation protocols are not described with sufficient specificity (e.g., no details on system sizes, interaction potentials, time steps, quench rate implementation, or how gradual destabilisation near criticality is realized). This prevents verification that the simulations capture the claimed thermodynamically undriven Brownian motion dynamics, which is load-bearing for the 4/15 exponent claim.
  2. [Results] Results section (and abstract): No error bars, fitting procedures, robustness checks against system size or multiple quench rates, or explicit comparison to the mean-field 1/6 scaling are provided for the extraction of the 4/15 exponent. The central claim that the length scale 'decreases with the 4/15 dynamic critical exponent' cannot be assessed for statistical support or independence from analysis choices.
  3. [Discussion] Discussion section: The extension to 3D systems with ordering liquid crystals or active particles is presented without supporting simulations, analytic arguments, or references to specific dynamic universality classes, making the claimed broader implications speculative rather than substantiated.
minor comments (1)
  1. [Abstract] Abstract: The phrase '4/15 dynamic critical exponent of the quench rate' would benefit from a brief parenthetical reference to the relevant dynamic universality class or prior literature for clarity.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their thoughtful review and constructive suggestions. We address each major comment below and will revise the manuscript to improve clarity and support for the central claims where appropriate.

read point-by-point responses

  1. Referee: [Methods] Methods section: The kinetic Monte Carlo and molecular dynamics simulation protocols are not described with sufficient specificity (e.g., no details on system sizes, interaction potentials, time steps, quench rate implementation, or how gradual destabilisation near criticality is realized). This prevents verification that the simulations capture the claimed thermodynamically undriven Brownian motion dynamics, which is load-bearing for the 4/15 exponent claim.

    Authors: We agree that the methods description requires more detail to enable independent verification. In the revised manuscript we will expand the Methods section with explicit information on system sizes, interaction potentials, time-stepping, the precise implementation of the gradual quench protocol, and confirmation that the dynamics remain in the undriven Brownian regime. revision: yes

  2. Referee: [Results] Results section (and abstract): No error bars, fitting procedures, robustness checks against system size or multiple quench rates, or explicit comparison to the mean-field 1/6 scaling are provided for the extraction of the 4/15 exponent. The central claim that the length scale 'decreases with the 4/15 dynamic critical exponent' cannot be assessed for statistical support or independence from analysis choices.

    Authors: We accept that the presentation of the scaling result needs strengthening. The revised version will include error bars on all data points, a clear description of the fitting procedure used to extract the exponent, finite-size robustness checks, data for additional quench rates, and a direct side-by-side comparison with the mean-field 1/6 prediction. revision: yes

  3. Referee: [Discussion] Discussion section: The extension to 3D systems with ordering liquid crystals or active particles is presented without supporting simulations, analytic arguments, or references to specific dynamic universality classes, making the claimed broader implications speculative rather than substantiated.

    Authors: The 3D implications are offered as an outlook based on the 2D scaling result and the known dynamic universality classes that govern undriven Brownian motion. While no new 3D simulations are performed, we will add explicit references to the relevant universality classes and a short analytic argument linking the 2D exponent to expected 3D behavior. If the referee considers this insufficient, we are prepared to remove or further qualify the paragraph. revision: partial

Circularity Check

0 steps flagged

No significant circularity; result emerges from independent simulations

full rationale

The paper's central claim—that the characteristic length scale scales as the 4/15 dynamic critical exponent rather than the mean-field 1/6—is obtained directly from kinetic Monte Carlo and molecular dynamics simulations of the binary surface fluid under gradual destabilisation. No load-bearing step reduces by construction to a fitted parameter, self-definition, or self-citation chain; the scaling is reported as an output of the numerical experiments. The derivation chain is therefore self-contained and externally falsifiable via the simulation protocols themselves.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The paper relies on standard simulation techniques and mean-field theory as a baseline comparison; no free parameters or invented entities are identifiable from the abstract alone.

axioms (1)
  • domain assumption Mean-field theory predicts a 1/6 power for the length scale dependence on quench rate
    Invoked in the abstract as the point of contrast for the simulation result.

pith-pipeline@v0.9.0 · 5703 in / 1268 out tokens · 34894 ms · 2026-05-25T01:33:08.607919+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

49 extracted references · 49 canonical work pages

  1. [1]

    5 3 −1 −0. 5 0 0 . 5 1 critical cooling 1.00 0.90 θ = 0. 75 ǫAB/k BT ψ 10 100 0.01 0 .1 1 10 increasing quench rate θ= 0. 95 ψ= 0 R∗ ΓTt Figure 4 Top: Phase diagrams for various surface coverages θ. εAB is the nearest-neighbour interaction energy between species A and B, kBT is the thermal energy, and ψ is the order parameter. The symbols are the binodal ...

    work page

  2. [2]

    All of those reported works stud- ied phase separation in response to a concentration rather than a temperature quench

    and ternary [17] mixtures. All of those reported works stud- ied phase separation in response to a concentration rather than a temperature quench. The ¯R values of the DDFT results, as well as the quench rates ¯Γ have been shifted for clarity . From Figure 7 we see that the microscopic modelling results us- ing kMC and MD show a much stronger dependence (...

    work page

  3. [3]

    Brangwynne

    Yongdae Shin and Clifford P. Brangwynne. Liquid phase con- densation in cell physiology and disease. Science, 357:1253, 2017

    work page 2017

  4. [4]

    Brabec, Srinivas Gowrisanker, Jonathan J

    Christoph J. Brabec, Srinivas Gowrisanker, Jonathan J. M. Halls, Darin Laird, Shijun Jia, and Shawn P. Williams. Polymer-fullerene bulk-heterojunction solar cells. Adv. Mater ., 22:3839–3856, 2010

    work page 2010

  5. [5]

    S. R. de Groot and P. Mazur. Non-equilibrium Thermodynam- ics. Dover Publications, Inc., New York, 1984

    work page 1984

  6. [6]

    Hansen and I

    J-P. Hansen and I. R. McDonald. Theory of simple liquids . Academic Press, Hoboken, NJ, 3 edition, 2006

    work page 2006

  7. [7]

    Dynamics of phase separation and critical phenomena in polymer mixtures

    Kurt Binder. Dynamics of phase separation and critical phenomena in polymer mixtures. Colloid & Polymer Sci. , 265:273–288, 1987

    work page 1987

  8. [8]

    Viscoelastic phase separation

    Hajime Tanaka. Viscoelastic phase separation. J. Phys.: Con- dens. Matter, 12:R207–R264, 2000

    work page 2000

  9. [9]

    Furtado and J

    K. Furtado and J. M. Yeomans. Lattice boltzmann simulations of phase separation in chemically reactive binary fluids.Phys. Rev. E, 73:066124, 2006

    work page 2006

  10. [10]

    Molecular dynamics study of phase separation in fluids with chemical reactions

    Raishma Krishnan and Sanjay Puri. Molecular dynamics study of phase separation in fluids with chemical reactions. Phys. Rev. E, 92:052316, 2015

    work page 2015

  11. [11]

    S. C. Glotzer, E. A. Di Marzio, and M. Muthukumar. Reaction- controlled morphology of phase-separating mixtures. Phys. Rev. Lett, 74:2034–2037, 1995

    work page 2034

  12. [12]

    Growth kinetics of nanoclusters in solution

    Awaneesh Singh, Sanjay Puri, and Chandan Dasgupta. Growth kinetics of nanoclusters in solution. J. Phys. Chem. B, 116:4519, 2012

    work page 2012

  13. [13]

    Prasad Perlekar, Roberto Benzi, Herman J. H. Clercx, David R. Nelson, and Federico Toschi. Spinodal decompo- sition in homogeneous and isotropic turbulence. Phys. Rev. Lett., 112:014502, 2014

    work page 2014

  14. [14]

    Benzi, M

    R. Benzi, M. Bernaschi, M. Sbragaglia, and S. Succi. Hetero- geneous diffuse interfaces: A new mechanism for arrested coarsening in binary mixtures. Eur . Phys. J. E, 34:93, 2011

    work page 2011

  15. [15]

    E. L. Huston, J. W. Cahn, and J. E. Hilliard. Spinodal decom- position during continuous cooling. Acta Metall. , 15:1053– 1062, 1966

    work page 1966

  16. [16]

    E. B. Nauman, M. V. Ariyapadi, N. P. Balsara, T. A. Gro- cela, J. S. Furno, S. H. Liu, and R. Mallikarjun. Composi- tional quenching: A process for forming polymer-in-polymer microdispersions and cocontinuous networks. Chem. Eng. Comm., 66:29–55, 1988

    work page 1988

  17. [17]

    Wodo and B

    O. Wodo and B. Ganapathysubramanian. Modeling morphol- ogy evolution during solvent-based fabrication of organic so- lar cells. Comp. Mater . Sci., 55:113–126, 2012

    work page 2012

  18. [18]

    Charley Schaefer, Paul van der Schoot, and Jasper J. Michels. Structuring of polymer solutions upon solvent evaporation. Phys. Rev. E, 91:022602, 2015

    work page 2015

  19. [19]

    Michels, and Paul van der Schoot

    Charley Schaefer, Jasper J. Michels, and Paul van der Schoot. Structuring of thin-film polymer mixtures upon solvent evap- oration. Macromolecules, 49:6858–6870, 2016

    work page 2016

  20. [20]

    Structuring of fluid adlayers upon ongoing unimolecular adsorption

    Charley Schaefer. Structuring of fluid adlayers upon ongoing unimolecular adsorption. Phys. Rev. Lett., 120:036001, 2018

    work page 2018

  21. [21]

    Mod- eling size controlled nanoparticle precipitation with the co- solvency method by spinodal decomposition

    Simon Keßler, Friederike Schmid, and Klaus Drese. Mod- eling size controlled nanoparticle precipitation with the co- solvency method by spinodal decomposition. Soft Matter , 12:7231–7240, 2016

    work page 2016

  22. [22]

    Dynamic density func- tional theories for inhomogeneous polymer systems com- pared to brownian dynamics simulations

    Shuanhu Qi and Friederike Schmid. Dynamic density func- tional theories for inhomogeneous polymer systems com- pared to brownian dynamics simulations. Macromolecules, 50:9831–9845, 2018

    work page 2018

  23. [23]

    Cirillo, Matteo Colangeli, Ellen Moons, Adrian Muntean, Stela-Andrea Muntean, and Jan van Stam

    Emilio N.M. Cirillo, Matteo Colangeli, Ellen Moons, Adrian Muntean, Stela-Andrea Muntean, and Jan van Stam. A lattice model approach to the morphology formation from ternary mixtures during the evaporation of one component. Eur . Phys. J. Special T opics, 228:55–68, 2019

    work page 2019

  24. [24]

    Carmesin, D

    H.-O. Carmesin, D. Heermann, and K. Binder. Influence of a continuous quenching procedure on the initial stages of spinodal decomposition. Z. Phys. B - Condens. Matter , 65:89– 102, 1986

    work page 1986

  25. [25]

    P. C. Hohenberg and B. I. Halperin. Theory of dynamic criti- cal phenomena. Rev. of Mod. Phys. , 49:435–479, 1977

    work page 1977

  26. [26]

    H. E. Cook. Brownian motion in spinodal decomposition. Acta Metall., 17:297–306, 1970

    work page 1970

  27. [27]

    J. S. Langer. Theory of spinodal decomposition in alloys. Ann. Phys., 65:53–86, 1971

    work page 1971

  28. [28]

    R. G. Petschek and H. Metiu. A computer simulation of the time-dependent ginzburg-landau model for spinodal decom- position. J. Chem. Phys. , 79:3443–3456, 1983

    work page 1983

  29. [29]

    James F. Lutsko. A dynamical theory of nucleation for col- loids and macromolecules. J. Chem. Phys. , 136:034509, 2012

    work page 2012

  30. [30]

    S. Buil, J. P. Delville, and A. Ducasse. Eur . Phys. J. E, 2:105, 2000

    work page 2000

  31. [31]

    Kashchiev

    D. Kashchiev. Nucleation: Basic theory with applications . Butterworth-Heinemann, Oxford, 2000

    work page 2000

  32. [32]

    Nucleation theory and dymamics of first-order phase transitions near a critical point

    Kurt Binder. Nucleation theory and dymamics of first-order phase transitions near a critical point. rnal de Physique Col- loques, 1, 41:C4–51–C4–62, 1980

    work page 1980

  33. [33]

    F. J. Alexander, D. A. Huse, and S. A. Janowsky . Dynamics of phase separation and critical phenomena in polymer mix- tures. Phys. Rev. B, 50:663–667, 1994

    work page 1994

  34. [34]

    Non-equilibrium critical dynamics of the fer- roomagnetic ising model with kawasaki dynamics

    Claude Godreche, Florent Krzakala, and Federico Ricci- Tersenghi. Non-equilibrium critical dynamics of the fer- roomagnetic ising model with kawasaki dynamics. J. Stat. Mech.: Theor . Exp., page P04007, 2004

    work page 2004

  35. [35]

    Cahn and John E

    John W. Cahn and John E. Hilliard. Free energy of a nonuni- form system. i. interfacial free energy .J. Chem. Phys., 28:258, 1958. 7

    work page 1958

  36. [36]

    Diffusion constants near the critical point for time-dependent ising models

    Kyozi Kawasaki. Diffusion constants near the critical point for time-dependent ising models. i. Phys. Rev., 145:224–230, 1965

    work page 1965

  37. [37]

    Segregation dynamics of binary mixtures with simple chemical reactions

    Puri and Frisch. Segregation dynamics of binary mixtures with simple chemical reactions. J. Phys. A: Math. Gen. , 27:6027–6038, 1994

    work page 1994

  38. [38]

    Phase separation in binary mixtures with chemical reactions

    Puri and Frisch. Phase separation in binary mixtures with chemical reactions. J. Mod. Phys. B , 12:1623–1641, 1998

    work page 1998

  39. [39]

    Statistical Mechanics

    Kerson Huang. Statistical Mechanics . John Wiley & Sons, Inc., Toronto, 2 edition, 1987

    work page 1987

  40. [40]

    C. G. W. Hermse and A. P. J. Jansen. Kinetics of surface re- actions with lateral interactions: Theory and simulations , vol- ume 19 of SPR - Catalysis . The Royal Society of Chemistry , 2006

    work page 2006

  41. [41]

    A. P. J. Jansen. An introduction to kinetic Monte Carlo simula- tions of surface reactions . Springer-Verlag, Berlin Heidelberg, 1 edition, 2012

    work page 2012

  42. [42]

    J. J. Lukkien, J. P. L. Segers, P. A. J. Hilbers, R. J. Gelten, and A. P. J. Jansen. Efficient monte carlo methods for the simu- lation of catalytic surface reactions. Phys. Rev. E , 58:2598– 2610, 1998

    work page 1998

  43. [43]

    Plimpton

    S. Plimpton. Fast parallel algorithms for short-range molec- ular dynamics. J. Comp. Phys. , 117:1–19, 1995

    work page 1995

  44. [44]

    van Franeker, Daniel Werhoff, Mathieu Turbiez, Martijn M

    Jacobus J. van Franeker, Daniel Werhoff, Mathieu Turbiez, Martijn M. Wienk, Volker Schmidt, and René A. J. Janssen. Controlling the dominant length scale of liquid-liquid phase separation in spin-coated organic semiconductor films. Adv. Func. Mater ., 25:855–863, 2015

    work page 2015

  45. [45]

    J. K. G. Dhont and W. J. Briels. Isotropic-nematic spinodal decomposition kinetics. Phys. Rev. E, 72:031404, 2005

    work page 2005

  46. [46]

    Paul Lettinga, Kyongok Kang, Peter Holmqvist, Arnout Imhof, Didi Derks, and Jan K

    M. Paul Lettinga, Kyongok Kang, Peter Holmqvist, Arnout Imhof, Didi Derks, and Jan K. G. Dhont. Nematic-isotropic spinodal decomposition kinetics of rodlike viruses. Phys. Rev. E, 73:011412, 2006

    work page 2006

  47. [47]

    Green, Robert A

    Micah J. Green, Robert A. Brown, and Robert C. Armstrong. Spinodal decomposition and nematic coarsening in a rigid- rod solution. J. Non-Newtonian Fluid Mech. , 161:30–36, 2009

    work page 2009

  48. [48]

    Dynamic scaling of growing interfaces

    Mehran Kardar and Giorgio Parisi Yi-Checg Zhang. Dynamic scaling of growing interfaces. Phys. Rev. Lett. , 56:889–892, 1986

    work page 1986

  49. [49]

    Allen, Davide Marenduzzo, and Michael E

    Raphael Wittkowski, Adriano Tiribocchi, Joakim Stenham- mar, Rosalind J. Allen, Davide Marenduzzo, and Michael E. Cates. Scalar φ 4 field theory for active-particle phase separa- tion. Nature Comm., 5:4351, 2014. 8

    work page 2014