2025-11-12T10:37:09.722784

Classical simulation of circuits with realistic odd-dimensional Gottesman-Kitaev-Preskill states

Calcluth, Hahn, Bermejo-Vega et al.
Classically simulating circuits with bosonic codes is challenging due to the prohibitive cost of simulating quantum systems with many, possibly infinite, energy levels. We propose an algorithm to simulate circuits with encoded Gottesman-Kitaev-Preskill (GKP) states, specifically for odd-dimensional encoded qudits. Our approach is tailored to be especially effective in the most challenging but practically relevant regime, where the codeword states exhibit high (but finite) squeezing. Our algorithm leverages the Zak-Gross Wigner function introduced by J. Davis et al. [arXiv:2407.18394], which represents infinitely squeezed encoded stabilizer states positively. The runtime of the algorithm scales with the negativity of the Wigner function, allowing for efficient simulation of certain large-scale circuits - namely, input stabilizer GKP states undergoing generalized GKP-encoded Clifford operations followed by modular measurements - with a high degree of squeezing. For stabilizer GKP states exhibiting 12 dB of squeezing, our algorithm can simulate circuits with up to 1,000 modes with less than double the number of samples required for a single input mode, in stark contrast to existing simulators. Therefore, this approach holds significant potential for benchmarking early implementations of quantum computing architectures utilizing bosonic codes.
academic

यथार्थवादी विषम-आयामी गॉटेसमैन-किटाएव-प्रेस्किल अवस्थाओं के साथ परिपथों का शास्त्रीय अनुकरण

मूल जानकारी

  • पेपर ID: 2412.13136
  • शीर्षक: यथार्थवादी विषम-आयामी गॉटेसमैन-किटाएव-प्रेस्किल अवस्थाओं के साथ परिपथों का शास्त्रीय अनुकरण
  • लेखक: कैमरून कैलक्लुथ, ओलिवर हान, जुआनी बर्मेजो-वेगा, एलेसेंड्रो फेरारो, गिउलिया फेरिनी
  • वर्गीकरण: quant-ph (क्वांटम भौतिकी)
  • प्रकाशन समय: 17 दिसंबर 2024 (arXiv प्रीप्रिंट)
  • पेपर लिंक: https://arxiv.org/abs/2412.13136

सारांश

यह पेपर गॉटेसमैन-किटाएव-प्रेस्किल (GKP) एन्कोडिंग पर आधारित विषम-आयामी क्वांटम परिपथों के अनुकरण के लिए एक शास्त्रीय अनुकरण एल्गोरिथ्म प्रस्तावित करता है। यह एल्गोरिथ्म विशेष रूप से उच्च संपीड़न के साथ परंतु सीमित संपीड़न वाले कोडवर्ड अवस्थाओं के सबसे चुनौतीपूर्ण किंतु व्यावहारिक रूप से प्रासंगिक परिदृश्य के लिए अनुकूलित है। एल्गोरिथ्म Zak-Gross Wigner फलन का उपयोग करता है, जो अनंत संपीड़न वाली एन्कोडित स्थिरांक अवस्थाओं का सकारात्मक प्रतिनिधित्व कर सकता है। एल्गोरिथ्म का चलन समय Wigner फलन की नकारात्मकता के समानुपाती है, जिससे कुछ बड़े परिपथों का कुशल अनुकरण संभव हो जाता है। 12dB संपीड़न वाली स्थिरांक GKP अवस्थाओं के लिए, यह एल्गोरिथ्म 1000 मोड तक के परिपथों का अनुकरण कर सकता है, जिसमें एकल-इनपुट मोड के लिए आवश्यक नमूनों से दो गुना से कम नमूनों की आवश्यकता होती है, जो मौजूदा अनुकारकों के साथ तीव्र विपरीतता दर्शाता है।

अनुसंधान पृष्ठभूमि और प्रेरणा

समस्या परिभाषा

बोसोनिक एन्कोडित क्वांटम परिपथों का शास्त्रीय अनुकरण विशाल चुनौतियों का सामना करता है, मुख्य कारण यह है कि कई (संभवतः अनंत) ऊर्जा स्तरों वाली क्वांटम प्रणालियों का अनुकरण करने की आवश्यकता है, जिसकी कम्प्यूटेशनल लागत अत्यंत अधिक है। मौजूदा अनुकरण एल्गोरिथ्म व्यावहारिक रूप से प्रासंगिक सीमित संपीड़न GKP अवस्थाओं को संभालने में महत्वपूर्ण सीमाएं रखते हैं।

