2025-11-22T03:34:16.302574

Analysis of quantities determining the critical inverse temperature in the annealed Potts model with Pareto vertex weights

Janssen

We consider in this work the crucial quantity $t_c$ that determines the critical inverse temperature $Î²_c$ in the $q$-state Potts model on sparse rank-1 random graphs where the vertices are equipped with a Pareto weight density $(Ï-1)\,w^{-Ï}\,{\cal X}_{[1,\infty)}(w)$. It is shown in \cite{ref1} that this $t_c$ is the unique positive zero of a function ${\cal K}$ that is obtained by an appropriate combination of the stationarity condition and the criticality condition for the case the external field $B$ equals 0 and that $q\geq3$ and $Ï\geq4$, see \cite{ref1}, Theorem~1.14 and Theorem ~1.21 and their proofs in \cite{ref1}, Section~7.1 and Section~7.3. From the proof of \cite{ref1}, Theorem~1.14, it is seen that ${\cal K}'$ and ${\cal K}''$ also have a unique positive zero, $t_c'$ and $t_c''$, respectively, and $t_c'=t_b$ and $t_c''=t_{\ast}$, where $t_b$ and $t_{\ast}$ are the unique positive zeros of ${\cal F}_0(t)-t\,{\cal F}_0'(t)$ and ${\cal F}_0''(t)$, respectively. Here, ${\cal F}_0(t)=E\,[W(e^{tW}-1)/(E\,[W]\,(e^{tW}+q-1))]$, and $t_c$, $t_b$ and $t_{\ast}$ play a key role in the graphical analysis of \cite{ref1}, Section~5.1 and Figure~1. Furthermore, $Î³_c=\exp(Î²_c)-1$ and $t_c$ are related according to $Î³_c=t_c/{\cal F}_0(t_c)$. We analyse $t_c$, $t_c'$ and $t_c''$ for general real $Ï\geq4$ and general real $q>2$ by an appropriate formulation of their defining equations ${\cal K}(t_c)={\cal K}'(t_c')={\cal K}''(t_c'')=0$. Thus we find, along with the inequality $0<t_c''<t_c'<t_c<\infty$, the simple upper bounds $t_c<2\,{\rm ln}(q-1)$, $t_c'<\frac32\,{\rm ln}(q-1)$, $t_c''<{\rm ln}(q-1)$, as well as certain sharpenings of these simple bounds and counterparts about the large-$q$ behaviour of $t_c$, $t_c$ and $t_c''$. We show that these bounds are sharp in the sense that they hold with equality for the limiting homogeneous case $Ï\to\infty$.

academic

Analysis of quantities determining the critical inverse temperature in the annealed Potts model with Pareto vertex weights

Basic Information

Paper ID: 2508.21409
Title: Analysis of quantities determining the critical inverse temperature in the annealed Potts model with Pareto vertex weights
Author: A.J.E.M. Janssen (Eindhoven University of Technology)
Classification: math-ph (Mathematical Physics), math.MP (Mathematical Physics), math.PR (Probability)
Publication Date: November 5, 2025 (arXiv v2)
Paper Link: https://arxiv.org/abs/2508.21409

Abstract

This paper provides an in-depth investigation of the critical quantity tc that determines the critical inverse temperature βc in the q-state Potts model with Pareto vertex weights on sparse rank-1 random graphs. The study demonstrates that tc is the unique positive zero of the function K, which is derived from the stationarity condition at zero external field B=0 and the criticality condition (q≥3, τ≥4). Through careful construction of defining equations K(tc)=K'(t'c)=K''(t''c)=0, the paper provides comprehensive analysis of tc, t'c, and t''c, establishing simple upper bounds tc<2ln(q-1), t'c<(3/2)ln(q-1), t''c<ln(q-1), along with refined versions and large-q asymptotics. The paper proves that these bounds are tight in the homogeneous limit τ→∞.

Research Background and Motivation

1. Research Problem

