Entropy of the cell fluid model with Curie-Weiss interaction

Basic Information

• Paper ID: 2510.22700
• Title: Entropy of the cell fluid model with Curie-Weiss interaction
• Authors: R. V. Romanik, O. A. Dobush, M. P. Kozlovskii, I. V. Pylyuk, M. A. Shpot
• Institution: Yukhnovskii Institute for Condensed Matter Physics, National Academy of Sciences of Ukraine
• Classification: cond-mat.stat-mech (Condensed Matter Physics - Statistical Mechanics)
• Submission Date: October 26, 2025
• Paper Link: https://arxiv.org/abs/2510.22700v1

Abstract

This paper derives analytical expressions for the entropy of a cell fluid model with Curie-Weiss interaction, expressed as functions of temperature and chemical potential. Parametric equations representing entropy as a function of density are derived. The characteristics of single-particle entropy and single-cell entropy are investigated in the entropy-density and entropy-chemical potential planes. The considered cell model is a multi-occupancy model that exhibits an infinite sequence of first-order phase transitions at sufficiently low temperatures. The study reveals that entropy exhibits pronounced minima near integer-valued particle densities, which may be a general characteristic of multi-occupancy models.

Research Background and Motivation

Research Problem

This paper investigates the entropy characteristics of multi-occupancy lattice gas models, particularly the cell fluid model with Curie-Weiss interaction. This model allows each cell to contain an arbitrary number of particles, in contrast to classical lattice gases where each lattice site contains at most one particle.

Importance of the Problem

1. Theoretical Value: Multi-occupancy models can reproduce fundamental thermodynamic characteristics of fluids with soft or penetrable interactions, representing important theoretical objects with non-trivial collective behavior in statistical mechanics
2. Central Role of Entropy: Entropy directly measures configurational degeneracy and is closely related to the occurrence of phase transitions
3. Research Gap: Although some studies have explored the thermodynamics and phase properties of multi-occupancy models, their entropy characteristics remain largely unexplored

Limitations of Existing Methods

• Research on multi-occupancy lattice gas models is far less extensive than on classical lattice gases
• Systematic analytical studies of entropy dependencies (particularly entropy-density and entropy-chemical potential relationships) for multi-occupancy models are lacking
• Existing literature primarily focuses on quantum systems (such as Hubbard models), with classical multi-occupancy systems' entropy studies nearly absent

Research Motivation

By exploiting the exact solvability of mean-field Curie-Weiss interactions, this work provides a convenient analytical reference for multi-occupancy systems, filling a research gap in this field and exploring universal entropy characteristics of multi-occupancy models.

Core Contributions

1. Analytical Entropy Expression: First derivation of exact analytical entropy expressions for multi-occupancy cell fluid models with Curie-Weiss interaction, including single-particle entropy S* and single-cell entropy S*_v
2. Parametric Representation: Establishment of parametric equations for entropy as a function of density, facilitating analysis of entropy-density relationships
3. Phase Transition Characteristic Analysis: Systematic investigation of discontinuous entropy jumps in infinite sequences of first-order phase transitions
4. Discovery of Integer Density Minima: Discovery that entropy exhibits pronounced minima near integer-valued densities, with argumentation that this is a universal characteristic of multi-occupancy models
5. Cross-Model Comparison: Through comparison with quantum systems such as Hubbard and Bose-Hubbard models, revelation of common features of multi-occupancy models

Detailed Methodology

Task Definition

Investigation of point particle entropy properties in open systems, where the system volume V is divided into N_v non-overlapping cubic cells, each with volume v. The objectives are to obtain:

• Single-particle entropy: S* = S/(k_B⟨N⟩)
• Single-cell entropy: S_v = S/(k_B N_v) as functions of temperature T, chemical potential μ*, and density ρ*

Model Architecture

Interaction Definition

Two-particle interaction energy is defined as: $\Phi_{N_v}(r_i, r_j) = -\frac{J_1}{N_v} + J_2 \sum_{l=1}^{N_v} I_{\Delta_l}(r_i)I_{\Delta_l}(r_j)$

• First term: Global Curie-Weiss (mean-field) attraction with strength parameter J_1 > 0
• Second term: Local repulsion between particles within the same cell with strength parameter J_2 > 0
• I_Δl(r): Indicator function of cell Δ_l

