انتقل إلى المحتوى الرئيسي
هذا المشروع في مرحلة تجريبية. وجدت خطأ؟ افتح مشكلة على GitHub. لديك سؤال أو فكرة؟ اسأل في Discussions — نحب سماع ملاحظاتك!
كيف جودة الترجمة في هذه الصفحة؟

تمييز الحالة الكمومية والتوموغرافيا

في الجزء الأخير من الدرس، سنتناول باختصار مهمتين مرتبطتين بالقياسات: تمييز الحالة الكمومية وتوموغرافيا الحالة الكمومية.

  1. تمييز الحالة الكمومية

    في تمييز الحالة الكمومية، لدينا مجموعة معروفة من الحالات الكمومية ρ0,…,ρm−1\rho_0,\ldots,\rho_{m-1}، مع احتمالات p0,…,pm−1p_0,\ldots,p_{m-1} مرتبطة بهذه الحالات. طريقة موجزة للتعبير عن ذلك هي القول إن لدينا مجموعة

    {(p0,ρ0),…,(pm−1,ρm−1)}\{(p_0,\rho_0),\ldots,(p_{m-1},\rho_{m-1})\}

    من الحالات الكمومية.

    يُختار عدد a∈{0,…,m−1}a\in\{0,\ldots,m-1\} بشكل عشوائي وفق الاحتمالات (p0,…,pm−1)(p_0,\ldots,p_{m-1}) ويُهيَّأ النظام X\mathsf{X} في الحالة ρa\rho_a. الهدف هو تحديد قيمة aa التي اختيرت عن طريق قياس X\mathsf{X} وحده.

    إذن، لدينا عدد محدود من البدائل مع توزيع مسبق — وهو معرفتنا باحتمالية اختيار كل aa — والهدف هو تحديد أي البدائل حدث فعلًا. قد يكون ذلك سهلًا لبعض خيارات الحالات والاحتمالات، وقد لا يكون ممكنًا لأخرى دون بعض احتمالية الخطأ.

  2. توموغرافيا الحالة الكمومية

    في توموغرافيا الحالة الكمومية، لدينا حالة كمومية مجهولة لنظام — لذا خلافًا لتمييز الحالة الكمومية، لا يوجد عادةً مسبق ولا معلومات عن البدائل الممكنة.

    لكن هذه المرة، لا تُتاح نسخة واحدة من الحالة بل نسخ مستقلة عديدة. أي أن NN نظامًا متطابقًا X1,…,XN\mathsf{X}_1,\ldots,\mathsf{X}_N يُهيَّأ كل منها بصورة مستقلة في الحالة ρ\rho لعدد (ربما كبير) NN. الهدف هو إيجاد تقريب للحالة المجهولة على هيئة مصفوفة كثافة عن طريق قياس الأنظمة.

التمييز بين حالتين​

أبسط حالة في تمييز الحالة الكمومية هي وجود حالتين، ρ0\rho_0 وρ1\rho_1، ينبغي التمييز بينهما.

تخيَّل موقفًا يُختار فيه بت aa بشكل عشوائي: a=0a = 0 باحتمالية pp وa=1a = 1 باحتمالية 1−p1 - p. يُهيَّأ النظام X\mathsf{X} في الحالة ρa\rho_a، أي ρ0\rho_0 أو ρ1\rho_1 اعتمادًا على قيمة aa، ويُسلَّم إلينا. هدفنا هو تخمين قيمة aa بشكل صحيح عن طريق قياس X\mathsf{X}. بدقة أكبر، سنسعى إلى تعظيم احتمالية صحة تخميننا.

قياس أمثل​

تبدأ الطريقة المثلى لحل هذه المسألة بالتحليل الطيفي للفرق الموزون بين ρ0\rho_0 وρ1\rho_1، حيث الأوزان هي الاحتمالات المقابلة.

