ข้ามไปยังเนื้อหาหลัก
โครงการนี้อยู่ในช่วงเบต้า พบข้อผิดพลาด? สร้าง GitHub Issue มีคำถามหรือไอเดีย? ถามใน Discussions — ยินดีรับฟังความคิดเห็น!
คุณภาพการแปลในหน้านี้เป็นอย่างไร?

การแทนค่าของช่อง

ต่อไปเราจะพูดถึงการแทนค่าทางคณิตศาสตร์ของช่อง

การส่งแบบเชิงเส้นจากเวกเตอร์ไปยังเวกเตอร์สามารถแทนด้วยเมทริกซ์ได้ในแบบที่คุ้นเคย โดยการกระทำของการส่งเชิงเส้นนั้นอธิบายด้วยการคูณเมทริกซ์กับเวกเตอร์ แต่ช่องเป็นการส่งเชิงเส้นจากเมทริกซ์ไปยังเมทริกซ์ ไม่ใช่จากเวกเตอร์ไปยังเวกเตอร์ ดังนั้น โดยทั่วไป เราจะแสดงช่องในเชิงคณิตศาสตร์ได้อย่างไร?

สำหรับช่องบางช่อง เราอาจมีสูตรง่าย ๆ ที่อธิบายมันได้ เช่น ตัวอย่างสามช่อง qubit ที่ไม่ใช่ยูนิทารีที่อธิบายก่อนหน้านี้ แต่ช่องโดยพลการอาจไม่มีสูตรที่สะอาดแบบนั้น จึงไม่ใช่เรื่องปฏิบัติได้ในทางทั่วไปที่จะแสดงช่องในลักษณะนี้

เพื่อเปรียบเทียบ ในการกำหนดสูตรควอนตัมแบบง่าย เราใช้ เมทริกซ์ยูนิทารี เพื่อแทนการกระทำบนเวกเตอร์สถานะควอนตัม โดยที่เมทริกซ์ยูนิทารีทุกตัวแทนการกระทำที่ถูกต้อง และทุกการกระทำที่ถูกต้องสามารถแสดงเป็นเมทริกซ์ยูนิทารีได้ โดยแก่นแท้แล้ว คำถามที่กำลังถามคือ: เราจะทำสิ่งที่คล้ายกันสำหรับช่องได้อย่างไร?

ในการตอบคำถามนี้ เราต้องการเครื่องมือทางคณิตศาสตร์เพิ่มเติม เราจะเห็นว่าช่องสามารถอธิบายทางคณิตศาสตร์ได้หลายวิธี รวมถึงการแทนค่าที่ตั้งชื่อเพื่อเป็นเกียรติแก่สามบุคคลที่มีบทบาทสำคัญในการพัฒนาสิ่งเหล่านี้: Stinespring, Kraus, และ Choi. วิธีต่าง ๆ เหล่านี้ในการอธิบายช่องร่วมกันเสนอมุมมองและการวิเคราะห์ที่แตกต่างกัน

การแทนค่าแบบ Stinespring​

การแทนค่าแบบ Stinespring มีพื้นฐานมาจากแนวคิดที่ว่าช่องทุกช่องสามารถนำไปใช้งานได้ในแบบมาตรฐาน โดยที่ระบบนำเข้าจะถูกรวมกับระบบพื้นที่ทำงานที่ถูกกำหนดค่าเริ่มต้นก่อน ก่อให้เกิดระบบรวม จากนั้นจึงดำเนินการยูนิทารีบนระบบรวม และสุดท้ายระบบพื้นที่ทำงานจะถูกทิ้ง (หรือ traced out) เหลือเพียงผลลัพธ์ของช่อง

ภาพต่อไปนี้แสดงการนำไปใช้งานดังกล่าวในรูปแบบ circuit diagram สำหรับช่องที่ระบบนำเข้าและระบบนำออกเป็นระบบเดียวกัน X\mathsf{X}

แผนภาพแสดงการแทนค่าแบบ Stinespring ของช่องที่ระบบนำเข้าและนำออกเป็นระบบเดียวกัน

ในแผนภาพนี้ สายไฟแทนระบบโดยพลการตามที่ระบุด้วยป้ายกำกับเหนือสายไฟ ไม่จำเป็นต้องเป็น qubit เดี่ยว นอกจากนี้ สัญลักษณ์ กราวด์ ที่ใช้กันทั่วไปในวิศวกรรมไฟฟ้าบ่งชี้อย่างชัดเจนว่า W\mathsf{W} ถูกทิ้ง

อธิบายเป็นคำพูด การนำไปใช้งานทำงานดังนี้ ระบบนำเข้า X\mathsf{X} เริ่มต้นในสถานะ ρ\rho บางสถานะ ในขณะที่ระบบพื้นที่ทำงาน W\mathsf{W} ถูกกำหนดค่าเริ่มต้นเป็นสถานะฐาน ∣0⟩\vert 0\rangle การดำเนินการยูนิทารี UU ถูกกระทำบนคู่ (W,X)(\mathsf{W},\mathsf{X}) และสุดท้ายระบบพื้นที่ทำงาน W\mathsf{W} ถูก trace out เหลือ X\mathsf{X} เป็นผลลัพธ์

สังเกตว่าเรากำลังสันนิษฐานว่า 00 เป็นสถานะคลาสสิกของ W\mathsf{W} และเราเลือกให้เป็นสถานะเริ่มต้นของระบบนี้ ซึ่งจะช่วยลดความซับซ้อนของคณิตศาสตร์ อย่างไรก็ตาม เราสามารถเลือกสถานะบริสุทธิ์คงที่ใด ๆ เพื่อแทนสถานะเริ่มต้นของ W\mathsf{W} โดยไม่เปลี่ยนแปลงคุณสมบัติพื้นฐานของการแทนค่า

นิพจน์ทางคณิตศาสตร์ของช่องผลลัพธ์ Φ\Phi มีดังนี้

Φ(ρ)=Tr⁡W(U(∣0⟩⟨0∣W⊗ρ)U†)\Phi(\rho) = \operatorname{Tr}_{\mathsf{W}} \bigl( U (\vert 0\rangle \langle 0 \vert_{\mathsf{W}} \otimes \rho) U^{\dagger} \bigr)

ตามที่เคยทำมา เราใช้ข้อตกลงการเรียงลำดับของ Qiskit: ระบบ X\mathsf{X} อยู่ด้านบนในแผนภาพ และดังนั้นสอดคล้องกับตัวประกอบ tensor ด้านขวาในสูตร

โดยทั่วไป ระบบนำเข้าและระบบนำออกของช่องไม่จำเป็นต้องเป็นระบบเดียวกัน นี่คือภาพแสดงการนำไปใช้งานของช่อง Φ\Phi ที่มีระบบนำเข้าเป็น X\mathsf{X} และระบบนำออกเป็น Y\mathsf{Y}

แผนภาพแสดงการแทนค่าแบบ Stinespring ของช่องที่ระบบนำเข้าและนำออกอาจต่างกัน

คราวนี้การดำเนินการยูนิทารีแปลง (W,X)(\mathsf{W},\mathsf{X}) เป็นคู่ (G,Y)(\mathsf{G},\mathsf{Y}) โดยที่ G\mathsf{G} เป็นระบบ "ขยะ" ที่ถูก trace out เหลือ Y\mathsf{Y} เป็นระบบผลลัพธ์ เพื่อให้ UU เป็นยูนิทารี มันต้องเป็นเมทริกซ์จัตุรัส ซึ่งกำหนดให้คู่ (G,Y)(\mathsf{G},\mathsf{Y}) มีจำนวนสถานะคลาสสิกเท่ากับคู่ (W,X)(\mathsf{W},\mathsf{X}) ดังนั้นระบบ W\mathsf{W} และ G\mathsf{G} ต้องถูกเลือกในลักษณะที่อนุญาตให้เป็นเช่นนี้ได้

เราได้นิพจน์ทางคณิตศาสตร์ของช่องผลลัพธ์ Φ\Phi ที่คล้ายกับที่เคยมีก่อนหน้า

Φ(ρ)=Tr⁡G(U(∣0⟩⟨0∣W⊗ρ)U†)\Phi(\rho) = \operatorname{Tr}_{\mathsf{G}} \bigl( U (\vert 0\rangle \langle 0 \vert_{\mathsf{W}} \otimes \rho) U^{\dagger} \bigr)