महत्व विश्लेषण

  1. व्यावहारिक मूल्य: GKP एन्कोडिंग माइक्रोवेव गुहा युग्मित सुपरकंडक्टिंग परिपथ, फंसे हुए आयन और फोटोनिक प्लेटफॉर्म जैसी प्रायोगिक प्रणालियों में महसूस की गई है, और यह सहनशील क्वांटम कंप्यूटिंग के लिए एक महत्वपूर्ण उम्मीदवार है
  2. सैद्धांतिक महत्व: प्रारंभिक बोसोनिक क्वांटम प्रोसेसर के कार्यान्वयन को सत्यापित करने के लिए बेंचमार्क परीक्षण उपकरण प्रदान करता है
  3. तकनीकी चुनौती: पारंपरिक विधियां उच्च संपीड़न GKP अवस्थाओं के बड़े परिपथ अनुकरण को संभाल नहीं सकती हैं

मौजूदा विधियों की सीमाएं

  1. Fock आधार विस्तार विधि: चूंकि GKP कोडवर्ड में अनंत कई ऊर्जा स्तर शामिल हैं, बल प्रयोग अनुकरण विधि कुछ GKP क्वांटम बिट्स के लिए ही अव्यावहारिक हो जाती है
  2. अर्ध-संभाव्यता वितरण विधि: एन्कोडित GKP अवस्थाओं का Wigner फलन बहुत नकारात्मक है, जिससे चलन समय में विस्फोट होता है
  3. गॉसियन अवस्था अपघटन विधि: व्यावहारिक रूप से प्रासंगिक उच्च संपीड़न GKP अवस्थाओं पर लागू नहीं है, क्योंकि इन अवस्थाओं में बड़ी गैर-गॉसियन विशेषताएं हैं

मुख्य योगदान

  1. बहु-मोड Zak-Gross Wigner (ZGW) फलन प्रस्तावित किया: एकल-मोड ZGW फलन को n-मोड स्थिति तक सामान्यीकृत किया और सिद्ध किया कि यह संशोधित Stratonovich-Weyl स्वयंसिद्धों को संतुष्ट करता है
  2. ZGW फलन और Gross Wigner फलन के बीच संबंध स्थापित किया: CV GKP अवस्थाओं के ZGW फलन और एन्कोडित तार्किक अवस्थाओं के Gross Wigner फलन के बीच पत्राचार सिद्ध किया
  3. कुशल अनुकरण एल्गोरिथ्म विकसित किया:
    • आदर्श स्थिरांक GKP अवस्थाओं के लिए दुर्बल अनुकरण एल्गोरिथ्म
    • सीमित संपीड़न GKP अवस्थाओं के लिए संभाव्यता अनुमान एल्गोरिथ्म
  4. उल्लेखनीय प्रदर्शन सुधार प्राप्त किया: 12dB संपीड़न वाली स्थिरांक GKP अवस्थाओं के लिए, 1000 मोड के परिपथ का अनुकरण एकल-मोड के लिए आवश्यक नमूनों से दो गुना से कम नमूनों की आवश्यकता है

विधि विवरण

कार्य परिभाषा

निम्नलिखित क्वांटम परिपथ के माप परिणाम संभाव्यता वितरण का अनुकरण करें:

  • इनपुट: मनमाना GKP एन्कोडित क्वांटम बिट, विपरीत संपीड़न पैरामीटर Δ के साथ
  • विकास: पूर्णांक सिम्पलेक्टिक मैट्रिक्स द्वारा वर्णित एकात्मक संचालन और मनमाना विस्थापन (सभी एन्कोडित Clifford संचालन सहित)
  • माप: स्थिति आधार में मोड माप MZ(s)=1dneisnT(0,n)M_Z(s) = \frac{1}{d\ell}\sum_n e^{-i\ell s \cdot n}T_{(0,n)}

ZGW फलन परिभाषा

परिभाषा 1: बहु-मोड विषम-आयामी ZGW फलन को इस प्रकार परिभाषित किया जाता है: Wρ^(η)=Tr(ρ^A^η)W_{\hat{\rho}}(\eta) = \text{Tr}(\hat{\rho}\hat{A}_\eta)

जहां चरण-बिंदु संचालक है: A^η=1(2π)naZ2nei[a,η]+iπaXTaZT^a\hat{A}_\eta = \frac{1}{(2\pi)^n}\sum_{a \in \mathbb{Z}^{2n}} e^{i\ell[a,\eta]+i\pi a_X^T a_Z}\hat{T}_a