pρ0−(1−p)ρ1=∑k=0n−1λk∣ψk⟩⟨ψk∣p \rho_0 - (1-p) \rho_1 = \sum_{k = 0}^{n-1} \lambda_k \vert \psi_k \rangle \langle \psi_k \vert

لاحظ أن لدينا إشارة طرح لا جمع في هذا التعبير: هذا فرق موزون لا مجموع موزون.

يمكننا تعظيم احتمالية التخمين الصحيح باختيار قياس إسقاطي {Π0,Π1}\{\Pi_0,\Pi_1\} كما يلي. أولًا لنُقسِّم عناصر {0,…,n−1}\{0,\ldots,n-1\} إلى مجموعتين منفصلتين S0S_0 وS1S_1 بحسب ما إذا كانت القيمة الذاتية المقابلة في الفرق الموزون غير سالبة أو سالبة.

S0={k∈{0,…,n−1}:λk≥0}S1={k∈{0,…,n−1}:λk<0}\begin{gathered} S_0 = \{k\in\{0,\ldots,n-1\} : \lambda_k \geq 0 \}\\[2mm] S_1 = \{k\in\{0,\ldots,n-1\} : \lambda_k < 0 \} \end{gathered}

يمكننا حينئذ اختيار قياس إسقاطي كما يلي.

Π0=∑k∈S0∣ψk⟩⟨ψk∣وΠ1=∑k∈S1∣ψk⟩⟨ψk∣\Pi_0 = \sum_{k \in S_0} \vert \psi_k \rangle \langle \psi_k \vert \quad\text{و}\quad \Pi_1 = \sum_{k \in S_1} \vert \psi_k \rangle \langle \psi_k \vert

(في الواقع، لا يهم في أي من S0S_0 أو S1S_1 ندرج قيم kk التي يكون λk=0\lambda_k = 0 فيها. هنا نختار بشكل تعسفي إدراج هذه القيم في S0S_0.)

هذا هو القياس الأمثل في الموقف المطروح، إذ يُقلِّل احتمالية التحديد الخاطئ للحالة المختارة إلى أدنى حد.

احتمالية الصحة​

الآن سنحدد احتمالية الصحة للقياس {Π0,Π1}\{\Pi_0,\Pi_1\}.

للبدء، لسنا بحاجة فعلًا إلى الاهتمام بالاختيار المحدد لـΠ0\Pi_0 وΠ1\Pi_1، وإن كان قد يكون مفيدًا الاحتفاظ به في الذهن. لأي قياس {P0,P1}\{P_0,P_1\} (ليس بالضرورة إسقاطيًا) يمكن كتابة احتمالية الصحة على النحو التالي.

pTr⁡(P0ρ0)+(1−p)Tr⁡(P1ρ1)p \operatorname{Tr}(P_0 \rho_0) + (1 - p) \operatorname{Tr}(P_1 \rho_1)

باستخدام حقيقة أن {P0,P1}\{P_0,P_1\} قياس، لذا P1=I−P0P_1 = \mathbb{I} - P_0، يمكن إعادة كتابة هذا التعبير على النحو التالي.

pTr⁡(P0ρ0)+(1−p)Tr⁡((I−P0)ρ1)=pTr⁡(P0ρ0)−(1−p)Tr⁡(P0ρ1)+(1−p)Tr⁡(ρ1)=Tr⁡(P0(pρ0−(1−p)ρ1))+1−pp \operatorname{Tr}(P_0 \rho_0) + (1 - p) \operatorname{Tr}((\mathbb{I} - P_0) \rho_1)\hspace*{3cm}\\[1mm] \begin{aligned} & = p \operatorname{Tr}(P_0 \rho_0) - (1 - p) \operatorname{Tr}(P_0 \rho_1) + (1-p) \operatorname{Tr}(\rho_1)\\[1mm] & = \operatorname{Tr}\bigl( P_0 (p \rho_0 - (1-p)\rho_1) \bigr) + 1 - p \end{aligned}