Reduced Variables

Standard dimensionless variables are employed:

• Reduced temperature: T* = k_B T/J_1
• Reduced chemical potential: μ* = μ/J_1
• Reduced density: ρ* = ⟨N⟩/(N_v)
• Reduced pressure: P* = Pv/J_1
• Parameter: a = J_2/J_1, v* = v/λ³ (λ related to thermal wavelength)

Core Formulas

Single-Particle Entropy

$S^*(T^*, \mu^*) = \left(\frac{3}{2} - \frac{\mu^*}{T^*}\right) - \frac{1}{T^*}\frac{K_1}{K_0} + \frac{K_0\ln K_0}{K_1} + \frac{a}{2T^*}\frac{K_2}{K_1}$

where special functions K_j are defined as: $K_j(T^*, \mu^*; \bar{y}_{max}) = \sum_{n=0}^{\infty} \frac{n^j(v^*T^{*3/2})^n}{n!} \exp\left[\left(\bar{y}_{max} + \frac{\mu^*}{T^*}\right)n - \frac{a}{2T^*}n^2\right]$

Single-Cell Entropy

$S_v^*(T^*, \mu^*) = \left(\frac{3}{2} - \frac{\mu^*}{T^*}\right)\frac{K_1}{K_0} - \frac{1}{T^*}\frac{K_1^2}{K_0^2} + \ln K_0 + \frac{a}{2T^*}\frac{K_2}{K_0}$

Parametric Representation

Introducing a new function: $\bar{z}_{max}(T^*, \mu^*) = \bar{y}_{max} + \frac{\mu^*}{T^*}$

Defining: $\tilde{K}_j(T^*; \bar{z}_{max}) = \sum_{n=0}^{\infty} \frac{n^j(v^*T^{*3/2})^n}{n!} \exp\left[\bar{z}_{max}n - \frac{a}{2T^*}n^2\right]$

The parametric representation of entropy becomes: $S^*(T^*; \bar{z}_{max}) = \frac{3}{2} - \bar{z}_{max} + \frac{\tilde{K}_0\ln\tilde{K}_0}{\tilde{K}_1} + \frac{a}{2T^*}\frac{\tilde{K}_2}{\tilde{K}_1}$

Density: $\rho^*(T^*; \bar{z}_{max}) = \frac{\tilde{K}_1(T^*; \bar{z}_{max})}{\tilde{K}_0(T^*; \bar{z}_{max})}$

Chemical potential: $\mu^* = T^*\bar{z}_{max} - \frac{\tilde{K}_1(T^*; \bar{z}_{max})}{\tilde{K}_0(T^*; \bar{z}_{max})}$

Technical Innovations

1. Exact Solvability: Exploitation of the mean-field nature of Curie-Weiss interactions to obtain exact solutions in the thermodynamic limit
2. Parametrization Technique: Transformation of complex implicit function relationships into parametric equations through introduction of z̄_max parameter, facilitating numerical computation and analysis
3. Extremum Conditions: Determination of ȳ_max through extremum conditions of functions E_1 and E_2, ensuring thermodynamic stability
4. Ideal Gas Contribution Separation: Clear identification of ideal gas contribution terms (3/2 - μ*/T*) in entropy expressions, facilitating understanding of interaction effects

Experimental Setup

Parameter Selection

• Interaction Parameter: a = J_2/J_1 = 1.2 (ratio of repulsive to attractive strength)
• Volume Parameter: v* = 5.0
• Temperature Range:
  • T* = 0.40 (well above critical temperature)
  • T* = 0.30 (slightly above critical temperature)
  • T* = 0.25 (slightly below critical temperature)
  • T* = 0.20 (significantly below critical temperature)

Critical Point Data

The paper provides detailed numerical data for the first four critical points (see Table 1):

• 1st critical point: T_c = 0.254567, ρ_c = 0.513896
• 2nd critical point: T_c = 0.261881, ρ_c = 1.50557
• 3rd critical point: T_c = 0.265254, ρ_c = 2.50303
• 4th critical point: T_c = 0.267242, ρ_c = 3.50191

Coexistence Line Data

At temperatures T* = 0.25 and T* = 0.20, coexistence line data for the first four phase transitions are provided (see Table 2), including density, entropy, chemical potential, and other physical quantities.

