pith. sign in

arxiv: 2606.22123 · v1 · pith:BKR3CDV3new · submitted 2026-06-20 · ❄️ cond-mat.other

Enhanced approach to calculation of cluster integrals for lattice models of matter

M. V. Ushcats , S. Yu. Ushcats This is my paper

Pith reviewed 2026-06-26 10:58 UTC · model grok-4.3

classification ❄️ cond-mat.other
keywords cluster integralsMayer expansionlattice modelsnumerical integrationsymmetry reductionhigh-order calculations
0
0 comments X

The pith

Optimizations enable new high-order cluster integral data for lattice models

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

The paper develops an enhanced method for computing Mayer's cluster integrals in lattice models of matter. It introduces two optimizations: simplifying the integrand at each integration point and reducing the number of points by eliminating physically identical configurations via symmetry. These changes make higher-order integrals computable for various two- and three-dimensional lattices. As a result, the authors report new numerical data on these integrals for multiple models.

Core claim

By simplifying integrand calculations and eliminating physically identical configurations through symmetry, the approach yields new data on high-order cluster integrals for multiple 2D and 3D lattice models.

What carries the argument

Symmetry-based reduction of integration points combined with simplified integrand evaluation in Mayer cluster integral calculations

If this is right

  • Higher-order cluster integrals become available for more lattice models.
  • Mayer series expansions can reach higher densities or more complex lattice systems.
  • Computational cost drops, allowing still larger integral orders to be tackled.

Where Pith is reading between the lines

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

  • The symmetry technique may transfer to integrals over other discrete or discretized spaces.
  • The new numerical values could serve as benchmarks for approximate theories of lattice gases.

Load-bearing premise

Eliminating configurations based on physical identity correctly captures all unique contributions without numerical error or bias.

What would settle it

Computing a known low-order cluster integral with the reduced point set and comparing it to the result from the full set or an exact analytic value.

Figures

Figures reproduced from arXiv: 2606.22123 by M. V. Ushcats, S. Yu. Ushcats.

