量子保真度

计算两个量子态保真度的公式

$$F(\rho ,\sigma )={\left(\text{Tr}\sqrt{\sqrt{\rho }\sigma \sqrt{\rho }}\right)}^2={\left(\text{Tr}\sqrt{\sqrt{\sigma }\rho \sqrt{\sigma }}\right)}^2=F(\sigma ,\rho )$$

纯态之间的保真度（两个状态向量）

if rho.ndim == 1 and sigma.ndim == 1:
    return np.abs(np.inner(rho.conj().T, sigma)) ** 2

$$F(\rho ,\sigma )={\left|⟨{\psi }_{\rho }|{\psi }_{\sigma }⟩\right|}^2$$

纯态与混态之间的保真度（状态向量与密度矩阵）

if rho.ndim == 1 and sigma.ndim == 2:
    return np.real(rho.conj().T @ sigma @ rho)

$$F(\rho ,\sigma )=\mathrm{Tr}(⟨{\psi }_{\rho }|\sigma |{\psi }_{\rho }⟩)$$

混态与纯态之间的保真度（密度矩阵与状态向量）

if rho.ndim == 2 and sigma.ndim == 1:
    return np.real(sigma.conj().T @ rho @ sigma)

$$F(\rho ,\sigma )=\mathrm{Tr}(⟨{\psi }_{\sigma }|\rho |{\psi }_{\sigma }⟩)$$

混态之间的保真度（两个密度矩阵）

return np.real(np.trace(rho @ sigma) + 2 * np.sqrt(det(rho) * det(sigma)))

$$F(\rho ,\sigma )=\text{Tr}(\rho \sigma )+2\ast \sqrt{|\rho ||\sigma |}$$
|*|表示矩阵表示的行列式的值