เมื่อช่องถูกอธิบายในลักษณะนี้ ในฐานะการดำเนินการยูนิทารีพร้อมกับข้อกำหนดว่าระบบพื้นที่ทำงานถูกกำหนดค่าเริ่มต้นอย่างไรและระบบนำออกถูกเลือกอย่างไร เราเรียกว่ามันถูกแสดงใน รูปแบบ Stinespring หรือเป็น การแทนค่าแบบ Stinespring ของช่อง

แม้จะไม่ใช่เรื่องที่เห็นได้ชัดเลย แต่ช่องทุกช่องมีการแทนค่าแบบ Stinespring จริง ๆ ดังที่เราจะเห็นเมื่อสิ้นสุดบทนี้ เราจะเห็นด้วยว่าการแทนค่าแบบ Stinespring ไม่ได้มีเพียงแบบเดียว จะมีวิธีต่าง ๆ ในการนำช่องเดียวกันไปใช้งานในลักษณะที่อธิบายไว้เสมอ

หมายเหตุ

ในบริบทของสารสนเทศเชิงควอนตัม คำว่า การแทนค่าแบบ Stinespring มักหมายถึงนิพจน์ทั่วไปกว่าเล็กน้อยของช่องในรูปแบบ

Φ(ρ)=Tr⁡G(AρA†)\Phi(\rho) = \operatorname{Tr}_{\mathsf{G}} \bigl( A \rho A^{\dagger} \bigr)

สำหรับ isometry AA ซึ่งเป็นเมทริกซ์ที่คอลัมน์ตั้งฉากกันตามบรรทัดฐาน แต่อาจไม่ใช่เมทริกซ์จัตุรัส สำหรับการแทนค่าแบบ Stinespring ที่มีรูปแบบที่เราใช้เป็นนิยาม เราสามารถได้นิพจน์ในรูปแบบอื่นนี้โดยการกำหนด

A=U(∣0⟩W⊗IX).A = U (\vert 0\rangle_{\mathsf{W}} \otimes \mathbb{I}_{\mathsf{X}}).

ช่อง dephasing แบบสมบูรณ์​

นี่คือการแทนค่าแบบ Stinespring ของช่อง dephasing ของ qubit Δ\Delta ในแผนภาพนี้ สายไฟทั้งสองแทน qubit เดี่ยว ดังนั้นนี่คือ circuit diagram ควอนตัมธรรมดา

แผนภาพ circuit ควอนตัมแทนช่อง dephasing แบบสมบูรณ์

เพื่อดูว่าผลของวงจรนี้ต่อ qubit นำเข้าถูกอธิบายโดยช่อง dephasing แบบสมบูรณ์จริง ๆ เราสามารถไล่ผ่านวงจรทีละขั้น โดยใช้การแทนด้วยเมทริกซ์ที่ชัดเจนของ partial trace ที่กล่าวถึงในบทที่แล้ว เราจะเรียก qubit บนสุดว่า X\mathsf{X} — นี่คือนำเข้าและนำออกของช่อง — และจะสมมติว่า X\mathsf{X} เริ่มต้นในสถานะโดยพลการ ρ\rho

ขั้นแรกคือการนำ qubit พื้นที่ทำงาน W\mathsf{W} เข้ามา ก่อนที่จะดำเนินการเกต controlled-NOT สถานะของคู่ (W,X)(\mathsf{W},\mathsf{X}) แทนด้วย density matrix ดังนี้

∣0⟩⟨0∣W⊗ρ=(1000)⊗(⟨0∣ρ∣0⟩⟨0∣ρ∣1⟩⟨1∣ρ∣0⟩⟨1∣ρ∣1⟩)=(⟨0∣ρ∣0⟩⟨0∣ρ∣1⟩00⟨1∣ρ∣0⟩⟨1∣ρ∣1⟩0000000000)\begin{aligned} \vert 0\rangle \langle 0 \vert_{\mathsf{W}} \otimes \rho & = \begin{pmatrix} 1 & 0\\ 0 & 0 \end{pmatrix} \otimes \begin{pmatrix} \langle 0 \vert \rho \vert 0 \rangle & \langle 0 \vert \rho \vert 1 \rangle\\[1mm] \langle 1 \vert \rho \vert 0 \rangle & \langle 1 \vert \rho \vert 1 \rangle \end{pmatrix}\\[4mm] & = \begin{pmatrix} \langle 0 \vert \rho \vert 0 \rangle & \langle 0 \vert \rho \vert 1 \rangle & 0 & 0 \\[1mm] \langle 1 \vert \rho \vert 0 \rangle & \langle 1 \vert \rho \vert 1 \rangle & 0 & 0 \\[1mm] 0 & 0 & 0 & 0 \\[1mm] 0 & 0 & 0 & 0 \end{pmatrix} \end{aligned}

ตามข้อตกลงการเรียงลำดับของ Qiskit qubit บนสุด X\mathsf{X} อยู่ทางขวาและ qubit ล่าง W\mathsf{W} อยู่ทางซ้าย เราใช้ density matrix แทน quantum state vector แต่มันถูก tensor เข้าด้วยกันในลักษณะที่คล้ายกับที่ทำในการกำหนดสูตรควอนตัมแบบง่าย

ขั้นต่อไปคือการดำเนินการ controlled-NOT โดยที่ X\mathsf{X} เป็น control และ W\mathsf{W} เป็น target ยังคำนึงถึงข้อตกลงการเรียงลำดับของ Qiskit การแทนด้วยเมทริกซ์ของเกตนี้เป็นดังนี้

(1000000100100100)\begin{pmatrix} 1 & 0 & 0 & 0\\ 0 & 0 & 0 & 1\\ 0 & 0 & 1 & 0\\ 0 & 1 & 0 & 0 \end{pmatrix}

นี่คือการดำเนินการยูนิทารี และในการนำไปใช้กับ density matrix เราจะ conjugate ด้วยเมทริกซ์ยูนิทารี conjugate-transpose ไม่เปลี่ยนแปลงเมทริกซ์นี้โดยเฉพาะ ดังนั้นผลลัพธ์เป็นดังนี้

(1000000100100100)(⟨0∣ρ∣0⟩⟨0∣ρ∣1⟩00⟨1∣ρ∣0⟩⟨1∣ρ∣1⟩0000000000)(1000000100100100)=(⟨0∣ρ∣0⟩00⟨0∣ρ∣1⟩00000000⟨1∣ρ∣0⟩00⟨1∣ρ∣1⟩)\begin{pmatrix} 1 & 0 & 0 & 0\\[1mm] 0 & 0 & 0 & 1\\[1mm] 0 & 0 & 1 & 0\\[1mm] 0 & 1 & 0 & 0 \end{pmatrix} \begin{pmatrix} \langle 0 \vert \rho \vert 0 \rangle & \langle 0 \vert \rho \vert 1 \rangle & 0 & 0 \\[1mm] \langle 1 \vert \rho \vert 0 \rangle & \langle 1 \vert \rho \vert 1 \rangle & 0 & 0 \\[1mm] 0 & 0 & 0 & 0 \\[1mm] 0 & 0 & 0 & 0 \end{pmatrix} \begin{pmatrix} 1 & 0 & 0 & 0\\[1mm] 0 & 0 & 0 & 1\\[1mm] 0 & 0 & 1 & 0\\[1mm] 0 & 1 & 0 & 0 \end{pmatrix}\\[3mm] = \begin{pmatrix} \langle 0 \vert \rho \vert 0 \rangle & 0 & 0 & \langle 0 \vert \rho \vert 1 \rangle\\[1mm] 0 & 0 & 0 & 0\\[1mm] 0 & 0 & 0 & 0\\[1mm] \langle 1 \vert \rho \vert 0 \rangle & 0 & 0 & \langle 1 \vert \rho \vert 1 \rangle \end{pmatrix}

สุดท้าย partial trace ถูกดำเนินการบน W\mathsf{W} จำการกระทำของการดำเนินการนี้บนเมทริกซ์ 4×44\times 4 ที่อธิบายในบทที่แล้ว เราได้ density matrix ผลลัพธ์ดังนี้

