Concentrated random variables are frequently used in representing deterministic delays in stochastic models. The squared coefficient of variation ($\mathrm {SCV}$) of the most concentrated phase-type distribution of order $N$ is $1/N$. To further reduce the $\mathrm {SCV}$, concentrated matrix exponential (CME) distributions with complex eigenvalues were investigated recently. It was obtained that the $\mathrm {SCV}$ of an order $N$ CME distribution can be less than $n^{-2.1}$ for odd $N=2n+1$ orders, and the matrix exponential distribution, which exhibits such a low $\mathrm {SCV}$ has complex eigenvalues. In this paper, we consider CME distributions with real eigenvalues (CME-R). We present efficient numerical methods for identifying a CME-R distribution with smallest SCV for a given order $n$. Our investigations show that the $\mathrm {SCV}$ of the most concentrated CME-R of order $N=2n+1$ is less than $n^{-1.85}$. We also discuss how CME-R can be used for numerical inverse Laplace transformation, which is beneficial when the Laplace transform function is impossible to evaluate at complex points.

A necessary condition for a rational Laplace–Stieltjes transform to correspond to a matrix exponential distribution is that the pole of maximal real part is real and negative. Given a rational Laplace–Stieltjes transform with such a pole, we present a method to determine whether or not the numerator polynomial admits a transform that corresponds to a matrix exponential distribution. The method relies on the minimization of a continuous function of one variable over the nonnegative real numbers. Using this approach, we give an alternative characterization for all matrix exponential distributions of order three.

In this paper we introduce certain Hankel matrices that can be used to study ME (matrix exponential) distributions, in particular to compute their ME orders. The Hankel matrices for a given ME probability distribution can be constructed if one of the following five types of information about the distribution is available: (i) an ME representation, (ii) its moments, (iii) the derivatives of its distribution function, (iv) its Laplace-Stieltjes transform, or (v) its distribution function. Using the Hankel matrices, a necessary and sufficient condition for a probability distribution to be an ME distribution is found and a method of computing the ME order of the ME distribution developed. Implications for the PH (phase-type) order of PH distributions are examined. The relationship between the ME order, the PH order, and a lower bound on the PH order given by Aldous and Shepp (1987) is discussed in numerical examples.