Pith. sign in

REVIEW 3 major objections 4 minor 27 references

Influence of microscopic parameters on phase behavior of a cell model with Curie-Weiss interaction

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

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

Pith's one-line read The paper establishes that cell volume leaves the pressure-temperature and temperature-density phase diagrams of the Curie-Weiss cell model unchanged, while the repulsion-to-attraction ratio produces quantitative shifts only, with no…

desk verdict A modest but clean parameter scan of an exactly solvable cell model; the v-independence is proven and the f>3 saturation is new, but Eq. (6) has a typo and the global stability check is not documented. read the letter →

arxiv 2505.14456 v2 pith:JNLP3FO6 submitted 2025-05-20 cond-mat.stat-mech

O.A. Dobush , M.P. Kozlovskii , I.V. Pylyuk , Yu.O. Plevachuk This is my paper
classification cond-mat.stat-mech
keywords cellmodelCurie-Weissinteractionfirst-orderphasetransitionsdiagramsgrandcanonicalensembleequationofstaterepulsion-to-attractionratiocriticalpoint
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

This paper asks which microscopic parameters actually change the phase behavior of an exactly solvable fluid model in which space is divided into cells and particles interact through a Curie-Weiss attraction together with a repulsion that acts only within a cell. The answer it argues for is that the cell volume $v$ matters only through an additive shift of the chemical potential, so the pressure-temperature and temperature-density phase diagrams are independent of cell size. The ratio $f = g_r/g_a$ of repulsive to attractive intensities, by contrast, does change the diagrams quantitatively: larger $f$ moves the coexistence lines to higher pressures, lowers the second and third critical temperatures toward the first, and eventually, for $f > 3$, makes all critical temperatures coincide near $T_c$ while the critical densities stay close to $\bar n = 0.5, 1.5, 2.5$. The central conclusion is qualitative stability: over the whole range $1 < f \le 10$, the model always shows the same sequence of three first-order transitions and four stable phases, with no triple point.

What carries the argument

The engine of the analysis is the family of special functions $K_m(z) = \sum_{n=0}^{\infty} \frac{1}{n!\,n^m} \exp(z n - f p_a n^2/2)$, which encode the competition between attraction and same-cell repulsion. The first two members, $K_0$ and $K_1$, give the density via $\bar n = K_1(\bar z)/K_0(\bar z)$, and together with the parameter $p_a = \beta g_a$ they give the pressure through equation (16). The cell volume never appears in these objects; it appears only as $-\ln v^*$ in the chemical potential, so it drops out of every phase-equilibrium calculation. Phase transitions are located by solving the coexistence equations $E_0(z_1) = E_0(z_2)$ and $\mu(z_1) = \mu(z_2)$, and critical points come from the conditions $d\mu/dz = 0$ and $d^2\mu/dz^2 = 0$.

What would settle it

Evaluate the same grand partition function directly, by summing the defining series for a large but finite number of cells or by Monte Carlo simulation of the Curie-Weiss cell Hamiltonian, at f = 1.2 and f = 2.0, and compare the predicted pressure plateaus and critical points with the reported values (e.g., at f = 1.2 the three coexistence pressures 0.198, 1.719, and 4.157, with critical densities near 0.51, 1.51, and 2.50); any extra or missing first-order transition would show that the saddle-point evaluation or the coexistence equations miss a phase branch.

Watch

Extended reading notes

Core claim

The paper studies a fluid model in which space is divided into cells of volume $v$, particles experience a Curie-Weiss attraction of strength $g_a$ with every other particle, and a repulsion $g_r$ acts between any two particles sharing a cell; the grand partition function of this model has an exact thermodynamic-limit evaluation. Using that solution, the paper writes the reduced density as $\bar n = K_1(\bar z)/K_0(\bar z)$, the chemical potential as $\beta\mu = \bar z - p_a K_1(\bar z)/K_0(\bar z) - \ln v^*$, and the pressure through $P v \beta_c = (\tau+1)[\ln K_0(\bar z) - (p_a/2)(K_1(\bar z)/K_0(\bar z))^2]$. The cell volume $v^*$ enters only through the additive $-\ln v^*$ term in the chemical potential, so it cancels from coexistence conditions and leaves the pressure and density untouched. The ratio $f = g_r/g_a$ enters the special functions $K_m(z)$ through the quadratic term $f p_a n^2/2$ in the exponent, and increasing $f$ raises the coexistence pressures of the second and third transitions, pulls the higher critical temperatures down toward $T_c$, and for $f > 3$ makes all critical temperatures equal to $T_c$ while the critical densities stay near $\bar n = 0.5, 1.5, 2.5$. The paper's conclusion is that the phase diagram is universal in $v^*$ and only quantitatively sensitive to $f$: the model always exhibits three first-order transitions and four stable phases, with no triple point, for $1 < f \le 10$.