Tr⁡W(⟨0∣ρ∣0⟩00⟨0∣ρ∣1⟩00000000⟨1∣ρ∣0⟩00⟨1∣ρ∣1⟩)=(⟨0∣ρ∣0⟩000)+(000⟨1∣ρ∣1⟩)=(⟨0∣ρ∣0⟩00⟨1∣ρ∣1⟩)=Δ(ρ)\begin{aligned} \operatorname{Tr}_{\mathsf{W}} \begin{pmatrix} \langle 0 \vert \rho \vert 0 \rangle & 0 & 0 & \langle 0 \vert \rho \vert 1 \rangle\\[1mm] 0 & 0 & 0 & 0\\[1mm] 0 & 0 & 0 & 0\\[1mm] \langle 1 \vert \rho \vert 0 \rangle & 0 & 0 & \langle 1 \vert \rho \vert 1 \rangle \end{pmatrix} & = \begin{pmatrix} \langle 0 \vert \rho \vert 0 \rangle & 0 \\[1mm] 0 & 0 \end{pmatrix} + \begin{pmatrix} 0 & 0 \\[1mm] 0 & \langle 1 \vert \rho \vert 1 \rangle \end{pmatrix}\\[3mm] & = \begin{pmatrix} \langle 0 \vert \rho \vert 0 \rangle & 0 \\[1mm] 0 & \langle 1 \vert \rho \vert 1 \rangle \end{pmatrix}\\[4mm] & = \Delta(\rho) \end{aligned}

เราสามารถคำนวณ partial trace ทางเลือกได้โดยแปลงเป็น Dirac notation ก่อน

(⟨0∣ρ∣0⟩00⟨0∣ρ∣1⟩00000000⟨1∣ρ∣0⟩00⟨1∣ρ∣1⟩)=⟨0∣ρ∣0⟩ ∣0⟩⟨0∣⊗∣0⟩⟨0∣+ ⟨0∣ρ∣1⟩ ∣0⟩⟨1∣⊗∣0⟩⟨1∣+ ⟨1∣ρ∣0⟩ ∣1⟩⟨0∣⊗∣1⟩⟨0∣+ ⟨1∣ρ∣1⟩ ∣1⟩⟨1∣⊗∣1⟩⟨1∣\begin{pmatrix} \langle 0\vert \rho \vert 0\rangle & 0 & 0 & \langle 0\vert \rho \vert 1\rangle\\[1mm] 0 & 0 & 0 & 0\\[1mm] 0 & 0 & 0 & 0\\[1mm] \langle 1\vert \rho \vert 0\rangle & 0 & 0 & \langle 1\vert \rho \vert 1\rangle \end{pmatrix} = \begin{array}{r} \langle 0\vert \rho \vert 0\rangle \, \vert 0\rangle\langle 0\vert \otimes \vert 0\rangle\langle 0\vert \\[1mm] +\, \langle 0\vert \rho \vert 1\rangle \, \vert 0\rangle\langle 1\vert \otimes \vert 0\rangle\langle 1\vert \\[1mm] +\, \langle 1\vert \rho \vert 0\rangle \, \vert 1\rangle\langle 0\vert \otimes \vert 1\rangle\langle 0\vert \\[1mm] +\, \langle 1\vert \rho \vert 1\rangle \, \vert 1\rangle\langle 1\vert \otimes \vert 1\rangle\langle 1\vert \end{array}

การ trace out qubit ด้านซ้ายให้คำตอบเดิม

⟨0∣ρ∣0⟩ ∣0⟩⟨0∣+ ⟨1∣ρ∣1⟩ ∣1⟩⟨1∣=Δ(ρ)\langle 0\vert \rho \vert 0\rangle \, \vert 0\rangle\langle 0\vert +\, \langle 1\vert \rho \vert 1\rangle \, \vert 1\rangle\langle 1\vert = \Delta(\rho)

วิธีที่เข้าใจง่ายในการคิดถึงวงจรนี้คือการดำเนินการ controlled-NOT โดยผลแล้วก็คือการคัดลอกสถานะคลาสสิกของ qubit นำเข้า และเมื่อสำเนาถูกทิ้งในถังขยะ qubit นำเข้าจะ "ยุบตัว" แบบสุ่มตามความน่าจะเป็นไปเป็นหนึ่งในสองสถานะคลาสสิกที่เป็นไปได้ ซึ่งเทียบเท่ากับ dephasing แบบสมบูรณ์

ช่อง dephasing แบบสมบูรณ์ (แบบทางเลือก)​

วงจรที่อธิบายข้างต้นไม่ใช่วิธีเดียวในการนำช่อง dephasing แบบสมบูรณ์ไปใช้งาน นี่คือวิธีที่ต่างออกไป

แผนภาพ circuit ควอนตัมทางเลือกแทนช่อง dephasing แบบสมบูรณ์

นี่คือการวิเคราะห์สั้น ๆ ที่แสดงว่าการนำไปใช้งานนี้ทำงาน หลังจากดำเนินการเกต Hadamard เราได้สถานะสองQubitนี้เป็น density matrix:

∣+⟩⟨+∣⊗ρ=12(1111)⊗(⟨0∣ρ∣0⟩⟨0∣ρ∣1⟩⟨1∣ρ∣0⟩⟨1∣ρ∣1⟩)=12(⟨0∣ρ∣0⟩⟨0∣ρ∣1⟩⟨0∣ρ∣0⟩⟨0∣ρ∣1⟩⟨1∣ρ∣0⟩⟨1∣ρ∣1⟩⟨1∣ρ∣0⟩⟨1∣ρ∣1⟩⟨0∣ρ∣0⟩⟨0∣ρ∣1⟩⟨0∣ρ∣0⟩⟨0∣ρ∣1⟩⟨1∣ρ∣0⟩⟨1∣ρ∣1⟩⟨1∣ρ∣0⟩⟨1∣ρ∣1⟩).\begin{aligned} \vert + \rangle\langle + \vert \otimes \rho & = \frac{1}{2}\begin{pmatrix} 1 & 1\\[1mm] 1 & 1 \end{pmatrix} \otimes \begin{pmatrix} \langle 0 \vert \rho \vert 0 \rangle & \langle 0 \vert \rho \vert 1 \rangle\\[1mm] \langle 1 \vert \rho \vert 0 \rangle & \langle 1 \vert \rho \vert 1 \rangle \end{pmatrix}\\[4mm] & = \frac{1}{2} \begin{pmatrix} \langle 0 \vert \rho \vert 0 \rangle & \langle 0 \vert \rho \vert 1 \rangle & \langle 0 \vert \rho \vert 0 \rangle & \langle 0 \vert \rho \vert 1 \rangle\\[1mm] \langle 1 \vert \rho \vert 0 \rangle & \langle 1 \vert \rho \vert 1 \rangle & \langle 1 \vert \rho \vert 0 \rangle & \langle 1 \vert \rho \vert 1 \rangle\\[1mm] \langle 0 \vert \rho \vert 0 \rangle & \langle 0 \vert \rho \vert 1 \rangle & \langle 0 \vert \rho \vert 0 \rangle & \langle 0 \vert \rho \vert 1 \rangle\\[1mm] \langle 1 \vert \rho \vert 0 \rangle & \langle 1 \vert \rho \vert 1 \rangle & \langle 1 \vert \rho \vert 0 \rangle & \langle 1 \vert \rho \vert 1 \rangle \end{pmatrix}. \end{aligned}

เกต controlled-σz\sigma_z ดำเนินการโดย conjugation ดังนี้

