خوارزمية تقدير الطور الكمومي

في الحوسبة الكمومية ، تُعدّ خوارزمية تقدير الطور الكمومي خوارزمية كمومية لتقدير الطور المقابل لقيمة ذاتية لمؤثر وحدوي مُعطى . ولأن القيم الذاتية للمؤثر الوحدوي لها دائمًا قيمة مطلقة تساوي واحدًا ، فإنها تُوصَف بطورها، وبالتالي يمكن وصف الخوارزمية بشكل مكافئ بأنها تسترجع إما الطور أو القيمة الذاتية نفسها. وقد طُرحت هذه الخوارزمية لأول مرة من قِبَل أليكسي كيتايف عام ١٩٩٥. [ ١ ] [ ٢ ] : ٢٤٦

يُستخدم تقدير الطور بشكل متكرر كإجراء فرعي في خوارزميات الكم الأخرى، مثل خوارزمية شور ، [ 2 ] : 131 خوارزمية الكم لأنظمة المعادلات الخطية ، وخوارزمية العد الكمي .

نظرة عامة على الخوارزمية

تعمل الخوارزمية على مجموعتين من الكيوبتات، يُشار إليهما في هذا السياق باسم المسجلات . تحتوي المسجلتان علىن{\displaystyle n}وم{\displaystyle m}الكيوبتات، على التوالي. لنفترضيو{\displaystyle U}يكون عاملًا وحدويًا يعمل علىم{\displaystyle m}- سجل الكيوبت . القيم الذاتية للمؤثر الوحدوي لها معيار يساوي واحدًا، وبالتالي تتميز بطورها. لذا إذا|ψ{\displaystyle |\psi \rangle }هو متجه ذاتي لـيو{\displaystyle U}، ثميو|ψ=هـ2πأناθ|ψ{\displaystyle U|\psi \rangle =e^{2\pi i\theta }\left|\psi \right\rangle }بالنسبة للبعضθR{\displaystyle \theta \in \mathbb {R} }بسبب دورية الدالة الأسية المركبة، يمكننا دائمًا افتراض0θ<1{\displaystyle 0\leq \theta <1}.

الهدف هو إنتاج تقريب جيد لـθ{\displaystyle \theta }باستخدام عدد قليل من البوابات واحتمالية نجاح عالية. تحقق خوارزمية تقدير الطور الكمومي هذا بافتراض الوصول عن بُعد إلىيو{\displaystyle U}وامتلاك|ψ{\displaystyle |\psi \rangle }متاح كحالة كمومية . هذا يعني أنه عند مناقشة كفاءة الخوارزمية، فإننا نهتم فقط بعدد المراتيو{\displaystyle U}يجب استخدامها، ولكن ليس فيما يتعلق بتكلفة التنفيذيو{\displaystyle U}نفسها.

وبشكل أدق، تُعيد الخوارزمية باحتمالية عالية تقريبًا لـθ{\displaystyle \theta }، ضمن الخطأ التراكميε{\displaystyle \varepsilon }، استخدامن=يا(سجل(1/ε)){\displaystyle n=O(\log(1/\varepsilon ))}الكيوبتات في السجل الأول، ويا(1/ε){\displaystyle O(1/\varepsilon )}عمليات U المتحكم بها . علاوة على ذلك، يمكننا تحسين احتمالية النجاح لـ1-Δ{\displaystyle 1-\Delta }لأيΔ>0{\displaystyle \Delta >0}باستخدام إجمالييا(سجل(1/Δ)/ε){\displaystyle O(\log(1/\Delta )/\varepsilon )}استخدامات U المتحكم بها، وهذا هو الأمثل. [ 3 ]

وصف تفصيلي للخوارزمية

الدائرة الخاصة بتقدير الطور الكمي.

استعداد الدولة

الحالة الابتدائية للنظام هي:

|Ψ0=|0ن|ψ،{\displaystyle |\Psi _{0}\rangle =|0\rangle ^{\otimes n}|\psi \rangle ,}

