The equal-time-scale conjecture for Mittag-Leffler renewal queues

At least 3 years old · documented by

Let α,β∈(0,1]\alpha,\beta\in(0,1], and let XαX_{\alpha} and XβX_{\beta} be independent Mittag-Leffler random variables with power indices α\alpha and β\beta, respectively, and equal time scales λ=μ\lambda=\mu. The embedded queue is governed by the arrival probability pα,βλ,μ=P{Xα<Xβ}p_{\alpha,\beta}^{\lambda,\mu}=\mathbb{P}\{X_{\alpha}<X_{\beta}\}. Equal-time-scale conjecture. One has

P{Xα<Xβ}=pα,βλ,λ=12.\mathbb{P}\{X_{\alpha}<X_{\beta}\}=p_{\alpha,\beta}^{\lambda,\lambda}=\frac{1}{2}.

Consequently, the embedded discrete queue is null recurrent if and only if λ=μ\lambda=\mu, and positive recurrent if and only if μ>λ\mu>\lambda; its behaviour depends only on the time scales, not on the power-tail indices. The preceding special cases and numerical simulations support the claim, but no proof or resolution is provided here.

References

Primary source

Jacob Butt, Nicos Georgiou and Enrico Scalas, “Queuing models with Mittag-Leffler inter-event times”, arXiv:2211.13127 (2022).

Progress summary

Refreshed
Claimed solved

A reader-written complete proof claims the conjecture and a stronger scale-threshold result, but neither has been independently verified.

The conjecture, posed by Jacob Butt, Nicos Georgiou, and Enrico Scalas in 2022, says that equal time scales make two independent Mittag-Leffler waiting times equally likely to occur first, regardless of their indices α,β∈(0,1]\alpha,\beta\in(0,1]. The original paper explicitly states that no rigorous proof was obtained.

Known results

  • The equality is proved for α=β\alpha=\beta by symmetry (Butt, Georgiou, and Scalas, 2022).
  • It is known when α=1\alpha=1 or β=1\beta=1; the source gives an explicit comparison formula in that case (Butt, Georgiou, and Scalas, 2022).
  • For general indices, a critical scale ratio exists, but identifying it with 11 is precisely the conjecture (Butt, Georgiou, and Scalas, 2022).

Posted attempt

A reader-written argument claims a complete proof: the Mellin transform makes log⁡Xα−log⁡Xβ\log X_\alpha-\log X_\beta symmetric, yielding probability 1/21/2 at equal scales and a universal threshold at scale ratio 11. The attempt has not been independently verified.

Current status (as of August 2026): A complete proof has been claimed, but the conjecture remains unverified; the classical special cases are settled and the general queue classification remains open.

Sources

Solutions 1

ProofThis solution needs a summarySee full solutionHide full solution

The conjecture holds for every pair of Mittag-Leffler indices. In fact, the complete arbitrary-scale comparison has a universal threshold independent of those indices.

Let 0<γ≤10<\gamma\le1, and let XγX_\gamma have the standard Mittag-Leffler waiting-time distribution

Ee−sXγ=11+sγ.\mathbb E e^{-sX_\gamma}=\frac1{1+s^\gamma}.

The source uses the scale convention

Xγ,λ=Xγλ.X_{\gamma,\lambda}=\frac{X_\gamma}{\lambda}.

For 0<q<γ0<q<\gamma, the fractional-moment identity and Euler's beta integral give

EXγq=qΓ(1−q)∫0∞sγ−q−11+sγ ds=Γ(1+q/γ)Γ(1−q/γ)Γ(1−q).\begin{aligned} \mathbb EX_\gamma^q &=\frac{q}{\Gamma(1-q)} \int_0^\infty \frac{s^{\gamma-q-1}}{1+s^\gamma}\,ds\\ &= \frac{\Gamma(1+q/\gamma)\Gamma(1-q/\gamma)} {\Gamma(1-q)}. \end{aligned}

Similarly,

EXγ−q=1Γ(q)∫0∞sq−11+sγ ds<∞.\mathbb EX_\gamma^{-q} = \frac1{\Gamma(q)} \int_0^\infty\frac{s^{q-1}}{1+s^\gamma}\,ds <\infty.

Hence the Mellin transform is holomorphic on

−γ<ℜz<γ,-\gamma<\Re z<\gamma,

and analytic continuation yields

EXγz=Γ(1+z/γ)Γ(1−z/γ)Γ(1−z).\mathbb EX_\gamma^z = \frac{\Gamma(1+z/\gamma)\Gamma(1-z/\gamma)} {\Gamma(1-z)}.

For independent Xα,XβX_\alpha,X_\beta, set

W=log⁡Xα−log⁡Xβ.W=\log X_\alpha-\log X_\beta.

Its characteristic function is

φW(t)=EXαit EXβ−it=∣Γ(1+it/α)∣2∣Γ(1+it/β)∣2∣Γ(1+it)∣2.\begin{aligned} \varphi_W(t) &=\mathbb EX_\alpha^{it}\,\mathbb EX_\beta^{-it}\\ &= \frac{ |\Gamma(1+it/\alpha)|^2 |\Gamma(1+it/\beta)|^2 }{ |\Gamma(1+it)|^2 }. \end{aligned}

This is real and even. Uniqueness of characteristic functions therefore gives

W=d−W.W\stackrel d=-W.

The waiting-time laws are continuous, so Pr⁡(W=0)=0\Pr(W=0)=0. Consequently

Pr⁡(Xα<Xβ)=12\Pr(X_\alpha<X_\beta)=\frac12

for every α,β∈(0,1]\alpha,\beta\in(0,1], and therefore

pα,βλ,λ=Pr⁡(Xα,λ<Xβ,λ)=12.p_{\alpha,\beta}^{\lambda,\lambda} = \Pr(X_{\alpha,\lambda}<X_{\beta,\lambda}) =\frac12.

More generally,

pα,βλ,μ=Pr⁡(W<log⁡λμ).p_{\alpha,\beta}^{\lambda,\mu} = \Pr\left(W<\log\frac{\lambda}{\mu}\right).

Both waiting-time densities are strictly positive on (0,∞)(0,\infty), so the distribution function of WW is strictly increasing. Its symmetry gives the complete classification

pα,βλ,μ{<1/2,λ<μ,=1/2,λ=μ,>1/2,λ>μ.p_{\alpha,\beta}^{\lambda,\mu} \begin{cases} <1/2,&\lambda<\mu,\\ =1/2,&\lambda=\mu,\\ >1/2,&\lambda>\mu. \end{cases}

Thus the critical scale ratio is exactly 11 for every pair of Mittag-Leffler indices.