12(100001000010000−1)(⟨0∣ρ∣0⟩⟨0∣ρ∣1⟩⟨0∣ρ∣0⟩⟨0∣ρ∣1⟩⟨1∣ρ∣0⟩⟨1∣ρ∣1⟩⟨1∣ρ∣0⟩⟨1∣ρ∣1⟩⟨0∣ρ∣0⟩⟨0∣ρ∣1⟩⟨0∣ρ∣0⟩⟨0∣ρ∣1⟩⟨1∣ρ∣0⟩⟨1∣ρ∣1⟩⟨1∣ρ∣0⟩⟨1∣ρ∣1⟩)(100001000010000−1)=12(⟨0∣ρ∣0⟩⟨0∣ρ∣1⟩⟨0∣ρ∣0⟩−⟨0∣ρ∣1⟩⟨1∣ρ∣0⟩⟨1∣ρ∣1⟩⟨1∣ρ∣0⟩−⟨1∣ρ∣1⟩⟨0∣ρ∣0⟩⟨0∣ρ∣1⟩⟨0∣ρ∣0⟩−⟨0∣ρ∣1⟩−⟨1∣ρ∣0⟩−⟨1∣ρ∣1⟩−⟨1∣ρ∣0⟩⟨1∣ρ∣1⟩)\frac{1}{2} \begin{pmatrix} 1 & 0 & 0 & 0\\[1mm] 0 & 1 & 0 & 0\\[1mm] 0 & 0 & 1 & 0\\[1mm] 0 & 0 & 0 & -1 \end{pmatrix} \begin{pmatrix} \langle 0 \vert \rho \vert 0 \rangle & \langle 0 \vert \rho \vert 1 \rangle & \langle 0 \vert \rho \vert 0 \rangle & \langle 0 \vert \rho \vert 1 \rangle\\[1mm] \langle 1 \vert \rho \vert 0 \rangle & \langle 1 \vert \rho \vert 1 \rangle & \langle 1 \vert \rho \vert 0 \rangle & \langle 1 \vert \rho \vert 1 \rangle\\[1mm] \langle 0 \vert \rho \vert 0 \rangle & \langle 0 \vert \rho \vert 1 \rangle & \langle 0 \vert \rho \vert 0 \rangle & \langle 0 \vert \rho \vert 1 \rangle\\[1mm] \langle 1 \vert \rho \vert 0 \rangle & \langle 1 \vert \rho \vert 1 \rangle & \langle 1 \vert \rho \vert 0 \rangle & \langle 1 \vert \rho \vert 1 \rangle \end{pmatrix} \begin{pmatrix} 1 & 0 & 0 & 0\\[1mm] 0 & 1 & 0 & 0\\[1mm] 0 & 0 & 1 & 0\\[1mm] 0 & 0 & 0 & -1 \end{pmatrix}\\[3mm] = \frac{1}{2} \begin{pmatrix} \langle 0 \vert \rho \vert 0 \rangle & \langle 0 \vert \rho \vert 1 \rangle & \langle 0 \vert \rho \vert 0 \rangle & -\langle 0 \vert \rho \vert 1 \rangle\\[1mm] \langle 1 \vert \rho \vert 0 \rangle & \langle 1 \vert \rho \vert 1 \rangle & \langle 1 \vert \rho \vert 0 \rangle & -\langle 1 \vert \rho \vert 1 \rangle\\[1mm] \langle 0 \vert \rho \vert 0 \rangle & \langle 0 \vert \rho \vert 1 \rangle & \langle 0 \vert \rho \vert 0 \rangle & -\langle 0 \vert \rho \vert 1 \rangle\\[1mm] -\langle 1 \vert \rho \vert 0 \rangle & -\langle 1 \vert \rho \vert 1 \rangle & -\langle 1 \vert \rho \vert 0 \rangle & \langle 1 \vert \rho \vert 1 \rangle \end{pmatrix}

สุดท้ายระบบพื้นที่ทำงาน W\mathsf{W} ถูก trace out

12Tr⁡W(⟨0∣ρ∣0⟩⟨0∣ρ∣1⟩⟨0∣ρ∣0⟩−⟨0∣ρ∣1⟩⟨1∣ρ∣0⟩⟨1∣ρ∣1⟩⟨1∣ρ∣0⟩−⟨1∣ρ∣1⟩⟨0∣ρ∣0⟩⟨0∣ρ∣1⟩⟨0∣ρ∣0⟩−⟨0∣ρ∣1⟩−⟨1∣ρ∣0⟩−⟨1∣ρ∣1⟩−⟨1∣ρ∣0⟩⟨1∣ρ∣1⟩)=12(⟨0∣ρ∣0⟩⟨0∣ρ∣1⟩⟨1∣ρ∣0⟩⟨1∣ρ∣1⟩)+12(⟨0∣ρ∣0⟩−⟨0∣ρ∣1⟩−⟨1∣ρ∣0⟩⟨1∣ρ∣1⟩)=(⟨0∣ρ∣0⟩00⟨1∣ρ∣1⟩)\frac{1}{2} \operatorname{Tr}_{\mathsf{W}} \begin{pmatrix} \langle 0 \vert \rho \vert 0 \rangle & \langle 0 \vert \rho \vert 1 \rangle & \langle 0 \vert \rho \vert 0 \rangle & -\langle 0 \vert \rho \vert 1 \rangle\\[1mm] \langle 1 \vert \rho \vert 0 \rangle & \langle 1 \vert \rho \vert 1 \rangle & \langle 1 \vert \rho \vert 0 \rangle & -\langle 1 \vert \rho \vert 1 \rangle\\[1mm] \langle 0 \vert \rho \vert 0 \rangle & \langle 0 \vert \rho \vert 1 \rangle & \langle 0 \vert \rho \vert 0 \rangle & -\langle 0 \vert \rho \vert 1 \rangle\\[1mm] -\langle 1 \vert \rho \vert 0 \rangle & -\langle 1 \vert \rho \vert 1 \rangle & -\langle 1 \vert \rho \vert 0 \rangle & \langle 1 \vert \rho \vert 1 \rangle \end{pmatrix}\\[3mm] \begin{aligned} & = \frac{1}{2} \begin{pmatrix} \langle 0 \vert \rho \vert 0 \rangle & \langle 0 \vert \rho \vert 1 \rangle\\[1mm] \langle 1 \vert \rho \vert 0 \rangle & \langle 1 \vert \rho \vert 1 \rangle \end{pmatrix} + \frac{1}{2} \begin{pmatrix} \langle 0 \vert \rho \vert 0 \rangle & -\langle 0 \vert \rho \vert 1 \rangle\\[1mm] -\langle 1 \vert \rho \vert 0 \rangle & \langle 1 \vert \rho \vert 1 \rangle \end{pmatrix}\\[4mm] & = \begin{pmatrix} \langle 0 \vert \rho \vert 0 \rangle & 0\\[2mm] 0 & \langle 1 \vert \rho \vert 1 \rangle \end{pmatrix} \end{aligned}

การนำไปใช้งานนี้มีพื้นฐานมาจากแนวคิดง่าย ๆ: dephasing เทียบเท่ากับการไม่ทำอะไร (นั่นคือ ใช้การดำเนินการ identity) หรือการใช้เกต σz\sigma_z แต่ละอย่างด้วยความน่าจะเป็น 1/21/2

12ρ+12σzρσz=12(⟨0∣ρ∣0⟩⟨0∣ρ∣1⟩⟨1∣ρ∣0⟩⟨1∣ρ∣1⟩)+12(⟨0∣ρ∣0⟩−⟨0∣ρ∣1⟩−⟨1∣ρ∣0⟩⟨1∣ρ∣1⟩)=(⟨0∣ρ∣0⟩00⟨1∣ρ∣1⟩)=Δ(ρ)\begin{aligned} \frac{1}{2} \rho + \frac{1}{2} \sigma_z \rho \sigma_z & = \frac{1}{2} \begin{pmatrix} \langle 0 \vert \rho \vert 0 \rangle & \langle 0 \vert \rho \vert 1 \rangle\\[1mm] \langle 1 \vert \rho \vert 0 \rangle & \langle 1 \vert \rho \vert 1 \rangle \end{pmatrix} + \frac{1}{2} \begin{pmatrix} \langle 0 \vert \rho \vert 0 \rangle & -\langle 0 \vert \rho \vert 1 \rangle\\[1mm] -\langle 1 \vert \rho \vert 0 \rangle & \langle 1 \vert \rho \vert 1 \rangle \end{pmatrix}\\[4mm] & = \begin{pmatrix} \langle 0 \vert \rho \vert 0 \rangle & 0\\[1mm] 0 & \langle 1 \vert \rho \vert 1 \rangle \end{pmatrix}\\[2mm] & = \Delta(\rho) \end{aligned}

นั่นคือ ช่อง dephasing แบบสมบูรณ์เป็นตัวอย่างของ mixed-unitary channel และโดยเฉพาะเจาะจงคือ Pauli channel

ช่อง reset ของ qubit​

ช่อง reset ของ qubit สามารถนำไปใช้งานได้ดังนี้