أين|ψ{\displaystyle |\psi \rangle }هوم{\displaystyle m}حالة الكيوبت التي تتطور من خلاليو{\displaystyle U}نطبق أولاً عملية بوابة هادامارد ذات n كيوبتحن{\displaystyle H^{\otimes n}}في السجل الأول، الذي ينتج الحالة:|Ψ1=(حنأنام)|Ψ0=12ن2(|0+|1)ن|ψ=12ن/2ج=02ن-1|ج|ψ.{\displaystyle |\Psi _{1}\rangle =(H^{\otimes n}\otimes I_{m})|\Psi _{0}\rangle ={\frac {1}{2^{\frac {n}{2}}}}(|0\rangle +|1\rangle )^{\otimes n}|\psi \rangle ={\frac {1}{2^{n/2}}}\sum _{j=0}^{2^{n}-1}|j\rangle |\psi \rangle .}لاحظ أننا هنا ننتقل بين النظام الثنائي ون{\displaystyle n}التمثيل -ary لـن{\displaystyle n}- سجل الكيوبت: الكيت|ج{\displaystyle |j\rangle }على الجانب الأيمن يوجد اختصار لـن{\displaystyle n}حالة الكيوبت|ج=0ن-1|ج{\displaystyle |j\rangle \equiv \bigotimes _{\ell =0}^{n-1}|j_{\ell }\rangle }، أينج==0ن-1ج2{\displaystyle j=\sum _{\ell =0}^{n-1}j_{\ell }2^{\ell }}هو التحلل الثنائي لـج{\displaystyle j}.

عمليات التحكم U

هذه الولاية|Ψ1{\displaystyle |\Psi _{1}\rangle }ثم يتطور من خلال التطور الوحدوي المتحكم فيهيوج{\displaystyle U_{C}}الذي يمكن كتابة فعله على النحو التالييوج(|ك|ψ)=|ك(يوك|ψ)،{\displaystyle U_{C}(|k\rangle \otimes |\psi \rangle )=|k\rangle \otimes (U^{k}|\psi \rangle ),}للجميعك=0،...،2ن-1{\displaystyle k=0,...,2^{n}-1}ويمكن كتابة هذا التطور بإيجاز على النحو التالي:يوج=ك=02ن-1|كك|يوك،{\displaystyle U_{C}=\sum _{k=0}^{2^{n}-1}|k\rangle \!\langle k|\otimes U^{k},}مما يسلط الضوء على طبيعته الخاضعة للرقابة: فهو ينطبقيوك{\displaystyle U^{k}}إلى السجل الثاني بشرط أن يكون السجل الأول|ك{\displaystyle |k\rangle }مع الأخذ في الاعتبار شرط القيمة الذاتية الذي ينطبق على|ψ{\displaystyle |\psi \rangle }تطبيقيوج{\displaystyle U_{C}}ل|Ψ1{\displaystyle |\Psi _{1}\rangle }وهذا يعطي|Ψ2يوج|Ψ1=(12ن/2ك=02ن-1هـ2πأناθك|ك)|ψ،{\displaystyle |\Psi _{2}\rangle \equiv U_{C}|\Psi _{1}\rangle =\left({\frac {1}{2^{n/2}}}\sum _{k=0}^{2^{n}-1}e^{2\pi i\theta k}|k\rangle \right)\otimes |\psi \rangle ,}حيث استخدمنايوك|ψ=هـ2πأناكθ|ψ{\displaystyle U^{k}|\psi \rangle =e^{2\pi ik\theta }|\psi \rangle }.