من ناحية أخرى، كان يمكننا تعويض P0=I−P1P_0 = \mathbb{I} - P_1 بدلًا من ذلك. لن يغير ذلك القيمة لكنه يمنحنا تعبيرًا بديلًا.

pTr⁡((I−P1)ρ0)+(1−p)Tr⁡(P1ρ1)=pTr⁡(ρ0)−pTr⁡(P1ρ0)+(1−p)Tr⁡(P1ρ1)=p−Tr⁡(P1(pρ0−(1−p)ρ1))p \operatorname{Tr}((\mathbb{I} - P_1) \rho_0) + (1 - p) \operatorname{Tr}(P_1 \rho_1)\hspace*{3cm}\\[1mm] \begin{aligned} & = p \operatorname{Tr}(\rho_0) - p \operatorname{Tr}(P_1 \rho_0) + (1 - p) \operatorname{Tr}(P_1 \rho_1)\\[1mm] & = p - \operatorname{Tr}\bigl( P_1 (p \rho_0 - (1-p)\rho_1) \bigr) \end{aligned}

التعبيران لهما القيمة ذاتها، لذا يمكننا أخذ متوسطهما للحصول على تعبير آخر لهذه القيمة. (أخذ المتوسط هو مجرد حيلة لتبسيط التعبير الناتج.)

12(Tr⁡(P0(pρ0−(1−p)ρ1))+1−p)+12(p−Tr⁡(P1(pρ0−(1−p)ρ1)))=12Tr⁡((P0−P1)(pρ0−(1−p)ρ1))+12\frac{1}{2} \bigl(\operatorname{Tr}\bigl( P_0 (p \rho_0 - (1-p)\rho_1) \bigr) + 1-p\bigr) + \frac{1}{2} \bigl(p - \operatorname{Tr}\bigl( P_1 (p \rho_0 - (1-p)\rho_1) \bigr)\bigr)\\ = \frac{1}{2} \operatorname{Tr}\bigl( (P_0-P_1) (p \rho_0 - (1-p)\rho_1)\bigr) + \frac{1}{2}

الآن يمكننا أن نرى لماذا من المنطقي اختيار الإسقاطات Π0\Pi_0 وΠ1\Pi_1 (كما حُدِّد أعلاه) من أجل P0P_0 وP1P_1 على التوالي — لأن ذلك يُمكِّننا من جعل الأثر في التعبير الأخير أكبر ما يمكن. بالتحديد،

(Π0−Π1)(pρ0−(1−p)ρ1)=∑k=0n−1∣λk∣⋅∣ψk⟩⟨ψk∣.(\Pi_0-\Pi_1) (p \rho_0 - (1-p)\rho_1) = \sum_{k = 0}^{n-1} \vert\lambda_k\vert \cdot \vert \psi_k \rangle \langle \psi_k \vert.

إذن حين نأخذ الأثر، نحصل على مجموع القيم المطلقة للقيم الذاتية — وهو ما يُعرف بـمعيار الأثر للفرق الموزون.

Tr⁡((Π0−Π1)(pρ0−(1−p)ρ1))=∑k=0n−1∣λk∣=∥pρ0−(1−p)ρ1∥1\operatorname{Tr}\bigl( (\Pi_0-\Pi_1) (p \rho_0 - (1-p)\rho_1)\bigr) = \sum_{k = 0}^{n-1} \vert\lambda_k\vert = \bigl\| p \rho_0 - (1-p)\rho_1 \bigr\|_1

وبالتالي، احتمالية أن يؤدي القياس {Π0,Π1}\{\Pi_0,\Pi_1\} إلى تمييز صحيح لـρ0\rho_0 وρ1\rho_1، المعطاة باحتمالية pp و1−p1-p على التوالي، هي كما يلي.

12+12∥pρ0−(1−p)ρ1∥1\frac{1}{2} + \frac{1}{2} \bigl\| p \rho_0 - (1-p)\rho_1 \bigr\|_1