แผนภาพ circuit ควอนตัมแทนช่อง reset ของ qubit

เกต swap เพียงแค่เลื่อนสถานะ ∣0⟩\vert 0\rangle ที่กำหนดค่าเริ่มต้นของ qubit พื้นที่ทำงานเพื่อให้มันได้รับการนำออก ในขณะที่สถานะนำเข้า ρ\rho ถูกย้ายไปยัง qubit ล่างแล้วถูก trace out

หรืออีกทางหนึ่ง ถ้าเราไม่ต้องการให้ผลลัพธ์ของช่องอยู่ด้านบน เราสามารถใช้วงจรที่ง่ายมากนี้เป็นการแทนค่าของเรา

แผนภาพ circuit ควอนตัมทางเลือกแทนช่อง reset ของ qubit

อธิบายเป็นคำพูด การ reset qubit ไปยังสถานะ ∣0⟩\vert 0\rangle เทียบเท่ากับการโยน qubit ในถังขยะแล้วรับ qubit ใหม่

การแทนค่าแบบ Kraus​

ต่อไปเราจะพูดถึง การแทนค่าแบบ Kraus ซึ่งเสนอวิธีสูตรที่สะดวกในการแสดงการกระทำของช่องผ่านการคูณและบวกเมทริกซ์ โดยเฉพาะ การแทนค่าแบบ Kraus คือข้อกำหนดของช่อง Φ\Phi ในรูปแบบดังนี้

Φ(ρ)=∑k=0N−1AkρAk†\Phi(\rho) = \sum_{k = 0}^{N-1} A_k \rho A_k^{\dagger}

ที่นี่ A0,…,AN−1A_0,\ldots,A_{N-1} เป็นเมทริกซ์ที่มีขนาดเหมือนกันทั้งหมด: คอลัมน์ของมันสอดคล้องกับสถานะคลาสสิกของระบบนำเข้า X\mathsf{X} และแถวของมันสอดคล้องกับสถานะคลาสสิกของระบบนำออก ไม่ว่าจะเป็น X\mathsf{X} หรือระบบอื่น Y\mathsf{Y} เพื่อให้ Φ\Phi เป็นช่องที่ถูกต้อง เมทริกซ์เหล่านี้ต้องตรงตามเงื่อนไขดังนี้

∑k=0N−1Ak†Ak=IX\sum_{k = 0}^{N-1} A_k^{\dagger} A_k = \mathbb{I}_{\mathsf{X}}

เงื่อนไขนี้เทียบเท่ากับเงื่อนไขที่ Φ\Phi รักษา trace คุณสมบัติอีกอย่างที่จำเป็นสำหรับช่อง ซึ่งก็คือ complete positivity นั้นตามมาจากรูปแบบทั่วไปของสมการสำหรับ Φ\Phi ในฐานะผลรวมของการ conjugate

บางครั้งการตั้งชื่อเมทริกซ์ A0,…,AN−1A_0,\ldots,A_{N-1} ต่างออกไปก็สะดวก ตัวอย่างเช่น เราอาจเริ่มนับจาก 11 หรือใช้สถานะในเซตสถานะคลาสสิกโดยพลการ Γ\Gamma แทนตัวเลขเป็น subscript:

Φ(ρ)=∑a∈ΓAaρAa†where∑a∈ΓAa†Aa=I.\Phi(\rho) = \sum_{a\in\Gamma} A_a \rho A_a^{\dagger} \quad \text{where} \quad \sum_{a\in\Gamma} A_a^{\dagger} A_a = \mathbb{I}.

วิธีต่าง ๆ ในการตั้งชื่อเมทริกซ์เหล่านี้ ซึ่งเรียกว่า เมทริกซ์ Kraus ล้วนเป็นที่นิยมและสะดวกในสถานการณ์ที่แตกต่างกัน แต่เราจะยึดชื่อ A0,…,AN−1A_0,\ldots,A_{N-1} ในบทนี้เพื่อความเรียบง่าย

จำนวน NN สามารถเป็นจำนวนเต็มบวกโดยพลการ แต่ไม่จำเป็นต้องมากเกินไป: ถ้าระบบนำเข้า X\mathsf{X} มี nn สถานะคลาสสิกและระบบนำออก Y\mathsf{Y} มี mm สถานะคลาสสิก แล้วช่องใด ๆ จาก X\mathsf{X} ไปยัง Y\mathsf{Y} จะมีการแทนค่าแบบ Kraus เสมอที่ NN มีค่าไม่เกินผลคูณ nmnm

ช่อง dephasing แบบสมบูรณ์​

เราได้การแทนค่าแบบ Kraus ของช่อง dephasing แบบสมบูรณ์โดยการกำหนด A0=∣0⟩⟨0∣A_0 = \vert 0\rangle\langle 0\vert และ A1=∣1⟩⟨1∣A_1 = \vert 1\rangle\langle 1\vert

∑k=01AkρAk†=∣0⟩⟨0∣ρ∣0⟩⟨0∣+∣1⟩⟨1∣ρ∣1⟩⟨1∣=⟨0∣ρ∣0⟩ ∣0⟩⟨0∣+⟨1∣ρ∣1⟩ ∣1⟩⟨1∣=(⟨0∣ρ∣0⟩00⟨1∣ρ∣1⟩)\begin{aligned} \sum_{k = 0}^1 A_k \rho A_k^{\dagger} & = \vert 0\rangle\langle 0 \vert \rho \vert 0\rangle\langle 0 \vert + \vert 1\rangle\langle 1 \vert \rho \vert 1\rangle\langle 1 \vert\\ & = \langle 0 \vert \rho \vert 0\rangle \, \vert 0\rangle\langle 0 \vert + \langle 1 \vert \rho \vert 1\rangle \, \vert 1\rangle\langle 1 \vert \\[2mm] & = \begin{pmatrix} \langle 0 \vert \rho \vert 0 \rangle & 0 \\[1mm] 0 & \langle 1 \vert \rho \vert 1 \rangle \end{pmatrix} \end{aligned}

เมทริกซ์เหล่านี้ตรงตามเงื่อนไขที่จำเป็น

∑k=01Ak†Ak=∣0⟩⟨0∣0⟩⟨0∣+∣1⟩⟨1∣1⟩⟨1∣=∣0⟩⟨0∣+∣1⟩⟨1∣=I\sum_{k = 0}^1 A_k^{\dagger} A_k = \vert 0\rangle\langle 0\vert 0\rangle\langle 0\vert + \vert 1\rangle\langle 1\vert 1\rangle\langle 1\vert = \vert 0\rangle\langle 0\vert + \vert 1\rangle\langle 1\vert = \mathbb{I}

หรืออีกทางหนึ่ง เราสามารถกำหนด A0=12IA_0 = \frac{1}{\sqrt{2}}\mathbb{I} และ A1=12σzA_1 = \frac{1}{\sqrt{2}}\sigma_z ดังนั้น

∑k=01AkρAk†=12ρ+12σzρσz=Δ(ρ),\sum_{k = 0}^1 A_k \rho A_k^{\dagger} = \frac{1}{2} \rho + \frac{1}{2} \sigma_z \rho \sigma_z = \Delta(\rho),

ดังที่คำนวณไว้ก่อนหน้า คราวนี้สามารถตรวจสอบเงื่อนไขที่จำเป็นได้ดังนี้

∑k=01Ak†Ak=12I+12σz2=12I+12I=I\sum_{k = 0}^1 A_k^{\dagger} A_k = \frac{1}{2} \mathbb{I} + \frac{1}{2} \sigma_z^2 = \frac{1}{2} \mathbb{I} + \frac{1}{2} \mathbb{I} = \mathbb{I}

ช่อง reset ของ qubit​

เราได้การแทนค่าแบบ Kraus ของช่อง reset ของ qubit โดยการกำหนด A0=∣0⟩⟨0∣A_0 = \vert 0\rangle\langle 0\vert และ A1=∣0⟩⟨1∣A_1 = \vert 0\rangle\langle 1\vert

