Pith. sign in

REVIEW 3 major objections 6 minor 29 references

Boundary-Driven Anisotropic Coarsening in Conserved Phase Separation

T0 review · 3 major / 6 minor · reviewed 2026-07-30 · grok-4.5

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Boundary forcing alone can break isotropic scaling in conserved phase separation, producing different coarsening laws parallel and perpendicular to the boundary.

desk verdict Boundary evaporation can produce direction-dependent effective coarsening in an otherwise isotropic conserved ternary system; the claim is real as a global phenomenon, with the main caveat already flagged by the authors themselves. read the letter →

arxiv 2607.26920 v1 pith:OOJL4IVM submitted 2026-07-29 cond-mat.stat-mech math-phmath.MP

Emilio N.M. Cirillo , Nicklas Jävergård , Adrian Muntean , Stela Andrea Muntean This is my paper · ORCID
classification cond-mat.stat-mechmath-phmath.MP
keywords phaseseparationternarymixtureevaporationanisotropiccoarseningModelBBlume–Capelboundaryforcingmorphologyformation
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

In conserved phase separation, domain growth is usually expected to be isotropic, with one universal law (Model B, size growing like t to the one-third) in every direction. This paper argues that forcing only at the boundary is enough to break that symmetry. In a ternary mixture where a passive component evaporates from the surface, mass loss builds macroscopic concentration gradients that make domains coarsen at different effective rates along and across the evaporation direction, while interior bulk regions still look like ordinary isotropic Model B. The claim matters because it separates boundary-driven anisotropy from any change in the intrinsic bulk dynamics, and it suggests that surface conditions alone can be used to steer morphology in nonequilibrium phase-separating materials.

What carries the argument

A ternary Blume–Capel lattice model (two conserved active species plus a passive evaporating component) and its continuum nonlocal drift–diffusion counterpart: Kawasaki bulk exchanges plus stochastic passive loss at the free surface generate macroscopic composition gradients that couple to active-species coarsening and yield direction-dependent two-point correlation lengths.

What would settle it

Restrict correlation analysis to many independent bulk subvolumes far from the evaporating faces across a range of evaporation rates: if those subvolumes still show isotropic t^{1/3} scaling while full-sample x/y versus z lengths keep splitting into the reported anisotropic exponents, the boundary-flux claim holds; if the split disappears once layers are synchronized or mass is held fixed, the global laws are composites rather than true anisotropic coarsening.

Watch

Extended reading notes

Core claim

Boundary forcing by surface evaporation of a passive species is sufficient to break dynamical scaling symmetry in an otherwise isotropic conserved system. Macroscopic concentration gradients produced by progressive mass loss drive anisotropic coarsening, with different effective global growth laws parallel and perpendicular to the boundary, while bulk regions retain standard Model B scaling—so the anisotropy is attributed to boundary fluxes rather than altered intrinsic dynamics.

Load-bearing premise

That full-sample correlation lengths along each axis remain clean direction-dependent coarsening laws even when the passive fraction is macroscopically uneven and different layers sit at different morphological stages, so global exponents are not mainly superposition artifacts.

Editorial extensions

If this is right

  • Surface evaporation or analogous boundary mass fluxes can be used as a control knob to set different coarsening rates along chosen axes without redesigning bulk interactions.
  • Morphology design in drying films and ternary coatings can exploit gradient-driven elongation rather than only composition and interfacial energy.
  • Claims of universal isotropic Model B scaling in conserved systems must be checked against boundary conditions whenever open or evaporating surfaces are present.
  • Layer-resolved measurements become necessary: global exponents alone can mix staggered local crossovers into apparent intermediate powers.