Load-bearing premise

The load-bearing premise is that the earlier exact solution, which evaluates the grand partition function in the thermodynamic limit and locates phase coexistence by equating pressure and chemical potential, correctly identifies all stable phases at every temperature studied here; the present paper relies on that solution rather than proving it anew.

Editorial extensions

If this is right

  • For any fixed $f$, the pressure-temperature and temperature-density phase diagrams hold for every cell volume; changing $v^*$ only shifts the chemical potential scale.
  • Between $f = 1$ and $f = 3$, the second and third transitions have critical temperatures above $T_c$ (e.g., $1.029\,T_c$ and $1.042\,T_c$ at $f = 1.2$); for $f > 3$ these temperatures coincide with $T_c$ to many decimal places.
  • Increasing $f$ stretches the phase diagram along the pressure axis: the second coexistence line moves from pressures around 3 at $f=1.5$ to around 17 at $f=5$, while the critical densities remain near $\bar n = 1.5$ and $2.5$.
  • Within the studied range $1 < f \le 10$, the model always has the same sequence of three first-order transitions and four stable phases, with no triple point.

Reading between the lines

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

  • If the exact solution is correct, the model shows that discretization scale can be a pure gauge parameter in a mean-field cell description: one can shift $v^*$ to fit chemical-potential scales without altering coexistence properties.
  • The saturation at $f > 3$, where every critical transition sits at $p_{ac} = 4.0$, suggests a universal strongly-repulsive regime in which the quadratic term in $K_m$ dominates; varying cell shape or adding a finite-range repulsion is a direct way to test whether this saturation is an artifact of the interaction form.
  • The absence of triple points across all studied $f$ is a structural prediction of the two-parameter interaction; adding a second attraction scale or an external field would be a natural way to see whether triple points can emerge in this class of cell models.
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 / 4 minor

Summary. The paper studies a cell model with Curie-Weiss interaction, using the exact grand-canonical solution from ref. [17] to examine how the cell volume v and the repulsion-to-attraction ratio f = gr/ga affect phase behavior. It reports that v shifts the chemical potential but leaves pressure and all coexistence and critical quantities unchanged, while increasing f changes critical parameters and coexistence pressures quantitatively without, within the studied range, altering the qualitative structure of three sequential first-order transitions and no triple points. Results are presented as chemical-potential-density and pressure-density isotherms, pressure-temperature and temperature-density phase diagrams, and tables of critical coordinates.

Significance. If the claims hold, the paper provides a systematic parametric map of a solvable continuum model with multiple first-order transitions, showing that a single dimensionless ratio f controls the entire phase diagram and that the cell volume is irrelevant except for a chemical-potential shift. This is a useful benchmark result for approximation methods and gives a concrete, falsifiable prediction within the model, most notably the saturation of all critical temperatures at Tc for f > 3. The manuscript does not introduce new theory and is heavily dependent on the previously published exact solution [17]; its main contribution is the quantitative survey and the documentation of the saturation behavior. The paper would be strengthened by a short description of the numerical method or by making the data available in machine-readable form.

major comments (3)
  1. [Section 2, Eq. (6)] Equation (6) as printed defines K_m(z) = sum_{n=0}^∞ 1/(n! n^m) exp(...), which is singular at n=0 for m≥1 and would force K_1/K_0 ≤ 1, contradicting the reduced densities n̄ ≈ 0.5, 1.5, 2.5 reported in Table 2 and Fig. 2. The intended definition must be K_m(z) = sum_{n=0}^∞ (n^m/n!) exp(...), which is consistent with Eq. (10) and with the reported values. Because every table and figure depends on K_m, this typo must be corrected; as printed, the equations are not reproducible.
  2. [Section 4, Fig. 2 and Eq. (15)] The coexistence conditions (15) are necessary but not sufficient for stable phase equilibrium: if a third stationary point z3 of E0(z; μ) has E0(z3) > E0(z1)=E0(z2), then neither z1 nor z2 is the stable phase, and the plotted 'stable' regions would in fact be metastable. The manuscript nowhere states that a global maximum over all stationary points of E0 was performed, nor does it explain how the regions labeled 'stable', 'metastable', and 'unstable' in Fig. 2 were assigned at each temperature. Since the central qualitative claims (order of the three coexistence lines, absence of triple points) are global statements about stability, this check must be documented explicitly.
  3. [Conclusions, last paragraph] The statement that changes in f 'do not lead to qualitative alterations' is a universal claim over f ≥ 1, but the evidence presented is a finite set of values (f = 1.01, 1.2, 1.5, 2.0, 5.0 in Fig. 2; f up to 10 in Table 2). No proof or exhaustive scan is given that triple points never appear or that the three-transition sequence persists for all f in the stability range. Please either restrict the claim to the studied interval or supply an argument, for example based on monotonicity of the coexistence curves, that no qualitative change occurs.