Figure 1
Figure 1. Figure 1: (Colour online) One of the many connected diagrams representing the 𝑏12 integrand (a) that, in turn, can be considered as a combination of six 1st-order irreducible integrands (b), one 2nd-order irreducible integrand (c), and one of the parts (d) of the 3rd-order irreducible integrand. 23503-2 [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: (Colour online) Integration volume per particle for 𝑛 = 3 in the 2D Lee–Yang lattice model. Numbers correspond to coordination spheres. the total number of integration points reaches [2𝑛 (𝑛 − 1) + 1] 𝑛−1 . At high orders (see the 1st and 2nd columns of table 1), this number makes it almost impossible to perform any real computations over all the integration points and, for other lattice models, the number … view at source ↗
Figure 3
Figure 3. Figure 3: (Colour online) Magnetization curves of the 2D square Ising model with nearest-neighbor interactions at 𝑇 = 0.8𝑇C (𝑇C is the Curie point). The circle corresponds to the exact spontaneous magnetization point of this Lee–Yang model [3, 5] at the same 𝑇. For each curve, Mayer’s expansion includes 10000 reducible cluster integrals calculated by using the scaling technique [29, 31, 32] on the basis of some limi… view at source ↗
read the original abstract

The study is devoted to enhancing the existing techniques of calculating Mayer's expansion cluster integrals for lattice models of matter. Two important optimizations are proposed: simplifying the calculation of the integrand at each integration point and reducing the number of such integration points due to eliminating physically identical configurations. Based on those optimizations, new data on high-order cluster integrals are obtained for a number of 2D and 3D lattice models.

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 proposes two optimizations to existing methods for computing Mayer cluster integrals on lattices: (1) simplification of the integrand evaluation at each integration point and (2) reduction of the number of integration points by eliminating physically identical configurations via symmetry. These changes are used to generate new numerical values for high-order cluster integrals on several 2D and 3D lattice models.

Significance. If the symmetry reduction is shown to be free of omissions or overcounting and the new integrals are validated, the work would supply useful high-order data for lattice models that could support improved equations of state or phase-transition studies. The optimizations themselves are standard in principle but their concrete implementation for high orders on 3D lattices would be the main contribution.

major comments (2)
  1. [Method description (abstract and main text)] The symmetry elimination of equivalent configurations is load-bearing for all claimed new high-order results, yet the manuscript provides no indication that the reduced integrals were cross-checked against independently published low-order values on the same lattices (e.g., orders 2–4 on square or simple-cubic lattices). Without such validation the numerical accuracy of the new data cannot be confirmed.
  2. [Results section] No error analysis, convergence tests, or comparison with existing tabulated cluster integrals appears in the provided description; this is required to establish that the integrand simplification preserves exactness while the symmetry reduction assigns correct multiplicities.
minor comments (2)
  1. [Abstract] The abstract states that 'new data' are obtained but does not specify the lattice types, the highest order reached, or the numerical precision; these details should be added for clarity.
  2. [Methods] Notation for the integrand and the symmetry group operations should be defined explicitly in the methods section to allow independent reproduction.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We agree that explicit validation of the symmetry reduction and integrand simplification is necessary to support the new high-order results. Below we address each major comment and describe the revisions that will be made.

read point-by-point responses

  1. Referee: [Method description (abstract and main text)] The symmetry elimination of equivalent configurations is load-bearing for all claimed new high-order results, yet the manuscript provides no indication that the reduced integrals were cross-checked against independently published low-order values on the same lattices (e.g., orders 2–4 on square or simple-cubic lattices). Without such validation the numerical accuracy of the new data cannot be confirmed.

    Authors: We agree that cross-validation against independently published low-order integrals is required to confirm the correctness of the symmetry factors and the absence of omissions or overcounting. Although the optimizations were verified during code development, the manuscript does not present these checks. In the revised version we will add a new subsection (or table) in the Results section that reports our computed values for orders 2–4 on the square and simple-cubic lattices together with the corresponding literature values, thereby demonstrating that the symmetry reduction assigns correct multiplicities and that the integrand simplification preserves exactness. revision: yes

  2. Referee: [Results section] No error analysis, convergence tests, or comparison with existing tabulated cluster integrals appears in the provided description; this is required to establish that the integrand simplification preserves exactness while the symmetry reduction assigns correct multiplicities.

    Authors: We acknowledge that the current manuscript lacks a dedicated error analysis, convergence tests, and systematic comparisons with tabulated integrals. The revised manuscript will incorporate an expanded Results section that includes (i) a quantitative error analysis of the numerical integration, (ii) convergence tests with respect to the number of integration points, and (iii) direct comparisons with existing tabulated cluster integrals for the lower orders already mentioned. These additions will establish that both proposed optimizations maintain the required numerical accuracy. revision: yes

Circularity Check

0 steps flagged

No circularity in optimization-based computation of cluster integrals

full rationale

The paper describes two computational optimizations (integrand simplification and symmetry-based elimination of equivalent configurations) to obtain new high-order cluster integral values for lattice models. No equations, fitted parameters, or self-citations are shown that would reduce any claimed result to its own inputs by construction. The symmetry reduction is presented as a methodological improvement without evidence of self-definition, renaming of known results, or load-bearing reliance on prior author work that itself lacks independent verification. The derivation chain for the new numerical data remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract supplies no information on free parameters, background axioms, or new postulated entities.

pith-pipeline@v0.9.1-grok · 5587 in / 871 out tokens · 13243 ms · 2026-06-26T10:58:46.333593+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

42 extracted references · 36 canonical work pages

  1. [1]

    E., Hemmer P

    Kac M., Uhlenbeck G. E., Hemmer P. C., J. Math. Phys., 1963,4, No. 2, 216–228, doi:10.1063/1.1703946

    work page doi:10.1063/1.1703946 1963

  2. [2]

    L., Penrose O., J

    Lebowitz J. L., Penrose O., J. Math. Phys., 1966,7, No. 1, 98–113, doi:10.1063/1.1704821

    work page doi:10.1063/1.1704821 1966

  3. [3]

    Crystal Statistics

    Onsager L., Phys. Rev., 1944,65, 117–149, doi:10.1103/PhysRev.65.117

    work page doi:10.1103/physrev.65.117 1944

  4. [4]

    Beitrag zur Theorie des Ferromagnetismus , year =

    Ising E., Z. Phys., 1925,31, 253–258, doi:10.1007/BF02980577

    work page doi:10.1007/bf02980577 1925

  5. [5]

    D., Yang C

    Lee T. D., Yang C. N., Phys. Rev., 1952,87, 410–419, doi:10.1103/PhysRev.87.410

    work page doi:10.1103/physrev.87.410 1952

  6. [6]

    E., Mayer M

    Mayer J. E., Mayer M. G., Statistical Mechanics, 2nd edn., John Wiley, New York, 1977

    1977

  7. [7]

    V., Bulavin L

    Ushcats M. V., Bulavin L. A., Sysoev V. M., Bardik V. Y., Alekseev A. N., J. Mol. Liq., 2016,224, 694–712, doi:10.1016/j.molliq.2016.09.100

    work page doi:10.1016/j.molliq.2016.09.100 2016

  8. [8]

    M., Physica A, 2015,419, 675–680, doi:10.1016/j.physa.2014.10.053

    Bannur V. M., Physica A, 2015,419, 675–680, doi:10.1016/j.physa.2014.10.053

    work page doi:10.1016/j.physa.2014.10.053 2015

  9. [9]

    V., Ushcats S

    Ushcats M. V., Ushcats S. Y., Bulavin L. A., Sysoev V. M., Phys. Rev. E, 2018,98, 032135, doi:10.1103/PhysRevE.98.032135

    work page doi:10.1103/physreve.98.032135 2018

  10. [10]

    L., Penrose O., J

    Lebowitz J. L., Penrose O., J. Math. Phys., 1964,5, No. 7, 841–847, doi:10.1063/1.1704186

    work page doi:10.1063/1.1704186 1964

  11. [11]

    V., Bulavin L

    Ushcats M. V., Bulavin L. A., Sysoev V. M., Ushcats S. Y., Phys. Rev. E, 2017,96, 062115, doi:10.1103/PhysRevE.96.062115

    work page doi:10.1103/physreve.96.062115 2017

  12. [12]

    UshcatsM.V.,BulavinL.A.,UshcatsS.Y.,Phys.Rev.E,2018,98,042127,doi:10.1103/PhysRevE.98.042127

    work page doi:10.1103/physreve.98.042127 2018

  13. [13]

    V., Bulavin L

    Ushcats M. V., Bulavin L. A., Sysoev V. M., Ushcats S. Y., Ukr. J. Phys., 2017,62, No. 6, 533–538, doi:10.15407/ujpe62.06.0533

    work page doi:10.15407/ujpe62.06.0533 2017

  14. [14]

    3, 31, doi:10.1007/s12043-018-1604-3

    Ushcats S., Ushcats M., Bulavin L., Svechnikova O., Mykheliev I., Pramana, 2018,91, No. 3, 31, doi:10.1007/s12043-018-1604-3

    work page doi:10.1007/s12043-018-1604-3 2018

  15. [15]

    V., Bulavin L

    Ushcats M. V., Bulavin L. A., Phys. Rev. E, 2020,101, 062128, doi:10.1103/PhysRevE.101.062128

    work page doi:10.1103/physreve.101.062128 2020

  16. [16]

    V., Bulavin L

    Ushcats M. V., Bulavin L. A., Ushcats S. Y., Markina L. M., Phys. Rev. E, 2020,102, 042130, doi:10.1103/PhysRevE.102.042130

    work page doi:10.1103/physreve.102.042130 2020

  17. [17]

    P., Udayanandan K

    Suresh T. P., Udayanandan K. M., Turk. J. Phys., 2018,42, No. 6, 668–674, doi:10.3906/fiz-1806-19

    work page doi:10.3906/fiz-1806-19 2018

  18. [18]

    P., Udayanandan K

    Suresh T. P., Udayanandan K. M., Pramana, 2020,94, No. 1, 26, doi:10.1007/s12043-019-1908-y

    work page doi:10.1007/s12043-019-1908-y 2020

  19. [19]

    Huang K., Statistical Mechanics, John Wiley, New York, 1963

    1963

  20. [20]

    L., Statistical Mechanics: Principles and Selected Applications, McGraw-Hill, New York, 1956

    Hill T. L., Statistical Mechanics: Principles and Selected Applications, McGraw-Hill, New York, 1956

    1956

  21. [21]

    Isihara A., Statistical Physics, Academic Press, 1971

    1971

  22. [22]

    K., Statistical Mechanics, Butterworth-Heinemann, Oxford, 1997

    Pathria R. K., Statistical Mechanics, Butterworth-Heinemann, Oxford, 1997

    1997

  23. [23]

    Gallavotti G., Statistical Mechanics, Springer, Berlin, Heidelberg, 1999, doi:10.1007/978-3-662-03952-6

    work page doi:10.1007/978-3-662-03952-6 1999

  24. [24]

    Baxter R., Exactly Solved Models in Statistical Mechanics, Dover Publications, 2007

    2007

  25. [25]

    Gallavotti G., Miracle-Sole S., Commun. Math. Phys., 1968,7, No. 4, 274–288, doi:10.1007/BF01646661

    work page doi:10.1007/bf01646661 1968

  26. [26]

    Matter Phys., 2011,14, 43002, doi:10.5488/CMP.14.43002

    Kozlovskii M., Romanik R., Condens. Matter Phys., 2011,14, 43002, doi:10.5488/CMP.14.43002

    work page doi:10.5488/cmp.14.43002 2011

  27. [27]

    Matter Phys., 2010,13, 43004, doi:10.5488/CMP.13.43004

    Kozlovskii M., Romanik R., Condens. Matter Phys., 2010,13, 43004, doi:10.5488/CMP.13.43004

    work page doi:10.5488/cmp.13.43004 2010

  28. [28]

    V., Bulavin L

    Ushcats M. V., Bulavin L. A., Sysoev V. M., Ushcats S. J., Phys. Rev. E, 2016,94, 012143, doi:10.1103/PhysRevE.94.012143

    work page doi:10.1103/physreve.94.012143 2016

  29. [29]

    Ushcats S., Ushcats M., Sysoev V., Gavryushenko D., Ukr. J. Phys., 2018,63, No. 12, 1066, doi:10.15407/ujpe63.12.1066

    work page doi:10.15407/ujpe63.12.1066 2018

  30. [30]

    Ushcats M., Bulavin L., Ushcats S., Lazarenko M., Labartkava A., Physica A, 2022,598, 127307, doi:10.1016/j.physa.2022.127307

    work page doi:10.1016/j.physa.2022.127307 2022

  31. [31]

    Ushcats M., Ushcats S., Kondratieva A., Koval S., Physica A, 2024,649, 129957, doi:10.1016/j.physa.2024.129957

    work page doi:10.1016/j.physa.2024.129957 2024

  32. [32]

    1, 13501, doi:10.5488/cmp.28.13501

    UshcatsM.,BulavinL.,UshcatsS.,BuruninaZ.,MaiborodaO.,RomanchukN.,ShapovalN.,Condens.Matter Phys., 2025,28, No. 1, 13501, doi:10.5488/cmp.28.13501

    work page doi:10.5488/cmp.28.13501 2025

  33. [33]

    J., Phys

    Wheatley R. J., Phys. Rev. Lett., 2013,110, 200601, doi:10.1103/PhysRevLett.110.200601

    work page doi:10.1103/physrevlett.110.200601 2013

  34. [34]

    K., Kofke D

    Singh J. K., Kofke D. A., Phys. Rev. Lett., 2004,92, 220601, doi:10.1103/PhysRevLett.92.220601

    work page doi:10.1103/physrevlett.92.220601 2004

  35. [35]

    J., Kofke D

    Schultz A. J., Kofke D. A., Mol. Phys., 2009,107, No. 21, 2309–2318, doi:10.1080/00268970903267053

    work page doi:10.1080/00268970903267053 2009

  36. [36]

    Ushcats M. V., Ukr. J. Phys., 2014,59, No. 2, 172–178, doi:10.15407/ujpe59.02.0172

    work page doi:10.15407/ujpe59.02.0172 2014

  37. [37]

    Ushcats M. V., Ukr. J. Phys., 2014,59, No. 7, 737–742, doi:10.15407/ujpe59.07.0737

    work page doi:10.15407/ujpe59.07.0737 2014

  38. [38]

    V., Ushcats S

    Ushcats M. V., Ushcats S. J., Mochalov A. A., Ukr. J. Phys., 2016,61, No. 2, 160–167, doi:10.15407/ujpe61.02.0160

    work page doi:10.15407/ujpe61.02.0160 2016

  39. [39]

    J., Chaudhary V., Kofke D

    Feng C., Schultz A. J., Chaudhary V., Kofke D. A., J. Chem. Phys., 2015,143, No. 4, 044504, doi:10.1063/1.4927339

    work page doi:10.1063/1.4927339 2015

  40. [40]

    TrokhymchukA.,MelnykR.,NezbedaI.,Condens.MatterPhys.,2015,18,13501,doi:10.5488/CMP.18.13501

    work page doi:10.5488/cmp.18.13501 2015

  41. [41]

    , year = 2012, month = oct, journal =

    Schultz A. J., Barlow N. S., Chaudhary V., Kofke D. A., Mol. Phys., 2013,111, No. 4, 535–543, doi:10.1080/00268976.2012.730642

    work page doi:10.1080/00268976.2012.730642 2013

  42. [42]

    Kofke D. A., J. Phys. Chem. B, 2023,127, No. 16, 3690–3700, doi:10.1021/acs.jpcb.3c00807. 23503-10 Enhanced approach to calculation of cluster integrals for lattice models of matter Вдосконалений пiдхiд до розрахункiв групових iнтегралiв для ґраткових моделей речовини М. В. Ушкац, С. Ю. Ушкац Нацiональний унiверситет кораблебудування iм. адмiрала Макарова...

    work page doi:10.1021/acs.jpcb.3c00807 2023