Two-qubit separability probabilities and Beta functions
Paul B. Slater
Phys. Rev. A 75, 032326 (2007) - Published 20 March, 2007
Due to recent important work of Życzkowski and Sommers [J. Phys. A 36, 10115 (2003); 36, 10083 (2003)], exact formulas are available (in terms of both the Hilbert-Schmidt and the Bures metrics) for the -dimensional and -dimensional volumes of the complex and real density matrices. However, no comparable formulas are available for the volumes (and, hence, probabilities) of various separable subsets of them. We seek to clarify this situation for the Hilbert-Schmidt metric for the simplest possible case of , that is, two-qubit systems. Making use of the density matrix parametrization of Bloore [J. Phys. A 9, 2059 (1976)], we are able to reduce each of the real and complex volume problems to the calculation of a one-dimensional integral, the single relevant variable being a certain ratio of diagonal entries, . The associated integrand in each case is the product of a known Jacobian (highly oscillatory near ) and a certain unknown univariate function, which our extensive numerical (quasi Monte Carlo) computations indicate is very closely proportional to an (incomplete) Beta function , with , in the real case, and in the complex case. Assuming the full applicability of these specific incomplete Beta functions, we undertake separable-volume calculations.