This paper investigates the critical inverse temperature in the annealed ferromagnetic q-state Potts model on sparse rank-1 random graphs equipped with Pareto weight density (τ-1)w^(-τ)X[1,∞)(w). The core focus is analyzing the critical quantity tc and related quantities t'c, t''c that determine the critical inverse temperature βc.

2. Problem Significance

Phase Transition Theory: The Potts model is fundamental in statistical mechanics for studying phase transitions; understanding critical temperatures is crucial for predicting system behavior
Complex Networks: Rank-1 random graphs are important mathematical models for studying real-world networks, with phase transition behavior having broad application value
Mathematical Physics: Analysis of behavior near critical points is fundamental to understanding macroscopic properties of systems

3. Limitations of Existing Methods

Previous work 1 established existence and uniqueness of tc but lacked precise bounds and asymptotic behavior for tc, t'c, t''c
For exponential weights 2, results are relatively complex and require extensive numerical computation
Systematic analysis across different parameter regions (τ and q) is lacking

4. Research Motivation

The author aims to:

Provide transparent, complete, and elegant analytical results for the Pareto weight case
Establish precise upper and lower bounds for tc, t'c, t''c
Reveal asymptotic behavior as q→∞ and q↓2
Compare finite τ and the homogeneous limit τ→∞

Core Contributions

Unified Integral Representation: Proves that K(t), K'(t), K''(t) can be expressed using a single integral D(t) (formulas 18-21), greatly simplifying analysis and numerical computation
Simple Upper Bounds: Establishes simple upper bounds for all τ≥4 and q>2:
- tc < 2ln(q-1)
- t'c < (3/2)ln(q-1)
- t''c < ln(q-1)
Refined Bounds: Provides more refined bounds:
- tc < 2(τ-2)/(τ-1)ln(q-1) (formula 24)
- t'c < T, where T is the unique positive solution of equation (27)
Large-q Asymptotics: Proves that as q→∞:
- tc/ln(q-1) → 2(τ-2)/(τ-1)
- t'c/ln(q-1) → 1
- t''c/ln(q-1) → 1
Small-q Behavior Classification: Systematically analyzes the decay behavior of t''c as q↓2, classifying into four cases based on τ (formulas 34-35):
- τ=4: exponential decay b·exp(-K₁/b)
- 4<τ<5: power-law decay K₂b^(1/(τ-4))
- τ=5: logarithmic correction K₃b/|ln b|
- τ>5: linear decay K₄b, where K₄=(τ-5)/(τ-4)
Exact Solutions in Homogeneous Limit: Proves that bounds are tight as τ→∞ and provides exact solutions for the homogeneous case

Methodology Details

Task Definition

Input:

Parameter τ≥4 (shape parameter of Pareto distribution)
Parameter q>2 (number of states in Potts model)
External field B=0

Output:

Precise bounds and asymptotic behavior of critical quantities tc, t'c, t''c
Critical inverse temperature βc = ln(1 + tc/F₀(tc))

Constraints:

F₀''(t) is first positive then negative, with unique inflection point t*>0
EW²<∞ (second moment of weights is finite)

Core Function Definitions

Basic Function F₀(t): $F_0(t) = E\left[\frac{W}{E[W]}\frac{e^{tW}-1}{e^{tW}+q-1}\right], \quad t\geq 0$

For Pareto weights, the explicit form is: $F_0(t) = 1 - \frac{q(τ-2)}{∞}\int_1^∞ \frac{w^{-τ+1}dw}{e^{tw}+q-1}$

K Function (formula 6): $K(t) = \frac{1}{E[W]}E\left[\ln\frac{e^{tW}+q-1}{q}\right] - \frac{q-1}{2q}tF_0(t) - \frac{t}{q}$

Key Technical Innovations

1. Unified Integral Representation (Section 3)

Through integration by parts and clever transformations, K, K', K'' are expressed as combinations of a single integral D(t):

$D(t) = \int_1^∞ \frac{w^{-τ+1}dw}{e^{tw}+q-1}$