∑k=01AkρAk†=∣0⟩⟨0∣ρ∣0⟩⟨0∣+∣0⟩⟨1∣ρ∣1⟩⟨0∣=⟨0∣ρ∣0⟩ ∣0⟩⟨0∣+⟨1∣ρ∣1⟩ ∣0⟩⟨0∣=Tr⁡(ρ)∣0⟩⟨0∣\begin{aligned} \sum_{k = 0}^1 A_k \rho A_k^{\dagger} & = \vert 0\rangle\langle 0 \vert \rho \vert 0\rangle\langle 0 \vert + \vert 0\rangle\langle 1 \vert \rho \vert 1\rangle\langle 0 \vert\\ & = \langle 0 \vert \rho \vert 0\rangle \, \vert 0\rangle\langle 0 \vert + \langle 1 \vert \rho \vert 1\rangle \, \vert 0\rangle\langle 0 \vert\\[2mm] & = \operatorname{Tr}(\rho) \vert 0\rangle \langle 0 \vert \end{aligned}

เมทริกซ์เหล่านี้ตรงตามเงื่อนไขที่จำเป็น

∑k=01Ak†Ak=∣0⟩⟨0∣0⟩⟨0∣+∣1⟩⟨0∣0⟩⟨1∣=∣0⟩⟨0∣+∣1⟩⟨1∣=I\sum_{k = 0}^1 A_k^{\dagger} A_k = \vert 0\rangle\langle 0\vert 0\rangle\langle 0\vert + \vert 1\rangle\langle 0\vert 0\rangle\langle 1\vert = \vert 0\rangle\langle 0\vert + \vert 1\rangle\langle 1\vert = \mathbb{I}

ช่อง depolarizing แบบสมบูรณ์​

วิธีหนึ่งในการได้การแทนค่าแบบ Kraus สำหรับช่อง depolarizing แบบสมบูรณ์คือการเลือกเมทริกซ์ Kraus A0,…,A3A_0,\ldots,A_3 ดังนี้

A0=∣0⟩⟨0∣2A1=∣0⟩⟨1∣2A2=∣1⟩⟨0∣2A3=∣1⟩⟨1∣2A_0 = \frac{\vert 0\rangle\langle 0\vert}{\sqrt{2}} \quad A_1 = \frac{\vert 0\rangle\langle 1\vert}{\sqrt{2}} \quad A_2 = \frac{\vert 1\rangle\langle 0\vert}{\sqrt{2}} \quad A_3 = \frac{\vert 1\rangle\langle 1\vert}{\sqrt{2}}

สำหรับ density matrix qubit ρ\rho ใด ๆ เรามี

∑k=03AkρAk†=12(∣0⟩⟨0∣ρ∣0⟩⟨0∣+∣0⟩⟨1∣ρ∣1⟩⟨0∣+∣1⟩⟨0∣ρ∣0⟩⟨1∣+∣1⟩⟨1∣ρ∣1⟩⟨1∣)=Tr⁡(ρ)I2=Ω(ρ).\begin{aligned} \sum_{k = 0}^3 A_k \rho A_k^{\dagger} & = \frac{1}{2} \bigl(\vert 0\rangle\langle 0\vert \rho \vert 0\rangle\langle 0\vert + \vert 0\rangle\langle 1\vert \rho \vert 1\rangle\langle 0\vert + \vert 1\rangle\langle 0\vert \rho \vert 0\rangle\langle 1\vert + \vert 1\rangle\langle 1\vert \rho \vert 1\rangle\langle 1\vert\bigr)\\ & = \operatorname{Tr}(\rho) \frac{\mathbb{I}}{2}\\[1mm] & = \Omega(\rho). \end{aligned}

การแทนค่าแบบ Kraus ทางเลือกได้จากการเลือกเมทริกซ์ Kraus ดังนี้

A0=I2A1=σx2A2=σy2A3=σz2A_0 = \frac{\mathbb{I}}{2} \quad A_1 = \frac{\sigma_x}{2} \quad A_2 = \frac{\sigma_y}{2} \quad A_3 = \frac{\sigma_z}{2}

เพื่อยืนยันความถูกต้องของการแทนค่าแบบ Kraus เหล่านี้ว่าแทนช่อง depolarizing แบบสมบูรณ์จริง ๆ มาสังเกตก่อนว่าการ conjugate เมทริกซ์ 2×22\times 2 โดยพลการด้วยเมทริกซ์ Pauli ทำงานดังนี้

σx(α0,0α0,1α1,0α1,1)σx=(α1,1α1,0α0,1α0,0)σy(α0,0α0,1α1,0α1,1)σy=(α1,1−α1,0−α0,1α0,0)σz(α0,0α0,1α1,0α1,1)σz=(α0,0−α0,1−α1,0α1,1)\begin{aligned} \sigma_x \begin{pmatrix} \alpha_{0,0} & \alpha_{0,1}\\[1mm] \alpha_{1,0} & \alpha_{1,1} \end{pmatrix} \sigma_x & = \begin{pmatrix} \alpha_{1,1} & \alpha_{1,0}\\[1mm] \alpha_{0,1} & \alpha_{0,0} \end{pmatrix}\\[5mm] \sigma_y \begin{pmatrix} \alpha_{0,0} & \alpha_{0,1}\\[1mm] \alpha_{1,0} & \alpha_{1,1} \end{pmatrix} \sigma_y & = \begin{pmatrix} \alpha_{1,1} & -\alpha_{1,0}\\[1mm] -\alpha_{0,1} & \alpha_{0,0} \end{pmatrix}\\[5mm] \sigma_z \begin{pmatrix} \alpha_{0,0} & \alpha_{0,1}\\[1mm] \alpha_{1,0} & \alpha_{1,1} \end{pmatrix} \sigma_z & = \begin{pmatrix} \alpha_{0,0} & -\alpha_{0,1}\\[1mm] -\alpha_{1,0} & \alpha_{1,1} \end{pmatrix} \end{aligned}

ซึ่งช่วยให้เราตรวจสอบความถูกต้องของการแทนค่าแบบ Kraus ของเรา

∑k=03AkρAk†=ρ+σxρσx+σyρσy+σzρσz4=14(⟨0∣ρ∣0⟩+⟨1∣ρ∣1⟩+⟨1∣ρ∣1⟩+⟨0∣ρ∣0⟩⟨0∣ρ∣1⟩+⟨1∣ρ∣0⟩−⟨1∣ρ∣0⟩−⟨0∣ρ∣1⟩⟨1∣ρ∣0⟩+⟨0∣ρ∣1⟩−⟨0∣ρ∣1⟩−⟨1∣ρ∣0⟩⟨1∣ρ∣1⟩+⟨0∣ρ∣0⟩+⟨0∣ρ∣0⟩+⟨1∣ρ∣1⟩)=Tr⁡(ρ)I2\begin{aligned} \sum_{k = 0}^3 A_k \rho A_k^{\dagger} & = \frac{\rho + \sigma_x \rho \sigma_x + \sigma_y \rho \sigma_y + \sigma_z \rho \sigma_z}{4} \\ & = \frac{1}{4} \begin{pmatrix} \langle 0\vert\rho\vert 0\rangle + \langle 1\vert\rho\vert 1\rangle + \langle 1\vert\rho\vert 1\rangle + \langle 0\vert\rho\vert 0\rangle & \langle 0\vert\rho\vert 1\rangle + \langle 1\vert\rho\vert 0\rangle - \langle 1\vert\rho\vert 0\rangle - \langle 0\vert\rho\vert 1\rangle \\[2mm] \langle 1\vert\rho\vert 0\rangle + \langle 0\vert\rho\vert 1\rangle - \langle 0\vert\rho\vert 1\rangle - \langle 1\vert\rho\vert 0\rangle & \langle 1\vert\rho\vert 1\rangle + \langle 0\vert\rho\vert 0\rangle + \langle 0\vert\rho\vert 0\rangle + \langle 1\vert\rho\vert 1\rangle \end{pmatrix} \\[4mm] & = \operatorname{Tr}(\rho) \frac{\mathbb{I}}{2} \end{aligned}

การแทนค่าแบบ Kraus นี้แสดงแนวคิดสำคัญว่าสถานะของ qubit สามารถสุ่มได้อย่างสมบูรณ์โดยการใช้เมทริกซ์ Pauli หนึ่งในสี่ตัว (รวมถึงเมทริกซ์ identity) ที่เลือกแบบ uniform at random ดังนั้น ช่อง depolarizing แบบสมบูรณ์จึงเป็นอีกตัวอย่างหนึ่งของ Pauli channel

