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

From papers

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.

Progress summary

Open

The conjecture that equal time scales make the two competing waiting times equally likely remains unproved, with only special cases and numerical support.

The conjecture asserts that independent Mittag-Leffler waiting times with equal scales satisfy P(Xα<Xβ)=1/2\mathbb{P}(X_{\alpha}<X_{\beta})=1/2 for all indices in (0,1](0,1]. The source presenting it states that this would make the queue’s recurrence classification depend only on the time scales, but gives no proof.

Known results

  • The equality holds when α=β\alpha=\beta, by symmetry.
  • It also holds when α=1\alpha=1 or β=1\beta=1.
  • High-precision numerical calculations support the conjecture.
  • The associated recurrence conclusions are established in the relevant special cases.

2026 scaling-limit update

A newer paper classifies queue scaling behavior according to whether pα,βλ,μp_{\alpha,\beta}^{\lambda,\mu} is above, equal to, or below 1/21/2, but does not prove or refute the equal-time-scale assertion.

Current status (as of August 2026): The conjecture is open; the cases α=β\alpha=\beta, α=1\alpha=1, or β=1\beta=1 and numerical evidence are known, while the general case remains unresolved.

Sources
Sources & referencesView supporting material

Primary source

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

Solutions 1

Proof

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

EesXγ=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Γ(1q)0sγq11+sγds=Γ(1+q/γ)Γ(1q/γ)Γ(1q).\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)0sq11+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/γ)Γ(1z/γ)Γ(1z).\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=logXαlogXβ.W=\log X_\alpha-\log X_\beta.

Its characteristic function is

φW(t)=EXαitEXβ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=dW.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.

0 endorsements
Shivam Patel ·