pith. sign in

arxiv: 2607.01009 · v1 · pith:7CO4FGY5new · submitted 2026-07-01 · ❄️ cond-mat.stat-mech · cond-mat.soft

Phase diagram of a double-occupancy cell model of a fluid with Curie-Weiss interaction

R. V. Romanik , O. A. Dobush , M. P. Kozlovskii , I. V. Pylyuk , M. A. Shpot This is my paper

Pith reviewed 2026-07-02 04:36 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech cond-mat.soft
keywords double-occupancy cell modelCurie-Weiss interactionphase diagramBlume-Capel modelcritical pointstricritical pointtriple pointfluid coexistence
0
0 comments X

The pith

A minimal double-occupancy fluid model with local repulsion and global attraction produces up to three coexisting phases and multiple critical points.

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

The paper maps a simple cell model allowing at most two particles per site, with short-range repulsion and infinite-range attraction, onto the Blume-Capel model on a complete graph. Analysis in the grand-canonical ensemble then shows that the ratio of repulsion to attraction strength controls the appearance of one or two critical points, a tricritical point, or a triple point. When repulsion is strong enough the system supports three distinct fluid densities, producing both ordinary gas-liquid and liquid-liquid coexistence. A sympathetic reader cares because the construction demonstrates that these complex features require only the competition between occupancy restriction and mean-field attraction, without additional interaction terms.

Core claim

The model is isomorphic to the Blume-Capel model on a complete graph. Depending on the ratio of local repulsive to global attractive interactions, its phase diagram contains regimes with a single critical point, two distinct critical points, tricritical behavior, or a triple point. Sufficiently strong repulsion yields three fluid phases of different densities and therefore both gas-liquid and liquid-liquid coexistence; the locations of all special points are located and the corresponding diagrams are constructed.

What carries the argument

Double-occupancy cell model with Curie-Weiss interaction, mapped by spin-to-occupancy transformation to the Blume-Capel model on a complete graph; the mapping converts the grand-canonical partition function into a solvable mean-field problem whose free energy is minimized with respect to density.

If this is right

  • The phase diagram topology changes qualitatively with the single control parameter that sets the ratio of local repulsion to long-range attraction.
  • Three distinct fluid phases appear once repulsion exceeds a threshold value, producing both gas-liquid and liquid-liquid first-order lines that meet at a triple point.
  • Critical and tricritical points can be located by solving the stationarity conditions of the mapped free-energy functional.
  • The same minimal ingredients suffice to generate both ordinary and liquid-liquid transitions without requiring explicit three-body forces.

Where Pith is reading between the lines

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

  • Because the complete-graph limit suppresses spatial fluctuations, finite-dimensional versions of the model may shift or eliminate some of the higher-order critical points.
  • The occupancy restriction used here is formally identical to a hard-core constraint on nearest-neighbor pairs, suggesting the phase structure could appear in lattice models of associating fluids.
  • The analytic location of the triple point supplies a concrete benchmark that could be checked by direct enumeration on small complete graphs.

Load-bearing premise

The model on the complete graph can be treated exactly by the mean-field free-energy minimization that follows from the Blume-Capel isomorphism.

What would settle it

Monte Carlo simulation of the same occupancy rules on a three-dimensional lattice that finds no liquid-liquid coexistence line for parameter values where the complete-graph calculation predicts a triple point.

read the original abstract

A double-occupancy cell model of a fluid with Curie-Weiss interaction is studied. First, we show that the model is isomorphic to the Blume-Capel model on a complete graph through a simple transformation from spin to occupancy variables. We then investigate its phase behavior within the grand-canonical ensemble using a combination of analytical and numerical methods. Despite its simplicity, the model exhibits a remarkably rich thermodynamic behavior depending on the ratio between the local repulsive and global attractive interactions. We identify regimes characterized by a single critical point, two distinct critical points, tricritical behavior, and triple-point formation. For sufficiently strong repulsion, the system possesses three fluid phases of different densities, leading to both gas-liquid and liquid-liquid coexistence. The locations of the critical, tricritical, and triple points are determined, and the corresponding phase diagrams are constructed. These results demonstrate that the competition between double-occupancy repulsion and long-range attraction is sufficient to generate complex phase behavior in a minimal multiple-occupancy lattice-gas model.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper introduces a double-occupancy cell model of a fluid with Curie-Weiss (infinite-range) interactions. It first establishes an isomorphism to the Blume-Capel model on the complete graph via a transformation between occupancy variables {0,1,2} and spin variables {-1,0,1}. Using analytic and numerical methods in the grand-canonical ensemble, the authors map out the phase diagram as a function of the ratio between local repulsion and global attraction, identifying regimes with one critical point, two critical points, tricritical points, triple points, and three coexisting fluid phases of different densities when repulsion is strong.

