Энтропия запутывания диагональных состояний колокола возвращает всегда один ⇐ Python

[Anonymous]

Я пытался вычислить энтропию запутывания для диагональных состояний колокола. Я пробовал много реализаций в Python, и все они дали мне, что все возможные государства имеют эту энтропию, равную 1. < /P>
Это совершенно неправильно, потому что единственными запутанными состояниями должны быть точки Тетраэдр, образованный диагональными состояниями колокола, как показано на рисунке 7 этой статьи. < /p>
Я понятия не имею, что я делаю не так. Это мой Код: < /p>

Код: Выделить всё

import numpy as np
import qutip as q

def tetrahedron(a1, a2, a3):

is_bellD_state = False

plano1 = a1 - a2 + a3
plano2 = a1 + a2 - a3
plano3 = a1 - a2 - a3
plano4 = a1 + a2 + a3

if plano1 >= -1 and plano2  >= -1 and plano3 

Подробнее здесь: [url]https://stackoverflow.com/questions/79452850/entanglement-entropy-of-bell-diagonal-states-returns-always-one[/url]