حقيقة أن هذه هي الاحتمالية المثلى للتمييز الصحيح لـρ0\rho_0 وρ1\rho_1 بالاحتمالات pp و1−p1-p، تُشار إليها عادةً بـمبرهنة هيلستروم-هوليفو (أو أحيانًا مبرهنة هيلستروم فقط).

التمييز بين ثلاث حالات أو أكثر​

في تمييز الحالة الكمومية حين توجد ثلاث حالات أو أكثر، لا يوجد حل مغلق معروف للقياس الأمثل، وإن كان بالإمكان صياغة المسألة كمسألة برمجة شبه محددة (SDP) — مما يتيح تقريبات عددية فعَّالة للقياسات المثلى بمساعدة الحاسوب.

يمكن أيضًا التحقق من (أو نفي) أمثلية قياس معيَّن في مهمة تمييز الحالة عبر شرط يُعرف بـشرط هوليفو-يون-كينيدي-لاكس. بالتحديد، من أجل مهمة تمييز الحالة المعرَّفة بالمجموعة

{(p0,ρ0),…,(pm−1,ρm−1)},\{(p_0,\rho_0),\ldots,(p_{m-1},\rho_{m-1})\},

يكون القياس {P0,…,Pm−1}\{P_0,\ldots,P_{m-1}\} أمثل إذا وفقط إذا كانت المصفوفة

Qa=∑b=0m−1pbρbPb−paρaQ_a = \sum_{b = 0}^{m-1} p_b \rho_b P_b - p_a \rho_a

شبه محددة الإيجابية لكل a∈{0,…,m−1}a\in\{0,\ldots,m-1\}.

على سبيل المثال، تأمَّل مهمة تمييز الحالة الكمومية حيث تُختار إحدى الحالات الرباعية السطوحية ∣ϕ0⟩,…,∣ϕ3⟩\vert\phi_0\rangle,\ldots,\vert\phi_3\rangle بشكل منتظم عشوائي. القياس الرباعي السطوحي {P0,P1,P2,P3}\{P_0,P_1,P_2,P_3\} ينجح باحتمالية

14Tr⁡(P0∣ϕ0⟩⟨ϕ0∣)+14Tr⁡(P1∣ϕ1⟩⟨ϕ1∣)+14Tr⁡(P2∣ϕ2⟩⟨ϕ2∣)+14Tr⁡(P3∣ϕ3⟩⟨ϕ3∣)=12.\frac{1}{4} \operatorname{Tr}(P_0 \vert\phi_0\rangle\langle \phi_0 \vert) + \frac{1}{4} \operatorname{Tr}(P_1 \vert\phi_1\rangle\langle \phi_1 \vert) + \frac{1}{4} \operatorname{Tr}(P_2 \vert\phi_2\rangle\langle \phi_2 \vert) + \frac{1}{4} \operatorname{Tr}(P_3 \vert\phi_3\rangle\langle \phi_3 \vert) = \frac{1}{2}.

هذا هو الأمثل وفقًا لشرط هوليفو-يون-كينيدي-لاكس، كما يكشف الحساب أن

Qa=14(I−∣ϕa⟩⟨ϕa∣)≥0Q_a = \frac{1}{4}(\mathbb{I} - \vert\phi_a\rangle\langle\phi_a\vert) \geq 0

for a=0,1,2,3a = 0,1,2,3.

توموغرافيا الحالة الكمومية​

أخيرًا، سنتناول باختصار مسألة توموغرافيا الحالة الكمومية. في هذه المسألة، يُمنح لنا عدد كبير NN من النسخ المستقلة لحالة كمومية مجهولة ρ\rho، والهدف هو إعادة بناء تقريب ρ~\tilde{\rho} لـρ\rho. للتوضيح، هذا يعني أننا نريد إيجاد وصف كلاسيكي لمصفوفة كثافة ρ~\tilde{\rho} أقرب ما يكون إلى ρ\rho.