यहां η[0,d)2n\eta \in [0,d\ell)^{2n}, dd एक सकारात्मक विषम संख्या है, T^a=eiπaXTaZ/dT^aXT^aZ\hat{T}_a = e^{i\pi a_X^T a_Z/d}\hat{T}_{a_X}\hat{T}_{a_Z}

मुख्य सैद्धांतिक परिणाम

लेम्मा 1: ZGW फलन और Gross Wigner फलन के बीच संबंध Wρ^(η)=Wˉρˉ(u)(t)W_{\hat{\rho}}(\eta) = \bar{W}_{\bar{\rho}(u)}(t) जहां η=(u+t)\eta = \ell(u+t), u=1ηmod1u = \frac{1}{\ell}\eta \bmod 1, t=1ηut = \frac{1}{\ell}\eta - u

प्रमेय 1: सिम्पलेक्टिक विकास के तहत ZGW फलन परिवर्तन पूर्णांक सिम्पलेक्टिक मैट्रिक्स SS द्वारा वर्णित गॉसियन एकात्मक विकास के लिए, विकसित अवस्था का ZGW फलन है: WU^Sρ^U^S(η)=Wρ^(Sηt)W_{\hat{U}_S\hat{\rho}\hat{U}_S^\dagger}(\eta) = W_{\hat{\rho}}(S\eta - t)

प्रमेय 2: माप संभाव्यता की गणना तार्किक संचालक Z^Ln\hat{Z}_L^{\otimes n} का माप संभाव्यता वितरण है: Tr(ρ^M^Z(s))=dηXWρ^((ηXs))\text{Tr}(\hat{\rho}\hat{M}_Z(s)) = \int d\eta_X W_{\hat{\rho}}\left(\begin{pmatrix}\eta_X\\s\end{pmatrix}\right)

अनुकरण एल्गोरिथ्म

आदर्श स्थिरांक GKP अवस्थाओं का दुर्बल अनुकरण

  1. प्रभावी संभाव्यता वितरण Wρ^0(η)W_{\hat{\rho}_0}(\eta) से वेक्टर η\eta का नमूना लें
  2. सिम्पलेक्टिक मैट्रिक्स SS और विस्थापन वेक्टर tt की कार्रवाई के तहत रूपांतरित करें: SηtS\eta - t
  3. रैखिक विस्थापन संचालक लागू करें: Sηt+cS\eta - t + c
  4. माप परिणाम वेक्टर के पश्च भाग द्वारा दिया जाता है

सीमित संपीड़न GKP अवस्थाओं के लिए संभाव्यता अनुमान

अर्ध-संभाव्यता वितरण नमूना विधि का उपयोग करें: Pr(η)=1MU^S,cρ^0U^S,cWU^S,cρ^0U^S,c(η)\text{Pr}(\eta) = \frac{1}{M_{\hat{U}_{S,c}\hat{\rho}_0\hat{U}_{S,c}^\dagger}}|W_{\hat{U}_{S,c}\hat{\rho}_0\hat{U}_{S,c}^\dagger}(\eta)|

आवश्यक नमूनों की संख्या है: N=2ϵ2Mρ^02log(2/δ)N = \frac{2}{\epsilon^2}M_{\hat{\rho}_0}^2\log(2/\delta)

जहां Mρ^M_{\hat{\rho}} ZGW फलन की नकारात्मकता (नकारात्मक आयतन) है।

प्रायोगिक सेटअप

सैद्धांतिक विश्लेषण ढांचा

यह पेपर मुख्य रूप से सैद्धांतिक विश्लेषण और संख्यात्मक गणना पर केंद्रित है, विशेष ध्यान:

  1. ZGW फलन नकारात्मकता की गणना: विभिन्न संपीड़न स्तरों पर GKP अवस्थाओं के लिए
  2. नमूना जटिलता विश्लेषण: नकारात्मकता के आधार पर स्केलिंग व्यवहार
  3. मौजूदा विधियों के साथ तुलना: सैद्धांतिक जटिलता तुलना

मूल्यांकन मेट्रिक्स

  • ZGW लॉगरिदमिक नकारात्मकता: logMρ^\log M_{\hat{\rho}}, अनुकरण जटिलता को मापता है
  • नमूना संख्या स्केलिंग: एकल-मोड स्थिति के सापेक्ष नमूना संख्या वृद्धि
  • मोड क्षमता: प्रभावी रूप से अनुकरण किए जा सकने वाली अधिकतम मोड संख्या

