The ultimate independence ratio of a graph $G$ is defined as $\mathscr{I}(G) = \lim_{k\rightarrow\infty } \frac{α(G^{\Box k})}{|V(G)|^k},$ where $α(G^{\Box k})$ is the independence number of the Cartesian product of $k$ copies of $G$. For all graphs $G$, Hahn, Hell, and Poljak (1995) proved that $\frac{1}{Ï(G)} \leq \mathscr{I}(G) \leq \frac{1}{Ï(G)}$ where $Ï(G)$ is the chromatic number, and $Ï(G)$ is the clique number of $G$. So all graphs $G$ with $Ï(G) = Ï(G)$ satisfy $\mathscr{I}(G) = \frac{1}{Ï(G)} = \frac{1}{Ï(G)}$. A construction of Zhu demonstrates that there exists a graph $G$ with $\frac{1}{Ï(G)} < \mathscr{I}(G) < \frac{1}{Ï(G)}$, so neither equality holds in general. In response, Hahn, Hell, and Poljak conjectured that all wheel graphs $W_n$ satisfy $\mathscr{I}(W_n) = \frac{1}{Ï(W_n)}$.
For even wheels $W_{2t}$ this follows from the fact $Ï(W_{2t}) = Ï(W_{2t}) = 3$. Odd wheels of length at least $5$ present a more challenging case, since $Ï(W_{2t+1}) = 4$ and $Ï(W_{2t+1}) = 3$. First, we prove that odd wheels of length at least $7$ satisfy $\mathscr{I}(W_{2t+1})\leq \frac{4t^2+6t}{3(2t+2)^2}<\frac{1}{3}$, which provides the best upper bound for large odd wheels. Next, we prove that $\mathscr{I}(W_5) \leq \frac{1019}{3888}$, improving an upper bound of Hahn, Hell, and Poljak that $\mathscr{I}(W_5) \leq \frac{11}{41}$. Our proofs combine counting arguments, recursive bounds on $α(W^{\Box k}_{2t+1})$, and computer-assisted calculation in the $W_5$ case.
Paper ID : 2511.18747Title : Improved Bounds for the Ultimate Independence Ratio of Odd WheelsAuthors : Alexander Clow, Hitesh Kumar, Shivaramakrishna Pragada (Simon Fraser University)Classification : math.CO (Combinatorics), math.OC (Optimization and Control)Submission Date : Submitted to arXiv on November 24, 2025Paper Link : https://arxiv.org/abs/2511.18747 This paper investigates the ultimate independence ratio of graphs. For a graph G G G , the ultimate independence ratio is defined as I ( G ) = lim k → ∞ α ( G □ k ) ∣ V ( G ) ∣ k \mathscr{I}(G) = \lim_{k\rightarrow\infty} \frac{\alpha(G^{\Box k})}{|V(G)|^k} I ( G ) = lim k → ∞ ∣ V ( G ) ∣ k α ( G □ k ) , where α ( G □ k ) \alpha(G^{\Box k}) α ( G □ k ) is the independence number of the k k k -fold Cartesian product of G G G . Hahn, Hell, and Poljak (1995) proved that 1 χ ( G ) ≤ I ( G ) ≤ 1 ω ( G ) \frac{1}{\chi(G)} \leq \mathscr{I}(G) \leq \frac{1}{\omega(G)} χ ( G ) 1 ≤ I ( G ) ≤ ω ( G ) 1 and conjectured that all wheel graphs W n W_n W n satisfy I ( W n ) = 1 χ ( W n ) \mathscr{I}(W_n) = \frac{1}{\chi(W_n)} I ( W n ) = χ ( W n ) 1 . This paper makes significant progress on this 30-year-old unsolved conjecture: it is proved that for odd wheels with t ≥ 3 t \geq 3 t ≥ 3 , I ( W 2 t + 1 ) ≤ 4 t 2 + 6 t 3 ( 2 t + 2 ) 2 < 1 3 \mathscr{I}(W_{2t+1}) \leq \frac{4t^2+6t}{3(2t+2)^2} < \frac{1}{3} I ( W 2 t + 1 ) ≤ 3 ( 2 t + 2 ) 2 4 t 2 + 6 t < 3 1 ; for the 5-wheel, the upper bound is improved to I ( W 5 ) ≤ 1019 3888 ≈ 0.262 \mathscr{I}(W_5) \leq \frac{1019}{3888} \approx 0.262 I ( W 5 ) ≤ 3888 1019 ≈ 0.262 (previous bound was 11 41 ≈ 0.268 \frac{11}{41} \approx 0.268 41 11 ≈ 0.268 ).
Definition and Significance of Ultimate Independence Ratio :The ultimate independence ratio characterizes the asymptotic behavior of maximum independent sets in Cartesian powers of graphs Closely related to Shannon capacity and graph homomorphism theory Hell, Yu, and Zhou (1994) proved that the sequence { i ( G □ k ) } \{i(G^{\Box k})\} { i ( G □ k )} is monotone decreasing and convergent Basic Theoretical Bounds :Hahn, Hell, and Poljak proved: 1 χ ( G ) ≤ I ( G ) ≤ 1 χ f ( G ) ≤ 1 ω ( G ) \frac{1}{\chi(G)} \leq \mathscr{I}(G) \leq \frac{1}{\chi_f(G)} \leq \frac{1}{\omega(G)} χ ( G ) 1 ≤ I ( G ) ≤ χ f ( G ) 1 ≤ ω ( G ) 1 For weakly perfect graphs (χ = ω \chi = \omega χ = ω ), the ultimate independence ratio can be determined Zhu (1996) constructed graphs satisfying 1 χ ( G ) < I ( G ) < 1 χ f ( G ) \frac{1}{\chi(G)} < \mathscr{I}(G) < \frac{1}{\chi_f(G)} χ ( G ) 1 < I ( G ) < χ f ( G ) 1 , showing that equality does not hold in general Special Properties of Wheel Graphs :Even wheels W 2 t W_{2t} W 2 t : χ ( W 2 t ) = ω ( W 2 t ) = 3 \chi(W_{2t}) = \omega(W_{2t}) = 3 χ ( W 2 t ) = ω ( W 2 t ) = 3 , thus I ( W 2 t ) = 1 3 \mathscr{I}(W_{2t}) = \frac{1}{3} I ( W 2 t ) = 3 1 Odd wheels W 2 t + 1 W_{2t+1} W 2 t + 1 : χ ( W 2 t + 1 ) = 4 \chi(W_{2t+1}) = 4 χ ( W 2 t + 1 ) = 4 , ω ( W 2 t + 1 ) = 3 \omega(W_{2t+1}) = 3 ω ( W 2 t + 1 ) = 3 , thus 1 4 ≤ I ( W 2 t + 1 ) ≤ 1 3 \frac{1}{4} \leq \mathscr{I}(W_{2t+1}) \leq \frac{1}{3} 4 1 ≤ I ( W 2 t + 1 ) ≤ 3 1 Core Conjecture (Conjecture 1.2): I ( W 2 t + 1 ) = 1 4 \mathscr{I}(W_{2t+1}) = \frac{1}{4} I ( W 2 t + 1 ) = 4 1 30-Year-Old Open Problem : Hahn revisited this problem at the Canadian Mathematical Society Winter Conference in 2024Smallest Unknown Case : The 5-wheel W 5 W_5 W 5 is the smallest graph with unknown ultimate independence ratioTheoretical Importance : May provide insights for the general theory of ultimate independence ratioComputational Challenge : Estimating I ( W 2 t + 1 ) \mathscr{I}(W_{2t+1}) I ( W 2 t + 1 ) to arbitrary precision using known methods is algorithmically infeasibleThe main contributions of this paper include:
New General Upper Bound Method (Theorem 1.3):For graphs G G G containing a k k k -clique, proves I ( G ) ≤ α ( G □ ℓ □ K k ) k ∣ V ( G ) ∣ ℓ \mathscr{I}(G) \leq \frac{\alpha(G^{\Box \ell} \Box K_k)}{k|V(G)|^\ell} I ( G ) ≤ k ∣ V ( G ) ∣ ℓ α ( G □ ℓ □ K k ) Generalizes results of Hahn, Hell, and Poljak on vertex-transitive graphs Improved Bounds for Large Odd Wheels (Theorem 1.5):For all t ≥ 3 t \geq 3 t ≥ 3 , proves I ( W 2 t + 1 ) ≤ 4 t 2 + 6 t 3 ( 2 t + 2 ) 2 \mathscr{I}(W_{2t+1}) \leq \frac{4t^2+6t}{3(2t+2)^2} I ( W 2 t + 1 ) ≤ 3 ( 2 t + 2 ) 2 4 t 2 + 6 t First non-trivial theoretical bound for large odd wheels (without computer assistance) Exact Computation of Critical Values (Theorem 1.6):Computer-assisted proof that α ( W 5 □ 3 □ K 3 ) = 170 \alpha(W_5^{\Box 3} \Box K_3) = 170 α ( W 5 □ 3 □ K 3 ) = 170 Combines structural arguments and integer programming Improved Bound for 5-Wheel (Theorem 1.7):Proves I ( W 5 ) ≤ 1019 3888 \mathscr{I}(W_5) \leq \frac{1019}{3888} I ( W 5 ) ≤ 3888 1019 First improvement of this bound in 30 years Computational Framework :Develops systematic methodology combining theoretical analysis and computational verification All code is publicly available and reproducible Core Task : Establish tighter upper bounds for the ultimate independence ratio I ( W 2 t + 1 ) \mathscr{I}(W_{2t+1}) I ( W 2 t + 1 ) of odd wheel graphs.
Technical Challenges :
Direct computation of α ( G □ k ) \alpha(G^{\Box k}) α ( G □ k ) is infeasible for large k k k Requires combining theoretical estimates with finite computation Standard methods fail for odd wheels where chromatic number differs from clique number The paper employs a three-level progressive approach:
Core Idea : Utilize clique structure in graphs to improve bounds.
Theorem Statement : If G G G contains a k k k -clique, then for any ℓ ≥ 1 \ell \geq 1 ℓ ≥ 1 :
I ( G ) ≤ α ( G □ ℓ □ K k ) k ∣ V ( G ) ∣ ℓ \mathscr{I}(G) \leq \frac{\alpha(G^{\Box \ell} \Box K_k)}{k|V(G)|^\ell} I ( G ) ≤ k ∣ V ( G ) ∣ ℓ α ( G □ ℓ □ K k )
and
I ( G ) = lim ℓ → ∞ α ( G □ ℓ □ K k ) k ∣ V ( G ) ∣ ℓ \mathscr{I}(G) = \lim_{\ell \to \infty} \frac{\alpha(G^{\Box \ell} \Box K_k)}{k|V(G)|^\ell} I ( G ) = lim ℓ → ∞ k ∣ V ( G ) ∣ ℓ α ( G □ ℓ □ K k )
Proof Technique :
Recurrence Relation : For maximum independent set I I I in G □ ( p + 1 ) G^{\Box (p+1)} G □ ( p + 1 ) , decompose by last coordinate:
α ( G □ ( p + 1 ) ) ≤ α ( G □ p □ K k ) + ( n − k ) α ( G □ p ) \alpha(G^{\Box (p+1)}) \leq \alpha(G^{\Box p} \Box K_k) + (n-k)\alpha(G^{\Box p}) α ( G □ ( p + 1 ) ) ≤ α ( G □ p □ K k ) + ( n − k ) α ( G □ p ) Limit Analysis :
i ( G □ ( p + 1 ) ) ≤ α ( G □ ℓ □ K k ) n ℓ + 1 ∑ i = 0 p − ℓ ( n − k n ) i + ( n − k ) p − ℓ + 1 α ( G □ ℓ ) n p + 1 i(G^{\Box (p+1)}) \leq \frac{\alpha(G^{\Box \ell} \Box K_k)}{n^{\ell+1}} \sum_{i=0}^{p-\ell} \left(\frac{n-k}{n}\right)^i + \frac{(n-k)^{p-\ell+1}\alpha(G^{\Box \ell})}{n^{p+1}} i ( G □ ( p + 1 ) ) ≤ n ℓ + 1 α ( G □ ℓ □ K k ) ∑ i = 0 p − ℓ ( n n − k ) i + n p + 1 ( n − k ) p − ℓ + 1 α ( G □ ℓ ) Geometric Series Summation : As p → ∞ p \to \infty p → ∞ , the second term vanishes and the first term converges to α ( G □ ℓ □ K k ) k n ℓ \frac{\alpha(G^{\Box \ell} \Box K_k)}{kn^\ell} k n ℓ α ( G □ ℓ □ K k ) Application to Odd Wheels (Corollary 1.4): Since W 2 t + 1 W_{2t+1} W 2 t + 1 contains K 3 K_3 K 3 , taking k = 3 k=3 k = 3 yields:
I ( W 2 t + 1 ) ≤ α ( W 2 t + 1 □ ℓ □ K 3 ) 3 ( 2 t + 2 ) ℓ \mathscr{I}(W_{2t+1}) \leq \frac{\alpha(W_{2t+1}^{\Box \ell} \Box K_3)}{3(2t+2)^\ell} I ( W 2 t + 1 ) ≤ 3 ( 2 t + 2 ) ℓ α ( W 2 t + 1 □ ℓ □ K 3 )
Key Lemma (Lemma 4.2): If I I I is an independent set in W 2 t + 1 □ 2 W_{2t+1}^{\Box 2} W 2 t + 1 □ 2 and I ∗ I_* I ∗ is the part involving the central vertex, then if ∣ I ∗ ∖ { ( w ∗ , w ∗ ) } ∣ = k |I_* \setminus \{(w_*, w_*)\}| = k ∣ I ∗ ∖ {( w ∗ , w ∗ )} ∣ = k :
∣ I ∣ ≤ t ( 2 t + 1 ) + 1 + ( 1 − t ) k |I| \leq t(2t+1) + 1 + (1-t)k ∣ I ∣ ≤ t ( 2 t + 1 ) + 1 + ( 1 − t ) k
Exact Value (Proposition 4.3):
α ( W 2 t + 1 □ 2 ) = ( 2 t + 1 ) t + 1 \alpha(W_{2t+1}^{\Box 2}) = (2t+1)t + 1 α ( W 2 t + 1 □ 2 ) = ( 2 t + 1 ) t + 1
Proof Strategy :
Constructive Lower Bound: Explicitly construct an independent set of size ( 2 t + 1 ) t + 1 (2t+1)t+1 ( 2 t + 1 ) t + 1 Upper Bound Proof: Using Lemma 4.2, if ∣ I ∣ > ( 2 t + 1 ) t + 1 |I| > (2t+1)t+1 ∣ I ∣ > ( 2 t + 1 ) t + 1 , then k ≥ 2 k \geq 2 k ≥ 2 , leading to contradiction Estimate for Ternary Product : For W 2 t + 1 □ 2 □ K 3 W_{2t+1}^{\Box 2} \Box K_3 W 2 t + 1 □ 2 □ K 3 , let S i S_i S i denote the part corresponding to the i i i -th vertex of K 3 K_3 K 3 . Through careful counting arguments:
α ( W 2 t + 1 □ 2 □ K 3 ) ≤ 4 t 2 + 6 t \alpha(W_{2t+1}^{\Box 2} \Box K_3) \leq 4t^2 + 6t α ( W 2 t + 1 □ 2 □ K 3 ) ≤ 4 t 2 + 6 t
This directly yields Theorem 1.5 .
Integer Programming Formulation :
Standard independent set IP:
max ∑ v ∈ V ( G ) x v s.t. B ( G ) T x ≤ 1 , x ∈ { 0 , 1 } ∣ V ( G ) ∣ \max \sum_{v \in V(G)} x_v \quad \text{s.t.} \quad B(G)^T x \leq \mathbf{1}, \quad x \in \{0,1\}^{|V(G)|} max ∑ v ∈ V ( G ) x v s.t. B ( G ) T x ≤ 1 , x ∈ { 0 , 1 } ∣ V ( G ) ∣
Improved Formulation for Cartesian Products (Program 7):
For G □ H G \Box H G □ H , introduce variable vectors x v x_v x v (v ∈ V ( H ) v \in V(H) v ∈ V ( H ) ):
max ∑ v ∈ V ( H ) 1 T x v \max \sum_{v \in V(H)} \mathbf{1}^T x_v max ∑ v ∈ V ( H ) 1 T x v s.t. B ( G ) T x v ≤ 1 ∀ v ∈ V ( H ) \text{s.t.} \quad B(G)^T x_v \leq \mathbf{1} \quad \forall v \in V(H) s.t. B ( G ) T x v ≤ 1 ∀ v ∈ V ( H ) x u + x v ≤ 1 ∀ u v ∈ E ( H ) x_u + x_v \leq \mathbf{1} \quad \forall uv \in E(H) x u + x v ≤ 1 ∀ uv ∈ E ( H )
Advantages : Allows convenient addition of structural constraints (e.g., 1 T x v ≥ k \mathbf{1}^T x_v \geq k 1 T x v ≥ k ) to study specific types of independent sets.
Branch-and-Bound Strategy (Lemma 6.2-6.4):
For W 5 □ 3 W_5^{\Box 3} W 5 □ 3 , systematically enumerate possible independent set size distributions:
Let I i I_i I i denote the part in the i i i -th coordinate layer Build decision tree by sizes ∣ I ∗ ∣ , ∣ I 0 ∣ , … , ∣ I 4 ∣ |I_*|, |I_0|, \ldots, |I_4| ∣ I ∗ ∣ , ∣ I 0 ∣ , … , ∣ I 4 ∣ Verify feasibility at each node using IP Key Computational Results :
Lemma 6.2(v): If ∣ I ∣ = 58 |I| = 58 ∣ I ∣ = 58 in W 5 □ 3 W_5^{\Box 3} W 5 □ 3 , then (WLOG) i ∈ { ( 9 , 11 , 9 , 11 , 9 , 9 ) , ( 8 , 11 , 9 , 10 , 10 , 10 ) } \mathbf{i} \in \{(9,11,9,11,9,9), (8,11,9,10,10,10)\} i ∈ {( 9 , 11 , 9 , 11 , 9 , 9 ) , ( 8 , 11 , 9 , 10 , 10 , 10 )} Lemma 6.4: Eliminates existence of independent sets of size 171 in W 5 □ 3 □ K 3 W_5^{\Box 3} \Box K_3 W 5 □ 3 □ K 3 Theoretical Innovation :Theorem 1.3 provides new methodology utilizing clique structure, independent of vertex-transitivity Limit equality I ( G ) = lim ℓ → ∞ α ( G □ ℓ □ K k ) k ∣ V ( G ) ∣ ℓ \mathscr{I}(G) = \lim_{\ell \to \infty} \frac{\alpha(G^{\Box \ell} \Box K_k)}{k|V(G)|^\ell} I ( G ) = lim ℓ → ∞ k ∣ V ( G ) ∣ ℓ α ( G □ ℓ □ K k ) provides new computational pathway Structural Insights :Lemma 4.2 establishes precise relationship between independent set size and central vertex usage Identifies key structural patterns of independent sets in W 2 t + 1 □ 2 □ K 3 W_{2t+1}^{\Box 2} \Box K_3 W 2 t + 1 □ 2 □ K 3 Computational Methodology :Organically combines theoretical bounds with finite computation Hybrid strategy of branch-and-bound + IP effectively handles exponential search space Utilizes automorphism group to simplify fractional chromatic number computation (Theorem 5.1) Reproducibility :All computational steps have corresponding code file links Provides detailed verification pathways Hardware : Authors' personal laptopSoftware :
Python with Gurobi optimization solver (integer programming) SageMath (basic graph theory computation) Open source code: https://github.com/Shivaramkratos/Ultimate_Independence_ratio_Python_Gurobi_code Independence Number Computation :α ( W 5 □ 2 ) = 11 \alpha(W_5^{\Box 2}) = 11 α ( W 5 □ 2 ) = 11 α ( W 5 □ 3 ) = 58 \alpha(W_5^{\Box 3}) = 58 α ( W 5 □ 3 ) = 58 α ( W 5 □ 2 □ K 3 ) = 29 \alpha(W_5^{\Box 2} \Box K_3) = 29 α ( W 5 □ 2 □ K 3 ) = 29 α ( W 5 □ 3 □ K 3 ) = 170 \alpha(W_5^{\Box 3} \Box K_3) = 170 α ( W 5 □ 3 □ K 3 ) = 170 (main contribution)Fractional Chromatic Number Computation :Use linear programming to compute χ f ( W 2 t + 1 □ 2 ) \chi_f(W_{2t+1}^{\Box 2}) χ f ( W 2 t + 1 □ 2 ) Simplify constraints through orbits of maximum independent sets Upper Bound Verification :α ( W 5 □ 4 ) ≤ 343 \alpha(W_5^{\Box 4}) \leq 343 α ( W 5 □ 4 ) ≤ 343 α ( W 5 □ 4 □ K 3 ) ≤ 1019 \alpha(W_5^{\Box 4} \Box K_3) \leq 1019 α ( W 5 □ 4 □ K 3 ) ≤ 1019 Branch-and-Bound Tree :
For example, verification of Lemma 6.2(v) involves 9 leaf nodes Each leaf node corresponds to an independent IP instance All infeasible cases have corresponding code file verification Case Analysis :
Symmetry Utilization: Reduce cases to check via automorphism group of W 2 t + 1 W_{2t+1} W 2 t + 1 Structural Pruning: Use α ( W 2 t + 1 □ 2 □ K 3 ) = 29 \alpha(W_{2t+1}^{\Box 2} \Box K_3) = 29 α ( W 2 t + 1 □ 2 □ K 3 ) = 29 to eliminate certain size combinations 2 t + 1 2t+1 2 t + 1 α ( W 2 t + 1 □ 3 ) \alpha(W_{2t+1}^{\Box 3}) α ( W 2 t + 1 □ 3 ) α ( W 2 t + 1 □ 2 □ K 3 ) \alpha(W_{2t+1}^{\Box 2} \Box K_3) α ( W 2 t + 1 □ 2 □ K 3 ) I ( W 2 t + 1 ) ≤ \mathscr{I}(W_{2t+1}) \leq I ( W 2 t + 1 ) ≤ Previous Bound 5 58 29 1019/3888 ≈ 0.262 11/41 ≈ 0.268 7 156 54 9/32 = 0.281 t/(3t+1) ≈ 0.304 9 336 87 29/100 = 0.29 0.310 11 620 128 8/27 ≈ 0.296 0.314 13 1032 177 59/196 ≈ 0.301 0.317
Key Observations :
All new bounds significantly improve previous bounds t 3 t + 1 \frac{t}{3t+1} 3 t + 1 t For t ≥ 3 t \geq 3 t ≥ 3 , general bound 4 t 2 + 6 t 3 ( 2 t + 2 ) 2 \frac{4t^2+6t}{3(2t+2)^2} 3 ( 2 t + 2 ) 2 4 t 2 + 6 t asymptotically approaches 1 3 \frac{1}{3} 3 1 from below Still has gap with conjectured value 1 4 \frac{1}{4} 4 1 Core Result : α ( W 5 □ 3 □ K 3 ) = 170 \alpha(W_5^{\Box 3} \Box K_3) = 170 α ( W 5 □ 3 □ K 3 ) = 170
Proof Structure :
Lower Bound : Figure 4 exhibits explicit independent set of size 170Upper Bound : Lemmas 6.2-6.5 systematically eliminate possibility of size ≥ 171Key Lemma Verification :
Lemma 6.2: 5 assertions involving approximately 15 IP instances Lemma 6.3: Handles size 57 case, 6 cases Lemma 6.4: Handles boundary case of size 170-171, approximately 15 IP instances Lemma 6.5: Final synthesis, confirms only two possible size distributions Theorem 6.6 : α ( W 5 □ 4 ) ≤ 343 \alpha(W_5^{\Box 4}) \leq 343 α ( W 5 □ 4 ) ≤ 343
Proof Strategy :
Assume ∣ I ∣ = 344 |I| = 344 ∣ I ∣ = 344 , decompose by coordinate layers Use α ( W 5 □ 3 □ K 3 ) = 170 \alpha(W_5^{\Box 3} \Box K_3) = 170 α ( W 5 □ 3 □ K 3 ) = 170 to establish constraints Derive contradiction: requires ∣ I ∗ ∣ = 54 |I_*| = 54 ∣ I ∗ ∣ = 54 and all ∣ I i ∣ = 58 |I_i| = 58 ∣ I i ∣ = 58 But this requires certain layers to satisfy impossible size combinations (Lemma 6.4) Theorem 6.7 : α ( W 5 □ 4 □ K 3 ) ≤ 1019 \alpha(W_5^{\Box 4} \Box K_3) \leq 1019 α ( W 5 □ 4 □ K 3 ) ≤ 1019
Proof Strategy :
Assume ∣ S ∣ = 1020 |S| = 1020 ∣ S ∣ = 1020 , implying all 6 coordinate layers achieve maximum value 170 Use Lemma 6.5, each layer must be specific size distribution By pigeonhole principle, at least one coordinate has two different K 3 K_3 K 3 parts not achieving maximum Results in total ≤ 1019 \leq 1019 ≤ 1019 Final Upper Bound :
I ( W 5 ) ≤ 1019 3888 ≈ 0.26217 \mathscr{I}(W_5) \leq \frac{1019}{3888} \approx 0.26217 I ( W 5 ) ≤ 3888 1019 ≈ 0.26217
This improves the 1995 bound 11 41 ≈ 0.26829 \frac{11}{41} \approx 0.26829 41 11 ≈ 0.26829 by approximately 2.3% .
Method (Section 5.2):
Compute 1 χ f ( G ) \frac{1}{\chi_f(G)} χ f ( G ) 1 via LP:
min z s.t. ∑ v x v = 1 , ∑ v ∈ I x v ≤ z ∀ I ∈ I max ( G ) \min z \quad \text{s.t.} \quad \sum_v x_v = 1, \quad \sum_{v \in I} x_v \leq z \quad \forall I \in \mathcal{I}_{\max}(G) min z s.t. ∑ v x v = 1 , ∑ v ∈ I x v ≤ z ∀ I ∈ I m a x ( G )
Automorphism Group Simplification (Theorem 5.1):
Optimal solution exists that is constant on vertex orbits, so only need to consider profiles of maximum independent sets.
Profiles for W₅² (from 7 ):
( 1 , 0 , 10 ) , ( 0 , 2 , 8 ) , ( 0 , 3 , 6 ) , ( 0 , 4 , 5 ) (1, 0, 10), (0, 2, 8), (0, 3, 6), (0, 4, 5) ( 1 , 0 , 10 ) , ( 0 , 2 , 8 ) , ( 0 , 3 , 6 ) , ( 0 , 4 , 5 )
where ( p 1 , p 2 , p 3 ) (p_1, p_2, p_3) ( p 1 , p 2 , p 3 ) denotes number of vertices in three orbits T 1 , T 2 , T 3 T_1, T_2, T_3 T 1 , T 2 , T 3 .
Computational Results :
χ f ( W 7 □ 2 ) = 39 11 \chi_f(W_7^{\Box 2}) = \frac{39}{11} χ f ( W 7 □ 2 ) = 11 39 χ f ( W 9 □ 2 ) = 127 37 \chi_f(W_9^{\Box 2}) = \frac{127}{37} χ f ( W 9 □ 2 ) = 37 127 χ f ( W 5 □ 3 ) = 722 189 \chi_f(W_5^{\Box 3}) = \frac{722}{189} χ f ( W 5 □ 3 ) = 189 722 (computationally intensive)Observed Pattern (Conjecture 7.3):
χ f ( W 2 t + 1 □ 2 ) = 6 t 2 + 7 t + 3 2 t 2 + t + 1 \chi_f(W_{2t+1}^{\Box 2}) = \frac{6t^2 + 7t + 3}{2t^2 + t + 1} χ f ( W 2 t + 1 □ 2 ) = 2 t 2 + t + 1 6 t 2 + 7 t + 3
This would yield I ( W 2 t + 1 ) ≥ 2 t 2 + t + 1 6 t 2 + 7 t + 3 \mathscr{I}(W_{2t+1}) \geq \frac{2t^2+t+1}{6t^2+7t+3} I ( W 2 t + 1 ) ≥ 6 t 2 + 7 t + 3 2 t 2 + t + 1 (asymptotically 1 3 \frac{1}{3} 3 1 ).
Appendix A : Shows maximum independent sets in W 2 t + 1 □ 2 □ K 3 W_{2t+1}^{\Box 2} \Box K_3 W 2 t + 1 □ 2 □ K 3 (as 3-colorings of W 2 t + 1 □ 2 W_{2t+1}^{\Box 2} W 2 t + 1 □ 2 ):
Figure 5: W 5 □ 2 □ K 3 W_5^{\Box 2} \Box K_3 W 5 □ 2 □ K 3 , size 29 Figure 6: W 7 □ 2 □ K 3 W_7^{\Box 2} \Box K_3 W 7 □ 2 □ K 3 , size 54 Figure 7: W 9 □ 2 □ K 3 W_9^{\Box 2} \Box K_3 W 9 □ 2 □ K 3 , size 87 Figure 8: W 11 □ 2 □ K 3 W_{11}^{\Box 2} \Box K_3 W 11 □ 2 □ K 3 , size 128 Structural Observation (Conjecture 7.1):
α ( W 2 t + 1 □ 2 □ K 3 ) = ( 2 t + 2 ) α ( W 2 t + 1 ) − ( t − 1 ) = 4 t 2 + 5 t + 3 \alpha(W_{2t+1}^{\Box 2} \Box K_3) = (2t+2)\alpha(W_{2t+1}) - (t-1) = 4t^2 + 5t + 3 α ( W 2 t + 1 □ 2 □ K 3 ) = ( 2 t + 2 ) α ( W 2 t + 1 ) − ( t − 1 ) = 4 t 2 + 5 t + 3
That is, the order of maximum 3-colorable subgraph is 4 t 2 + 5 t + 3 4t^2 + 5t + 3 4 t 2 + 5 t + 3 .
Foundational Work :Hell, Yu, Zhou (1994): Formalize definition, prove limit existence Hahn, Hell, Poljak (1995): Establish basic bounds 1 χ ≤ I ≤ 1 ω \frac{1}{\chi} \leq \mathscr{I} \leq \frac{1}{\omega} χ 1 ≤ I ≤ ω 1 Tightness of General Bounds :Zhu (1996): Constructs examples with 1 χ < I < 1 χ f \frac{1}{\chi} < \mathscr{I} < \frac{1}{\chi_f} χ 1 < I < χ f 1 Introduces star chromatic number for new bounds Related Limit Problems :Shannon Capacity: lim k → ∞ α ( G ⊠ k ) k \lim_{k \to \infty} \sqrt[k]{\alpha(G^{\boxtimes k})} lim k → ∞ k α ( G ⊠ k ) (strong product) Fractional Independence Ratio (Albert et al. 2001) Tensor Graph Powers (Alon & Lubetzky 2007) Chromatic and Clique Numbers :Even wheels: χ = ω = 3 \chi = \omega = 3 χ = ω = 3 (perfect graphs) Odd wheels: χ = 4 \chi = 4 χ = 4 , ω = 3 \omega = 3 ω = 3 (4-critical graphs) Fractional Chromatic Number :χ f ( W 2 t + 1 ) = 3 + 1 t \chi_f(W_{2t+1}) = 3 + \frac{1}{t} χ f ( W 2 t + 1 ) = 3 + t 1 (by connectivity properties)χ f ( C 2 t + 1 ) = 2 + 1 t \chi_f(C_{2t+1}) = 2 + \frac{1}{t} χ f ( C 2 t + 1 ) = 2 + t 1 (known)Independence Numbers of Cartesian Products :Proposition 2.1: α ( G □ H ) ≤ min { ∣ V ( G ) ∣ α ( H ) , ∣ V ( H ) ∣ α ( G ) } \alpha(G \Box H) \leq \min\{|V(G)|\alpha(H), |V(H)|\alpha(G)\} α ( G □ H ) ≤ min { ∣ V ( G ) ∣ α ( H ) , ∣ V ( H ) ∣ α ( G )} Integer Programming :Standard independent set IP (Program 6) Improved formulation for Cartesian products (Program 7) Fractional Chromatic Number Computation :LP formulation (Program 8) Symmetry utilization (Theorem 5.1) Graph Homomorphisms :Theorem 1.1: If homomorphism G → H G \to H G → H exists, then I ( H ) ≤ I ( G ) \mathscr{I}(H) \leq \mathscr{I}(G) I ( H ) ≤ I ( G ) Provides proof of basic bounds Theoretical Contributions :Proposes new upper bound method utilizing clique structure (Theorem 1.3) Establishes non-trivial theoretical bounds 4 t 2 + 6 t 3 ( 2 t + 2 ) 2 \frac{4t^2+6t}{3(2t+2)^2} 3 ( 2 t + 2 ) 2 4 t 2 + 6 t for all t ≥ 3 t \geq 3 t ≥ 3 Computational Breakthrough :Exactly determines α ( W 5 □ 3 □ K 3 ) = 170 \alpha(W_5^{\Box 3} \Box K_3) = 170 α ( W 5 □ 3 □ K 3 ) = 170 Improves 30-year-old bound for 5-wheel Methodology :Demonstrates effective combination of theoretical analysis and computational verification Provides complete reproducible framework Gap with Conjecture :New bound 1019 3888 ≈ 0.262 \frac{1019}{3888} \approx 0.262 3888 1019 ≈ 0.262 still has ~5% gap from conjectured value 1 4 = 0.25 \frac{1}{4} = 0.25 4 1 = 0.25 Bound for large odd wheels 4 t 2 + 6 t 3 ( 2 t + 2 ) 2 \frac{4t^2+6t}{3(2t+2)^2} 3 ( 2 t + 2 ) 2 4 t 2 + 6 t asymptotically approaches 1 3 \frac{1}{3} 3 1 , not 1 4 \frac{1}{4} 4 1 Computational Limitations :Cannot directly compute α ( W 5 □ 4 ) \alpha(W_5^{\Box 4}) α ( W 5 □ 4 ) (estimated at 338) Higher-order computations (e.g., χ f ( W 7 □ 3 ) \chi_f(W_7^{\Box 3}) χ f ( W 7 □ 3 ) ) are extremely intensive Theoretical Tools :Bound in Lemma 4.2 may not be tight Requires deeper structural understanding Generalization :Methods highly dependent on special structure of wheel graphs Applicability to other non-perfect graphs unknown Authors propose multiple conjectures in Section 7:
Conjecture 7.1 (Structural Conjecture):
α ( W 2 t + 1 □ 2 □ K 3 ) = 4 t 2 + 5 t + 3 \alpha(W_{2t+1}^{\Box 2} \Box K_3) = 4t^2 + 5t + 3 α ( W 2 t + 1 □ 2 □ K 3 ) = 4 t 2 + 5 t + 3 If true, would yield I ( W 2 t + 1 ) ≤ 4 t 2 + 5 t + 3 3 ( 2 t + 2 ) 2 \mathscr{I}(W_{2t+1}) \leq \frac{4t^2+5t+3}{3(2t+2)^2} I ( W 2 t + 1 ) ≤ 3 ( 2 t + 2 ) 2 4 t 2 + 5 t + 3 (slight improvement).Conjecture 7.2 : α ( W 5 □ 4 ) = 338 \alpha(W_5^{\Box 4}) = 338 α ( W 5 □ 4 ) = 338 Heuristic search supports this value. If correct, could further improve bound for I ( W 5 ) \mathscr{I}(W_5) I ( W 5 ) .Conjecture 7.3 (Fractional Chromatic Number Pattern):
χ f ( W 2 t + 1 □ 2 ) = 6 t 2 + 7 t + 3 2 t 2 + t + 1 \chi_f(W_{2t+1}^{\Box 2}) = \frac{6t^2 + 7t + 3}{2t^2 + t + 1} χ f ( W 2 t + 1 □ 2 ) = 2 t 2 + t + 1 6 t 2 + 7 t + 3 Would yield lower bound I ( W 2 t + 1 ) ≥ 2 t 2 + t + 1 6 t 2 + 7 t + 3 \mathscr{I}(W_{2t+1}) \geq \frac{2t^2+t+1}{6t^2+7t+3} I ( W 2 t + 1 ) ≥ 6 t 2 + 7 t + 3 2 t 2 + t + 1 .Conjecture 7.4 (Method Advantage):
For all t ≥ 3 t \geq 3 t ≥ 3 and ℓ ≥ 1 \ell \geq 1 ℓ ≥ 1 ,
α ( W 2 t + 1 □ ℓ □ K 3 ) 3 ( 2 t + 2 ) ℓ < 1 χ f ( W 2 t + 1 □ ℓ ) \frac{\alpha(W_{2t+1}^{\Box \ell} \Box K_3)}{3(2t+2)^\ell} < \frac{1}{\chi_f(W_{2t+1}^{\Box \ell})} 3 ( 2 t + 2 ) ℓ α ( W 2 t + 1 □ ℓ □ K 3 ) < χ f ( W 2 t + 1 □ ℓ ) 1 Conjecture 7.5 (Generalization):
For graphs G G G with χ > ω \chi > \omega χ > ω , if H H H is maximum induced subgraph with χ ( H ) ≤ ω ( G ) \chi(H) \leq \omega(G) χ ( H ) ≤ ω ( G ) , then constant c c c exists such that
χ f ( G ) < c ⋅ ω ( G ) ∣ V ( G ) ∣ ∣ V ( H ) ∣ \chi_f(G) < c \cdot \frac{\omega(G)|V(G)|}{|V(H)|} χ f ( G ) < c ⋅ ∣ V ( H ) ∣ ω ( G ) ∣ V ( G ) ∣ Theoretical Novelty :Theorem 1.3 provides new theoretical tool independent of special graph class assumptions Limit equality establishes computational pathway Lemma 4.2 reveals deep structure of odd wheel Cartesian products Methodological Rigor :Theoretical proofs clear and complete Computational part has complete verification pathways Each computational assertion has corresponding code file Practical Breakthrough :First improvement of I ( W 5 ) \mathscr{I}(W_5) I ( W 5 ) bound in 30 years Provides unified theoretical bound for all large odd wheels Reproducibility :Code completely open source Detailed computational framework explanation (Section 5) Visualization aids understanding (Appendix A) Writing Quality :Clear structure, rigorous logic Appropriate background introduction Good balance between technical details and intuitive explanation Distance from Conjecture :New bound insufficient to prove or disprove Conjecture 1.2 Asymptotic behavior of theoretical bound (approaching 1/3) mismatches conjectured value (1/4) Computational Bottlenecks :Depends on personal computer computational power Cannot handle higher-order cases (e.g., W 5 □ 5 W_5^{\Box 5} W 5 □ 5 ) Fractional chromatic number computation extremely difficult for large graphs Limitations of Theoretical Tools :Bound in Lemma 4.2 may not be tight (gap approximately t t t ) No exact formula provided for α ( W 2 t + 1 □ 2 □ K 3 ) \alpha(W_{2t+1}^{\Box 2} \Box K_3) α ( W 2 t + 1 □ 2 □ K 3 ) Insufficient Generalization :Methods highly customized for wheel graphs Applicability to other χ > ω \chi > \omega χ > ω graphs unclear Limited Lower Bound Work :Primarily focuses on upper bounds Less investigation of lower bounds (e.g., via construction) Contribution to Field :Reactivates 30-year-old open problem Provides new theoretical tool (Theorem 1.3) May inspire research on other non-perfect graphs Practical Value :Computational framework applicable to other graphs' ultimate independence ratio research Integer programming method has general applicability Theoretical Significance :Reveals role of clique structure in ultimate independence ratio Connects independence number, chromatic number, and fractional chromatic number Openness :Proposes 5 specific conjectures Provides clear research directions Direct Applications :Non-homomorphism proofs in graph homomorphism theory Shannon capacity-related problems in network coding Methodological Reference :Hybrid method combining theoretical bounds and finite computation Branch-and-bound + IP strategy Symmetry utilization for computational simplification Generalization Directions :Ultimate independence ratio of other critical graphs Similar problems for other graph products (strong product, lexicographic product) Hahn, Hell, Poljak (1995) : Foundational paper establishing basic theoretical frameworkZhu (1996) : Proves tightness of general boundsHell, Yu, Zhou (1994) : Formal definition of ultimate independence ratioScheinerman & Ullman (2011) : Fractional graph theory textbookHammack et al. (2011) : Graph products handbookThis paper makes substantial progress on the 30-year-old problem of ultimate independence ratio of odd wheels. Through innovative theoretical tools (Theorem 1.3), deep structural analysis (Lemma 4.2), and systematic computational verification, the authors improve bounds for all odd wheels, particularly improving the 5-wheel bound from 0.268 to 0.262. Although still distant from the conjectured value 0.25, this represents important progress in the field. The paper is methodologically rigorous, highly reproducible, and provides solid foundation for subsequent research. The main challenge lies in further narrowing the gap between theoretical bounds and conjectured value, which likely requires deeper understanding of independent set structure in odd wheel Cartesian products.