Higher-Order Approximations of Sojourn Times in M/G/1 Queues via Stein's Method
Bihan Chatterjee, Siva Theja Maguluri, Debankur Mukherjee
TL;DR
This work characterizes the pre-asymptotic behavior of the stationary sojourn time in an M/G/1 queue under heavy traffic by developing a higher-order Stein's method framework. The authors derive explicit $O(\varepsilon^k)$ bounds in the Zolotarev distance between the scaled sojourn time and an exponential limit when the first $k+1$ moments of the service-time distribution match those of the exponential distribution, using a generator-expansion approach and tight derivative bounds of the Stein solution. They establish a hierarchy of accuracy controlled by moment matching, with corollaries translating to Wasserstein bounds and moment convergence. This provides a systematic refinement of the classical exponential heavy-traffic approximation and has practical implications for more accurate performance evaluation in service systems, especially for modest orders of $k$ where a small number of moment-matches yields substantial gains.
Abstract
We study the stationary sojourn time distribution in an M/G/1 queue operating under heavy traffic. It is known that the sojourn time converges to an exponential distribution in the limit. Our focus is on obtaining pre-asymptotic, higher-order approximations that go beyond the classical exponential limit. Using Stein's method, we develop an approach based on higher-order expansions of the generator of the underlying Markov process. The key technical step is to represent higher-order derivatives in terms of lower-order ones and control the resulting error via derivative bounds of the Stein equation. Under suitable moment-matching conditions on the service distribution, we show that the approximation error decays as a high-order power of the slack parameter $\varepsilon=1-ρ$. Error bounds are established in the Zolotarev metric, which further imply bounds on the Wasserstein distance as well as the moments. Our results demonstrate that the accuracy of the exponential approximation can be systematically improved by matching progressively more moments of the service distribution.