يمكننا وصف الإعداد بصيغة بديلة على النحو التالي. تُختار مصفوفة كثافة مجهولة ρ\rho، ويُمنح لنا الوصول إلى NN نظامًا كموميًا X1,…,XN\mathsf{X}_1,\ldots,\mathsf{X}_N، كل منها مُهيَّأ بصورة مستقلة في الحالة ρ\rho. وبالتالي، حالة النظام المركَّب (X1,…,XN)(\mathsf{X}_1,\ldots,\mathsf{X}_N) هي

ρ⊗N=ρ⊗ρ⊗⋯⊗ρ(N مرة)\rho^{\otimes N} = \rho \otimes \rho \otimes \cdots \otimes \rho \quad \text{($N$ مرة)}

الهدف هو إجراء قياسات على الأنظمة X1,…,XN\mathsf{X}_1,\ldots,\mathsf{X}_N وبناءً على نتائج تلك القياسات، حساب مصفوفة كثافة ρ~\tilde{\rho} تُقرِّب ρ\rho بشكل جيد. هذه مسألة رائعة وثمة أبحاث جارية عليها.

يمكن النظر في أنواع مختلفة من الاستراتيجيات لمعالجة المسألة. على سبيل المثال، يمكن تصوُّر استراتيجية يُقاس فيها كل نظام من X1,…,XN\mathsf{X}_1,\ldots,\mathsf{X}_N بشكل منفصل تباعًا، منتِجًا سلسلة من نتائج القياسات. يمكن اتخاذ خيارات محددة مختلفة لأي القياسات تُجرى، بما يشمل الاختيارات التكيُّفية وغير التكيُّفية. بمعنى آخر، قد يعتمد اختيار القياس الذي يُجرى على نظام معيَّن أو لا يعتمد على نتائج القياسات السابقة. بناءً على سلسلة نتائج القياسات، يُستنتج تخمين ρ~\tilde{\rho} للحالة ρ\rho — وهنا أيضًا توجد منهجيات مختلفة للقيام بذلك.

نهج بديل هو إجراء قياس مشترك واحد للمجموعة بأكملها، حيث ننظر إلى (X1,…,XN)(\mathsf{X}_1,\ldots,\mathsf{X}_N) كنظام واحد ونختار قياسًا واحدًا مخرجه تخمين ρ~\tilde{\rho} للحالة ρ\rho. قد يؤدي ذلك إلى تقدير أفضل مما هو ممكن بقياسات منفصلة للأنظمة الفردية، وإن كان القياس المشترك لجميع الأنظمة معًا أصعب تطبيقًا على الأرجح.

توموغرافيا الكيوبت باستخدام قياسات باولي​

سنتناول الآن توموغرافيا الحالة الكمومية في الحالة البسيطة التي تكون فيها ρ\rho مصفوفة كثافة كيوبت. نفترض أننا مُعطَون كيوبتات X1,…,XN\mathsf{X}_1,\ldots,\mathsf{X}_N كل منها في الحالة ρ\rho بصورة مستقلة، وهدفنا حساب تقريب ρ~\tilde{\rho} قريب من ρ\rho.