Experimental Results

Main Results

1. Entropy-Chemical Potential Relationship (Figure 1a)

• Supercritical Region (red T*=0.40 and green T*=0.30 curves): Single-particle entropy S* decreases monotonically with chemical potential μ*
• Subcritical Region (black T*=0.25 and blue T*=0.20 curves):
  • Entropy of phase II is lower than phase I, phase III lower than phase II, decreasing sequentially
  • At phase transition points, entropy jumps from higher to lower values
  • Except in phase I, entropy dependence on chemical potential within individual phases may be non-monotonic

2. Entropy-Density Relationship (Figure 1b)

• Supercritical Temperature: Single-particle entropy S* is a decreasing function of density ρ*
• Subcritical Temperature:
  • Entropy remains decreasing in phase I
  • Key Finding: In subsequent phases, entropy exhibits pronounced minima near integer density values (ρ* ≈ 1, 2, 3...)
  • Both entropy and density undergo discontinuous jumps at first-order phase transitions
  • Entropy becomes negative in high-density phases (phase V and higher) - an artifact of classical treatment

Specific numerical values (T* = 0.20, phase II):

• At ρ* = 1.12438, S* = 0.996714
• At ρ* = 1.88058, S* = 0.601031
• Minimum exists near ρ* ≈ 1.5

3. Single-Cell Entropy-Chemical Potential Relationship (Figure 1c)

• High Temperature (T*=0.40): S*_v first increases then decreases, with a maximum
• Near Critical Temperature (T*=0.30): Complex behavior with pronounced curvature near the critical point
• Low Temperature:
  • S_v increases with μ in phase I
  • Phases II and III exhibit highest single-cell entropy
  • Subsequent phases jump to lower values at phase transitions

4. Single-Cell Entropy-Density Relationship (Figure 1d)

• High Temperature: S*_v first increases then decreases
• Low Temperature:
  • Increasing in phase I
  • Non-monotonic behavior within each subsequent phase
  • Minima exist near integer density values (ρ* ≈ 1, 2, 3...)

Important Findings

Universality of Integer Density Minima

The paper argues through comparison with multiple quantum models that integer density entropy minima are universal characteristics of multi-occupancy models:

1. Three-Dimensional Fermionic Hubbard Model16 (maximum occupancy=2):
   • Single-site entropy shows pronounced dip at ρ*=1.0
   • Qualitatively similar to results at T*=0.3 in this work
2. One-Dimensional Repulsive Hubbard Model17,18:
   • Single maximum in single-site entropy at high temperature
   • Pronounced minimum at ρ*=1.0 at low temperature
3. Bose-Hubbard Model19 (reported to filling number n=3.0):
   • Two depressions observed at n=1.0 and n=2.0
   • Corresponding to minima at ρ*≈1.0 and ρ*≈2.0 in this work

Unified Pattern

• Hubbard model (ρ_max=2.0) → minimum at ρ=1.0
• Bose-Hubbard model (reported to ρ*≤3.0) → minima at ρ*=1.0 and 2.0
• Curie-Weiss model in this work (no occupancy limit) → infinite sequence of minima at ρ*≈1, 2, 3...

Negative Entropy Phenomenon

At high-density phases (phase V and higher) or small v* values, entropy becomes negative. This is an artifact of classical treatment, similar to Stillinger's observation in Gaussian core models. Increasing v* pushes this behavior to higher densities, but the qualitative pattern remains unchanged.

Related Work

Multi-Occupancy Lattice Gas Models

1. Lee-Yang Theory1: Statistical theory foundation for classical lattice gases
2. Soft Interaction Fluids3-8:
   • Gaussian core systems3
   • Penetrable spheres4
   • Double-occupancy lattice gases7,8
3. Mean-Field Theory9,10: Curie-Weiss interaction models on which this work is based
4. Multi-Layer Adsorption Models11: Lattice gases allowing multi-occupancy
5. Cluster Crystal Models12,13: Generalized models of exponent 4

Quantum Multi-Occupancy Models

1. Hubbard Model16-18:
   • Dynamical mean-field theory studies of three-dimensional fermionic Hubbard model
   • Density functional theory of one-dimensional repulsive Hubbard model
2. Bose-Hubbard Model19: Quantum Monte Carlo method studies

Comparison with Van der Waals Fluid