لإثبات ذلكيوج{\displaystyle U_{C}}ويمكن أيضًا تنفيذه بكفاءة، لاحظ أنه يمكننا كتابةيوج==0ن-1ج(يو2){\displaystyle U_{C}=\prod _{\ell =0}^{n-1}C_{\ell }(U^{2^{\ell }})}، أينج(يو2){\displaystyle C_{\ell }(U^{2^{\ell }})}يشير إلى عملية التطبيقيو2{\displaystyle U^{2^{\ell }}}إلى السجل الثاني بشرط أن{\displaystyle \ell }الكيوبت رقم - من السجل الأول هو|1{\displaystyle |1\rangle }. رسميًا، يمكن وصف هذه البوابات من خلال عملها على النحو التالي:ج(يوك)(|ج|ψ)=|ج(يوجك|ψ).{\displaystyle C_{\ell }(U^{k})(|j\rangle \otimes |\psi \rangle )=|j\rangle \otimes (U^{j_{\ell }k}|\psi \rangle ).}يمكن تفسير هذه المعادلة على أنها تعني أن الحالة تبقى دون تغيير عندماج=0{\displaystyle j_{\ell }=0}أي عندما{\displaystyle \ell }الكيوبت رقم - هو|0{\displaystyle |0\rangle }بينما البوابةيوك{\displaystyle U^{k}}يتم تطبيق ذلك على السجل الثاني عندما{\displaystyle \ell }الكيوبت رقم - هو|1{\displaystyle |1\rangle }وبالتالي فإن تركيب هذه البوابات المتحكم بها يعطي=0ن-1ج(يو2)(|ج|ψ)=|ج(يو=0ن-1ج2|ψ)=يوج(|ج|ψ)،{\displaystyle \prod _{\ell =0}^{n-1}C_{\ell }(U^{2^{\ell }})(|j\rangle \otimes |\psi \rangle )=|j\rangle \otimes \left(U^{\sum _{\ell =0}^{n-1}j_{\ell }2^{\ell }}|\psi \rangle \right)=U_{C}\left(|j\rangle \otimes |\psi \rangle \right),}وتأتي الخطوة الأخيرة مباشرة بعد التفكيك الثنائيج==0ن-1ج2{\displaystyle j=\sum _{\ell =0}^{n-1}j_{\ell }2^{\ell }}.

من هذه النقطة فصاعدًا، يُترك السجل الثاني دون تغيير، وبالتالي يصبح من الملائم كتابة|Ψ2=|Ψ~2|ψ{\displaystyle |\Psi _{2}\rangle =|{\tilde {\Psi }}_{2}\rangle \otimes |\psi \rangle }، مع|Ψ~2{\displaystyle |{\tilde {\Psi }}_{2}\rangle }حالةن{\displaystyle n}سجل الكيوبت، وهو السجل الوحيد الذي نحتاج إلى مراعاته لبقية الخوارزمية.

تطبيق تحويل فورييه الكمي العكسي

يتضمن الجزء الأخير من الدائرة تطبيق تحويل فورييه الكمي العكسي (QFT).سؤالFتي{\displaystyle {\mathcal {QFT}}}في السجل الأول لـ|Ψ2{\displaystyle |\Psi _{2}\rangle }:|Ψ~3=سؤالFتي2ن-1|Ψ~2.{\displaystyle |{\tilde {\Psi }}_{3}\rangle ={\mathcal {QFT}}_{2^{n}}^{-1}|{\tilde {\Psi }}_{2}\rangle .}تتميز نظرية الحقل الكمومي وعكسها بتأثيرهما على حالات الأساس كما يلي:سؤالFتيشمال|ك=شمال-1/2ج=0شمال-1هـ2πأناشمالجك|ج،سؤالFتيشمال-1|ك=شمال-1/2ج=0شمال-1هـ-2πأناشمالجك|ج.{\displaystyle {\begin{aligned}{\mathcal {QFT}}_{N}|k\rangle &=N^{-1/2}\sum _{j=0}^{N-1}e^{{\frac {2\pi i}{N}}jk}|j\rangle ,\\{\mathcal {QFT}}_{N}^{-1}|k\rangle &=N^{-1/2}\sum _{j=0}^{N-1}e^{-{\frac {2\pi i}{N}}jk}|j\rangle .\end{aligned}}}ويترتب على ذلك أن