Reading between the lines

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

  • The same boundary-gradient mechanism should appear in other conserved multi-component models (not only Blume–Capel) whenever one species is removed at a face and the others remain locally conserved.
  • If evaporation rate and sample thickness can be tuned so the composition front sweeps uniformly, one might switch continuously between isotropic Model B and strongly anisotropic global growth.
  • Thin-film device processing that already uses solvent evaporation may already be operating in this anisotropic regime without recognizing the scaling split.
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 / 6 minor

Summary. The manuscript argues that boundary forcing alone can break the usual isotropic dynamical scaling of conserved phase separation. Using a three-state Blume–Capel lattice model with Kawasaki dynamics and stochastic evaporation of the passive species at one boundary, together with a continuum nonlocal Model-B-like system with a Robin flux for the passive component, the authors report direction-dependent effective growth of characteristic lengths extracted from two-point correlations: different power laws parallel and perpendicular to the evaporation direction. Layer-wise sections and passive-fraction profiles show macroscopic z-gradients and staggered morphological stages, while bulk sub-volumes away from the boundaries recover isotropic t^{1/3} Model B scaling. The anisotropy is therefore attributed to boundary-induced concentration gradients and fluxes rather than to a change of the intrinsic bulk dynamics. A simple diffusion-length argument is offered for the longitudinal t^{1/2}-like regime at weak evaporation.

Significance. If the anisotropic global scaling is cleanly established and not an averaging artifact, the result is a useful conceptual contribution: it identifies purely external boundary mass loss as a generic route to direction-dependent coarsening in otherwise isotropic conserved systems, without built-in anisotropic interactions or shear. That bulk cubes retain standard Model B while the full sample does not is a clear and falsifiable distinction. The dual lattice/continuum evidence, the layer-resolved passive profiles, and the explicit no-evaporation baseline from prior work strengthen the case. The finding is relevant to evaporative morphology control in polymer blends and thin films. The main scientific value hinges on whether the reported global exponents are genuine direction-dependent coarsening laws or composites of staggered layers; resolving that point determines how far the claim travels.

major comments (3)
  1. [Figs. 2–4 and discussion of effective 1/4 exponents] The central claim of distinct effective global growth laws parallel vs. perpendicular to the boundary is undercut by the paper’s own layer-wise analysis (Figs. 3–4 and surrounding text). Passive fraction is strongly inhomogeneous in z; successive planes sit at different morphological stages and cross over to binary-like Model B growth at different times. The authors explicitly state that this staggered progression “can produce an effective transverse growth exponent smaller than the local Model B value,” which is how the intermediate ~1/4 exponents for α=10^{-2} and 10^{-1} are explained. Full-sample correlation lengths along x/y versus z are therefore at risk of being superposition measures of an inhomogeneous sample rather than directionally distinct dynamical processes. Bulk cubes recovering isotropic t^{1/3} reinforce this reading. The manuscript needs either (i) a quantitative decom
  2. [Fig. 2 and Fig. 7] Power-law identification rests on visual guide lines and selected collapse windows (Fig. 2 top/bottom; Fig. 7) rather than documented fitting ranges, uncertainties, or sensitivity to the definition of the characteristic length from G. Prefactors are free (0.12 t^{1/4}, 0.45 t^{1/3}, 0.12 t^{1/2}, etc.), and the continuum exponents shift with α' in a way that is only loosely tied to the lattice cases. For the weak-evaporation longitudinal t^{1/2} claim and the fast-evaporation transverse t^{1/4} claim to support the abstract’s “different effective global growth laws,” the paper should report how the exponents are extracted (fit windows, alternative length definitions such as first zero or half-height of G, run-to-run variation) and whether the collapses remain stable under those choices.
  3. [Eqs. (2)–(3), Figs. 6–7] The continuum counterpart (Eqs. 2–3, Figs. 6–7) is used to argue that the anisotropy is not a lattice artifact, but the analysis is thinner: fewer diagnostics, no layer-wise passive profiles or bulk-cube checks analogous to the lattice, and evaporation imposed at both top and bottom. A parallel bulk-versus-global comparison in the continuum model would substantially strengthen the robustness claim; without it, the continuum evidence mainly shows morphological elongation, not the same scaling dichotomy.
