graphics – Adopt Mathematica codes for circumscribed and inscribed ellipsoids to related sets

Users Daniel Huber and Dominic have provided an interesting answer each to the question JohnEllipsoids of constructing circumscribing and inscribing ellipsoids for a certain three-dimensional convex set of quantum-information-theoretic interest. (The issue posed of minimality and maximality of volume has, however, not yet been addressed. Since the inscribed approach involves numerical fitting it seems unlikely to be strictly maximal, but the minimality of the circumscribed one seems less clear. Dominic has noted the availability of a python code “Inner and outer Löwner-John Ellipsoids” Mosekpythoncode for these tasks. However, it appears to require the convex set in question to be a polytope, which is not the case at hand, although the code might be adoptable to a non-polytope input (see comments in Adamaszek1 Adamaszek2).)

The defining constraint for the convex set (the “ordered spectra of two-qubit density matrices”) under study was

 1 > x && x >= y && y >= z && z >= 1 - x - y - z >= 0 && 
 x - z < 2 Sqrt(y (1 - x - y - z)

Further, it was also noted that
there are two other sets of associated interest Adhikari, also circumscribing and inscribed within the convex body under examination. These are given by the constraints

 1 > x && x >= y && y >= z && z >= 1 - x - y - z >= 0 && 
  x^2 + y^2 + (1 - x - y - z)^2 + z^2 < 3/8)

and

 1 > x && x >= y && y >= z && z >= 1 - x - y - z >= 0 && 
  x^2 + y^2 + (1 - x - y - z)^2 + z^2 < 1/3)  .

(Certainly, at least, the latter set is convex–corresponding to separable density matrices–so the question of the John ellipsoids for it seem valid ones to ask.)

We have attempted (clearly not fully successfully) to adopt the answer of Huber to the construction of a circumscribing ellipsoid for the set defined by the latter constraint and the answer of Dominic for the former case. (The volume of the former set is $frac{left(14-3 sqrt{6}right) pi }{3456 sqrt{3}} approx 0.0034909$ and of the latter set, $frac{pi }{864 sqrt{3}} approx 0.0020993$, while that of the original main set of interest is
$frac{1}{576} left(8-6 sqrt{2}-9 sqrt{2} pi +24 sqrt{2} cos^{-1}left(frac{1}{3}right)right) approx 0.00227243$.)

In the circumscribed case, Huber employed the four extremal points

pts={{1/3, 1/3, 1/3}, {1/4, 1/4, 1/4}, {1/2, 1/6, 1/6}, {1/8 (2 + Sqrt(2)), 1/8 (2 + Sqrt(2)), 1/2 (1 + 1/4 (-2 - Sqrt(2)))}}

It appears possible to pursue the modified problem with three of the same points, only replacing

{1/8 (2 + Sqrt(2)), 1/8 (2 + Sqrt(2)), 1/2 (1 + 1/4 (-2 - Sqrt(2)))}

by

{1/12 (3 + Sqrt(3)), 1/12 (3 + Sqrt(3)), 1/12 (3 - Sqrt(3))}

(The second and fourth points remain the furthest apart, as Huber requires, with Norm(pts(2) – pts(4)) equalling $frac{1}{2 sqrt{3}}$.)

My modification of the Huber code and its output is
ModifiedHuber

However, this led to a set of semiaxes

{1/3 I Sqrt(2/949 (81 + 56 Sqrt(3))), Sqrt( 2/429 (-17 + 12 Sqrt(3))), 1/(4 Sqrt(3))}

which when inputted to the volume formula

4/3 a1 a2 a3 (Pi)

gives a clearly unacceptable (imaginary!) volume of

2/351 I Sqrt(1/803 (639 + 20 Sqrt(3))) (Pi)

approx 0. + 0.0163957 I .

Now, here is our adoption of the inscribed ellipsoid code of Dominic (only replacing the original defining constraint by the one with the 3/8 bound).

ModifiedDominic

The output certainly seems less problematical than our modification of the circumscribed ellipsoid code (no imaginary volume,…).

Incorporation into the modified Huber and Dominic two-set plots of a third–the original set (involving the x – z < 2 Sqrt(y (1 – x – y – z) inequality)–should be of interest.

Computation of the volumes of new ellipsoids should be of interest. (For thoroughness, although not of immediate relevance to the problems at hand, let us note that the volumes of the three sets of interest have also been computed when Hilbert-Schmidt measure is attached to them ThreeVolumes. Listed increasingly, these volumes are $frac{35 pi }{23328 sqrt{3}} approx 0.00272132$, $frac{29902415923}{497664}-frac{50274109}{512 sqrt{2}}-frac{3072529845 pi }{32768
sqrt{2}}+frac{1024176615 cos ^{-1}left(frac{1}{3}right)}{4096 sqrt{2}} approx 0.00365826$
and $frac{35 sqrt{frac{1}{3} left(2692167889921345-919847607929856 sqrt{6}right)} pi}{27518828544} approx 0.0483353$. These calculations are of relevance when the sets are viewed in the 15-dimensional context of the two-qubit density matrices, rather than strictly a three-dimensional one.)