प्रायोगिक परिणाम

सीमित संपीड़न GKP अवस्थाओं का ZGW फलन

स्पष्ट गणना के माध्यम से, वास्तविक 0-तार्किक GKP अवस्था का ZGW फलन प्राप्त किया जाता है: WGKP((uv))ϑ(Γ;z)W_{\text{GKP}}\left(\begin{pmatrix}u\\v\end{pmatrix}\right) \propto \vartheta(\Gamma; z)

जहां z=(v/(d),u/(d),0,0)Tz = (v/(d\ell), -u/(d\ell), 0, 0)^T, Γ\Gamma संपीड़न पैरामीटर Δ\Delta द्वारा निर्धारित एक जटिल मैट्रिक्स है।

मुख्य संख्यात्मक परिणाम

Δ=0.25\Delta = 0.25 के लिए (12dB संपीड़न के अनुरूप):

  • एकल-मोड नकारात्मकता: Mρ^0Δe3×104M_{\hat{\rho}_0^\Delta} \approx e^{3 \times 10^{-4}}
  • हजार-मोड अनुकरण: नकारात्मकता की गुणात्मक प्रकृति के कारण, 1000 इनपुट मोड का अनुकरण करने के लिए एकल-मोड के दो गुना से कम नमूनों की आवश्यकता है
  • प्रदर्शन सुधार: मौजूदा विधियों की तुलना में कई परिमाण के क्रम में सुधार

जादू अवस्था विश्लेषण

π अवस्था ψπ=13(0L+1L2L)|\psi_\pi\rangle = \frac{1}{\sqrt{3}}(|0_L\rangle + |1_L\rangle - |2_L\rangle) के लिए:

  • सभी Δ\Delta मानों पर उल्लेखनीय नकारात्मकता प्रदर्शित करता है
  • बड़े परिपथ अनुकरण के समय महत्वपूर्ण ओवरहेड, जैसा कि अपेक्षित है

संबंधित कार्य

शास्त्रीय अनुकरण एल्गोरिथ्म विकास

  1. गॉटेसमैन-नील प्रमेय: स्थिरांक परिपथ का बहुपद समय अनुकरण
  2. अर्ध-संभाव्यता विधि: Gross Wigner फलन पर आधारित नमूना एल्गोरिथ्म
  3. बोसोनिक प्रणाली अनुकरण: अर्ध-संभाव्यता वितरण सकारात्मकता, टेंसर नेटवर्क, तारकीय प्रतिनिधित्व आदि का उपयोग

GKP एन्कोडिंग अनुकरण की वर्तमान स्थिति

मौजूदा विधियां मुख्य रूप से सीमित हैं:

  • आदर्श (अनंत संपीड़न) GKP अवस्थाएं
  • गॉसियन संचालन और शून्य-विषमता माप
  • व्यावहारिक रूप से प्रासंगिक सीमित संपीड़न स्थिति को संभाल नहीं सकते

इस पेपर की नवीनता

मौजूदा कार्य की तुलना में, यह पेपर पहली बार सीमित संपीड़न GKP अवस्थाओं का कुशल अनुकरण प्राप्त करता है, विशेष रूप से उच्च संपीड़न व्यवस्था में उल्लेखनीय प्रदर्शन सुधार के साथ।

निष्कर्ष और चर्चा

मुख्य निष्कर्ष

  1. सैद्धांतिक योगदान: ZGW फलन का एक संपूर्ण सैद्धांतिक ढांचा स्थापित किया, जिसमें बहु-मोड सामान्यीकरण और विकास नियम शामिल हैं
  2. एल्गोरिथ्म नवीनता: वास्तविक GKP अवस्थाओं के लिए उपयुक्त कुशल अनुकरण एल्गोरिथ्म प्रस्तावित किया
  3. व्यावहारिक मूल्य: बोसोनिक क्वांटम प्रोसेसर को सत्यापित करने के लिए एक शक्तिशाली बेंचमार्क परीक्षण उपकरण प्रदान करता है