Detailed analysis of van der Waals fluid entropy in Appendix A reveals:

• Only the phase I to phase II transition can be qualitatively compared with van der Waals fluid
• Indicating this first-order transition belongs to the gas-liquid type

Advantages of This Work

1. Analytical Exactness: Provides completely analytical entropy expressions independent of numerical approximations
2. Systematic Completeness: Simultaneous investigation of single-particle and single-cell entropy across multiple thermodynamic planes
3. Universality Arguments: Through cross-model comparison, reveals common features of multi-occupancy systems
4. Infinite Phase Sequence: Model exhibits infinite sequence of first-order phase transitions, richer than finite-occupancy models

Conclusions and Discussion

Main Conclusions

1. Successful Derivation of Analytical Expressions: Exact analytical entropy expressions for multi-occupancy cell fluid models with Curie-Weiss interaction have been obtained, including single-particle and single-cell entropy as functions of temperature and chemical potential
2. Entropy Characteristics in Phase Transitions: Below critical temperature, entropy exhibits discontinuous jumps at first-order phase transitions, completely corresponding to infinite cascades of coexistence regions
3. Integer Density Minima: Entropy exhibits pronounced minima near integer-valued densities, directly related to discrete cell occupancy characteristics, representing a universal thermodynamic feature of multi-occupancy lattice systems
4. Cross-Model Universality: Through qualitative comparison with lattice models of different microscopic origins (Hubbard, Bose-Hubbard), confirmation that all systems allowing more than single occupancy develop entropy minima at integer densities
5. Phase Transition Type Identification: Qualitative similarity of the phase I to phase II transition with van der Waals fluid indicates this transition belongs to the gas-liquid type

Limitations

1. Classical Treatment Limitations:
   • Entropy becomes negative at high densities or small v* values
   • This is an inherent artifact of classical statistical treatment
   • Quantum corrections are needed for complete accuracy
2. Mean-Field Approximation:
   • Curie-Weiss interaction is mean-field type
   • Neglects short-range fluctuations and spatial correlations
   • Critical exponents are mean-field values
3. Parameter Dependence:
   • Results are sensitive to v* and a parameters
   • Limited number of positive entropy phases when v* is small
4. Lack of Experimental Verification:
   • As theoretical research, lacks direct comparison with experimental data
   • Experimental realization of multi-occupancy classical systems is rare

Future Directions

1. Verification with Other Multi-Occupancy Models:
   • Double-occupancy lattice gases7,8 should show entropy dip at ρ*=1.0
   • Entropy characteristics of multi-layer adsorption models11
   • Integer density minima in cluster crystal models12,13
2. Beyond Mean-Field:
   • Consideration of short-range fluctuation corrections
   • Investigation of finite-size effects
   • Renormalization group analysis
3. Quantum Extensions:
   • Development of quantum versions of cell models
   • Quantitative comparison with Bose-Hubbard models
4. Practical Applications:
   • Soft matter systems (colloids, polymers)
   • Cluster-forming interaction systems
   • Theoretical guidance for cold atom optical lattice experiments

In-Depth Evaluation

Strengths

1. Methodological Innovation

• Completely Analytical Solution: Obtaining exact analytical entropy expressions in multi-occupancy models is extremely rare, providing important theoretical benchmarks for the field
• Parametrization Technique: Elegant handling of complex implicit function relationships through introduction of z̄_max parameter
• Systematicity: Simultaneous investigation of single-particle and single-cell entropy, providing complete picture across multiple thermodynamic planes

2. Scientific Discoveries

• Universality of Integer Density Minima: The most important finding of this work, revealing fundamental characteristics of multi-occupancy models
• Cross-Disciplinary Connections: Linking classical statistical mechanics models with quantum many-body systems, demonstrating profound physical insights
• Infinite Phase Sequence: Systematic investigation of models with infinite sequences of first-order phase transitions, theoretically very challenging

3. Presentation Quality

• Clear Figures: Four panels in Figure 1 systematically present different aspects of entropy with clear color coding
• Detailed Numerics: Tables 1 and 2 provide precise critical point and coexistence line data for convenient reference in subsequent research
• Logical Rigor: Clear structure from model definition, analytical derivation, numerical results, to physical discussion

4. Theoretical Value

