Implementasi Quantum Uhlmann Fidelity di Mathematica

Question

Saya ingin menerapkan kesetiaan Uhlmann kuantum $$F(\rho, \sigma) := \mathrm{tr} \left[\sqrt{\sqrt{\rho} \sigma \sqrt{\rho}} \right]^2.$$ di Mathematica sebagai ukuran "kedekatan" antara dua matriks kepadatan.

Namun, menurut saya penerapan saya salah.

$$fidelity[\rho\_, \sigma\_] := (\mathrm{Tr} \left[\sqrt{(\sqrt{\rho}).\sigma.(\sqrt{\rho})} \right])^2$$

Saat saya mengujinya di tes$\rho$ = {{1/9, 2/3, 2/9}, {1/3, 4/9, 4/9}, {2/9, 4/9, 4/9}};

uji$\sigma$ = {{1/9, 1/3, 2/9}, {1/3, 4/9, 4/9}, {2/9, 4/9, 4/9}};

Saya mendapatkan kesetiaan 1, yang jelas tidak benar, karena matriks ini tidak sama. Saya berharap bahwa kesetiaannya adalah 1 jika dan hanya jika$\rho=\sigma$.

Answer 1

Berikut adalah definisi yang tidak sensitif terhadap normalisasi matriks kerapatan:

F[ρ_, σ_] := Tr[MatrixPower[MatrixPower[ρ, 1/2].σ.MatrixPower[ρ, 1/2], 1/2]]^2/(Tr[ρ]*Tr[σ])

Tes dengan semi pasti positif $3\times3$ matriks:

testρ = {{1/5, 1/5, 0}, {1/5, 3/5, 1/5}, {0, 1/5, 1/5}};
testσ = {{1/3, -1/6, 1/6}, {-1/6, 1/3, -1/6}, {1/6, -1/6, 1/3}};

F[testρ, testρ] // FullSimplify
(*    1    *)

F[testρ, testσ] // FullSimplify
(*    8/15    *)

F[testσ, testρ] // FullSimplify
(*    8/15    *)

F[testσ, testσ] // FullSimplify
(*    1    *)