minor comments (4)
  1. [Section 2, Eqs. (4) and (8)] The prefactor (N_v/(2π p_a))^{1/2} in Eq. (4) is omitted in Eq. (8); since it is subdominant in ln Ξ this is harmless, but the omission should be stated explicitly.
  2. [Section 3, Fig. 1b caption] The phrase 'at arbitrary values of v*' is imprecise; the curve is identical for all tested values v* = 1, 2, 5, 10, 15, but that is not the same as an arbitrary continuous value.
  3. [Section 4, Fig. 2 (b,d,f,h,j)] The spinodal curves (thin red lines) are not defined in the text; the authors should specify the condition used (e.g., dμ/dn̄ = 0) and the numerical procedure by which they were generated.
  4. [Section 4, Table 2] The row for f = 1.00001 is presented as representative of the f → 1 limit; it would help to state explicitly that f = 1 is excluded by the stability condition gr > ga.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: computed parameter study, no fitted predictions.

full rationale

The chain of derivation is a numerical parameter study built on the explicit closed-form expressions (9), (10), (16) and the coexistence conditions (15). The two varied parameters v and f are inputs of the model, not fitting parameters tuned to any target data. The statement that the cell volume does not affect the pressure-temperature-density phase diagram is a direct, stated consequence of the equation of state (16), which contains no v, while the only v-dependence is the additive -ln v shift in the chemical potential (9); this is a mathematical implication of the cited exact solution, not a prediction forced by a fit. The f-study likewise evaluates the same equations for fixed values f = 1.01, 1.2, 1.5, 2.0, 5.0 without adjusting any parameter to reproduce the reported saturation of critical temperatures or the absence of triple points. The only external input is the saddle-point evaluation of the grand partition function (4)-(8), attributed to the authors' earlier work [17]; that citation supplies the starting equations, and the manuscript explicitly lists all subsequent relations and computes the phase diagrams from them, so the central claims are not imported as conclusions. Concerns about whether condition (15) selects globally stable phases are a mathematical-rigor issue, not a circularity. No fitted parameter is renamed as a prediction, no definition is made in terms of the target result, and no self-citation is used to forbid alternatives.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

No parameters are fitted to data; v and f are model inputs varied across chosen values. The load-bearing assumptions are the validity of the self-cited exact solution from ref. [17], the Laplace evaluation, the stability condition f>=1, and the completeness of the coexistence equations. No new physical entities are introduced.

assumptions (4)
  • standard math The grand partition function can be evaluated asymptotically by the Laplace method, yielding Eq. (8) as the thermodynamic-limit free energy.
    Section 2, Eqs. (4)-(8) invoke the Laplace method and cite ref. [27]; the paper uses this result without proving convergence or branch selection.
  • domain assumption The exact solution of the cell model with Curie-Weiss interaction in ref. [17] is correct and applies at all temperatures considered, including T/Tc = 0.02.
    All numerical results inherit the validity of ref. [17]; the current paper does not rederive or independently test this solution.
  • domain assumption Stability requires gr > ga, i.e., f >= 1, per Ruelle [26]; the study restricts to 1 < f <= 10.
    Section 4 and Conclusions cite [26]; if the stability bound were different, the phase diagrams beyond f=1 would be invalid.
  • domain assumption The coexistence conditions (15), E0(z1)=E0(z2) and μ(z1)=μ(z2), identify all stable first-order transitions in the sequence.
    Section 2, Eq. (15); the paper assumes no missing transitions or triple points for the parameter ranges explored.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Influence of microscopic parameters on phase behavior of a cell model with Curie-Weiss interaction." pith.science (2026). https://pith.science/paper/JNLP3FO6

@misc{pith2026250514456,
  author       = {Pith},
  title        = {Pith review of: Influence of microscopic parameters on phase behavior of a cell model with Curie-Weiss interaction},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/JNLP3FO6}},
  note         = {Machine review of arXiv:2505.14456}
}
read the original abstract