Significance. If the central mapping is exact and the subsequent calculations are correct, the work shows that a minimal lattice-gas model with only local double-occupancy repulsion and long-range attraction is sufficient to produce multiple fluid phases and multicritical behavior. This is of interest for understanding liquid-liquid coexistence in simple fluids and provides a solvable benchmark for more complex models.

major comments (2)
  1. [mapping to Blume-Capel model] The isomorphism to the Blume-Capel model (introduced immediately after the model definition) is load-bearing for every subsequent result. The transformation must be shown to map the grand-canonical weights exactly, including the fugacity factor e^{βμ n_i} for n_i = 0,1,2, without omitted Jacobians or residual constraints. Any mismatch would shift or eliminate the reported locations of the two critical points, tricritical point, and triple point.
  2. [phase behavior analysis] The phase-diagram construction (analytic mean-field plus numerics) relies on the mapped Blume-Capel Hamiltonian being free of extra terms. The paper should explicitly verify that the chemical-potential term becomes the correct linear field in the spin representation and that the partition function is identical (up to an overall factor independent of the order parameters).
minor comments (2)
  1. Notation for the ratio of repulsive to attractive couplings should be defined once and used consistently in all figures and equations.
  2. Figure captions for the phase diagrams should state the ensemble (grand-canonical) and the numerical method used to locate the coexistence lines.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and for highlighting the importance of an explicit verification of the mapping. We agree that the isomorphism is central and will revise the manuscript to provide a fully detailed derivation of the grand-canonical equivalence, including the fugacity mapping and partition-function identity.

read point-by-point responses

  1. Referee: [mapping to Blume-Capel model] The isomorphism to the Blume-Capel model (introduced immediately after the model definition) is load-bearing for every subsequent result. The transformation must be shown to map the grand-canonical weights exactly, including the fugacity factor e^{βμ n_i} for n_i = 0,1,2, without omitted Jacobians or residual constraints. Any mismatch would shift or eliminate the reported locations of the two critical points, tricritical point, and triple point.

    Authors: We agree that the mapping requires an explicit, step-by-step demonstration. The manuscript introduces the bijection n_i = s_i + 1 (s_i ∈ {-1,0,1} ↔ n_i ∈ {0,1,2}) and states the resulting Blume-Capel equivalence. In the revision we will add a dedicated subsection (or appendix) that writes the grand partition function in the occupancy representation, substitutes the transformation, and shows that the term ∏_i e^{βμ n_i} becomes exactly ∏_i e^{β h s_i} times an overall constant factor independent of the order parameters. Because the map is a discrete bijection on a finite state space per site, no Jacobian appears. The interaction term likewise maps without remainder. This will confirm that the locations of the critical, tricritical, and triple points are unaffected. revision: yes

  2. Referee: [phase behavior analysis] The phase-diagram construction (analytic mean-field plus numerics) relies on the mapped Blume-Capel Hamiltonian being free of extra terms. The paper should explicitly verify that the chemical-potential term becomes the correct linear field in the spin representation and that the partition function is identical (up to an overall factor independent of the order parameters).

    Authors: This point is addressed by the same explicit derivation described above. We will insert the verification that the chemical-potential contribution maps precisely to the linear field h s_i (plus a constant) and that the full grand partition function in the two representations differs only by a multiplicative factor independent of the magnetization and quadrupole moment. With this addition the analytic mean-field equations and the numerical results remain rigorously justified on the mapped model. revision: yes