minor comments (6)
  1. Notation for the initial passive fraction switches between c_0 (lattice) and 1−ϕ_0 (continuum); a single convention or an explicit dictionary would help.
  2. [Methods paragraph and Eq. (3)] Evaporation is said to act only at the “top” boundary in the lattice description, yet Fig. 4 and the continuum BC (3) indicate loss at both ends; clarify the lattice implementation (one face vs two).
  3. Inverse temperature is given as 0.80 for the lattice and β=1 for continuum without relating the two scales; a brief remark on corresponding reduced temperatures would aid comparison.
  4. [Fig. 1] Fig. 1 caption plane labels “(x,0,z)”, “(L−1,y,z)”, “(x,y,L−1)” are slightly ambiguous; stating which face is the evaporating boundary in each panel would improve readability.
  5. Typos/style: “J¨ averg ˚ ard”, “att=” spacing, “conponent”, “scillating”, and repeated “in this case” in the abstract; standard copy-edit pass needed.
  6. [Scaling argument paragraph] The simple ℓ∼(Dt)^{1/2} argument is plausible for weak evaporation but is not checked against the measured passive-gradient width versus time; a single panel of gradient width vs t would make the scaling argument quantitative.

Circularity Check

1 steps flagged · score 1.0 of 10

No meaningful circularity: anisotropic exponents are measured outputs of forward simulations, not quantities forced by definition or self-citation.

  1. self citation load bearing [Main text, paragraphs on no-evaporation baseline (refs [21,22]) and continuum well-posedness (ref [29])]
    "In the absence of evaporation, we previously reported morphology formation in both two- and three-dimensional systems [21, 22]. In three dimensions, the characteristic domain size follows a t^{1/3} growth law. Coarsening remains isotropic and is consistent with conserved Model B dynamics. ... The inequality |m| ≤ ϕ ≤ 1 is preserved during the evolution [29]."

    The isotropic Model B baseline and continuum well-posedness are justified by overlapping-author citations rather than re-derived here. This is ordinary background self-citation and is not load-bearing for the new anisotropic exponents under evaporation, which come from fresh simulations; it does not force the claimed direction-dependent growth laws.

full rationale

The central claim—that surface evaporation of a passive species induces direction-dependent effective coarsening exponents while bulk regions retain Model B t^{1/3}—is obtained by running conserved Kawasaki (lattice) and nonlocal Model-B-type (continuum) dynamics with an explicit evaporative boundary rule, then extracting characteristic lengths from two-point correlations along x/y/z. Those lengths and the fitted power laws are simulation outputs, not inputs renamed as predictions. Prior self-citations ([21,22] no-evaporation morphology and isotropic t^{1/3}; [29] continuum well-posedness) supply only the equilibrium/baseline comparison and mathematical background; they do not define or force the anisotropic exponents under evaporation. The paper’s own caveat that fast-evaporation transverse ~1/4 exponents may be composites of staggered layer crossovers is an interpretive limitation on what the global lengths mean, not a circular reduction of the derivation. No self-definitional loop, fitted-input-as-prediction, uniqueness import, or ansatz smuggling is present. Score 1 reflects only routine non-load-bearing self-citation of the authors’ prior baseline work.

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

The paper rests on standard conserved phase-separation modeling (Blume–Capel + Kawasaki; nonlocal Model B-like hydrodynamics) plus an imposed evaporative boundary that removes only the passive species. No new physical entities are postulated. Load-bearing choices are simulation parameters (α, c0, β, system size) and the operational definition of directional domain size from two-point correlations on an inhomogeneous sample. The t^{1/3} bulk benchmark and diffusive t^{1/2} gradient-relaxation sketch are taken from established scaling lore rather than re-derived.

