--- title: Curie–Weiss Cell Model Overview url: https://www.emergentmind.com/topics/curie-weiss-cell-model type: topic --- # Curie–Weiss Cell Model Overview Searching arXiv for recent and foundational papers on Curie–Weiss cell models and related multi-group/cell-fluid formulations. A Curie–Weiss cell model is a mean-field interacting system in which the relevant microscopic variables are attached to “cells,” groups, or coarse-grained occupancy sites, while the dominant collective interaction has Curie–Weiss form: each degree of freedom couples to a global average rather than to a finite-range neighborhood. In the literature, this designation covers at least three closely related constructions: multi-group spin models with group magnetizations as order parameters, conditional or partially frozen Curie–Weiss populations, and cell-fluid models in which the volume is partitioned into cells with multiple occupancy, global attraction, and local repulsion [1806.06708], [1608.07363], [2112.00808]. Across these variants, the common mathematical structure is a reduction of the thermodynamic or dynamical problem to a small set of collective variables—typically magnetizations or mean occupancies—together with phase transitions characterized by law-of-large-numbers limits, Gaussian or non-Gaussian fluctuation regimes, and coexistence between symmetry-broken or density-distinct phases [1712.08477], [2607.01009]. ## 1. Mean-field definition and principal variants In spin formulations, the Curie–Weiss cell idea appears when the full system is decomposed into interacting populations or “cells,” each represented by its own empirical magnetization. A basic heterogeneous two-cell construction has spins \(X_1,\dots,X_{N_1}\) in group 1 and \(Y_1,\dots,Y_{N_2}\) in group 2, with \(X_i,Y_j\in\{-1,+1\}\), proportions \(\alpha_1=\lim N_1/N\), \(\alpha_2=\lim N_2/N\), and a symmetric coupling matrix \[ J=\begin{pmatrix}J_1 & J\\ J & J_2\end{pmatrix}, \qquad A:=J_1J_2-J^2>0, \] where \(J_1>0\) and \(J_2>0\) are intra-group couplings and \(J>0\) is the inter-group coupling [1806.06708]. The associated Hamiltonian is quadratic in the group magnetizations \(S_1=\sum_{i=1}^{N_1}X_i\) and \(S_2=\sum_{j=1}^{N_2}Y_j\), which is the defining Curie–Weiss feature: macroscopic observables enter through a mean-field quadratic form rather than through geometry-dependent pair distances [1806.06708]. A simpler two-group model keeps the classical homogeneous Curie–Weiss coupling \(J_{ij}=1/N\) for all spins and only distinguishes two subsets at the level of observation. In that case, the two group magnetizations \[ M_1=\frac{1}{N_1}\sum_{i=1}^{N_1}X_i, \qquad M_2=\frac{1}{N_2}\sum_{j=1}^{N_2}Y_j \] are monitored inside a single homogeneous mean-field environment [1712.08477]. This produces a “two-cell” interpretation without changing the microscopic interaction law. A different construction is the conditional Curie–Weiss model, where one part of the population is fixed and only the remaining spins fluctuate. The population is partitioned into a group fixed to \(+1\), a group fixed to \(-1\), and an undecided group of proportion \(1-s-r\), with the conditional Gibbs measure obtained by restricting the standard Curie–Weiss model to configurations respecting the fixed spins [1608.07363]. In this setting, the dynamic part is “a Curie-Weiss model on the sites in \(V_{N,t_N}\) with \(N\)-dependent coupling strength \(1-s_N-r_N\) and \(N\)-dependent external field \(h+s_N-r_N\)” [1608.07363]. In fluid and lattice-gas formulations, the cells are actual spatial boxes. The volume \(V\subset\mathbb{R}^3\) is partitioned into \(N_v\) congruent cubic cells \(\Delta_l\) of volume \(v\), each cell carries an occupation number, and the interaction combines a global Curie–Weiss attraction with a cell-local repulsion [1610.01845], [2112.00808], [2505.14456]. One version allows unbounded occupancy \(n_l=0,1,2,\dots\) [2112.00808], [2505.14456], while another imposes double occupancy \(n_l\in\{0,1,2\}\) and is isomorphic to the Blume–Capel model on a complete graph [2607.01009]. ## 2. Spin-cell formulations: group magnetizations and coupling structure For the two-group heterogeneous spin model, the empirical magnetizations are \[ m_1=\frac{1}{N_1}\sum_{i=1}^{N_1}X_i, \qquad m_2=\frac{1}{N_2}\sum_{j=1}^{N_2}Y_j, \qquad \bar m=\alpha_1m_1+\alpha_2m_2. \] The Gibbs measure is \[ P_J(A)=Z^{-1}\sum_{(X,Y)\in A}e^{-H_J(X,Y)}, \] with coupling constants absorbing inverse temperature via \(J_1=\beta \tilde J_1\), \(J_2=\beta \tilde J_2\), \(J=\beta \tilde J\) [1806.06708]. The resulting phase structure is governed by the matrix \[ H=J^{-1}-\alpha, \qquad \alpha=\begin{pmatrix}\alpha_1&0\\0&\alpha_2\end{pmatrix}, \] so the distinction between high temperature, criticality, and low temperature is encoded by the spectral properties of \(J^{-1}-\alpha\) rather than by a single scalar threshold [1806.06708]. The high-temperature regime is characterized by \(H\) positive definite, equivalently by \[ J_1<\frac{1}{\alpha_1},\qquad J_2<\frac{1}{\alpha_2},\qquad J^2<\left(\frac{1}{\alpha_1}-J_1\right)\left(\frac{1}{\alpha_2}-J_2\right). \] Under these conditions, \((m_1,m_2)\to(0,0)\) in distribution and the scaled magnetizations satisfy a bivariate central limit theorem with covariance matrix \[ C= \begin{pmatrix} \frac{1-\alpha_2J_2}{(1-\alpha_1J_1)(1-\alpha_2J_2)-\alpha_1\alpha_2J^2} & \frac{\alpha_1\alpha_2J}{(1-\alpha_1J_1)(1-\alpha_2J_2)-\alpha_1\alpha_2J^2} \\[6pt] \frac{\alpha_1\alpha_2J}{(1-\alpha_1J_1)(1-\alpha_2J_2)-\alpha_1\alpha_2J^2} & \frac{1-\alpha_1J_1}{(1-\alpha_1J_1)(1-\alpha_2J_2)-\alpha_1\alpha_2J^2} \end{pmatrix}, \] so cross-group couplings induce correlated Gaussian fluctuations even when the macroscopic magnetization vanishes [1806.06708]. At criticality, \[ J_1<\frac{1}{\alpha_1},\qquad J_2<\frac{1}{\alpha_2},\qquad J^2=\left(\frac{1}{\alpha_1}-J_1\right)\left(\frac{1}{\alpha_2}-J_2\right), \] the law of large numbers still yields \((m_1,m_2)\to(0,0)\), but the Gaussian central limit theorem fails and the correct scaling is \(N^{3/4}\), with a non-Gaussian limit determined by quartic-order fluctuations [1806.06708]. In the low-temperature regime, the magnetization vector converges in distribution to \[ \frac12\delta_{m^*}+\frac12\delta_{-m^*}, \] where \(m^*=(m_1^*,m_2^*)\) is the positive solution of the mean-field equations \[ m_1=\tanh\big(J_1\alpha_1m_1-J\alpha_2m_2\big), \qquad m_2=\tanh\big(J_2\alpha_2m_2-J\alpha_1m_1\big), \] so the two cells exhibit spontaneous symmetry breaking into two symmetry-related ferromagnetic phases [1806.06708]. The homogeneous two-group model yields a more restrictive outcome. Since all spins share the same global Curie–Weiss interaction, the two group magnetizations converge jointly to \[ \tfrac12\delta_{(-m(\beta),-m(\beta))}+\tfrac12\delta_{(m(\beta),m(\beta))}, \] where \(m(\beta)\) solves \(m(\beta)=\tanh(\beta m(\beta))\) [1712.08477]. This excludes distinct-sign limiting phases for the two groups and shows that any macroscopic subsets are locked to the same global phase in the homogeneous model [1712.08477]. ## 3. Conditional and exchangeable formulations The conditional Curie–Weiss model introduces a frozen-core architecture. The population is partitioned into \(V_{N,s_N}\) fixed to \(+1\), \(V_{N,r_N}\) fixed to \(-1\), and \(V_{N,t_N}\) of undecided spins, with \(s_N\to s\), \(r_N\to r\), and \(t_N\to 1-s-r\) [1608.07363]. The specific magnetization satisfies \[ m(\beta,s,r,h)=s-r+(1-s-r)\,z(\beta,s,r,h), \] where \(z\) minimizes \[ i_{\beta,s,r,h}(z) = -\frac12\beta(1-s-r)z^2 -\beta(s-r+h)z +\frac{1-z}{2}\log(1-z) +\frac{1+z}{2}\log(1+z) \] on \([-1,1]\) [1608.07363]. The threshold \[ \beta^*=(1-s-r)^{-1} \] marks a first-order phase transition in magnetization at \(h=r-s\): for \(\beta\le \beta^*\), the magnetization is continuous at that field, whereas for \(\beta>\beta^*\) it jumps discontinuously [1608.07363]. The paper explicitly distinguishes a discontinuous jump in the order parameter from a genuine change of phase, because the sign of the total magnetization can remain fixed when one frozen opinion has a sufficiently large majority [1608.07363]. A different probabilistic viewpoint treats the Curie–Weiss model through exchangeability. The Gibbs law is represented as a De Finetti mixture \[ \mathbb{P}_n^{(\beta)}(X_1,\dots,X_n)=\int_{[0,1]}\mathbb{P}_p^{\otimes n}(X_1,\dots,X_n)\,\nu_{n,\beta}(dp), \] so conditional on the latent parameter \(p\), the spins are i.i.d. Bernoulli and the dependence is carried by the random mixing measure \(\nu_{n,\beta}\) [2305.06872]. In that framework, the magnetization admits the surrogate \[ M_n^{\mathrm{sur}(\beta)}=\sqrt{n}\,G\sqrt{1-T_{n,\beta}^2}+nT_{n,\beta}, \] with \(G\sim N(0,1)\) independent of the De Finetti random variable \(T_{n,\beta}\) [2305.06872]. This suggests that the phase transition can be interpreted as a competition between “the randomness of the i.i.d. spin variables given the latent environment” and “the randomness of the latent environment \(T_{n,\beta}\) itself” [2305.06872]. That decoupling was sharpened through almost-sure Laplace inversion. Using the De Finetti representation, the uniform coupling of Bernoulli variables, and the inverse Laplace transform on a complex line, the magnetization is expressed as a functional of an i.i.d. field and an independent randomisation field; the limiting Gaussian objects are identified as a Gaussian analytic function, a Brownian bridge, and a Brownian sheet [2507.17459]. The same analysis reproduces the subcritical Gaussian regime, the \(N^{3/4}\) critical scaling with non-Gaussian limit, and supercritical fluctuations around \(t_\beta\), where \(t_\beta\) solves \[ t_\beta=\tanh(\beta t_\beta) \] [2507.17459]. A plausible implication is that the “cell” interpretation of Curie–Weiss models extends naturally to a hierarchy of probabilistic fields: microscopic i.i.d. cells, a global randomisation field, and Gaussian functional limits governing fluctuations. ## 4. Cell-fluid models with multiple occupancy In continuum cell-fluid formulations, the volume \(V\) is partitioned into \(N_v\) cells \(\Delta_l\) of equal volume \(v\), and the occupation number of cell \(l\) is \[ \rho_l(\eta)=\sum_{x\in\eta}I_{\Delta_l}(x). \] The pair potential is \[ \Psi_{N_v}(x,y) = -\frac{g_a}{N_v} + g_r\sum_{l=1}^{N_v}I_{\Delta_l}(x)I_{\Delta_l}(y), \] with \(g_r>g_a>0\) or equivalently \(f=g_r/g_a\ge 1\) for stability [2505.14456], [2112.00808]. The inter-cell attraction is uniform and of Curie–Weiss type, while the intra-cell repulsion penalizes multiple occupancy [2112.00808]. In the exactly solved grand-canonical reduction, the partition function takes the one-dimensional form \[ \Xi = \left(\frac{N_v}{2\pi p_a}\right)^{1/2} \int dz\,\exp\!\left[N_vE_0(z,\mu)\right], \] where \[ E_0(z,\mu) = -\frac{1}{2p_a}(z-\beta\mu-\ln v^*)^2+\ln K_0(z), \] \[ K_m(z) = \sum_{n=0}^{\infty} \frac{1}{n!}n^m\exp\!\left(zn-\frac{fp_a}{2}n^2\right), \qquad p_a=\beta g_a,\quad p_r=\beta g_r \] [2505.14456]. The saddle-point equation gives \[ \beta\mu=\bar z-p_a\frac{K_1(\bar z)}{K_0(\bar z)}-\ln v^*, \qquad \bar n=\frac{K_1(\bar z)}{K_0(\bar z)}, \] so the reduced density is the average number of particles per cell [2505.14456]. The pressure is \[ \tilde P = (\tau+1)\left[\ln K_0(\bar z)-\frac{p_a}{2}\left(\frac{K_1(\bar z)}{K_0(\bar z)}\right)^2\right] \] in the notation of that work [2505.14456]. An earlier formulation on \(\mathbb{R}^d\) writes the occupation-space partition function as \[ \Xi_N(p,\mu) = \sum_{\varrho\in\mathbb{N}_0^N} \exp\left[ \frac{p}{2N}\left(\sum_{\ell=1}^N\varrho_\ell\right)^2 \right] \prod_{\ell=1}^N \pi(\varrho_\ell,\mu), \qquad \pi(n,\mu)=\frac{\upsilon^n}{n!}\exp\left(\mu n-\frac{ap}{2}n^2\right), \] with \(p=\beta J_1\), \(a=J_2/J_1>1\), and \(\upsilon=c^d\) [1610.01845]. There the thermodynamic phase is a product measure over cell occupation numbers determined by a Curie–Weiss self-consistency equation for the scalar parameter \(y\), \[ -\frac{y}{p}+\frac{K_1(y,p,\mu)}{K(y,p,\mu)}=0, \] and the pressure in the thermodynamic limit is \(\upsilon^{-1}E(\bar y,p,\mu)\), where \(E(y,p,\mu)=-y^2/(2p)+\ln K(y,p,\mu)\) [1610.01845]. For sufficiently small attraction there is a single phase, while for large enough \(p\) the model exhibits phase-coexistence points with two thermodynamic phases of distinct densities [1610.01845]. A double-occupancy version restricts \(n_l\in\{0,1,2\}\) and isomorphicly maps to the Blume–Capel model on a complete graph via \(n_l=S_l+1\), \(S_l\in\{-1,0,+1\}\) [2607.01009]. Its Hamiltonian contains a Curie–Weiss attraction term \(-\frac{J}{2\mathscr N}\sum_{l\ne m}n_ln_m\), a local double-occupancy repulsion \(\Delta\sum_l n_l(n_l-1)\), and an effective linear term, with the grand partition function reduced by a Hubbard–Stratonovich transformation to a saddle-point problem over a single auxiliary variable [2607.01009]. ## 5. Phase structure, criticality, and coexistence A characteristic feature of the unbounded-occupancy Curie–Weiss cell fluid is a sequence of first-order phase transitions below the first critical temperature. As the chemical potential or density increases, the system passes through phases I, II, III, and IV with density ranges approximately \([0,0.5]\), \([0.5,1.5]\), \([1.5,2.5]\), and \([2.5,3.5]\), and the transitions I–II, II–III, and III–IV occur near \(\bar n\approx 0.5\), \(1.5\), and \(2.5\), respectively [2505.14456]. The coexistence points are obtained microscopically from \[ E_0(z_1)-E_0(z_2)=0, \qquad \mu(z_1)-\mu(z_2)=0, \] which is the mean-field analogue of Maxwell construction but derived directly from the saddle-point grand potential [2505.14456]. The explicit equation of state derived for the same class of models is \[ Pv = k_BT\ln K_0(\eta)-\frac{g_av^2}{2}\eta^2, \] with \[ \eta=\frac{1}{v}\frac{K_1(z)}{K_0(z)}, \] where \(z\) is the saddle variable and \(\eta\) is the density [2112.00808]. In that analysis, critical points are defined microscopically by \[ \frac{d\mu}{dz}=0, \qquad \frac{d^2\mu}{dz^2}=0, \] which yields a sequence of critical temperatures \[ k_BT_c^{(n)}=\frac{g_a}{p_{ac}^{(n)}}, \qquad T_c^{(n+1)}>T_c^{(n)} \] [2112.00808]. For \(f=1.2\) and \(v=1\), the first three critical densities are approximately \(0.5139\), \(1.5056\), and \(2.5030\), demonstrating multiple critical points associated with successive density transitions [2112.00808]. The double-occupancy model displays a different but equally rich phase diagram. With \(a=\Delta/J\), it shows one critical point for \(a\le a_{\rm tc}=\frac12\ln2\), two critical points for \(a>a_{\rm tc}\), a tricritical point at \[ a_{\rm tc}=\frac12\ln2, \qquad T^*_{\rm tc}=\frac13, \qquad \rho^*_{\rm tc}=1, \] and a line of triple points for \(a_{\rm tc}