|Ψ~3=12ن2ك=02ن-1هـ2πأناθك(12ن2x=02ن-1هـ-2πأناكx2ن|x)=12نx=02ن-1ك=02ن-1هـ-2πأناك2ن(x-2نθ)|x.{\displaystyle |{\tilde {\Psi }}_{3}\rangle ={\frac {1}{2^{\frac {n}{2}}}}\sum _{k=0}^{2^{n}-1}e^{2\pi i\theta k}\left({\frac {1}{2^{\frac {n}{2}}}}\sum _{x=0}^{2^{n}-1}e^{\frac {-2\pi ikx}{2^{n}}}|x\rangle \right)={\frac {1}{2^{n}}}\sum _{x=0}^{2^{n}-1}\sum _{k=0}^{2^{n}-1}e^{-{\frac {2\pi ik}{2^{n}}}\left(x-2^{n}\theta \right)}|x\rangle .}

تحليل الحالة في الأساس الحسابي على النحو التالي:|Ψ~3=x=02ن-1جx|x،{\textstyle |{\tilde {\Psi }}_{3}\rangle =\sum _{x=0}^{2^{n}-1}c_{x}|x\rangle ,}وبالتالي فإن المعاملات تساويجx12نك=02ن-1هـ-2πأناك2ن(x-2نθ)=12نك=02ن-1هـ-2πأناك2ن(x-أ)هـ2πأنادلتاك،{\displaystyle c_{x}\equiv {\frac {1}{2^{n}}}\sum _{k=0}^{2^{n}-1}e^{-{\frac {2\pi ik}{2^{n}}}(x-2^{n}\theta )}={\frac {1}{2^{n}}}\sum _{k=0}^{2^{n}-1}e^{-{\frac {2\pi ik}{2^{n}}}\left(x-a\right)}e^{2\pi i\delta k},}حيث كتبنا2نθ=أ+2ندلتا،{\displaystyle 2^{n}\theta =a+2^{n}\delta ,}معأ{\displaystyle a}هو أقرب عدد صحيح إلى2نθ{\displaystyle 2^{n}\theta }الفرق2ندلتا{\displaystyle 2^{n}\delta } يجب أن يفي بالتعريف0|2ندلتا|12{\displaystyle 0\leqslant |2^{n}\delta |\leqslant {\tfrac {1}{2}}}وهذا يعادل تقريبًا قيمةθ[0،1]{\displaystyle \theta \in [0,1]}عن طريق التقريب2نθ{\displaystyle 2^{n}\theta }إلى أقرب عدد صحيح.

قياس