free parameters (5)
  • evaporation probability α (lattice) = 10^{-3}, 10^{-2}, 10^{-1}
    Hand-chosen rates 10^{-3}, 10^{-2}, 10^{-1} that select which effective exponent pair is observed; central phenomenology depends on this control parameter.
  • continuum evaporation coefficient α' = 1, 10, 100
    Boundary flux strength in the continuum model; scanned as 1, 10, 100 (and 100 in morphology snapshots) without first-principles fixation.
  • initial passive fraction c0 / 1-ϕ0 = 0.1, 0.4, 0.8
    Composition control (0.1, 0.4, 0.8) that gates how visible anisotropy is; not predicted, imposed.
  • inverse temperature β (and lattice 0.80) = lattice 0.80; continuum β=1
    Sets quench depth/interfacial structure; fixed by authors (lattice inverse temperature 0.80; continuum β=1 in shown runs).
  • power-law prefactors on domain-size guides = e.g. 0.12, 0.45, 0.8, 1.2 (figure guides)
    Solid-line references such as 0.12 t^{1/4}, 0.45 t^{1/3}, 0.12 t^{1/2} are amplitude choices guiding the eye; exponents are the scientific claim but amplitudes are fit/reference parameters.
assumptions (6)
  • domain assumption Conserved Model B bulk coarsening has characteristic domain growth ~ t^{1/3} in the late stage for diffusive dynamics.
    Used as the isotropic baseline and bulk benchmark throughout (abstract; no-evaporation paragraph; Fig. 3 references).
  • domain assumption Blume–Capel Hamiltonian with Kawasaki exchanges correctly represents two conserved active species plus a passive screener on the lattice.
    Eq. (1) and methods paragraphs define the microscopic model inherited from authors’ prior ternary studies.
  • domain assumption Hydrodynamic/Kac limit yields the stated coupled nonlocal equations for m and ϕ with |m|≤ϕ≤1 preserved.
    Eqs. (2)–(3) and citations to Marra–Mourragui, Giacomin–Lebowitz, Lyons et al.; continuum evidence stream depends on this limit.
  • ad hoc to paper Evaporation may be modeled as stochastic conversion of passive sites only at the top boundary (lattice) or as the Robin-like flux (3) (continuum), without direct creation/removal of active species.
    Core modeling choice that isolates boundary forcing; stated explicitly in the evaporation implementation paragraphs.
  • ad hoc to paper Directional characteristic lengths from normalized two-point correlations G along x,y,z are valid probes of anisotropic coarsening even under macroscopic z-inhomogeneity.
    Operational basis of Figs. 2, 5, 7 and all exponent claims.
  • domain assumption Relaxation of the boundary-induced passive gradient is diffusion-dominated, motivating longitudinal ℓ ~ (Dt)^{1/2} in the weak-evaporation regime.
    Simple scaling argument near the end; explains one observed exponent pair but is not derived from the full free-boundary problem.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Boundary-Driven Anisotropic Coarsening in Conserved Phase Separation." pith.science (2026). https://pith.science/paper/OOJL4IVM

@misc{pith2026260726920,
  author       = {Pith},
  title        = {Pith review of: Boundary-Driven Anisotropic Coarsening in Conserved Phase Separation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/OOJL4IVM}},
  note         = {Machine review of arXiv:2607.26920}
}
read the original abstract

Universal scaling in phase separation is typically assumed to be isotropic in systems with conserved dynamics. Here we show that boundary forcing alone can break this dynamical scaling symmetry, leading to anisotropic coarsening, with different effective global growth laws observed parallel and perpendicular to the boundary. We consider a ternary mixture with two conserved components and a passive species undergoing surface evaporation, which provides a simple setting to investigate this effect. In this case, evaporation leads to a progressive mass loss and to the formation of macroscopic concentration gradients, which, in turn, drive anisotropic coarsening, with different effective growth laws observed parallel and perpendicular to the boundary. At the same time, bulk regions appear to retain the standard Model B scaling, suggesting that the observed anisotropy is mainly induced by boundary fluxes rather than by changes in the intrinsic dynamics. Our results indicate that boundary conditions can play an important role in breaking scaling symmetry and may offer a way to influence coarsening behavior in nonequilibrium phase separation.

