Formulas for Generalized Two-Qubit Separability Probabilities

To begin, we find certain formulas $Q(k,\alpha)= G_1^k(\alpha) G_2^k(\alpha)$,for$k = -1, 0, 1,…,9$. These yield that part of the total separability probability, $P(k,\alpha )$, for generalized (real, complex, quaternionic,\ldots) two-qubit states endowed with random induced measure, for which the determinantal inequality $\mathrm{\mid }{\rho }^{PT}\mathrm{\mid }>\mathrm{\mid }\rho \mathrm{\mid }$ holds. Here $\rho$ denotes a $4\times 4$ density matrix, obtained by tracing over the pure states in $4\times (4+k)$-dimensions, and ${\rho }^{PT}$, its partial transpose. Further, $\alpha$ is a Dyson-index-like parameter with $\alpha =1$ for the standard (15-dimensional) convex set of (complex) two-qubit states. For $k=0$, we obtain the previously reported Hilbert-Schmidt formulas, with (the real case) $Q(0,\frac{1}{2})=\frac{29}{128}$, (the standard complex case) $Q(0,1)=\frac{4}{33}$, and (the quaternionic case) $Q(0,2)= \frac{13}{323}$---thethreesimplyequalling$ P(0,\alpha)/2$\mathrm{.}Thefactors$G_2^k(\alpha)$ are sums of polynomial-weighted generalized hypergeometric functions ${}_p{F}_{p-1}$, $p\ge 7$, all with argument $z=\frac{27}{64} =(\frac{3}{4})^3$. We find number-theoretic-based formulas for the upper ($u_{ik}$) and lower ($b_{ik}$) parameter sets of these functions and, then, equivalently express $G_2^k(\alpha )$ in terms of first-order difference equations. Applications of Zeilberger’s algorithm yield “concise” forms, parallel to the one obtained previously for $P(0,\alpha )=2Q(0,\alpha )$. For nonnegative half-integer and integer values of $\alpha$, $Q(k,\alpha )$ has descending roots starting at $k=-\alpha -1$. Then, we (C. Dunkl and I) construct a remarkably compact (hypergeometric) form for $Q(k,\alpha )$ itself. The possibility of an analogous “master” formula for $P(k,\alpha )$ is, then, investigated, and a number of interesting results found.