ไม่สามารถหาการแทนค่าแบบ Kraus สำหรับช่อง depolarizing แบบสมบูรณ์ Ω\Omega ที่มีเมทริกซ์ Kraus สามตัวหรือน้อยกว่าได้ ต้องมีอย่างน้อยสี่ตัวสำหรับช่องนี้

ช่องยูนิทารี​

ถ้าเรามีเมทริกซ์ยูนิทารี UU แทนการดำเนินการบนระบบ X\mathsf{X} เราสามารถแสดงการกระทำของการดำเนินการยูนิทารีนี้เป็นช่องได้:

Φ(ρ)=UρU†.\Phi(\rho) = U \rho U^{\dagger}.

นิพจน์นี้เป็นการแทนค่าแบบ Kraus ที่ถูกต้องแล้วของช่อง Φ\Phi ที่เราบังเอิญมีเมทริกซ์ Kraus เพียงตัวเดียว A0=UA_0 = U ในกรณีนี้ เงื่อนไขที่จำเป็น

∑k=0N−1Ak†Ak=IX\sum_{k = 0}^{N-1} A_k^{\dagger} A_k = \mathbb{I}_{\mathsf{X}}

กลายเป็นรูปแบบที่ง่ายกว่ามากคือ U†U=IXU^{\dagger} U = \mathbb{I}_{\mathsf{X}} ซึ่งเรารู้ว่าเป็นจริงเพราะ UU เป็นยูนิทารี

การแทนค่าแบบ Choi​

ต่อไปเราจะพูดถึงวิธีที่สามที่ช่องสามารถอธิบายได้ ผ่าน การแทนค่าแบบ Choi วิธีที่มันทำงานคือแต่ละช่องถูกแทนด้วยเมทริกซ์เดี่ยวที่รู้จักกันในชื่อ เมทริกซ์ Choi ถ้าระบบนำเข้ามี nn สถานะคลาสสิกและระบบนำออกมี mm สถานะคลาสสิก แล้วเมทริกซ์ Choi ของช่องจะมี nmnm แถวและ nmnm คอลัมน์

เมทริกซ์ Choi ให้การแทนค่าแบบ faithful ของช่อง หมายความว่าสองช่องเหมือนกันก็ต่อเมื่อมีเมทริกซ์ Choi เหมือนกัน เหตุผลหนึ่งที่สำคัญคือมันให้วิธีในการพิจารณาว่าคำอธิบายที่แตกต่างกันสองแบบสอดคล้องกับช่องเดียวกันหรือช่องที่ต่างกัน: เราเพียงแค่คำนวณเมทริกซ์ Choi แล้วเปรียบเทียบว่าเท่ากันหรือไม่ ในทางตรงกันข้าม การแทนค่าแบบ Stinespring และ Kraus ไม่ได้มีเพียงแบบเดียวในลักษณะนี้ ดังที่เราเห็นมาแล้ว

เมทริกซ์ Choi ยังมีประโยชน์ในการเปิดเผยคุณสมบัติทางคณิตศาสตร์ต่าง ๆ ของช่องอีกด้วย

นิยาม​

ให้ Φ\Phi เป็นช่องจากระบบ X\mathsf{X} ไปยังระบบ Y\mathsf{Y} และสมมติว่าเซตสถานะคลาสสิกของระบบนำเข้า X\mathsf{X} คือ Σ\Sigma การแทนค่าแบบ Choi ของ Φ\Phi ซึ่งแทนด้วย J(Φ)J(\Phi) นิยามด้วยสมการดังนี้

J(Φ)=∑a,b∈Σ∣a⟩⟨b∣⊗Φ(∣a⟩⟨b∣)J(\Phi) = \sum_{a,b\in\Sigma} \vert a\rangle\langle b \vert \otimes \Phi\bigl( \vert a\rangle\langle b \vert\bigr)

ถ้าเราสมมติว่า Σ={0,…,n−1}\Sigma = \{0,\ldots, n-1\} สำหรับจำนวนเต็มบวก nn บางตัว แล้วเราสามารถแสดง J(Φ)J(\Phi) เป็น block matrix แทนได้:

J(Φ)=(Φ(∣0⟩⟨0∣)Φ(∣0⟩⟨1∣)⋯Φ(∣0⟩⟨n−1∣)Φ(∣1⟩⟨0∣)Φ(∣1⟩⟨1∣)⋯Φ(∣1⟩⟨n−1∣)⋮⋮⋱⋮Φ(∣n−1⟩⟨0∣)Φ(∣n−1⟩⟨1∣)⋯Φ(∣n−1⟩⟨n−1∣))J(\Phi) = \begin{pmatrix} \Phi\bigl(\vert 0\rangle\langle 0\vert\bigr) & \Phi\bigl(\vert 0\rangle\langle 1\vert\bigr) & \cdots & \Phi\bigl(\vert 0\rangle\langle n-1\vert\bigr) \\[1mm] \Phi\bigl(\vert 1\rangle\langle 0\vert\bigr) & \Phi\bigl(\vert 1\rangle\langle 1\vert\bigr) & \cdots & \Phi\bigl(\vert 1\rangle\langle n-1\vert\bigr) \\[1mm] \vdots & \vdots & \ddots & \vdots\\[1mm] \Phi\bigl(\vert n-1\rangle\langle 0\vert\bigr) & \Phi\bigl(\vert n-1\rangle\langle 1\vert\bigr) & \cdots & \Phi\bigl(\vert n-1\rangle\langle n-1\vert\bigr) \end{pmatrix}

นั่นคือ เมื่อแสดงเป็น block matrix เมทริกซ์ Choi ของช่องมีหนึ่ง block Φ(∣a⟩⟨b∣)\Phi(\vert a\rangle\langle b\vert) สำหรับแต่ละคู่ (a,b)(a,b) ของสถานะคลาสสิกของระบบนำเข้า โดย block เรียงในลักษณะที่เป็นธรรมชาติ

สังเกตว่าเซต {∣a⟩⟨b∣ : 0≤a,b<n}\{\vert a\rangle\langle b\vert \,:\, 0\leq a,b < n\} เป็นฐานสำหรับปริภูมิของเมทริกซ์ n×nn\times n ทั้งหมด เนื่องจาก Φ\Phi เป็นเชิงเส้น จึงตามมาว่าการกระทำของมันสามารถกู้คืนได้จากเมทริกซ์ Choi โดยการรวม block เชิงเส้น

สถานะ Choi ของช่อง​

อีกวิธีในการคิดเกี่ยวกับเมทริกซ์ Choi ของช่องคือมันเป็น density matrix ถ้าเราหารด้วย n=∣Σ∣n = \vert\Sigma\vert มาโฟกัสที่สถานการณ์ที่ Σ={0,…,n−1}\Sigma = \{0,\ldots,n-1\} เพื่อความเรียบง่าย และจินตนาการว่าเรามีสำเนาของ X\mathsf{X} สองชุดที่เหมือนกันซึ่งอยู่ในสถานะ entangled ร่วมกัน

∣ψ⟩=1n∑a=0n−1∣a⟩⊗∣a⟩.\vert \psi \rangle = \frac{1}{\sqrt{n}} \sum_{a = 0}^{n-1} \vert a \rangle \otimes \vert a \rangle.

เมื่อแสดงเป็น density matrix สถานะนี้เป็นดังนี้

∣ψ⟩⟨ψ∣=1n∑a,b=0n−1∣a⟩⟨b∣⊗∣a⟩⟨b∣\vert \psi \rangle \langle \psi \vert = \frac{1}{n} \sum_{a,b = 0}^{n-1} \vert a\rangle\langle b \vert \otimes \vert a\rangle\langle b \vert

ถ้าเราใช้ Φ\Phi กับสำเนา X\mathsf{X} ด้านขวา เราได้เมทริกซ์ Choi หารด้วย nn

(Id⁡⊗ Φ)(∣ψ⟩⟨ψ∣)=1n∑a,b=0n−1∣a⟩⟨b∣⊗Φ(∣a⟩⟨b∣)=J(Φ)n(\operatorname{Id}\otimes \,\Phi) \bigl(\vert \psi \rangle \langle \psi \vert\bigr) = \frac{1}{n} \sum_{a,b = 0}^{n-1} \vert a\rangle\langle b \vert \otimes \Phi\bigl(\vert a\rangle\langle b \vert\bigr) = \frac{J(\Phi)}{n}

