Aqua and Fisher (Phys. Rev. E 79, 011118 (2009)) exactly solved a multicomponent generalization of the Berlin–Kac spherical model — an S-species lattice gas with species-resolved interlaced sublattices encoding hard cores and one spherical constraint per species — and found that for binary neutral fluids with short-range overall-attractive interactions the near-critical equation of state acquires an unphysical linear "demagnetization" term that prevents the thermodynamic susceptibility from diverging at criticality for all fixed compositions except special submanifolds. The open problem is to construct an exactly soluble multicomponent spherical-type model (or comparable exactly soluble model in the n→∞n\to\infty class) of a neutral fluid mixture whose susceptibility diverges with the standard spherical critical singularities at every fixed composition, or to prove that such a linear term cannot be eliminated identically over the physical composition domain in exactly soluble spherical-type models with per-species spherical constraints and vectorial (n→∞n\to\infty) order parameters. The authors' own 2019 ionic companion shows the artifact disappears for binary ionic fluids, but only because electroneutrality confines those systems to the composition submanifold on which the artifact is already absent; for neutral mixtures at imposed nonzero composition the artifact stands. Companion facets arising from the source's explicit directions and acknowledged limitations — behavior for d≤2d\le 2 and for interactions violating the assumed asymmetry bounds, extensions to ternary and higher mixtures, and a principled prescription for the mixing parameters — also remain open.

Background

The spherical model of Berlin and Kac replaces the fixed-length spins of the Ising model by real spin variables subject to a single global constraint ∑isi2=N\sum_i s_i^2 = N, which makes the partition function exactly evaluable in any spatial dimension dd while retaining a genuine critical point for d>2d>2. Stanley showed that it is the n→∞n\to\infty limit of the O(n)O(n) family of vector-spin models (order-parameter dimensionality nn, with n=1,2,3n=1,2,3 corresponding to Ising, XY, and Heisenberg models), which made it a standard exactly solvable laboratory for critical phenomena — the "poor man's Ising model". Through the usual lattice-gas–spin correspondence (site occupancy ↔\leftrightarrow spin direction; chemical potential ↔\leftrightarrow magnetic field), the same mathematics describes fluids, where the magnetic susceptibility becomes the fluid compressibility.

To model a mixture of SS distinct molecular species one must distinguish the species statistically while keeping the model solvable. Aqua and Fisher (2009) achieved this by placing each species σ\sigma on its own interlaced sublattice of a common lattice — particles of different species can then never sit on the same site region, which encodes hard-core exclusion exactly — and by imposing one spherical constraint per species, enforced by its own Lagrange multiplier λσ\lambda_\sigma. The exact free energy and all pair correlation functions are then expressed through the eigenvalues of the S×SS\times S matrix of Fourier-transformed pair interactions. For binary systems (S=2S=2) with short-range, overall-attractive interactions and d>2d>2, criticality occurs on a locus Tc(m^†)T_c(\hat m^\dagger) in the space of temperature and the two thermodynamic fields, where m^†\hat m^\dagger is the mixed variable characterizing composition (essentially the density difference between the species) and m^\hat m the mixed order-parameter density. Exposing this standard structure requires linear recombinations of the species fields and densities with state-dependent mixing parameters θ\theta, θm\theta_m, θh\theta_h.

The unexpected outcome — the reason for "cautions" in the paper's title — is that the near-critical equation of state acquires a linear term n m^\mathfrak n\,\hat m: schematically, h^−h^c−n(m^,m^†) m^∼m^{ctt+cmm^2}γ\hat h - \hat h_c - \mathfrak n(\hat m,\hat m^\dagger)\,\hat m \sim \hat m\{c_t t + c_m \hat m^2\}^{\gamma} with reduced temperature t=(T−Tc)/Tct=(T-T_c)/T_c and the usual spherical exponent γ\gamma. The name is borrowed from magnetism, where the internal field seen by the spins differs from the applied field by a demagnetization correction Hint=Hext−NMH_{\rm int} = H_{\rm ext} - N M with a geometry-dependent coefficient NN; here the coefficient n\mathfrak n is generated internally by the model itself. Consequently the inverse susceptibility 1/χ=(∂h^/∂m^)t,m^†1/\chi = (\partial\hat h/\partial\hat m)_{t,\hat m^\dagger} approaches the finite limit nc(m^†)≠0\mathfrak n_c(\hat m^\dagger)\neq 0 at criticality: the susceptibility fails to diverge for every fixed composition except special submanifolds (m^†=0\hat m^\dagger = 0 in symmetric systems; two isolated loci in asymmetric ones), and the chemical potential and pressure become linear in the density through the two-phase region — both features impossible for a real fluid near criticality. A mean-field analysis of the corresponding asymmetric XY model (two-component spins, the minimal representative of the vectorial order parameter) reveals the mechanism: below TcT_c the two sublattice magnetizations cant away from the field direction, and the transverse degrees of freedom absorb the applied field. This identifies the artifact with the n→∞n\to\infty vectorial character of the order parameter combined with compositional asymmetry, suggesting it is intrinsic to the model class — but the 2009 paper does not establish that, and no modified construction avoiding the artifact is known.