تتضمن الخطوة الأخيرة إجراء قياس في قاعدة البيانات الحسابية على السجل الأول. وهذا ينتج عنه النتيجة|y{\displaystyle |y\rangle }باحتمالبرو(y)=|جy|2=|12نك=02ن-1هـ-2πأناك2ن(y-أ)هـ2πأنادلتاك|2.{\displaystyle \Pr(y)=|c_{y}|^{2}=\left|{\frac {1}{2^{n}}}\sum _{k=0}^{2^{n}-1}e^{{\frac {-2\pi ik}{2^{n}}}(y-a)}e^{2\pi i\delta k}\right|^{2}.}ويترتب على ذلك أنبرو(أ)=1{\displaystyle \operatorname {Pr} (a)=1}لودلتا=0{\displaystyle \delta =0}أي عندماθ{\displaystyle \theta }يمكن كتابتها على النحو التاليθ=أ/2ن{\displaystyle \theta =a/2^{n}}يجد المرء دائماً النتيجةy=أ{\displaystyle y=a}من ناحية أخرى، إذادلتا0{\displaystyle \delta \neq 0}، والاحتمالية هيبرو(أ)=122ن|ك=02ن-1هـ2πأنادلتاك|2=122ن|1-هـ2πأنا2ندلتا1-هـ2πأنادلتا|2.{\displaystyle \operatorname {Pr} (a)={\frac {1}{2^{2n}}}\left|\sum _{k=0}^{2^{n}-1}e^{2\pi i\delta k}\right|^{2}={\frac {1}{2^{2n}}}\left|{\frac {1-{e^{2\pi i2^{n}\delta }}}{1-{e^{2\pi i\delta }}}}\right|^{2}.}من هذا التعبير يمكننا أن نرى أنبرو(أ)4π20.405{\displaystyle \Pr(a)\geqslant {\frac {4}{\pi ^{2}}}\approx 0.405}متىدلتا0{\displaystyle \delta \neq 0}ولتوضيح ذلك، نلاحظ أنه من تعريفدلتا{\displaystyle \delta }لدينا عدم المساواة|دلتا|12ن+1{\displaystyle |\delta |\leqslant {\tfrac {1}{2^{n+1}}}}وبالتالي: [ 4 ] : ​​157 [ 5 ] : 348برو(أ)=122ن|1-هـ2πأنا2ندلتا1-هـ2πأنادلتا|2ل دلتا0=122ن|2الخطيئة(π2ندلتا)2الخطيئة(πدلتا)|2|1-هـ2أناx|2=4|الخطيئة(x)|2=122ن|الخطيئة(π2ندلتا)|2|الخطيئة(πدلتا)|2122ن|الخطيئة(π2ندلتا)|2|πدلتا|2|الخطيئة(πدلتا)||πدلتا|122ن|22ندلتا|2|πدلتا|2|22ندلتا||الخطيئة(π2ندلتا)| ل |دلتا|12ن+14π2.{\displaystyle {\begin{aligned}\Pr(a)&={\frac {1}{2^{2n}}}\left|{\frac {1-{e^{2\pi i2^{n}\delta }}}{1-{e^{2\pi i\delta }}}}\right|^{2}&&{\text{for }}\delta \neq 0\\&={\frac {1}{2^{2n}}}\left|{\frac {2\sin \left(\pi 2^{n}\delta \right)}{2\sin(\pi \delta )}}\right|^{2}&&\left|1-e^{2ix}\right|^{2}=4\left|\sin(x)\right|^{2}\\&={\frac {1}{2^{2n}}}{\frac {\left|\sin \left(\pi 2^{n}\delta \right)\right|^{2}}{|\sin(\pi \delta )|^{2}}}\\&\geqslant {\frac {1}{2^{2n}}}{\frac {\left|\sin \left(\pi 2^{n}\delta \right)\right|^{2}}{|\pi \delta |^{2}}}&&|\sin(\pi \delta )|\leqslant |\pi \delta |\\&\geqslant {\frac {1}{2^{2n}}}{\frac {|2\cdot 2^{n}\delta |^{2}}{|\pi \delta |^{2}}}&&|2\cdot 2^{n}\delta |\leqslant |\sin(\pi 2^{n}\delta )|{\text{ for }}|\delta |\leqslant {\frac {1}{2^{n+1}}}\\&\geqslant {\frac {4}{\pi ^{2}}}.\end{aligned}}}