सीमाएं

  1. विषम-आयामी प्रतिबंध: वर्तमान विधि केवल विषम-आयामी GKP क्वांटम बिट्स पर लागू होती है, सम-आयामी (जैसे क्वांटम बिट्स) तक विस्तार अभी भी एक खुली समस्या है
  2. संचालन प्रतिबंध: पूर्णांक सिम्पलेक्टिक मैट्रिक्स द्वारा वर्णित संचालन तक सीमित, मनमाना एकल-मोड संपीड़न को संभाल नहीं सकते
  3. शोर मॉडल: नुकसान जैसे व्यावहारिक शोर के ZGW फलन नकारात्मकता पर प्रभाव पर विचार नहीं किया गया है

भविष्य की दिशाएं

  1. GKP क्वांटम बिट्स तक विस्तार: द्विआयामी स्थिति के लिए उपयुक्त विधियां विकसित करें
  2. अधिक व्यापक बोसोनिक एन्कोडिंग: ढांचे को GKP से परे अन्य बोसोनिक एन्कोडिंग तक सामान्यीकृत करें
  3. अर्ध-संभाव्यता ढांचा: बोसोनिक एन्कोडिंग के लिए विशेष व्यापक अर्ध-संभाव्यता सिद्धांत विकसित करें

गहन मूल्यांकन

शक्तियां

  1. सैद्धांतिक कठोरता: गणितीय व्युत्पत्ति संपूर्ण, प्रमाण विस्तृत, सैद्धांतिक आधार मजबूत
  2. व्यावहारिक शक्ति: सीमित संपीड़न GKP अवस्थाओं के अनुकरण की व्यावहारिक समस्या को हल करता है
  3. उत्कृष्ट प्रदर्शन: उच्च संपीड़न व्यवस्था में परिमाण के क्रम में प्रदर्शन सुधार प्राप्त करता है
  4. विधि नवीनता: ZGW फलन का बहु-मोड सामान्यीकरण और अनुप्रयोग महत्वपूर्ण सैद्धांतिक योगदान है

कमियां

  1. सीमित लागू क्षेत्र: केवल विषम-आयामी प्रणालियों पर लागू, सार्वभौमिकता को सीमित करता है
  2. प्रायोगिक सत्यापन की कमी: मुख्य रूप से सैद्धांतिक कार्य, वास्तविक क्वांटम प्रणालियों के साथ तुलना की कमी
  3. जटिलता विश्लेषण: बड़े परिपथों के लिए वास्तविक कम्प्यूटेशनल जटिलता विश्लेषण पर्याप्त विस्तृत नहीं है

प्रभाव

  1. शैक्षणिक मूल्य: क्वांटम अनुकरण सिद्धांत के लिए नए उपकरण और दृष्टिकोण प्रदान करता है
  2. व्यावहारिक संभावनाएं: प्रारंभिक सहनशील क्वांटम कंप्यूटिंग प्रयोगों के लिए महत्वपूर्ण मार्गदर्शन
  3. तकनीकी प्रेरणा: अधिक बोसोनिक एन्कोडिंग अनुकरण विधियों के विकास को प्रेरित कर सकता है

लागू परिदृश्य

  • बोसोनिक क्वांटम प्रोसेसर का बेंचमार्क परीक्षण और सत्यापन
  • GKP एन्कोडिंग सहनशील क्वांटम कंप्यूटिंग योजना का सैद्धांतिक विश्लेषण
  • निरंतर चर क्वांटम कंप्यूटिंग का शास्त्रीय अनुकरण अनुसंधान

संदर्भ

यह पेपर 68 संबंधित संदर्भों का हवाला देता है, मुख्य रूप से शामिल हैं:

  • गॉटेसमैन-किटाएव-प्रेस्किल एन्कोडिंग का मूल कार्य
  • क्वांटम अनुकरण एल्गोरिथ्म की शास्त्रीय साहित्य
  • बोसोनिक क्वांटम कंप्यूटिंग की नवीनतम प्रगति
  • Wigner फलन और अर्ध-संभाव्यता वितरण का सैद्धांतिक आधार

समग्र मूल्यांकन: यह GKP एन्कोडित क्वांटम परिपथों के शास्त्रीय अनुकरण में एक उच्च गुणवत्ता वाला सैद्धांतिक भौतिकी पेपर है जो महत्वपूर्ण सफलता प्राप्त करता है। हालांकि कुछ सीमाएं हैं, इसके सैद्धांतिक योगदान और व्यावहारिक मूल्य दोनों बहुत महत्वपूर्ण हैं, और क्वांटम कंप्यूटिंग क्षेत्र के लिए महत्वपूर्ण अर्थ रखते हैं।