Figures

Figures reproduced from arXiv: 2607.26920 by the authors.

Figure 2
Figure 2. FIG. 2. Top row: characteristic domain size as a function [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 4
Figure 4. FIG. 4. Fraction of passive conponent as function of the [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Configuration at time 10000 on the planes [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figures from the paper (2 more)
Figure 6
Figure 6. Figure 6: FIG. 6. Typical configurations of the continuum model in [PITH_FULL_IMAGE:figures/full_fig_p005_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7. Characteristic domain size as a function of time for [PITH_FULL_IMAGE:figures/full_fig_p005_7.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

29 extracted references

  1. [1]

    Mouritsen, and Søren Tox- vaerd

    Mohamed Laradji, Ole G. Mouritsen, and Søren Tox- vaerd. Spinodal decomposition in multicomponent fluid mixtures: A molecular dynamics study.Physical Review E, 53(4):3673–3681, 1996

    1996

  2. [2]

    Hoppe and N

    H. Hoppe and N. S. Sariciftci. Organic solar cells: An overview.Journal of Materials Research, 19(7):1924– 1945, 2004

    1924

  3. [3]

    E. Moons. Conjugated polymer blends: linking film morphology to performance of lightemitting diodes and photodiodes.Journal of Physics: Condensed Matter, 14(47):12235–12260, 2002

    2002

  4. [4]

    Nilsson, A

    S. Nilsson, A. Bernasik, A. Budkowski, and E. Moons. Morphology and phase segregation of spin-casted films of polyfluorene/pcbm blends.Macromolecules, 40(23):8291–8301, 2007

    2007

  5. [5]

    J. W. Cahn and J. E. Hilliard. Free energy of a nonuni- form system. i. interfacial free energy.Journal of Chem- ical Physics, 28:258–267, 1958

    1958

  6. [6]

    A. J. Bray. Theory of phase-ordering kinetics.Advances in Physics, 43(3):357–459, 1994

    1994

  7. [7]

    P. C. Hohenberg and B. I. Halperin. Theory of dynamic critical phenomena.Reviews of Modern Physics, 49:435– 479, 1977

    1977

  8. [8]

    A. J. Bray. Universal scaling function for domain growth in the Glauber-Ising chain.Journal of Physics A: Math- ematical and General, 23(2):L67, 1990

    1990

Show all 29 references
  1. [9]

    Chauhan, D

    A. Chauhan, D. Gogoi, S. Puri, and A. Singh. Ef- fect of amphiphilic polymers on phase separating binary mixtures: A dpd simulation study.J. Chem. Phys., 159:204901, 2023

    2023

  2. [10]

    E. N. M. Cirillo, G. Gonnella, and S. Stramaglia. Anisotropic dynamical scaling in a spin model with com- peting interactions.Physical Review E, 56:5065–5074, 1997

    1997

  3. [11]

    Corberi, G

    F. Corberi, G. Gonnella, and A. Lamura. Two-scale com- petition in phase separation with shear.Phys. Rev. Lett., 83:4057, 1999

    1999

  4. [12]

    Parameshwaran, and Bhaskar Sen Gupta

    Daniya Davis, A. Parameshwaran, and Bhaskar Sen Gupta. Phase separation and rheology of segre- gating binary fluid under shear.Physical Review E, 111(6):065414, June 2025

    2025

  5. [13]

    Onuki.Phase Transition Dynamics

    A. Onuki.Phase Transition Dynamics. Cambridge Uni- versity Press, 2002

    2002

  6. [14]

    R. D. Deegan et al. Capillary flow as the cause of ring stains from dried liquid drops.Nature, 389:827–829, 1997

    1997

  7. [15]

    R. G. Larson. Transport and deposition patterns in dry- ing sessile droplets.AIChE Journal, 60:1538–1571, 2014

    2014

  8. [16]

    Hopkinson and M

    I. Hopkinson and M. Myatt. Phase separation in ternary polymer solutions induced by solvent loss.Macro- molecules, 35(13), 2002

    2002

  9. [17]

    Bishnu, R

    R. Bishnu, R. Mukherjee, N. Bhandaru, and A. Dutta. Morphology prediction for polymer blend thin films using machine learning.Soft Matter, 21:5824, 2025

    2025

  10. [18]

    U. Thiele. Thin film evolution equations from evaporat- ing dewetting liquid layers to epitaxial growth.Journal of Physics: Condensed Matter, 22:084019, 2010

    2010

  11. [19]

    Puri and K

    S. Puri and K. Binder. Surface-directed spinodal decom- position: Phenomenology and numerical results.Physical Review A, 46(8):R4487–R4489, 1992

    1992

  12. [20]

    Kinetics and morphology of surface-directed spinodal decomposition.Physica A: Sta- tistical Mechanics and its Applications, 695:131612, 2026

    Rounak Bhattacharyya. Kinetics and morphology of surface-directed spinodal decomposition.Physica A: Sta- tistical Mechanics and its Applications, 695:131612, 2026

    2026

  13. [21]

    Cirillo, N

    E.N.M. Cirillo, N. J¨ averg ˚ ard, R. Lyons, A. Muntean, and S. A. Muntean. 3D morphology formation in a mixture of three differently averse components.Modelling Simul. Mater. Sci. Eng., 33:055014, 2025. 7

    2025

  14. [22]

    S. A. Muntean, V. C. E. Kronberg, M. Colangeli, A. Muntean, J. van Stam, E. Moons, and E. N. M. Cir- illo. Quantitative analysis of phase formation and growth in ternary mixtures upon evaporation of one component. Physical Review E, 106:025306, 2022

    2022

  15. [23]

    M. Blume. Theory of the First–Order Magnetic Phase Change in UO 2.Physical Review, 141:517, 1966

    1966

  16. [24]

    Cirillo, V

    E.N.M. Cirillo, V. Jacquier, and C. Spitoni. Homoge- neous and heterogeneous nucleation in the three-state Blume–Capel model.Physica D, 461:134125, 2024

    2024

  17. [25]

    Mendes, G

    C. Mendes, G. M. Buendia, and P. A. Rikvold. Numerical simulation of a two–dimensional Blume–Capel ferromag- net in an scillating magnetic field with a constant bias. Physical Review E, 110:044133, 2024

    2024

  18. [26]

    Marra and M

    R. Marra and M. Mourragui. Phase segregation dynam- ics for the Blume–Capel model with Kac interaction. Stochastic Processes and their Applications, 88(1):79– 124, 2000

    2000

  19. [27]

    Lebowitz and O

    J.L. Lebowitz and O. Penrose. Rigorous treatment of the van der Walls-Maxwell theory of the liquid-vapor transi- tion.Journal of Mathematical Physics, 7:98–113, 1966

    1966

  20. [28]

    Giacomin and J

    G. Giacomin and J. L. Lebowitz. Phase segregation dy- namics in particle systems with long range interactions. Journal of Statistical Physics, 87:37–61, 1997

    1997

  21. [29]

    Lyons, E

    R. Lyons, E. N. M. Cirillo, and A. Muntean. Phase sep- aration and morphology formation in interacting ternary mixtures under evaporation: Well-posedness and numer- ical simulation of a non-local evolution system.Nonlinear Analysis: Real World Applications, 77:104039, 2024

    2024

Pith tools

Reviewed July 30, 2026 · model on record in the stance chip above.