نستنتج أن الخوارزمية توفر الأفضلن{\displaystyle n}تقدير بت واحد (أي، واحد ضمن1/2ن{\displaystyle 1/2^{n}}(من الإجابة الصحيحة) منθ{\displaystyle \theta }باحتمالية لا تقل عن4/π2{\displaystyle 4/\pi ^{2}}بإضافة عدد من الكيوبتات الإضافية من رتبةيا(سجل(1/ϵ)){\displaystyle O(\log(1/\epsilon ))}وبحذف الكيوبتات الإضافية، يمكن أن تزداد الاحتمالية إلى1-ϵ{\displaystyle 1-\epsilon }[ 5 ]

أمثلة على الألعاب

لنفترض أبسط مثال ممكن للخوارزمية، حيث فقطن=1{\displaystyle n=1}الكيوبت، بالإضافة إلى الكيوبتات المطلوبة للترميز|ψ{\displaystyle |\psi \rangle }، متضمنة. لنفترض أن القيمة الذاتية لـ|ψ{\displaystyle |\psi \rangle }يقرأλ=هـ2πأناθ{\displaystyle \lambda =e^{2\pi i\theta }}،θ[0،1){\displaystyle \theta \in [0,1)}يُولّد الجزء الأول من الخوارزمية حالة الكيوبت الواحد|ϕ12(|0+λ|1){\textstyle |\phi \rangle \equiv {\frac {1}{\sqrt {2}}}(|0\rangle +\lambda |1\rangle )}بتطبيق نظرية الحقل الكمومي العكسية، يُعادل ذلك في هذه الحالة تطبيق بوابة هادامارد . وبالتالي، تكون احتمالات النتيجة النهائية كما يلي:ص±=|±|ϕ|2{\displaystyle p_{\pm }=|\langle \pm |\phi \rangle |^{2}}أين|±12(|0±|1){\textstyle |\pm \rangle \equiv {\frac {1}{\sqrt {2}}}(|0\rangle \pm |1\rangle )}أو بشكل أكثر وضوحاً،ص±=|1±λ|24=1±كوس(2πθ)2.{\displaystyle p_{\pm }={\frac {|1\pm \lambda |^{2}}{4}}={\frac {1\pm \cos(2\pi \theta )}{2}}.}يفترضλ=1{\displaystyle \lambda =1}، معنى|ϕ=|+{\displaystyle |\phi \rangle =|+\rangle }. ثمص+=1{\displaystyle p_{+}=1}،ص-=0{\displaystyle p_{-}=0}ونستعيد بشكل حتمي القيمة الدقيقة لـλ{\displaystyle \lambda }من نتائج القياس. وينطبق الشيء نفسه إذاλ=-1{\displaystyle \lambda =-1}.

أما من ناحية أخرىλ=هـ2πأنا/3{\displaystyle \lambda =e^{2\pi i/3}}، ثمص±=[1±كوس(2π/3)]/2{\displaystyle p_{\pm }=[1\pm \cos(2\pi /3)]/2}، إنه،ص+=1/4{\displaystyle p_{+}=1/4}وص-=3/4{\displaystyle p_{-}=3/4}في هذه الحالة، لا تكون النتيجة حتمية، لكننا مع ذلك نجد النتيجة.|-{\displaystyle |-\rangle }على الأرجح، بما يتوافق مع حقيقة أن2/3{\displaystyle 2/3}أقرب إلى 1 منه إلى 0.

وبشكل أعم، إذاλ=هـ2πأناθ{\displaystyle \lambda =e^{2\pi i\theta }}، ثمص+1/2{\displaystyle p_{+}\geq 1/2}إذا وفقط إذا|θ|1/4{\displaystyle |\theta |\leq 1/4}وهذا يتوافق مع النتائج المذكورة أعلاه لأنه في الحالاتλ=±1{\displaystyle \lambda =\pm 1}، بما يتوافق معθ=0،1/2{\displaystyle \theta =0,1/2}يتم استرجاع الطور بشكل حتمي، ويتم استرجاع الأطوار الأخرى بدقة أعلى كلما كانت أقرب إلى هذين الطورين.

انظر أيضاً

مراجع

  1. كيتايف، أ. يو (20-11-1995). "القياسات الكمومية ومسألة المثبت الأبلي". arXiv : quant-ph/9511026 .
  2. 1 2 نيلسن، مايكل أ. وإسحاق ل. تشوانغ (2001). الحوسبة الكمومية والمعلومات الكمومية (طبعة مُعاد طباعتها ). كامبريدج [ua]: مطبعة جامعة كامبريدج. ISBN  978-0521635035.
  3. ماندي، نيخيل س.؛ رونالد دي وولف (2023). "حدود دقيقة لتقدير الطور الكمومي والمشاكل ذات الصلة". arXiv : 2305.04908 [ quant-ph ].
  4. ^ بينينتي، جويليانو؛ كاساتي، جوليو؛ ستريني، جوليانو (2004). مبادئ الحساب الكمي والمعلومات (أعيد طبعه. إد.). نيو جيرسي [وا]: العالم العلمي. رقم ISBN  978-9812388582.
  5. 1 2 كليف، ر.؛ إيكرت، أ.؛ ماكيافيلو، س.؛ موسكا، م. (8 يناير 1998). "إعادة النظر في الخوارزميات الكمومية". وقائع الجمعية الملكية أ: العلوم الرياضية والفيزيائية والهندسية . 454 (1969): 339-354 . arXiv : quant-ph/9708016 . Bibcode : 1998RSPSA.454..339C . doi : 10.1098/rspa.1998.0164 . S2CID 16128238 .