Representation Forms (formulas 18-20): $K(t) = \frac{τ-2}{τ-1}\ln\frac{e^t+q-1}{q} + \left(\frac{1}{τ-1}-\frac{q+1}{2q}\right)t + \frac{(τ-2)(τ-3)}{2(τ-1)}t(q-1)D$

$K'(t) = \frac{1}{2}(q-1)\left(\frac{1}{q} - \frac{τ-2}{e^t+q-1} + (τ-2)(τ-3)D\right)$

$K''(t) = \frac{(q-1)(τ-2)}{2t}\left(\frac{te^t}{(e^t+q-1)^2} - \frac{τ-3}{e^t+q-1} + (τ-2)(τ-3)D\right)$

Innovation Significance:

Avoids complex multiple integrals
Facilitates numerical computation (Newton iteration)
Unifies the analytical framework

2. Convexity Arguments (Section 4)

Utilizes strict convexity of function f(w) = 1/((q-1)^(αw-1)+1) to establish upper bounds.

Key Lemma: For α≥1, function f(w) is strictly convex on w>1, therefore: $f(w) > f(1) + f'(1)(w-1), \quad w>1$

Through integral inequalities: $(q-1)D > \frac{f(1)}{τ-2} + \frac{f'(1)}{(τ-2)(τ-3)}$

Applications:

α=2: proves tc < 2ln(q-1)
α=3/2: proves t'c < (3/2)ln(q-1)
α=1: proves t''c < ln(q-1)

3. Newton Iteration Refinement (Section 5)

Utilizes convexity of K(t) for t≥tc to refine bounds through the first step of Newton iteration:

$t_c < T - \frac{K(T)}{K'(T)} < T$

Key Observation (formula 66): $K(2\ln(q-1)) = \frac{2\ln(q-1)}{τ-1}K'(2\ln(q-1))$

This leads to the refined bound: $t_c < 2\frac{τ-2}{τ-1}\ln(q-1)$

For t'c, a special upper bound T is constructed satisfying equation (27): $\frac{te^t}{(e^t+q-1)^2} + \frac{1}{e^t+q-1} - \frac{1}{q} = 0$

such that K'(T)/K''(T) has a simple form.

4. Asymptotic Analysis Techniques

Large-q Behavior (Section 6): Uses dominated convergence theorem and refined remainder estimates. For α∈(1, 2(τ-2)/(τ-1)):

$K(α\ln(q-1)) = \frac{1}{2}\left(α - 2\frac{τ-2}{τ-1} - \frac{α}{q} + O\left(\frac{1}{(q-1)^{α-1}}\right)\right)\ln(q-1)$

Small-q Behavior (Section 7): Utilizes that t''c is the unique positive zero of Φ(t), where: $Φ(t) = \int_t^∞ x^{-τ+3}a(x)dx, \quad a(x) = \frac{(q-1)e^x - e^{2x}}{(e^x+q-1)^3}$

The key is analyzing the balance of integrals T₁ and T₂ (formula 116).

Geometric Interpretation

tc has geometric significance (formula 13): $\int_0^{t_c} F_0(s)ds = \frac{1}{2}t_c F_0(t_c)$

This means the area under the curve y=F₀(s) equals the area under the line through the origin and (tc, F₀(tc)).

Experimental Setup

Numerical Verification Cases

This paper is primarily theoretical, but provides concrete numerical examples to verify theoretical results:

Homogeneous Case (W≡1, q=100):

F₀(t) = 1 - 100/(e^t + 99)
t''c = ln(q-1) = 4.5951...
tc = 2ln(q-1) = 9.1902...
t'c = 1.3648...ln(q-1) = 6.2716...

Numerical Examples for Pareto Case (Section 5):

q=20, τ=6: t'c=3.1829, T=4.1914
q=20, τ=11: t'c=3.7205, T=4.194
q=20, τ=18: t'c=3.9245, T=4.1914

Computational Methods

Uses Newton iteration (formula 39): $t^{(j+1)} = t^{(j)} - \frac{L(t^{(j)})}{L'(t^{(j)})}, \quad L=K \text{ or } K'$

Advantages:

Representations of K and K' require computing only the same integral D
Simpler than formula (C.10) in reference 1
Rapid convergence (guaranteed by convexity)

Experimental Results

Main Theoretical Results

1. Simple Upper Bounds (Theorem Form)

Theorem 4.1 (formula 22): For all τ≥4 and q>2, $t_c < 2\ln(q-1), \quad t'_c < \frac{3}{2}\ln(q-1), \quad t''_c < \ln(q-1)$

Proof Outline:

Utilizes convexity inequality f(w) > f(1) + f'(1)(w-1)
Sign determination of K, K', K'' at specific points
Rigorous inequality chain derivation

2. Refined Bounds

Theorem 5.1 (formula 24): For τ≥4 and q>2, $t_c < 2\frac{τ-2}{τ-1}\ln(q-1)$

Improvement Degree:

τ=4: bound is (4/3)ln(q-1), 67% tighter than 2ln(q-1)
τ→∞: bound approaches 2ln(q-1), consistent with simple bound

Theorem 5.2 (formula 28): There exists T satisfying equation (27) such that $t'_c < T < \frac{3}{2}\ln(q-1), \quad \ln(q-1) < T < \frac{3}{2}\ln(q-1)$

3. Asymptotic Limits

Theorem 6.1 (formula 31): For any τ≥4, $\lim_{q→∞}\frac{t_c}{\ln(q-1)} = 2\frac{τ-2}{τ-1}, \quad \lim_{q→∞}\frac{t'_c}{\ln(q-1)} = 1, \quad \lim_{q→∞}\frac{t''_c}{\ln(q-1)} = 1$

Physical Significance:

For large q, the leading term of tc is 2(τ-2)/(τ-1)ln(q-1)
Leading terms of t'c and t''c are both ln(q-1)
τ's influence persists in tc but diminishes in t'c, t''c

Corollary (formula 93): For γc = exp(βc)-1, $γ_c = t_c\left(1 + O\left(\frac{1}{q-1}\right)^α\right), \quad q→∞$

4. Small-q Classification Theorem

Theorem 7.1 (formulas 34-35): As b=ln(q-1)↓0, the asymptotic behavior of t''c is classified as follows:

τ Range	Asymptotic Behavior of t''c	Decay Type
τ=4	b·exp(-K₁/b)	Super-exponential decay
4<τ<5	K₂b^(1/(τ-4))	Power-law decay
τ=5	K₃b/\|ln b\|	Logarithmic-corrected linear
τ>5	((τ-5)/(τ-4))b	Linear decay

Physical Interpretation:

Larger τ brings system closer to homogeneous case with more regular decay
τ=5 is a critical transition point
Coefficient K₄=(τ-5)/(τ-4) has clear physical meaning

Exact Solutions in Homogeneous Limit

Theorem 8.1: As τ→∞,

t_{H,c} = 2ln(q-1) (exact)
t'_{H,c} = T (solution of equation 27)
t''_{H,c} = ln(q-1) (exact)

Tightness of Bounds: Inequalities for finite τ all become equalities as τ→∞, proving optimality of bounds.

Refined Analysis of T(q) (Section 9):

Small-q Expansion (formula 142): $T(q) = \frac{3}{2}\left(q-2 - \frac{1}{2}(q-2)^2 + \frac{13}{40}(q-2)^3 + O((q-2)^4)\right)$

Large-q Asymptotics (formula 144): $T(q) = \ln q + \ln(\ln B) + O\left(\frac{\ln(\ln B)}{\ln B}\right), \quad B=\frac{q}{e}$

Monotonicity: T(q) is strictly increasing on q>2, while T(q)/ln(q-1) is strictly decreasing on q>2, decreasing from 3/2 to 1.

Numerical Verification

The paper verifies theoretical predictions through concrete calculations:

Verification Case 1 (q=20, τ=6):

Theoretical prediction: t'c < T = 4.1914
Numerical computation: t'c = 3.1829 ✓
(τ-3)/(τ-2)·T = 3.1436 < t'c (indicates K' is non-convex on t'c,T)

Verification Case 2 (q=20, τ=18):

t'c = 3.9245 < T = 4.1914 ✓
K''(t) is increasing on t'c, 4.25 (convexity holds)

1. Foundational Theoretical Framework

Reference 1 (Giardinà et al., 2025):

Establishes complete theory of annealed Potts model on rank-1 inhomogeneous graphs
Proves existence and uniqueness of tc (Theorem 1.14)
Provides properties of F₀ and general conditions for phase transitions
Present paper provides in-depth analysis of key quantities in 1

2. Exponential Weight Case

Reference 2 (Janssen, 2025):

Studies exponential density exp(-w)X[0,∞)(w) case
Results are relatively complex, requiring extensive numerical computation
Pareto case in present paper yields more transparent, complete results

Comparison:

Pareto case integral D(t) is more tractable
Bound forms are simpler and more explicit
Asymptotic behavior classification is more systematic

Lehmer Inequality (Reference 3):

Used for conjecture (86): μ₃/μ₄ < μ₀/μ₁
Provides conjecture for bounds with more general weight distributions

Newton Iteration:

Classical numerical method, used here for theoretical bound refinement
Convexity guarantees convergence

Advantages of Present Work

Compared to existing work:

Completeness: Systematically analyzes all critical quantities (tc, t'c, t''c)
Precision: Provides tight upper/lower bounds and precise asymptotics
Transparency: Results have simple forms with clear physical meaning
Universality: Covers all parameter regions (τ≥4, q>2)

Conclusions and Discussion

Main Conclusions

Unified Theoretical Framework: Establishes complete analytical system based on single integral D(t)
Hierarchical Structure of Precise Bounds:
- Simple bounds: tc<2ln(q-1), t'c<(3/2)ln(q-1), t''c<ln(q-1)
- Refined bounds: tc<2(τ-2)/(τ-1)ln(q-1), t'c<T
- Asymptotic bounds: match large-q and small-q behavior
Complete Picture of Parameter Dependence:
- τ influence: persists in tc, weakens in t'c, t''c
- q influence: unified large-q behavior, classified small-q behavior
- Homogeneous limit (τ→∞) establishes bound optimality
Physical Significance: Connects to critical inverse temperature βc through γc=tc/F₀(tc)

Limitations

Parameter Restrictions:
- Requires τ≥4 to ensure F₀'' has unique inflection point
- Requires q>2 to ensure phase transition exists
- Assumes zero external field B=0
Lower Bound Issues:
- Lower bound (32) for tc remains conjecture, proven only for large q
- Systematic lower bounds for t'c and t''c are lacking
Non-convex Regions:
- K' is non-convex in some parameter ranges, affecting Newton iteration
- (τ-3)/(τ-2)·T is not always an upper bound for t'c
Numerical Aspects:
- Integral D(t) still requires numerical computation
- T(q) has no closed-form solution

Future Directions

Generalization to Arbitrary Weights:
- Verify conjecture (86) for other distributions
- Establish unified theory based on moments
Non-zero External Field:
- Study behavior for B≠0
- Analyze external field's effect on critical points
Other Graph Models:
- Extend to more general random graphs
- Study higher-rank cases
Computational Optimization:
- Develop more efficient numerical algorithms
- Find better approximations for T(q)
Physical Applications:
- Connect to phase transition phenomena in real networks
- Apply to social networks, biological networks, etc.

In-Depth Evaluation

Strengths

1. Mathematical Rigor

Complete Proofs: All main results have rigorous proofs
Technical Innovation: Convexity arguments, Newton iteration refinement, etc. are clever
Sufficient Detail: Logic from lemmas to theorems is clear

2. Depth of Results

Multi-level Bounds: From simple to refined, well-structured
Asymptotic Completeness: Both large-q and small-q behavior fully characterized
Systematic Classification: Behavior as q↓2 systematically classified by τ values

3. Clarity of Presentation

Geometric Intuition: Figure 1 clearly shows geometric meaning of tc, t'c, t''c
Simple Formulas: Unified integral representation greatly simplifies forms
Logical Structure: Progresses from basics to applications systematically

4. Theoretical Value

Tight Bounds Proven: Homogeneous limit proves bound optimality
Universal Framework: Methods generalizable to other distributions
Physical Insights: Reveals mechanisms of parameter influence on critical behavior

Weaknesses

1. Theoretical Limitations

Missing Lower Bounds: Lower bound (32) for tc not fully proven
Convexity Assumptions: Some results depend on K or K' convexity, which doesn't always hold
Parameter Restrictions: Requirement τ≥4 limits applicability

2. Computational Complexity

Integral D(t): While unified, still requires numerical integration
Equation (27): T's defining equation has no closed-form solution
Small-q Case: Constants K₁, K₂, K₃ determination requires additional computation

3. Practical Utility

Pure Theory: Lacks connection to real systems
Limited Numerical Verification: Only few parameter combinations tested
No Software Implementation: No reproducible code provided

4. Presentation Issues

Symbol Proliferation: Symbols tc, t'c, t''c, tb, t* easily confused
High Technical Density: Not accessible to non-specialists
Insufficient Figures: Only one figure, lacks more visualizations

Impact

1. Academic Contribution

Methodology: Unified integral representation and convexity arguments applicable to other models
Theoretical Completeness: Fills analytical gaps left by 1
Benchmark Results: Provides precise comparison baseline for future research

2. Practical Value

Indirect Applications: Provides theoretical foundation for phase transitions on complex networks
Algorithm Guidance: Newton iteration refinement useful for practical computation
Parameter Understanding: Helps understand τ and q's influence on system behavior

3. Reproducibility

Strengths:
- Complete proofs, verifiable step-by-step
- Numerical examples provide verification points
- Explicit formulas, easy to implement
Challenges:
- Requires strong mathematical background
- Some proof techniques difficult to reproduce
- Lacks standardized implementation

4. Follow-up Research

Direct Extensions: Verify conjectures (32) and (86)
Generalization Directions: Other weight distributions, non-zero field, higher dimensions
Application Exploration: Connect to real network data

Applicable Scenarios

1. Theoretical Research

Statistical Physics: Phase transition theory in Potts models
Random Graph Theory: Rank-1 random graph properties
Probability Theory: Large deviations theory, extreme value theory

2. Computational Applications

Critical Point Calculation: Provides efficient Newton iteration scheme
Parameter Estimation: Bounds enable quick critical temperature range estimation
Numerical Verification: Provides theoretical guarantees for numerical algorithms

Complex Networks: Phase transitions in social, biological networks
Machine Learning: Theoretical analysis of graph neural networks
Optimization: Phase transition phenomena in combinatorial optimization

4. Inapplicable Scenarios

Not Applicable to:
- τ<4 cases (F₀'' inflection point property fails)
- Strong external field (B≠0)
- Non-Pareto distributions (requires reanalysis)

References

Main references cited in this paper:

1 C. Giardinà, C. Giberti, R. van der Hofstad, A.J.E.M. Janssen, and N. Maitra, Annealed Potts models on rank-1 inhomogeneous graphs, arXiv:2502.10553v1, 2025.

Establishes foundational theoretical framework

2 A.J.E.M. Janssen, The critical temperature in the annealed Potts model with exponential vertex weights, Eurandom preprint series, 2025-08.

Comparative study for exponential weight case

3 P.S. Bullen, Handbook of means and their inequalities, Springer, 1987.

Source for Lehmer inequality

Overall Assessment: This is a high-quality mathematical physics theory paper with rigorous technical depth and systematic complete results. While primarily theoretical, it provides important mathematical foundations for understanding critical behavior in Potts models, with methodology having generalization value. Recommended future work should strengthen numerical verification and practical application connections.