• Provides theoretical tools for soft matter physics
• Offers new perspective for understanding thermodynamics of cluster-forming systems
• Establishes unified understanding framework for classical and quantum multi-occupancy systems

Weaknesses

1. Depth of Physical Interpretation

• Microscopic mechanism of integer density minima not sufficiently explained
• Why does multi-occupancy lead to this universal behavior? Lacks more fundamental statistical mechanics arguments
• Physical pictures of various phases (phases I, II, III, etc.) insufficiently described

2. Limited Numerical Range

• Only first 4-5 phases shown in results
• Insufficient discussion of behavior in higher phases (phase VI and above)
• Choice of v*=5.0 lacks sufficient justification

3. Connection with Experiment

• Purely theoretical work lacking comparison with experiments
• No discussion of how to observe predicted entropy features in actual systems
• Correspondence between parameters a=1.2, v*=5.0 and actual materials unclear

4. Model Limitation Discussion

• Insufficient discussion of limitations from long-range nature of Curie-Weiss interactions
• Why three-dimensional space? (Why not one or two dimensions?)
• Applicability range of point particle approximation

5. Mathematical Rigor

• Convergence conditions of special functions K_j not discussed in detail
• Mathematical rigor of thermodynamic limit treatment relatively brief
• Uniqueness and globality of extremum conditions (11) not sufficiently established

Impact Assessment

Contribution to Field

• Theoretical Benchmark: Provides important analytical reference for multi-occupancy model research
• Universal Principle: Integer density minima may become standard characteristic for identifying and classifying multi-occupancy systems
• Cross-Disciplinary Bridge: Connects classical statistical mechanics, soft matter physics, and quantum many-body theory

Practical Value

• Moderate: As fundamental theoretical research, direct applications are limited
• May guide:
  • Understanding phase behavior of colloidal systems
  • Theoretical predictions for cold atom experiments
  • Design of novel soft matter materials

Reproducibility

• High: Complete analytical formulas and parameters provided
• Numerical calculations directly reproducible
• Recommendation for authors to release computational code to promote applications

Applicable Scenarios

Theoretical Research

1. Benchmark testing for multi-occupancy lattice models
2. Teaching case for mean-field theory
3. Research object for phase transition theory

Potential Applications

1. Soft Matter Systems:
   • Penetrable colloidal particles
   • Star polymer solutions
   • Dendritic macromolecules
2. Cold Atom Systems:
   • Bosonic gases in optical lattices
   • Multi-component atomic mixtures
3. Adsorption Phenomena:
   • Multi-layer gas adsorption
   • Surface cluster formation
4. Theoretical Method Development:
   • Testing beyond mean-field methods
   • Development of new numerical algorithms

Reference Literature (Key References)

1. Lee & Yang (1952): Statistical theory foundation for lattice gases and Ising model
2. Stanley (1971): Classical textbook on phase transitions and critical phenomena
3. Stillinger (1976): Phase transition research on Gaussian core systems
4. Kozlovskii et al. (2020, 2022): Prior work on Curie-Weiss cell model on which this work is based
5. De Leo et al. (2011): Thermodynamic research on three-dimensional Hubbard model
6. Pollet et al. (2008): Quantum Monte Carlo research on Bose-Hubbard model
7. Hill (1956): Classical statistical mechanics textbook, reference for ideal gas entropy

Overall Evaluation

This is a high-quality theoretical statistical mechanics paper making important contributions to entropy research in multi-occupancy lattice models. The paper's main value lies in: (1) providing completely analytical entropy expressions; (2) discovering and arguing the universality of integer density entropy minima; (3) establishing unified understanding of different types of multi-occupancy systems. Despite some limitations (primarily mean-field approximation and lack of experimental verification), as fundamental theoretical research, the paper demonstrates profound physical insights and rigorous mathematical treatment. This work will become an important reference for multi-occupancy model research and may inspire further investigation in soft matter and cold atom fields.

Recommendation Index: ⭐⭐⭐⭐☆ (4/5)

• Theoretical Innovation: ⭐⭐⭐⭐⭐
• Experimental Relevance: ⭐⭐⭐☆☆
• Methodological Rigor: ⭐⭐⭐⭐☆
• Presentation Clarity: ⭐⭐⭐⭐⭐
• Impact Potential: ⭐⭐⭐⭐☆