استراتيجيتنا هي تقسيم الكيوبتات NN وهي X1,…,XN\mathsf{X}_1,\ldots,\mathsf{X}_N إلى ثلاث مجموعات متقاربة الحجم، واحدة لكل من مصفوفات باولي الثلاث σx,\sigma_x, σy,\sigma_y, وσz\sigma_z. يُقاس كل كيوبت بصورة مستقلة على النحو التالي.

  1. لكل كيوبت في المجموعة المرتبطة بـσx\sigma_x نجري قياس σx\sigma_x. هذا يعني قياس الكيوبت بالنسبة للأساس {∣+⟩,∣−⟩}\{\vert + \rangle, \vert -\rangle\}، وهو أساس متعامد معياري من متجهات ذاتية لـσx\sigma_x، ونتائج القياس المقابلة هي القيم الذاتية المرتبطة بالمتجهين الذاتيين: +1+1 للحالة ∣+⟩\vert + \rangle و−1-1 للحالة ∣−⟩\vert -\rangle. بأخذ متوسط النتائج على جميع الحالات في المجموعة المرتبطة بـσx\sigma_x، نحصل على تقريب لقيمة التوقع

    ⟨+∣ρ∣+⟩−⟨−∣ρ∣−⟩=Tr⁡(σxρ).\langle + \vert \rho \vert + \rangle - \langle - \vert \rho \vert - \rangle = \operatorname{Tr}(\sigma_x \rho).
  2. لكل كيوبت في المجموعة المرتبطة بـσy\sigma_y نجري قياس σy\sigma_y. هذا القياس مشابه لقياس σx\sigma_x، غير أن أساس القياس هو {∣ ⁣+ ⁣i⟩,∣ ⁣− ⁣i⟩}\{\vert\! +\!i \rangle, \vert\! -\!i \rangle\}، وهو المتجهات الذاتية لـσy\sigma_y. بأخذ متوسط النتائج على جميع الحالات في المجموعة المرتبطة بـσy\sigma_y، نحصل على تقريب لقيمة التوقع

    ⟨+i∣ρ∣ ⁣+ ⁣i⟩−⟨−i∣ρ∣ ⁣− ⁣i⟩=Tr⁡(σyρ).\langle +i \vert \rho \vert \!+\!i \rangle - \langle -i \vert \rho \vert \!-\!i \rangle = \operatorname{Tr}(\sigma_y \rho).
  3. لكل كيوبت في المجموعة المرتبطة بـσz\sigma_z نجري قياس σz\sigma_z. هنا أساس القياس هو الأساس القياسي {∣0⟩,∣1⟩}\{\vert 0\rangle, \vert 1 \rangle\}، وهو المتجهات الذاتية لـσz\sigma_z. بأخذ متوسط النتائج على جميع الحالات في المجموعة المرتبطة بـσz\sigma_z، نحصل على تقريب لقيمة التوقع

    ⟨0∣ρ∣0⟩−⟨1∣ρ∣1⟩=Tr⁡(σzρ).\langle 0 \vert \rho \vert 0 \rangle - \langle 1 \vert \rho \vert 1 \rangle = \operatorname{Tr}(\sigma_z \rho).

بمجرد حصولنا على التقريبات

αx≈Tr⁡(σxρ),  αy≈Tr⁡(σyρ),  αz≈Tr⁡(σzρ)\alpha_x \approx \operatorname{Tr}(\sigma_x \rho),\; \alpha_y \approx \operatorname{Tr}(\sigma_y \rho),\; \alpha_z \approx \operatorname{Tr}(\sigma_z \rho)

بأخذ متوسط نتائج القياسات لكل مجموعة، يمكننا تقريب ρ\rho بـ

ρ~=I+αxσx+αyσy+αzσz2≈I+Tr⁡(σxρ)σx+Tr⁡(σyρ)σy+Tr⁡(σzρ)σz2=ρ.\tilde{\rho} = \frac{\mathbb{I} + \alpha_x \sigma_x + \alpha_y \sigma_y + \alpha_z \sigma_z}{2} \approx \frac{\mathbb{I} + \operatorname{Tr}(\sigma_x \rho) \sigma_x + \operatorname{Tr}(\sigma_y \rho) \sigma_y + \operatorname{Tr}(\sigma_z \rho) \sigma_z}{2} = \rho.

في الحد عندما يؤول NN إلى اللانهاية، يتقارب هذا التقريب احتماليًا (تقاربًا في الاحتمال) إلى مصفوفة الكثافة الحقيقية ρ\rho بموجب قانون الأعداد الكبيرة، ويمكن استخدام حدود إحصائية معروفة (كـمتراجحة هوفدينغ) لتقييد احتمالية انحراف التقريب ρ~\tilde{\rho} عن ρ\rho بمقادير متنوعة.