We investigate how varying two microscopic parameters - cell volume and the ratio between repulsion and attraction intensities - affect the phase behavior of a cell model with a Curie-Weiss-type interaction. The analysis is based on an exact solution previously derived for this model in the grand canonical ensemble. At sufficiently low temperatures, the cell model exhibits multiple first-order phase transitions. By varying the cell volume and the repulsion-to-attraction ratio, we represent a quantitative comparison of the chemical potential and pressure isotherms, along with the pressure-temperature and temperature-density phase diagrams. Our results demonstrate that altering these microscopic parameters induces quantitative changes in the phase diagrams of the cell model.

Figures

Figures reproduced from arXiv: 2505.14456 by the authors.

Figure 1
Figure 1. Part a: influence of varying the cell size, [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Phase diagrams in the pressure-temperature (plots a, c, e, g, i) and temperature [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

27 extracted references · 23 canonical work pages

  1. [17]

    KozlovskiiI, O.A

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

    work page arXiv 2022

  2. [1]

    Hansen, I.R

    J.-P. Hansen, I.R. McDonald, Theory of Simple Liquids: With Applications to Soft Matter, Academic Press, Oxford, 2013

    work page 2013

  3. [2]

    Hagen, D

    M.H.J. Hagen, D. Frenkel, Determination of phase diagrams for the hard-core attractive Yukawa system, J. Chem. Phys. 101 (1994) 4093–4097, https://doi.org/10.1063/1.467526

    work page doi:10.1063/1.467526 1994

  4. [3]

    Dijkstra, Phase behavior of hard spheres with a short-range Yukawa attraction, Phys

    M. Dijkstra, Phase behavior of hard spheres with a short-range Yukawa attraction, Phys. Rev. E 66 (2002) 021402, https://doi.org/10.1103/PhysRevE.66.021402

    work page doi:10.1103/physreve.66.021402 2002

  5. [4]

    X. Wang, S. Ram´irez-Hinestrosa, J. Dobnikar, D. Frenkel, The Lennard-Jones poten- tial: when (not) to use it, Physical Chemistry Chemical Physics 22 (2020) 10624, https://doi.org/10.1039/C9CP05445F

    work page doi:10.1039/c9cp05445f 2020

  6. [5]

    Phase transitions and cluster structures of the new finite range Lennard-Jones like model

    O.-F. Adesida, S. Havens, L.B. P´ artay, Phase transitions and cluster struc- tures of the new finite range Lennard-Jones like model, arXiv:2407.14688 (2024), https://doi.org/10.48550/arXiv.2407.14688

    work page Pith review arXiv doi:10.48550/arxiv.2407.14688 2024

  7. [6]

    Gonz´ alez-Melchor, A

    M. Gonz´ alez-Melchor, A. Trokhymchuk, J. Alejandre, Surface tension at the vapor/liquid interface in an attractive hard-core Yukawa fluid, J. Chem. Phys. 155 (2001) 3862–3872, https://doi.org/10.1063/1.1384553

    work page doi:10.1063/1.1384553 2001

  8. [7]

    Galicia-Pimentel, J

    U.F. Galicia-Pimentel, J. L´ opez-Lemus, P. Orea, Liquid–vapor interfacial prop- erties of attractive Yukawa fluids, Fluid Phase Equilibria 265 (2008) 205–208, https://doi.org/10.1016/j.fluid.2007.11.021

    work page doi:10.1016/j.fluid.2007.11.021 2008

Show all 27 references
  1. [8]

    J.K. Singh, Surface tension and vapour–liquid phase coexistence of variable- range hard-core attractive Yukawa fluids, Molecular Simulation 35 (2009) 880–887, https://doi.org/10.1080/08927020902787796

    2009 doi

  2. [9]

    Lomba, N.G

    E. Lomba, N.G. Almarza, Role of the interaction range in the shaping of phase diagrams in simple fluids. The hard sphere Yukawa fluid as a case study, J. Chem. Phys. 100 (1994) 8367–8372, https://doi.org/10.1063/1.466781

    1994 doi

  3. [10]

    Soares, A.G

    E.D.A. Soares, A.G. Barreto Jr., F.W. Tavares, A self-consistent perturba- tive density functional theory for hard-core fluids: phase diagrams, struc- tural and interfacial properties, Fluid Phase Equilibria 542–543 (2021) 113095, https://doi.org/10.1016/j.fluid.2021.113095

    2021

  4. [11]

    L. Vega, E. de Miguel, L.F. Rull, G. Jackson, I.A. McLure, Phase equilibria and critical behavior of square-well fluids of variable width by Gibbs ensemble Monte Carlo simulation, J. Chem. Phys. 96 (1992) 2296–2305, https://doi.org/10.1063/1.462080

    1992 doi

  5. [12]

    Fuentes-Herrera, J.A

    M. Fuentes-Herrera, J.A. Moreno-Razo, O. Guzm´ an, J. L´ opez-Lemus, B. Ibarra-Tandi, Separating the effects of repulsive and attractive forces on the phase diagram, in- terfacial, and critical properties of simple fluids, J. Chem. Phys. 144 (2016) 214502, https://doi.org/10.1...

    2016 doi

  6. [13]

    Werth, K

    S. Werth, K. St¨ obener, M. Horsch, H. Hasse, Simultaneous description of bulk and inter- facial properties of fluids by the Mie potential, Molecular Physics 115 (2017) 1017–1030, https://doi.org/10.1080/00268976.2016.1206218

    2017

  7. [14]

    Stephan, M

    S. Stephan, M. Urschel, Characteristic curves of the Mie fluid, J. Mol. Liq. 383 (2023) 122088, https://doi.org/10.1016/j.molliq.2023.122088. 10

    2023

  8. [15]

    Errington, A.Z

    Je.R. Errington, A.Z. Panagiotopoulos, Phase Equilibria of the Modified Buckingham Exponential-6PotentialfromHamiltonianScalingGrandCanonicalMonteCarlo, J.Chem. Phys. 109 (1998) 1093–1100, https://doi.org/10.1063/1.476652

    1998 doi

  9. [16]

    Panagiotopoulos, Monte Carlo methods for phase equilibria of fluids, J

    A.Z. Panagiotopoulos, Monte Carlo methods for phase equilibria of fluids, J. Phys.: Con- dens. Matter 12 (2000) R25–R52, https://doi.org/10.1088/0953-8984/12/3/201

    2000 doi

  10. [18]

    Rebenko, Cell gas model of classical statistical systems, Rev

    A.L. Rebenko, Cell gas model of classical statistical systems, Rev. Math. Phys. 25 (2013) 1330006, https://doi.org/10.1142/S0129055X13300069

    2013 doi

  11. [19]

    Boluh, A.L

    V.A. Boluh, A.L. Rebenko. Cell gas free energy as an approximation of the continuous model, J. Mod. Phys. 6 (2015) 168–175, https://doi.org/10.4236/jmp.2015.62022

    2015

  12. [20]

    Kozitsky, M

    Y. Kozitsky, M. Kozlovskii, O. Dobush, Phase transitions in a continuum Curie-Weiss system: A quantitative analysis, in: L.A. Bulavin, A.V. Chalyi (Eds.), Modern Problems of Molecular Physics, Springer, Cham, 2018, pp. 229–251, https://doi.org/10.1007/978-3- 319-61109-9_11

    2018 doi

  13. [21]

    Kozitsky, M.P

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

    2020 doi

  14. [22]

    Kozlovskii, I.V

    M.P. Kozlovskii, I.V. Pylyuk, O.A. Dobush, The equation of state of a cell fluid model in the supercritical region, Condens. Matter Phys. 21 (2018) 43502, https://doi.org/10.5488/CMP.21.43502

    2018 doi

  15. [23]

    Pylyuk, Fluid critical behavior at liquid–gas phase transition: Ana- lytic method for microscopic description, J

    I.V. Pylyuk, Fluid critical behavior at liquid–gas phase transition: Ana- lytic method for microscopic description, J. Mol. Liq. 310 (2020) 112933, https://doi.org/10.1016/j.molliq.2020.112933

    2020

  16. [24]

    Pylyuk, M.P

    I.V. Pylyuk, M.P. Kozlovskii, O.A. Dobush, M.V. Dufanets, Morse fluids in the immediate vicinity of the critical point: Calculation of thermodynamic coefficients, J. Mol. Liq. 385 (2023) 122322, https://doi.org/10.1016/j.molliq.2023.122322

    2023

  17. [25]

    Hill, Statistical Mechanics, McGraw-Hill, New York, 1956

    T.L. Hill, Statistical Mechanics, McGraw-Hill, New York, 1956

    1956

  18. [26]

    18 (1970) 127–159, https://doi.org/10.1007/BF01646091

    D.Ruelle, Superstableinteractionsinclassicalstatisticalmechanics, Commun.Math.Phys. 18 (1970) 127–159, https://doi.org/10.1007/BF01646091

    1970 doi

  19. [27]

    Fedoryuk, Asymptotic methods in analysis, in: R.V

    M.V. Fedoryuk, Asymptotic methods in analysis, in: R.V. Gamkrelidze (Ed.), Analy- sis I: Integral Representations and Asymptotic Methods, Encyclopaedia of Mathematical Sciences, vol. 13, Springer-Verlag, Berlin Heidelberg, 1989, pp. 83–191. 11

    1989

Pith tools

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