The companion application to charged fluids was carried out by the authors a decade later (2019): in a binary ionic spherical model with Coulomb plus short-range nonionic interactions, electroneutrality forces the composition variable m†m^\dagger to vanish, and the artifact is absent — the susceptibility diverges normally, the universality class is the standard spherical one, and the charge correlations acquire their own new critical phenomenology (coupling of charge and density fluctuations, common divergence of charge and density correlation lengths under asymmetry, violation of the Stillinger–Lovett sum rule — the exact second-moment condition on the small-wavenumber charge structure factor that characterizes a conducting fluid). That escape route, however, is available only when a conservation law pins the composition; for neutral mixtures at imposed nonzero composition, the artifact stands exactly as found in 2009.

Problem Statement

The problem, as left open by Aqua and Fisher (2009) and unresolved by subsequent literature, has one core question and three companion facets that arise from the source's explicit directions and acknowledged limitations as part of the same objective: making exactly soluble multicomponent models physically reliable for mixture criticality. Resolving the core question is by itself sufficient to resolve the task; the companion facets (i)–(iii) are additional open directions within this objective, not separately required deliverables.

Core question. For neutral binary fluid mixtures with hard cores and short-range overall-attractive interactions in spatial dimension d>2d>2, either (a) construct an exactly soluble model in the multicomponent spherical class — a lattice gas with species-resolved interlaced sublattices and per-species spherical constraints, allowing whatever modifications of the constraints, sublattice geometry, interaction structure, or ensemble preserve exact solubility, or a comparable alternative exactly soluble model with an n→∞n\to\infty vector order parameter and per-species constraints — in which the appropriately defined thermodynamic susceptibility diverges at criticality with the standard spherical-model singularities for every fixed composition, i.e. in which the demagnetization-like linear field–density term n m^\mathfrak n\,\hat m is absent or vanishes identically over the physical composition domain; or (b) prove that such a construction is impossible — that, within the class of exactly soluble spherical-type models described in (a), the linear coupling between the ordering field and the order-parameter density at criticality cannot be eliminated identically over the physical composition domain, so that every model in the class retains nonzero composition asymmetries m^†\hat m^\dagger at which the artifact is present — and characterize exactly the compositions and submanifolds that escape the artifact.

Companion facets, scientifically linked to the core question: (i) the critical behavior of the multicomponent spherical model in the regimes excluded by the 2009 analysis — spatial dimension d≤2d\le2, where the basic lattice integrals diverge, and interactions whose intra-species asymmetry violates the "not too extreme" bounds assumed there to locate the singularities; (ii) the extension to ternary and higher mixtures, for example a neutral solvent plus positive and negative ions, and whether additional species generate or mitigate demagnetization-like effects; (iii) a principled prescription for the mixing parameters θ\theta, θm\theta_m, θh\theta_h (or alternative thermodynamic variables) that renders the standard spherical critical singularities manifest without unphysical linear terms in general asymmetric multicomponent settings.

The binary ionic case with Coulomb interactions, in which electroneutrality pins m†=0m^\dagger=0 so that no artifact arises, was analyzed by the authors themselves in 2019 and is not part of what remains open; the neutral-mixture artifact at imposed nonzero composition is the load-bearing open target.

References

  1. Jean-Noël Aqua and Michael E. Fisher, "Criticality in multicomponent spherical models: Results and cautions", Physical Review E 79, 011118 (2009). DOI: 10.1103/physreve.79.011118; arXiv:0810.5539. https://doi.org/10.1103/physreve.79.011118
  2. Jean-Noël Aqua and Michael E. Fisher, "Critical charge and density coupling in ionic spherical models", Physical Review E 100, 052145 (2019). DOI: 10.1103/physreve.100.052145. https://doi.org/10.1103/physreve.100.052145
  3. William de Castilho and S. R. Salinas, "Spherical model with Dzyaloshinskii–Moriya interactions", Journal of Statistical Mechanics: Theory and Experiment 2022, 073106 (2022). DOI: 10.1088/1742-5468/ac7e3e; arXiv:2203.09970. https://doi.org/10.1088/1742-5468/ac7e3e
  4. T. H. Berlin and M. Kac, "The Spherical Model of a Ferromagnet", Physical Review 86, 821–835 (1952). DOI: 10.1103/PhysRev.86.821. https://doi.org/10.1103/PhysRev.86.821
  5. Jean-Noël Aqua and Michael E. Fisher, "Ionic Criticality: An Exactly Soluble Model", Physical Review Letters 92, 135702 (2004). DOI: 10.1103/physrevlett.92.135702. https://doi.org/10.1103/physrevlett.92.135702
  6. Jean-Noël Aqua and Michael E. Fisher, "Charge and density fluctuations lock horns: ionic criticality with power-law forces", Journal of Physics A: Mathematical and General 37, L241–L248 (2004). DOI: 10.1088/0305-4470/37/24/l02. https://doi.org/10.1088/0305-4470/37/24/l02
  7. H. E. Stanley, "Spherical Model as the Limit of Infinite Spin Dimensionality", Physical Review 176, 718–722 (1968). DOI: 10.1103/PhysRev.176.718. https://doi.org/10.1103/PhysRev.176.718
  8. Angel Alastuey and Subir K. Das, "Algebraic infection of charge correlations at the ionic critical point", Contributions to Plasma Physics 63, e202300016 (2023). DOI: 10.1002/ctpp.202300016. https://doi.org/10.1002/ctpp.202300016