غير أن شيئًا مهمًا ينبغي الإشارة إليه: المصفوفة ρ~\tilde{\rho} التي نحصل عليها بهذه الطريقة قد لا تكون مصفوفة كثافة. بالتحديد، رغم أن أثرها سيساوي دائمًا 11، فقد تفشل في كونها شبه محددة الإيجابية. ثمة استراتيجيات معروفة مختلفة "لتدوير" مثل هذا التقريب ρ~\tilde{\rho} (أي إسقاطه) إلى مصفوفة كثافة، إحداها حساب التحليل الطيفي واستبدال أي قيم ذاتية سالبة بـ00 ثم إعادة التطبيع (بقسمة المصفوفة الناتجة على أثرها).

توموغرافيا الكيوبت باستخدام القياس الرباعي السطوحي​

خيار آخر لإجراء توموغرافيا الكيوبت هو قياس كل كيوبت X1,…,XN\mathsf{X}_1,\ldots,\mathsf{X}_N باستخدام القياس الرباعي السطوحي {P0,P1,P2,P3}\{P_0,P_1,P_2,P_3\} المذكور سابقًا. أي

P0=∣ϕ0⟩⟨ϕ0∣2,P1=∣ϕ1⟩⟨ϕ1∣2,P2=∣ϕ2⟩⟨ϕ2∣2,P3=∣ϕ3⟩⟨ϕ3∣2P_0 = \frac{\vert \phi_0 \rangle \langle \phi_0 \vert}{2}, \quad P_1 = \frac{\vert \phi_1 \rangle \langle \phi_1 \vert}{2}, \quad P_2 = \frac{\vert \phi_2 \rangle \langle \phi_2 \vert}{2}, \quad P_3 = \frac{\vert \phi_3 \rangle \langle \phi_3 \vert}{2}

من أجل

∣ϕ0⟩=∣0⟩∣ϕ1⟩=13∣0⟩+23∣1⟩∣ϕ2⟩=13∣0⟩+23e2πi/3∣1⟩∣ϕ3⟩=13∣0⟩+23e−2πi/3∣1⟩.\begin{aligned} \vert \phi_0 \rangle & = \vert 0 \rangle\\ \vert \phi_1 \rangle & = \frac{1}{\sqrt{3}} \vert 0 \rangle + \sqrt{\frac{2}{3}} \vert 1 \rangle\\ \vert \phi_2 \rangle & = \frac{1}{\sqrt{3}} \vert 0 \rangle + \sqrt{\frac{2}{3}} e^{2\pi i/3} \vert 1 \rangle\\ \vert \phi_3 \rangle & = \frac{1}{\sqrt{3}} \vert 0 \rangle + \sqrt{\frac{2}{3}} e^{-2\pi i/3} \vert 1 \rangle. \end{aligned}

كل نتيجة تُحصَّل عليها عدد معيَّن من المرات، سنرمز إليه بـnan_a لكل a∈{0,1,2,3}a\in\{0,1,2,3\}، بحيث n0+n1+n2+n3=Nn_0 + n_1 + n_2 + n_3 = N. نسبة هذه الأعداد إلى NN تُقدِّم تقديرًا للاحتمالية المرتبطة بكل نتيجة ممكنة:

naN≈Tr⁡(Paρ).\frac{n_a}{N} \approx \operatorname{Tr}(P_a \rho).

أخيرًا، سنستخدم الصيغة الرائعة التالية:

ρ=∑a=03(3Tr⁡(Paρ)−12)∣ϕa⟩⟨ϕa∣.\rho = \sum_{a=0}^3 \Bigl( 3 \operatorname{Tr}(P_a \rho) - \frac{1}{2}\Bigr) \vert \phi_a \rangle \langle \phi_a \vert.

لإثبات هذه الصيغة، يمكننا استخدام المعادلة التالية للقيم المطلقة المربَّعة للجداءات الداخلية للحالات الرباعية السطوحية، والتي يمكن التحقق منها بحسابات مباشرة.