อธิบายเป็นคำพูด หากไม่นับตัวประกอบ normalization 1/n1/n เมทริกซ์ Choi ของ Φ\Phi คือ density matrix ที่เราได้จากการประเมิน Φ\Phi บนครึ่งหนึ่งของคู่ระบบนำเข้าที่ entangle แบบสูงสุด ดังที่ภาพต่อไปนี้แสดง

แผนภาพแสดงสถานะ Choi ของช่อง

สังเกตโดยเฉพาะว่าสิ่งนี้หมายความว่าเมทริกซ์ Choi ของช่องต้องเป็น positive semidefinite เสมอ

เราเห็นด้วยว่า เนื่องจากช่อง Φ\Phi ถูกใช้กับระบบขวา/บนเพียงอย่างเดียว มันไม่สามารถส่งผลต่อ reduced state (สถานะลดรูป) ของระบบซ้าย/ล่างได้ ในกรณีนี้ สถานะนั้นคือสถานะ mixed แบบสมบูรณ์ IX/n\mathbb{I}_{\mathsf{X}}/n ดังนั้น

Tr⁡Y(J(Φ)n)=IXn.\operatorname{Tr}_{\mathsf{Y}} \biggl(\frac{J(\Phi)}{n}\biggr) = \frac{\mathbb{I}_{\mathsf{X}}}{n}.

การตัด nn จากสองฝั่งให้ Tr⁡Y(J(Φ))=IX\operatorname{Tr}_{\mathsf{Y}} (J(\Phi)) = \mathbb{I}_{\mathsf{X}}

เราสามารถสรุปข้อสรุปเดียวกันได้โดยใช้ข้อเท็จจริงที่ว่าช่องต้องรักษา trace เสมอ ดังนั้น

Tr⁡Y(J(Φ))=∑a,b∈ΣTr⁡(Φ(∣a⟩⟨b∣)) ∣a⟩⟨b∣=∑a,b∈ΣTr⁡(∣a⟩⟨b∣) ∣a⟩⟨b∣=∑a∈Σ∣a⟩⟨a∣=IX.\begin{aligned} \operatorname{Tr}_{\mathsf{Y}} (J(\Phi)) & = \sum_{a,b\in\Sigma} \operatorname{Tr}\bigl(\Phi( \vert a\rangle\langle b \vert)\bigr) \, \vert a\rangle\langle b \vert \\ & = \sum_{a,b\in\Sigma} \operatorname{Tr}\bigl(\vert a\rangle\langle b \vert\bigr) \, \vert a\rangle\langle b \vert \\ & = \sum_{a\in\Sigma} \vert a\rangle\langle a \vert \\ & = \mathbb{I}_{\mathsf{X}}. \end{aligned}

โดยสรุป การแทนค่าแบบ Choi J(Φ)J(\Phi) สำหรับช่อง Φ\Phi ใด ๆ ต้องเป็น positive semidefinite และต้องตรงตาม

Tr⁡Y(J(Φ))=IX.\operatorname{Tr}_{\mathsf{Y}} (J(\Phi)) = \mathbb{I}_{\mathsf{X}}.

ดังที่เราจะเห็นเมื่อสิ้นสุดบทนี้ เงื่อนไขทั้งสองนี้ไม่เพียงแค่จำเป็นแต่ยังเพียงพอด้วย หมายความว่าการส่งเชิงเส้น Φ\Phi ใด ๆ จากเมทริกซ์ไปยังเมทริกซ์ที่ตรงตามข้อกำหนดเหล่านี้จะต้องเป็นช่องจริง ๆ

ช่อง dephasing แบบสมบูรณ์​

การแทนค่าแบบ Choi ของช่อง dephasing แบบสมบูรณ์ Δ\Delta คือ

J(Δ)=∑a,b=01∣a⟩⟨b∣⊗Δ(∣a⟩⟨b∣)=∑a=01∣a⟩⟨a∣⊗∣a⟩⟨a∣=(1000000000000001).\begin{aligned} J(\Delta) & = \sum_{a,b = 0}^{1} \vert a\rangle\langle b \vert \otimes \Delta\bigl(\vert a\rangle\langle b \vert\bigr) \\ & = \sum_{a = 0}^{1} \vert a\rangle\langle a \vert \otimes \vert a\rangle\langle a \vert \\[4mm] & = \begin{pmatrix} 1 & 0 & 0 & 0\\[1mm] 0 & 0 & 0 & 0\\[1mm] 0 & 0 & 0 & 0\\[1mm] 0 & 0 & 0 & 1 \end{pmatrix}. \end{aligned}

ช่อง depolarizing แบบสมบูรณ์​

การแทนค่าแบบ Choi ของช่อง depolarizing แบบสมบูรณ์คือ

J(Ω)=∑a,b=01∣a⟩⟨b∣⊗Ω(∣a⟩⟨b∣)=∑a=01∣a⟩⟨a∣⊗12I=12I⊗I=(12000012000012000012).\begin{aligned} J(\Omega) & = \sum_{a,b = 0}^{1} \vert a\rangle\langle b \vert \otimes \Omega\bigl(\vert a\rangle\langle b \vert\bigr)\\ & = \sum_{a = 0}^{1} \vert a\rangle\langle a \vert \otimes \frac{1}{2} \mathbb{I} \\[4mm] & = \frac{1}{2} \mathbb{I} \otimes \mathbb{I}\\[3mm] & = \begin{pmatrix} \frac{1}{2} & 0 & 0 & 0\\[1mm] 0 & \frac{1}{2} & 0 & 0\\[1mm] 0 & 0 & \frac{1}{2} & 0\\[1mm] 0 & 0 & 0 & \frac{1}{2} \end{pmatrix}. \end{aligned}

ช่อง reset ของ qubit​

การแทนค่าแบบ Choi ของช่อง reset ของ qubit Φ\Phi คือ

J(Λ)=∑a,b=01∣a⟩⟨b∣⊗Λ(∣a⟩⟨b∣)=∑a=01∣a⟩⟨a∣⊗∣0⟩⟨0∣=I⊗∣0⟩⟨0∣=(1000000000100000).\begin{aligned} J(\Lambda) & = \sum_{a,b = 0}^{1} \vert a\rangle\langle b \vert \otimes \Lambda\bigl(\vert a\rangle\langle b \vert\bigr)\\ & = \sum_{a = 0}^{1} \vert a\rangle\langle a \vert \otimes \vert 0\rangle \langle 0\vert\\[4mm] & = \mathbb{I} \otimes \vert 0\rangle \langle 0\vert\\[3mm] & = \begin{pmatrix} 1 & 0 & 0 & 0\\[1mm] 0 & 0 & 0 & 0\\[1mm] 0 & 0 & 1 & 0\\[1mm] 0 & 0 & 0 & 0 \end{pmatrix}. \end{aligned}

ช่อง identity​

การแทนค่าแบบ Choi ของช่อง identity ของ qubit Id⁡\operatorname{Id} คือ

J(Id⁡)=∑a,b=01∣a⟩⟨b∣⊗Id⁡(∣a⟩⟨b∣)=∑a,b=01∣a⟩⟨b∣⊗∣a⟩⟨b∣=(1001000000001001).\begin{aligned} J(\operatorname{Id}) & = \sum_{a,b = 0}^{1} \vert a\rangle\langle b \vert \otimes \operatorname{Id}\bigl(\vert a\rangle\langle b \vert\bigr) \\ & = \sum_{a,b = 0}^{1} \vert a \rangle\langle b \vert \otimes \vert a\rangle \langle b \vert \\ & = \begin{pmatrix} 1 & 0 & 0 & 1\\[1mm] 0 & 0 & 0 & 0\\[1mm] 0 & 0 & 0 & 0\\[1mm] 1 & 0 & 0 & 1 \end{pmatrix}. \end{aligned}

สังเกตโดยเฉพาะว่า J(Id⁡)J(\operatorname{Id}) ไม่ใช่เมทริกซ์ identity การแทนค่าแบบ Choi ไม่ได้อธิบายการกระทำของช่องโดยตรงในแบบที่เมทริกซ์แทนการส่งเชิงเส้นตามปกติ