Circularity Check

0 steps flagged

No circularity; isomorphism presented as independently shown before phase analysis

full rationale

The paper states it first shows the isomorphism to the Blume-Capel model via a transformation from occupancy to spin variables, then proceeds to analyze phase behavior with analytical and numerical methods on the mapped model. No step reduces a claimed prediction or central result to a fitted parameter, self-citation chain, or definitional renaming. The mapping is asserted as derivable rather than smuggled via prior self-work, and the reported multicritical points follow from subsequent calculations on the transformed Hamiltonian. The derivation chain is therefore self-contained against external benchmarks with no load-bearing reduction to its own inputs.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of the isomorphism to the Blume-Capel model and the assumption that varying the single interaction ratio is sufficient to classify all thermodynamic regimes.

free parameters (1)
  • ratio of local repulsive to global attractive interactions
    Phase behavior is classified by varying this ratio; no specific fitted values are given in the abstract.
axioms (1)
  • domain assumption The double-occupancy cell model is isomorphic to the Blume-Capel model on a complete graph via a simple spin-to-occupancy transformation
    This mapping is stated as the starting point that enables the phase-diagram analysis.

pith-pipeline@v0.9.1-grok · 5737 in / 1390 out tokens · 32802 ms · 2026-07-02T04:36:55.029255+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