∣⟨ϕa∣ϕb⟩∣2={1a=b13a≠b.\bigl\vert \langle \phi_a \vert \phi_b \rangle \bigr\vert^2 = \begin{cases} 1 & a=b\\ \frac{1}{3} & a\neq b. \end{cases}

المصفوفات الأربع

∣ϕ0⟩⟨ϕ0∣=(1000)∣ϕ1⟩⟨ϕ1∣=(13232323)∣ϕ2⟩⟨ϕ2∣=(1323e−2πi/323e2πi/323)∣ϕ3⟩⟨ϕ3∣=(1323e2πi/323e−2πi/323)\begin{aligned} \vert\phi_0\rangle \langle \phi_0 \vert & = \begin{pmatrix} 1 & 0\\[2mm] 0 & 0\end{pmatrix}\\[3mm] \vert\phi_1\rangle \langle \phi_1 \vert & = \begin{pmatrix} \frac{1}{3} & \frac{\sqrt{2}}{3}\\[2mm] \frac{\sqrt{2}}{3} & \frac{2}{3}\end{pmatrix}\\[3mm] \vert\phi_2\rangle \langle \phi_2 \vert & = \begin{pmatrix} \frac{1}{3} & \frac{\sqrt{2}}{3}e^{-2\pi i/3}\\[2mm] \frac{\sqrt{2}}{3}e^{2\pi i/3} & \frac{2}{3}\end{pmatrix}\\[3mm] \vert\phi_3\rangle \langle \phi_3 \vert & = \begin{pmatrix} \frac{1}{3} & \frac{\sqrt{2}}{3}e^{2\pi i/3}\\[2mm] \frac{\sqrt{2}}{3}e^{-2\pi i/3} & \frac{2}{3}\end{pmatrix} \end{aligned}

مستقلة خطيًا، لذا يكفي إثبات صحة الصيغة حين ρ=∣ϕb⟩⟨ϕb∣\rho = \vert\phi_b\rangle\langle\phi_b\vert من أجل b=0,1,2,3b = 0,1,2,3. بالتحديد،

3Tr⁡(Pa∣ϕb⟩⟨ϕb∣)−12=32∣⟨ϕa∣ϕb⟩∣2−12={1a=b0a≠b3 \operatorname{Tr}(P_a \vert\phi_b\rangle\langle\phi_b\vert) - \frac{1}{2} = \frac{3}{2} \vert \langle \phi_a \vert \phi_b \rangle \vert^2 - \frac{1}{2} = \begin{cases} 1 & a=b\\ 0 & a\neq b \end{cases}

وبالتالي

∑a=03(3Tr⁡(Pa∣ϕb⟩⟨ϕb∣)−Tr⁡(∣ϕb⟩⟨ϕb∣)2)∣ϕa⟩⟨ϕa∣=∣ϕb⟩⟨ϕb∣.\sum_{a=0}^3 \biggl( 3 \operatorname{Tr}(P_a \vert\phi_b\rangle\langle\phi_b\vert) - \frac{\operatorname{Tr}(\vert\phi_b\rangle\langle\phi_b\vert)}{2}\biggr) \vert \phi_a \rangle \langle \phi_a \vert = \vert \phi_b\rangle\langle \phi_b \vert.

نصل إلى تقريب لـρ\rho:

ρ~=∑a=03(3naN−12)∣ϕa⟩⟨ϕa∣.\tilde{\rho} = \sum_{a=0}^3 \Bigl( \frac{3 n_a}{N} - \frac{1}{2}\Bigr) \vert \phi_a \rangle \langle \phi_a \vert.

سيكون هذا التقريب دائمًا مصفوفة إرميتية ذات أثر يساوي واحدًا، لكنه قد يفشل في كونه شبه محدد الإيجابية. في هذه الحالة، يجب "تدوير" التقريب إلى مصفوفة كثافة، بصورة مشابهة للاستراتيجية المستخدمة في قياسات باولي.