46 extracted references · 39 canonical work pages · 1 internal anchor

  1. [1]

    Likos, C.N.: Effective interactions in soft condensed matter physics. Phys. Rep.348(4), 267–439 (2001) https://doi.org/10.1016/S0370-1573(00)00141-1

    work page doi:10.1016/s0370-1573(00)00141-1 2001

  2. [2]

    soft colloids

    Louis, A.A., Bolhuis, P.G., Hansen, J.P., Meijer, E.J.: Can polymer coils be modeled as “soft colloids”? Phys. Rev. Lett.85, 2522–2525 (2000) https://doi.org/10.1103/PhysRevLett.85.2522

    work page doi:10.1103/physrevlett.85.2522 2000

  3. [3]

    Bolhuis, P.G., Louis, A.A., Hansen, J.P., Meijer, E.J.: Accurate effective pair potentials for polymer solutions. J. Chem. Phys.114(9), 4296–4311 (2001) https://doi.org/10.1063/1.1344606

    work page doi:10.1063/1.1344606 2001

  4. [4]

    Macromolecules34, 2914–2920 (2001) https://doi

    Likos, C.N., Schmidt, M., L¨ owen, H., Ballauff, M., P¨ otschke, D., Lindner, P.: Soft interaction between dissolved flexible dendrimers: Theory and experiment. Macromolecules34, 2914–2920 (2001) https://doi. org/10.1021/ma001346x

    work page doi:10.1021/ma001346x 2001

  5. [5]

    Mladek, B.M., Kahl, G., Likos, C.N.: Computer assembly of cluster-forming amphiphilic dendrimers. Phys. Rev. Lett.100, 028301 (2008) https://doi.org/10.1103/PhysRevLett.100.028301

    work page doi:10.1103/physrevlett.100.028301 2008

  6. [6]

    Likos, C.N., L¨ owen, H., Watzlawek, M., Abbas, B., Jucknischke, O., Allgaier, J., Richter, D.: Star polymers viewed as ultrasoft colloidal particles. Phys. Rev. Lett.80, 4450–4453 (1998) https://doi.org/10.1103/ PhysRevLett.80.4450

    1998

  7. [7]

    Soft Matter2, 478–498 (2006) https://doi.org/10.1039/ B601916C 14

    Likos, C.N.: Soft matter with soft particles. Soft Matter2, 478–498 (2006) https://doi.org/10.1039/ B601916C 14

    2006

  8. [8]

    Likos, C.N., Watzlawek, M., L¨ owen, H.: Freezing and clustering transitions for penetrable spheres. Phys. Rev. E58, 3135–3144 (1998) https://doi.org/10.1103/PhysRevE.58.3135

    work page doi:10.1103/physreve.58.3135 1998

  9. [9]

    Lang, A., Likos, C.N., Watzlawek, M., L¨ owen, H.: Fluid and solid phases of the Gaussian core model. J. Phys.: Condens. Matter12(24), 5087 (2000) https://doi.org/10.1088/0953-8984/12/24/302

    work page doi:10.1088/0953-8984/12/24/302 2000

  10. [10]

    Mladek, B.M., Gottwald, D., Kahl, G., Neumann, M., Likos, C.N.: Formation of polymorphic cluster phases for a class of models of purely repulsive soft spheres. Phys. Rev. Lett.96, 045701 (2006) https: //doi.org/10.1103/PhysRevLett.96.045701

    work page doi:10.1103/physrevlett.96.045701 2006

  11. [11]

    Mladek, B.M., Charbonneau, P., Frenkel, D.: Phase coexistence of cluster crystals: Beyond the Gibbs phase rule. Phys. Rev. Lett.99, 235702 (2007) https://doi.org/10.1103/PhysRevLett.99.235702

    work page doi:10.1103/physrevlett.99.235702 2007

  12. [12]

    Lenz, D.A., Mladek, B.M., Likos, C.N., Kahl, G., Blaak, R.: Monomer-resolved simulations of cluster- forming dendrimers. J. Phys. Chem. B115(22), 7218–7226 (2010) https://doi.org/10.1021/jp109132m

    work page doi:10.1021/jp109132m 2010

  13. [13]

    Neuhaus, T., Likos, C.N.: Phonon dispersions of cluster crystals. J. Phys.: Condens. Matter23(23), 234112 (2011) https://doi.org/10.1088/0953-8984/23/23/234112

    work page doi:10.1088/0953-8984/23/23/234112 2011

  14. [14]

    Europhys

    Wilding, N.B., Sollich, P.: A Monte Carlo method for chemical potential determination in single and multiple occupancy crystals. Europhys. Lett.101(1), 10004 (2013) https://doi.org/10.1209/0295-5075/ 101/10004

    work page doi:10.1209/0295-5075/ 2013

  15. [15]

    Wilding, N.B., Sollich, P.: Demixing cascades in cluster crystals. J. Chem. Phys.141(9), 094903 (2014) https://doi.org/10.1063/1.4894374

    work page doi:10.1063/1.4894374 2014

  16. [16]

    Prestipino, S.: Cluster phases of penetrable rods on a line. Phys. Rev. E90, 042306 (2014) https://doi. org/10.1103/PhysRevE.90.042306

    work page doi:10.1103/physreve.90.042306 2014

  17. [17]

    Prestipino, S., Gazzillo, D., Tasinato, N.: Probing the existence of phase transitions in one-dimensional fluids of penetrable particles. Phys. Rev. E92, 022138 (2015) https://doi.org/10.1103/PhysRevE.92.022138

    work page doi:10.1103/physreve.92.022138 2015

  18. [18]

    Entropy 25(2) (2023) https://doi.org/10.3390/e25020356

    de Mello, M., D´ ıaz-M´ endez, R., Mendoza-Coto, A.: Ultrasoft classical systems at zero temperature. Entropy 25(2) (2023) https://doi.org/10.3390/e25020356

    work page doi:10.3390/e25020356 2023

  19. [19]

    Baxter, R.J.: Exactly Solved Models in Statistical Mechanics, p. 486. Academic Press, London (1982)

    1982

  20. [20]

    A Concrete Mathematical Introduction

    Friedli, S., Velenik, Y.: Statistical Mechanics of Lattice Systems. A Concrete Mathematical Introduction. Cambridge University Press, Cambridge (2017). https://doi.org/10.1017/9781316882603

    work page doi:10.1017/9781316882603 2017

  21. [21]

    Finken, R., Hansen, J.-P., Louis, A.A.: Phase separation of a multiple occupancy lattice gas. J. Phys. A: Math. Gen.37(3), 577–590 (2004) https://doi.org/10.1088/0305-4470/37/3/003

    work page doi:10.1088/0305-4470/37/3/003 2004

  22. [22]

    Frydel, D., Levin, Y.: Soft-particle lattice gas in one dimension: One- and two-component cases. Phys. Rev. E98, 062123 (2018) https://doi.org/10.1103/PhysRevE.98.062123

    work page doi:10.1103/physreve.98.062123 2018

  23. [23]

    Frydel, D., Levin, Y.: Thermodynamic collapse in a lattice-gas model for a two-component system of penetrable particles. Phys. Rev. E102, 032101 (2020) https://doi.org/10.1103/PhysRevE.102.032101

    work page doi:10.1103/physreve.102.032101 2020

  24. [24]

    Kozitsky, Yu.V., Kozlovskii, M.P., Dobush, O.A.: A phase transition in a Curie-Weiss system with binary interactions. Condens. Matter Phys.23(2), 23502 (2020) https://doi.org/10.5488/CMP.23.23502

    work page doi:10.5488/cmp.23.23502 2020

  25. [25]

    Kozlovskii, M.P., Dobush, O.A.: Phase behavior of a cell model with Curie-Weiss interaction. J. Mol. Liq. 353, 118843 (2022) https://doi.org/10.1016/j.molliq.2022.118843

    work page doi:10.1016/j.molliq.2022.118843 2022

  26. [26]

    Romanik, R.V., Dobush, O.A., Kozlovskii, M.P., Pylyuk, I.V., Shpot, M.A.: A multiple occupancy cell fluid model with competing attraction and repulsion interactions. J. Stat. Phys.193(14) (2026) https: 15 //doi.org/10.1007/s10955-026-03568-4

    work page doi:10.1007/s10955-026-03568-4 2026

  27. [27]

    Liu, W., Yan, Z., Zhou, G.: Phase transitions of a double occupancy lattice gas. Eur. Phys. J. B94(6), 133 (2021) https://doi.org/10.1140/epjb/s10051-021-00142-7

    work page doi:10.1140/epjb/s10051-021-00142-7 2021

  28. [28]

    Blume, M.: Theory of the first-order magnetic phase change in UO 2. Phys. Rev.141, 517–524 (1966) https://doi.org/10.1103/PhysRev.141.517

    work page doi:10.1103/physrev.141.517 1966

  29. [29]

    Physica32(5), 966–988 (1966) https://doi.org/10.1016/0031-8914(66)90027-9

    Capel, H.W.: On the possibility of first-order phase transitions in Ising systems of triplet ions with zero-field splitting. Physica32(5), 966–988 (1966) https://doi.org/10.1016/0031-8914(66)90027-9

    work page doi:10.1016/0031-8914(66)90027-9 1966

  30. [30]

    Lee, T.D., Yang, C.N.: Statistical theory of equations of state and phase transitions. II. Lattice gas and Ising model. Phys. Rev.87, 410–419 (1952) https://doi.org/10.1103/PhysRev.87.410

    work page doi:10.1103/physrev.87.410 1952

  31. [31]

    Low Temp

    Honchar, Y., Krasnytska, M., Berche, B., Holovatch, Y., Kenna, R.: Partition function zeros for the Blume–Capel model on a complete graph. Low Temp. Phys.51(5), 567–577 (2025) https://doi.org/10. 1063/10.0036500

    2025

  32. [32]

    Mancini, F.: The extended Hubbard model in the ionic limit. Eur. Phys. J. B47(4), 527–540 (2005) https://doi.org/10.1140/epjb/e2005-00357-2

    work page doi:10.1140/epjb/e2005-00357-2 2005

  33. [33]

    Fisher, M.E., Ruelle, D.: The stability of many-particle systems. J. Math. Phys.7(2), 260–270 (1966) https://doi.org/10.1063/1.1704928

    work page doi:10.1063/1.1704928 1966

  34. [34]

    Heyes, D.M., Rickayzen, G.: The stability of many-body systems. J. Phys.: Condens. Matter19(41), 416101 (2007) https://doi.org/10.1088/0953-8984/19/41/416101

    work page doi:10.1088/0953-8984/19/41/416101 2007

  35. [35]

    Bender, C.M., Orszag, S.A.: Advanced Mathematical Methods for Scientists and Engineers I: Asymptotic Methods and Perturbation Theory, p. 593. Springer, New York (1999)

    1999

  36. [36]

    The cell fluid model with Curie-Weiss interactions: special cases and analytical results

    Dobush, O.A., Kozlovskii, M.P., Pylyuk, I.V., Romanik, R.V., Shpot, M.A.: The cell fluid model with Curie-Weiss interactions: special cases and analytical results. arXiv:2605.22520 (2026). https://doi.org/10. 48550/arXiv.2605.22520

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  37. [37]

    Physica A681, 131133 (2026) https://doi

    Dobush, O.A., Kozlovskii, M.P., Pylyuk, I.V., Plevachuk, Y.O.: Influence of microscopic parameters on phase behavior of a cell model with Curie–Weiss interaction. Physica A681, 131133 (2026) https://doi. org/10.1016/j.physa.2025.131133

    work page doi:10.1016/j.physa.2025.131133 2026

  38. [38]

    Physica A692, 131510 (2026) https://doi.org/10

    Kozlovskii, M.P., Dobush, O.A., Pylyuk, I.V., Romanik, R.V., Shpot, M.A.: Triple point in a cell fluid model with effective temperature-dependent attraction. Physica A692, 131510 (2026) https://doi.org/10. 1016/j.physa.2026.131510

    work page arXiv 2026

  39. [39]

    Franzese, G.: Differences between discontinuous and continuous soft-core attractive potentials: The appear- ance of density anomaly. J. Mol. Liq.136(3), 267–273 (2007) https://doi.org/10.1016/j.molliq.2007.08. 021

    work page doi:10.1016/j.molliq.2007.08 2007

  40. [40]

    Louis, A.A., Bolhuis, P.G., Hansen, J.P.: Mean-field fluid behavior of the Gaussian core model. Phys. Rev. E62, 7961–7972 (2000) https://doi.org/10.1103/PhysRevE.62.7961

    work page doi:10.1103/physreve.62.7961 2000

  41. [41]

    Nature403, 170–173 (2000) https://doi.org/10.1038/ 35003143

    Katayama, Y., Mizutani, T., Utsumi, W., Shimomura, O., Yamakata, M., Funakoshi, K.-i.: A first- order liquid–liquid phase transition in phosphorus. Nature403, 170–173 (2000) https://doi.org/10.1038/ 35003143

    2000

  42. [42]

    Monaco, G., Falconi, S., Crichton, W.A., Mezouar, M.: Nature of the first-order phase transition in fluid phosphorus at high temperature and pressure. Phys. Rev. Lett.90, 255701 (2003) https://doi.org/10.1103/ PhysRevLett.90.255701 16

    2003

  43. [43]

    Yang, M., Karmakar, T., Parrinello, M.: Liquid-liquid critical point in phosphorus. Phys. Rev. Lett.127, 080603 (2021) https://doi.org/10.1103/PhysRevLett.127.080603

    work page doi:10.1103/physrevlett.127.080603 2021

  44. [44]

    Science306(5697), 845–848 (2004) https://doi.org/10.1126/science.1103073

    Kurita, R., Tanaka, H.: Critical-like phenomena associated with liquid-liquid transition in a molecular liquid. Science306(5697), 845–848 (2004) https://doi.org/10.1126/science.1103073

    work page doi:10.1126/science.1103073 2004

  45. [45]

    de Oliveira, A.B., Netz, P.A., Colla, T., Barbosa, M.C.: Thermodynamic and dynamic anomalies for a three-dimensional isotropic core-softened potential. J. Chem. Phys.124(8), 084505 (2006) https://doi.org/ 10.1063/1.2168458

    work page doi:10.1063/1.2168458 2006

  46. [46]

    arXiv.2511.15362 (To be published in Ukr

    Romanik, R.V., Dobush, O.A., Kozlovskii, M.P., Pylyuk, I.V., Shpot, M.A.: Pair distribution function of the cell fluid model with Curie-Weiss interaction. arXiv.2511.15362 (To be published in Ukr. J. Phys.,71, 8 (2026)) (2025). https://doi.org/10.48550/arXiv.2511.15362 17

    work page doi:10.48550/arxiv.2511.15362 2026