The study of real time dynamics of nuclear systems is of great importance to provide theoretical predictions of cross sections relevant for both terrestrial experiments as well as applications in astrophysics. First principles simulations of these dynamical processes is however hindered by an exponential cost in classical resources and the possibility of performing scalable simulations using quantum computers is currently an active field of research. In this work we provide the first complete characterization of the resource requirements for studying nuclear dynamics with the full Leading Order (LO) pionless EFT Hamiltonian in first quantization employing simulation strategies using both product formulas as well as Quantum Signal Processing. In particular, we show that time evolution of such an Hamiltonian can be performed with polynomial resources in the number of particles, and logarithmic resources in the number of single-particle basis states. This result provides an exponential improvement compared with previous work on the same Hamiltonian model in second quantization. We find that interesting simulations for low energy nuclear scattering could be achievable with tens of millions of T gates and few hundred logical qubits suggesting that the study of simple nuclear reactions could be amenable for early fault tolerant quantum platforms.
पेपर ID : 2507.22814शीर्षक : प्रथम परिमाणीकरण में परमाणु गतिशीलता का क्वांटम सिमुलेशनलेखक : लुका स्पैग्नोली, चियारा लिसोनी, एलेसांड्रो रोजेरो (ट्रेंटो विश्वविद्यालय और INFN-TIFPA)वर्गीकरण : quant-ph (क्वांटम भौतिकी), nucl-th (परमाणु सिद्धांत)प्रकाशन तिथि : 10 नवंबर 2025 (arXiv v2)पेपर लिंक : https://arxiv.org/abs/2507.22814 परमाणु प्रणालियों की वास्तविक समय गतिशीलता का अध्ययन जमीनी प्रयोगों और खगोल भौतिकी अनुप्रयोगों में बिखरने के अनुप्रस्थ काट की भविष्यवाणी के लिए महत्वपूर्ण है। हालांकि, इन गतिशील प्रक्रियाओं का प्रथम-सिद्धांत सिमुलेशन शास्त्रीय कंप्यूटिंग संसाधनों की घातीय लागत से सीमित है। यह पेपर पहली बार प्रथम परिमाणीकरण का उपयोग करके पूर्ण अग्रणी-क्रम (LO) पायनलेस प्रभावी क्षेत्र सिद्धांत (EFT) हैमिल्टनियन परमाणु गतिशीलता के संसाधन आवश्यकताओं को पूरी तरह से चिन्हित करता है, उत्पाद सूत्र (product formulas) और क्वांटम सिग्नल प्रोसेसिंग (QSP) दोनों रणनीतियों का उपयोग करके। अध्ययन से पता चलता है कि ऐसे हैमिल्टनियन का समय विकास कणों की संख्या में बहुपद संसाधनों और एकल-कण आधार अवस्थाओं की संख्या में लॉगरिदमिक संसाधनों के साथ प्राप्त किया जा सकता है, जो द्वितीय परिमाणीकरण विधियों की तुलना में घातीय सुधार प्रदान करता है। कम ऊर्जा परमाणु बिखरने का दिलचस्प सिमुलेशन संभवतः केवल दसियों मिलियन T गेट्स और कुछ सौ तार्किक क्यूबिट्स की आवश्यकता होगी, जो सरल परमाणु प्रतिक्रियाओं के अध्ययन को प्रारंभिक दोष-सहिष्णु क्वांटम प्लेटफॉर्म के लिए व्यावहारिक बनाता है।
यह अनुसंधान परमाणु प्रणालियों की गतिशील प्रक्रियाओं की क्वांटम सिमुलेशन संसाधन आवश्यकताओं की समस्या को हल करने का लक्ष्य रखता है, विशेष रूप से:
कम ऊर्जा परमाणु बिखरने और परमाणु प्रतिक्रिया प्रक्रियाओं को कुशलतापूर्वक कैसे सिमुलेट करें परमाणु गतिशीलता सिमुलेशन में शास्त्रीय कंप्यूटिंग की घातीय संसाधन बाधा को कैसे दूर करें प्रारंभिक दोष-सहिष्णु क्वांटम कंप्यूटर के लिए परमाणु भौतिकी अनुप्रयोग कैसे डिजाइन करें वैज्ञानिक महत्व : परमाणु गतिशीलता सिमुलेशन तारकीय परमाणु संश्लेषण, न्यूट्रिनो गुण प्रयोगों, न्यूट्रिनो-रहित दोहरे बीटा क्षय जैसे मौलिक भौतिकी प्रश्नों को समझने के लिए महत्वपूर्ण हैकंप्यूटिंग चुनौती : हालांकि परमाणु संरचना के प्रथम-सिद्धांत तरीकों में महत्वपूर्ण प्रगति हुई है, परमाणु प्रतिक्रियाओं का पूर्ण गतिशील विवरण अभी भी सार्वभौमिक, नियंत्रणीय और स्केलेबल सिमुलेशन रणनीति की कमी हैक्वांटम लाभ की संभावना : क्वांटम कंप्यूटिंग गतिशील प्रक्रिया सिमुलेशन में निकट भविष्य का सबसे महत्वपूर्ण प्रभाव प्रदान करने की उम्मीद हैद्वितीय परिमाणीकरण विधि की समस्याएं (जैसे Watson et al. 40 का काम):
उच्च मेमोरी आवश्यकता : क्यूबिट्स की संख्या एकल-कण स्पेस आकार Ω के साथ रैखिक रूप से बढ़ती है, O(Ω)उच्च गेट लागत : 8×8×8 जाली के लिए, कम से कम 2048-3072 क्यूबिट्स की आवश्यकता हैमजबूत आयतन निर्भरता : T गेट्स की संख्या कम से कम आयतन पर रैखिक रूप से निर्भर है, O(Ω), परमाणु बिखरने की समस्याओं के लिए (विरल सीमा η/Ω ≪ 1) अक्षम हैविशेष रूप से, द्वितीय-क्रम उत्पाद सूत्र की संसाधन आवश्यकता है:
C 2 n d = O ( t 3 / 2 η 1 / 2 ϵ 1 / 2 Ω log ( t η Ω ϵ ) ) C^{2nd} = O\left(\frac{t^{3/2}\eta^{1/2}}{\epsilon^{1/2}}\Omega \log\left(\frac{t\eta\Omega}{\epsilon}\right)\right) C 2 n d = O ( ϵ 1/2 t 3/2 η 1/2 Ω log ( ϵ t η Ω ) )
विरलता का उपयोग : परमाणु प्रतिक्रिया सिमुलेशन में कणों की संख्या η निश्चित है और एकल-कण अवस्थाओं की संख्या Ω से बहुत कम है, प्रथम परिमाणीकरण इस विशेषता का उपयोग कर सकता हैसंसाधन आवश्यकता में कमी : क्यूबिट्स की संख्या O(η log Ω) और गेट लागत की लॉगरिदमिक निर्भरता प्राप्त करने की अपेक्षा करेंप्रारंभिक अनुप्रयोग व्यावहारिकता : आने वाले दोष-सहिष्णु क्वांटम कंप्यूटर के लिए व्यावहारिक परमाणु भौतिकी अनुप्रयोग प्रदान करेंइस पेपर के मुख्य योगदान में शामिल हैं:
पहली बार पूर्ण संसाधन चिन्हन : प्रथम परिमाणीकरण के तहत पूर्ण LO पायनलेस EFT हैमिल्टनियन (तीन-शरीर इंटरैक्शन सहित) के क्वांटम सिमुलेशन के लिए पूर्ण संसाधन अनुमान प्रदान करता हैघातीय सुधार : द्वितीय परिमाणीकरण विधि की तुलना में, आयतन Ω पर घातीय सुधार:क्यूबिट्स की संख्या: O(Ω) → O(η log Ω) T गेट्स की संख्या (द्वितीय-क्रम Trotter): O(Ω) → Õ(η^{3/2} log Ω) कई एल्गोरिथम कार्यान्वयन : तीन प्रकार के एल्गोरिदम के विस्तृत कार्यान्वयन और संसाधन विश्लेषण प्रदान करता है:प्रथम-क्रम, द्वितीय-क्रम, चतुर्थ-क्रम उत्पाद सूत्र (Product Formulas) क्वांटम सिग्नल प्रोसेसिंग (QSP) सामान्यीकृत क्वांटम सिग्नल प्रोसेसिंग (GQSP) नवीन ब्लॉक एन्कोडिंग योजना : संपर्क इंटरैक्शन (contact interactions) के लिए नई ब्लॉक एन्कोडिंग विधि डिजाइन की, फर्मियन सांख्यिकी विशेषताओं का उपयोग करके O(η) गेट लागत प्राप्त कीफर्मियन सेमिनॉर्म सामान्यीकरण : द्वितीय परिमाणीकरण में फर्मियन सेमिनॉर्म की अवधारणा को प्रथम परिमाणीकरण तक सामान्यीकृत किया, Trotter त्रुटि सीमा के सटीक अनुमान के लिएव्यावहारिक पैरामीटर मूल्यांकन : वास्तविक परमाणु भौतिकी पैरामीटर का उपयोग करके संसाधन अनुमान, 16 परमाणु नाभिकों का सिमुलेशन लगभग 10^7 T गेट्स और ~200 तार्किक क्यूबिट्स के साथ पूरा किया जा सकता हैइनपुट :
η परमाणु नाभिकों की प्रारंभिक क्वांटम अवस्था पायनलेस EFT हैमिल्टनियन H = T + V₂ + V₃ विकास समय t और लक्ष्य त्रुटि ϵ आउटपुट :
अनुमानित समय विकास ऑपरेटर U(t) जैसे कि ‖U(t) - e^{-iHt}‖ < ϵ बाधा शर्तें :
स्पेस जाली: d-आयामी, प्रत्येक आयाम में M = 2^m जाली बिंदु, आवधिक सीमा शर्तें एकल-कण अवस्थाओं की संख्या: Ω = 4·M^d (4 स्पिन और आइसोस्पिन से आता है) विरल सीमा: η/Ω ≪ 1 प्रथम परिमाणीकरण एन्कोडिंग योजना :
प्रत्येक कण को log₂(Ω) क्यूबिट्स द्वारा प्रतिनिधित्व किया जाता है स्थिति एन्कोडिंग: d·log₂(M) = dm क्यूबिट्स स्पिन/आइसोस्पिन: 2 क्यूबिट्स कुल क्यूबिट्स की संख्या: Q = (dm + 2)η गतिज ऊर्जा ऑपरेटर (गति स्पेस में विकर्ण):
T = K ( Q F T ⊗ d η ) † ∑ i = 0 η − 1 ∑ w = 0 d − 1 ∑ p = 0 2 m − 1 q w 2 Π w , i ( p ) ( Q F T ⊗ d η ) T = K(QFT^{\otimes d\eta})^\dagger \sum_{i=0}^{\eta-1}\sum_{w=0}^{d-1}\sum_{p=0}^{2^m-1} q_w^2 \Pi_{w,i}(p) (QFT^{\otimes d\eta}) T = K ( QF T ⊗ d η ) † ∑ i = 0 η − 1 ∑ w = 0 d − 1 ∑ p = 0 2 m − 1 q w 2 Π w , i ( p ) ( QF T ⊗ d η )
जहां K = ℏ²/(2μa²)(2π/2^m)², q_w आवधिक सीमा शर्तों के बाद गति घटक है।
संभावित ऊर्जा ऑपरेटर (स्थिति स्पेस में विकर्ण):
V 2 = C 2 ∑ i , j ≠ i ∑ r ⃗ i , r ⃗ j δ r ⃗ i , r ⃗ j Π i ( r ⃗ i ) Π j ( r ⃗ j ) V_2 = \frac{C}{2}\sum_{i,j\neq i}\sum_{\vec{r}_i,\vec{r}_j} \delta_{\vec{r}_i,\vec{r}_j}\Pi_i(\vec{r}_i)\Pi_j(\vec{r}_j) V 2 = 2 C ∑ i , j = i ∑ r i , r j δ r i , r j Π i ( r i ) Π j ( r j ) V 3 = G 6 ∑ i , j ≠ i , k ≠ i , j ∑ r ⃗ i , r ⃗ j , r ⃗ k δ r ⃗ i , r ⃗ j δ r ⃗ i , r ⃗ k Π i ( r ⃗ i ) Π j ( r ⃗ j ) Π k ( r ⃗ k ) V_3 = \frac{G}{6}\sum_{i,j\neq i, k\neq i,j}\sum_{\vec{r}_i,\vec{r}_j,\vec{r}_k} \delta_{\vec{r}_i,\vec{r}_j}\delta_{\vec{r}_i,\vec{r}_k}\Pi_i(\vec{r}_i)\Pi_j(\vec{r}_j)\Pi_k(\vec{r}_k) V 3 = 6 G ∑ i , j = i , k = i , j ∑ r i , r j , r k δ r i , r j δ r i , r k Π i ( r i ) Π j ( r j ) Π k ( r k )
विकास ऑपरेटर को गतिज और संभावित ऊर्जा विकास के उत्पाद में विघटित करें:
e − i H t ≈ ( e − i T τ e − i V τ ) r , τ = t / r e^{-iHt} \approx (e^{-iT\tau}e^{-iV\tau})^r, \quad \tau = t/r e − i H t ≈ ( e − i T τ e − iV τ ) r , τ = t / r
रणनीति : गति स्पेस में गतिज ऊर्जा को विकर्ण करने के लिए QFT का उपयोग करें
चरण :
सभी कणों के स्थान रजिस्टर पर QFT लागू करें प्रत्येक स्थान आयाम और कण के लिए:
p² की गणना करें (वर्ग सर्किट की आवश्यकता है, लागत 4m² - 4m T गेट्स) phase kickback विधि का उपयोग करके विकर्ण यूनिटरी e^{-iKq²t} लागू करें व्युत्क्रम QFT लागू करें संसाधन लागत :
T T ( ϵ , t ) = O [ η m ( m + log ( t η λ T ϵ ) ) ] T_T(\epsilon, t) = O\left[\eta m\left(m + \log\left(\frac{t\eta\lambda_T}{\epsilon}\right)\right)\right] T T ( ϵ , t ) = O [ η m ( m + log ( ϵ t η λ T ) ) ]
जहां λ_T = dKη2^{2m-2} = (ℏ²dηπ²)/(2μa²) गतिज ऊर्जा मानदंड की ऊपरी सीमा है।
मुख्य नवाचार : V₂ और V₃ के कार्यान्वयन को मर्ज करें, पहचान का उपयोग करें:
V 3 ( r ⃗ 0 , . . . , r ⃗ η − 1 ) = G ∑ i < j δ r ⃗ i , r ⃗ j ∑ k > j δ r ⃗ j , r ⃗ k V_3(\vec{r}_0,...,\vec{r}_{\eta-1}) = G\sum_{i<j}\delta_{\vec{r}_i,\vec{r}_j}\sum_{k>j}\delta_{\vec{r}_j,\vec{r}_k} V 3 ( r 0 , ... , r η − 1 ) = G ∑ i < j δ r i , r j ∑ k > j δ r j , r k
कार्यान्वयन रणनीति :
प्रत्येक कण जोड़ी (i,j) के लिए:
dm CNOT और 1 C^{dm}X गेट का उपयोग करके यह जांचें कि क्या स्थिति समान है परिणाम को सहायक क्यूबिट में संग्रहीत करें RZ(2Ct) गेट लागू करें (दो-शरीर संभावना) प्रत्येक तीन-गुना (i,j,k), k>j के लिए:
r⃗_j = r⃗_k की जांच करें Toffoli गेट और सहायक बिट्स का उपयोग करके RZ(2Gt) लागू करें संसाधन लागत :
T V ( ϵ ) = O [ η 3 ( m + log ( η ϵ ) ) ] T_V(\epsilon) = O\left[\eta^3\left(m + \log\left(\frac{\eta}{\epsilon}\right)\right)\right] T V ( ϵ ) = O [ η 3 ( m + log ( ϵ η ) ) ]
dm+1 सहायक क्यूबिट्स की आवश्यकता है।
विघटन रूप :
U ( t ) = ( e − i T τ / 2 e − i V τ e − i T τ / 2 ) r U(t) = \left(e^{-iT\tau/2}e^{-iV\tau}e^{-iT\tau/2}\right)^r U ( t ) = ( e − i T τ /2 e − iV τ e − i T τ /2 ) r
त्रुटि सीमा : फर्मियन सेमिनॉर्म का उपयोग करते हुए, त्रुटि है
ϵ T r o t t e r ( 2 ) ≤ t 3 r 2 α 2 \epsilon_{Trotter}^{(2)} \leq \frac{t^3}{r^2}\alpha_2 ϵ T ro tt er ( 2 ) ≤ r 2 t 3 α 2
जहां
α 2 = η 24 ( d K 2 2 m − 2 ) 2 M ( 4 ) + η 12 ( d K 2 2 m − 2 ) M ( 3 ) 2 \alpha_2 = \frac{\eta}{24}(dK2^{2m-2})^2 M(4) + \frac{\eta}{12}(dK2^{2m-2})M(3)^2 α 2 = 24 η ( d K 2 2 m − 2 ) 2 M ( 4 ) + 12 η ( d K 2 2 m − 2 ) M ( 3 ) 2
M(s) संशोधित युग्मन स्थिरांक के बाद संभावित ऊर्जा की ऊपरी सीमा है।
चरण चयन : r = ⌈ 4 t 3 α 2 / ϵ ⌉ r = \lceil\sqrt{4t^3\alpha_2/\epsilon}\rceil r = ⌈ 4 t 3 α 2 / ϵ ⌉
कुल T गेट्स की संख्या :
T t r ( 2 ) ( ϵ , t ) = O ~ ( t 3 / 2 η 3 / 2 ϵ 1 / 2 log ( t η Ω ϵ ) ( η 2 + log Ω ) ) T_{tr}^{(2)}(\epsilon,t) = \tilde{O}\left(\frac{t^{3/2}\eta^{3/2}}{\epsilon^{1/2}}\log\left(\frac{t\eta\Omega}{\epsilon}\right)(\eta^2 + \log\Omega)\right) T t r ( 2 ) ( ϵ , t ) = O ~ ( ϵ 1/2 t 3/2 η 3/2 log ( ϵ t η Ω ) ( η 2 + log Ω ) )
लक्ष्य : यूनिटरी ऑपरेटर U_H का निर्माण करें जैसे कि a⟨0|U_H|0⟩_a = H/λ_H
गतिज ऊर्जा ब्लॉक एन्कोडिंग (Lemma 4):
Su et al. 51 की रणनीति का पालन करें PREPARE: अवस्था तैयार करें |+⟩|i⟩_f|ω⟩_g∑_r 2^{r/2}|r⟩∑_s 2^{s/2}|s⟩ SELECT: b रजिस्टर को सशर्त रूप से फ्लिप करें, जब तक गति बिट्स मेल न खाएं एक-मानदंड: λ_T = dKη2^{2(m-1)} लागत: O(ηm(m + log(m/ϵ))) T गेट्स संभावित ऊर्जा ब्लॉक एन्कोडिंग (Lemma 5):
नवाचार बिंदु : पहचान का उपयोग करें
V = ( C 2 − G 6 ) ∑ i , j ≠ i δ r ⃗ i , r ⃗ j + G 6 ( ∑ i , j ≠ i δ r ⃗ i , r ⃗ j ) 2 V = \left(\frac{C}{2} - \frac{G}{6}\right)\sum_{i,j\neq i}\delta_{\vec{r}_i,\vec{r}_j} + \frac{G}{6}\left(\sum_{i,j\neq i}\delta_{\vec{r}_i,\vec{r}_j}\right)^2 V = ( 2 C − 6 G ) ∑ i , j = i δ r i , r j + 6 G ( ∑ i , j = i δ r i , r j ) 2
कार्यान्वयन :
ऑरेकल U_match को परिभाषित करें: |0⟩S|r⃗₀⟩···|r⃗ {η-1}⟩ → |Λ⟩S|r⃗₀⟩···|r⃗ {η-1}⟩
Λ पहले कण के समान स्थिति में कणों की संख्या गिनता है (अधिकतम 3, 2 क्यूबिट्स की आवश्यकता) PREPARE: अवस्था तैयार करें |+⟩_b|i⟩_f|φ⟩_l|κ⟩_p|κ⟩_q
|φ⟩_l C और G के सापेक्ष वजन को एन्कोड करता है |κ⟩ = (|0⟩ + √2|1⟩)/√3 वर्ग पद के लिए उपयोग किया जाता है SELECT: U_match और सशर्त Z गेट का उपयोग करके कार्यान्वयन करें एक-मानदंड : λ_V = η(3|C| + 4G)/2
लागत : O(ηm + log(1/ϵ)) T गेट्स, dm+4 सहायक क्यूबिट्स
मूल निर्माण :
W = ((2|P⟩⟨P| - 1)⊗1)U_S को परिभाषित करें V_φ = (e^{-iφZ/2}⊗1)V_0(e^{iφZ/2}⊗1}) को परिभाषित करें V_φ के बार-बार अनुप्रयोग के माध्यम से बहुपद P_Q(t) ≈ e^{it} को कार्यान्वित करें बहुपद क्रम : Jacobi-Anger विस्तार का उपयोग करते हुए, की आवश्यकता है
Q ( t , ϵ t ) = ⌈ 2 λ H ∣ t ∣ + 3 log ( 6 / ϵ t ) ⌉ Q(t,\epsilon_t) = \lceil 2\lambda_H|t| + 3\log(6/\epsilon_t)\rceil Q ( t , ϵ t ) = ⌈ 2 λ H ∣ t ∣ + 3 log ( 6/ ϵ t )⌉
पुनरावृत्तियां
कुल संसाधन :
T Q S P ( ϵ , t ) = O ( ( η t + log ( 1 / ϵ ) ) log ( Ω ) ( η + log ( m η t / ϵ ) ) ) T_{QSP}(\epsilon,t) = O\left((ηt + \log(1/\epsilon))\log(\Omega)(η + \log(mηt/\epsilon))\right) T QSP ( ϵ , t ) = O ( ( η t + log ( 1/ ϵ )) log ( Ω ) ( η + log ( m η t / ϵ )) )
क्यूबिट्स की संख्या : O(η + m + log(t/ϵ))
अनुकूलन : W†W = U_S के संबंध का उपयोग करते हुए, लगभग आधी SELECT कॉल को कम करें
संसाधन : QSP की तुलना में लगभग 2 गुना गेट लागत में कमी, क्यूबिट्स की आवश्यकता समान है
फर्मियन सेमिनॉर्म सामान्यीकरण :‖O‖_A = ‖Π_A O Π_A‖ को परिभाषित करें, Π_A एंटीसिमेट्रिक सबस्पेस प्रक्षेपण है पाउली अपवर्जन सिद्धांत का उपयोग करें: प्रत्येक जाली बिंदु पर अधिकतम 4 फर्मियन कम्यूटेटर मानदंड को सटीकता से अनुमानित करें, Trotter चरणों को अनुकूलित करें संभावित ऊर्जा कुशल कार्यान्वयन :V₂ और V₃ को मर्ज करें दोहराई गई गणना से बचने के लिए सहायक बिट्स पर रोटेशन लागू करें कारक 2 बचाएं फर्मियन व्यस्तता संख्या सीमा का उपयोग करके सर्किट को सरल बनाएं ब्लॉक एन्कोडिंग नवाचार :संभावित ऊर्जा ब्लॉक एन्कोडिंग V₃ की वर्ग संरचना का उपयोग करता है U_match ऑरेकल को केवल O(η) लागत की आवश्यकता है QSP के साथ संगतता के लिए SELECT स्व-व्युत्क्रमता बनाए रखें चरण ढाल अवस्था पुनः उपयोग :QFT और विकर्ण यूनिटरी के लिए समान चरण रजिस्टर का उपयोग करें सहायक क्यूबिट्स ओवरहेड को कम करें पिछले काम 39,40 के साथ सुसंगत पैरामीटर का उपयोग करें:
जाली अंतराल: a = 1.4 fm परमाणु नाभिक द्रव्यमान: μ = 939 MeV युग्मन स्थिरांक (SU(4) सममिति):
ℏ²/(2μa²) = 10.58 MeV C = -98.23 MeV G = 127.84 MeV स्पेस आयाम: d = 3 परिदृश्य 1: पार करने का समय (crossing time) t c r o s s ( L ) = a L ℏ c μ 2 E t_{cross}(L) = \frac{aL}{\hbar c}\sqrt{\frac{\mu}{2E}} t cross ( L ) = ℏ c a L 2 E μ
जहां L = 2^m जाली रैखिक आकार है, E = 10 MeV परमाणु नाभिक ऊर्जा है।
उद्देश्य: परमाणु टकराव या अर्ध-एक्सक्लूसिव प्रक्रिया सिमुलेशन लागत का मूल्यांकन करें परिदृश्य 2: प्रतिक्रिया संकल्प समय t r ( Δ ω ) = ( ⌈ Δ H Δ ω ⌉ − 1 ) 2 π Δ H t_r(\Delta\omega) = \left(\left\lceil\frac{\Delta H}{\Delta\omega}\right\rceil - 1\right)\frac{2\pi}{\Delta H} t r ( Δ ω ) = ( ⌈ Δ ω Δ H ⌉ − 1 ) Δ H 2 π
जहां ΔH = ‖T‖ + ‖V‖ + 18η, Δω = 100 MeV।
उद्देश्य: रैखिक प्रतिक्रिया विधि परमाणु अनुप्रस्थ काट गणना की लागत का मूल्यांकन करें द्वितीय परिमाणीकरण Trotter 40 : Watson et al. का कार्यान्वयनप्रथम परिमाणीकरण द्वितीय-क्रम Trotter : इस पेपर का कार्यान्वयनप्रथम परिमाणीकरण QSP : इस पेपर का कार्यान्वयनप्रथम परिमाणीकरण GQSP : इस पेपर का कार्यान्वयनT गेट्स की संख्या : गैर-Clifford गेट लागत (मुख्य बाधा)तार्किक क्यूबिट्स की संख्या : क्वांटम मेमोरी आवश्यकतालक्ष्य त्रुटि : ϵ ∈ {10^{-1}, 10^{-3}}क्यूबिट एन्कोडिंग: स्थिति dm क्यूबिट्स, स्पिन/आइसोस्पिन 2 क्यूबिट्स QFT कार्यान्वयन: अनुमानित QFT 85 का उपयोग करें, सटीकता O(m log(m/ϵ)) रोटेशन संश्लेषण: Kliuchnikov et al. 86 विधि का उपयोग करें बहु-नियंत्रण गेट्स: Gidney निर्माण 61 का उपयोग करें, प्रत्येक Toffoli 4 T गेट्स तालिका IV सारांश (t_, ϵ=0.1):
विधि परमाणु नाभिकों की संख्या T गेट्स क्यूबिट्स 2nd-quant Trotter 16 3.64×10⁹ 3072 2nd-quant Trotter 40 5.89×10⁹ 3072 1st-quant Trotter 16 9.3×10⁸ 206 1st-quant Trotter 40 2.34×10¹⁰ 422 1st-quant QSP 16 9.74×10⁷ 235 1st-quant QSP 40 5.81×10⁸ 500 1st-quant GQSP 16 3.75×10⁷ 235 1st-quant GQSP 40 2.11×10⁸ 500
मुख्य निष्कर्ष :
क्यूबिट लाभ स्पष्ट : प्रथम परिमाणीकरण सभी कण संख्याओं में द्वितीय परिमाणीकरण की तुलना में क्यूबिट आवश्यकता में बहुत कम है (~15 गुना अंतर)GQSP सर्वोत्तम : सभी परीक्षण परिदृश्यों में, GQSP सबसे कम गेट लागत प्रदान करता हैक्रॉसओवर बिंदु : प्रथम-क्रम Trotter η≈10-20 के बाद द्वितीय परिमाणीकरण Trotter की तुलना में अधिक महंगा हो जाता है (η^{3/2} निर्भरता के कारण)पार करने का समय सिमुलेशन (η=16, ϵ=10^{-2}):
जाली आकार: 8³ से 4096³ (Ω 2048 से 2.7×10¹¹ तक) द्वितीय परिमाणीकरण : T गेट्स 10⁸ से 10²² तक बढ़ते हैं (Ω^{3/2} वृद्धि)प्रथम-क्रम Trotter : T गेट्स 10⁸ से 10¹² तक बढ़ते हैं (Ω^{1/2} वृद्धि)GQSP : T गेट्स 10⁷ से 10¹¹ तक बढ़ते हैं (Ω^{1/3} वृद्धि)प्रतिक्रिया संकल्प सिमुलेशन (η=16, Δω=100 MeV):
द्वितीय परिमाणीकरण : रैखिक वृद्धि O(Ω)प्रथम-क्रम विधि : लॉगरिदमिक वृद्धि O(log²Ω)GQSP लाभ : जब Ω 8 परिमाण क्रम बढ़ता है, तो गेट लागत 5 गुना से कम बढ़ता हैचित्र 10 (8³ जाली, t_):
ϵ=10^{-1} पर: GQSP द्वितीय परिमाणीकरण Trotter से ~100 गुना तेज है ϵ=10^{-3} पर: GQSP लाभ अधिक स्पष्ट है (लॉगरिदमिक निर्भरता बनाम वर्गमूल निर्भरता) प्रथम-क्रम Trotter ϵ=10^{-1} और η=40 पर द्वितीय परिमाणीकरण से थोड़ा बेहतर है क्रॉसओवर बिंदु विश्लेषण :
प्रथम-क्रम बनाम द्वितीय-क्रम Trotter : η≈10-20 (छोटी जाली)QSP बनाम द्वितीय-क्रम Trotter : η≈200-255GQSP : सभी परीक्षण श्रेणियों (η≤294) में लाभ बनाए रखता हैस्केलिंग नियम सत्यापन :
प्रथम-क्रम Trotter: T ∝ η^{3/2} (सिद्धांत के अनुरूप) QSP/GQSP: T ∝ η² (सिद्धांत के अनुरूप) द्वितीय परिमाणीकरण: T ∝ η^{1/2} (सिद्धांत के अनुरूप) न्यूनतम सिमुलेशन उदाहरण (¹⁶O, 8³ जाली, ϵ=0.1):
GQSP : 3.75×10⁷ T गेट्स, 235 क्यूबिट्सतुलना : Fermi-Hubbard QPE की तुलना में ~10 गुना अधिक महंगा, आणविक आधार अवस्था सिमुलेशन के समाननिष्कर्ष : प्रारंभिक दोष-सहिष्णु क्वांटम कंप्यूटर व्यावहारिक हैद्वितीय परिमाणीकरण विधि :Roggero et al. 39 : पहली बार पायनलेस EFT क्वांटम सिमुलेशन, O(Ω) संसाधन Watson et al. 40 : सुधारी गई Trotter सूत्र, O(η^{1/2}Ω) गेट लागत Jordan-Wigner या Verstraete-Cirac एन्कोडिंग का उपयोग करें प्रथम परिमाणीकरण अन्वेषण :Roggero et al. 39 : triton खिलौना मॉडल का पहला प्रदर्शन Weiss et al. 53 : तीन-शरीर इंटरैक्शन के बिना पायनलेस EFT, O(η² log Ω) एकल-चरण लागत यह पेपर: पहली बार पूर्ण संसाधन चिन्हन, तीन-शरीर इंटरैक्शन सहित Abrams & Lloyd 47 : फर्मियन प्रणाली सिमुलेशन का प्रारंभिक काम Kassal et al. 48 : रासायनिक गतिशीलता बहुपद समय एल्गोरिथम Su et al. 51 : दोष-सहिष्णु रसायन विज्ञान सिमुलेशन, O(n log Ω) क्यूबिट्स Georges et al. 52 : मनमाने आधार समूह के साथ रसायन विज्ञान सिमुलेशन Suzuki 54 , Lloyd 55 : मूल सिद्धांत Childs et al. 63 : कम्यूटेटर स्केलिंग सिद्धांत Campbell 74 : Hubbard मॉडल प्रारंभिक दोष-सहिष्णु सिमुलेशन Su et al. 62 : फर्मियन सेमिनॉर्म Low & Chuang 56,57 : QSP और qubitization Gilyén et al. 65,66 : QSVT Motlagh & Wiebe 58 : GQSP Berry et al. 67 : GQSP दक्षता दोगुनी घातीय सुधार : प्रथम परिमाणीकरण आयतन निर्भरता में घातीय सुधार प्राप्त करता हैक्यूबिट्स: O(Ω) → O(η log Ω) गेट लागत: O(Ω) → O(log² Ω) (निश्चित समय के लिए) सर्वोत्तम एल्गोरिथम : GQSP सभी परीक्षण परिदृश्यों में सर्वोत्तम प्रदर्शन प्रदान करता हैगेट लागत: O((ηt + log(1/ϵ))log(Ω)(η + log(mηt/ϵ))) क्यूबिट्स: O(η log Ω + log(t/ϵ)) प्रारंभिक अनुप्रयोग व्यावहारिकता : सरल परमाणु प्रतिक्रिया सिमुलेशन (जैसे ¹⁶O) प्रारंभिक दोष-सहिष्णु क्वांटम कंप्यूटर पर संभव हो सकता है~10⁷ T गेट्स ~200 तार्किक क्यूबिट्स प्रयोज्यता श्रेणी : प्रथम परिमाणीकरण विरल सीमा (η/Ω ≪ 1) में द्वितीय परिमाणीकरण से बेहतर हैकण संख्या निर्भरता :प्रथम-क्रम विधि η पर अधिक निर्भर है (η^{3/2} बनाम η^{1/2}) बहुत बड़े η (>200) के लिए, द्वितीय परिमाणीकरण अधिक अनुकूल हो सकता है अधूरी लागत शामिल :प्रारंभिक अवस्था तैयारी लागत विस्तार से विश्लेषण नहीं की गई (एंटीसिमेट्रिक O(η log η log M)) अवलोकन योग्य माप और सर्किट पुनरावृत्ति संख्या का अनुमान नहीं लगाया गया केवल SU(4) सममिति इंटरैक्शन पर विचार करें पैरामीटर निर्भरता :संसाधन अनुमान युग्मन स्थिरांक C, G के प्रति संवेदनशील है जाली अंतराल a की पसंद सटीकता और लागत को प्रभावित करती है एल्गोरिथम अनुकूलन स्पेस :सीमित अंतर विधि अन्वेषण नहीं की गई (Trotter त्रुटि में सुधार हो सकता है) संभावित ऊर्जा घातीय O(η³) से O(η²) तक अनुकूलित किया जा सकता है गहराई अनुकूलन पर्याप्त रूप से विचार नहीं किया गया अंत-से-अंत संसाधन अनुमान :पूर्ण परमाणु बिखरने सिमुलेशन प्रवाह प्रारंभिक अवस्था तैयारी रणनीति अवलोकन योग्य निष्कर्षण विधि अधिक सटीक हैमिल्टनियन :SU(4) ब्रेकिंग शर्तें (सिंगलेट-ट्रिपलेट भेद) उच्च-क्रम EFT इंटरैक्शन उच्च ऊर्जा टकराव एल्गोरिथम सुधार :सीमित अंतर गतिज ऊर्जा प्रतिनिधित्व समानांतर रोटेशन गेट्स (Hamming Weight Phasing) लॉगरिदमिक गहराई बहु-नियंत्रण गेट्स 81 गहराई अनुकूलन बनाम गेट संख्या अनुकूलन ट्रेडऑफ अन्य एन्कोडिंग योजनाएं :सतह कोड के अलावा QEC कोड 82-84 विशिष्ट हार्डवेयर के लिए अनुकूलन भौतिकी अनुप्रयोग :न्यूट्रिनो-परमाणु बिखरना हल्के परमाणु संलयन प्रतिक्रिया परमाणु प्रतिक्रिया कार्य गणना सैद्धांतिक कठोरता :पूर्ण गणितीय व्युत्पत्ति और त्रुटि विश्लेषण फर्मियन सेमिनॉर्म का औपचारिक सामान्यीकरण सभी प्रमेय विस्तृत प्रमाण के साथ (परिशिष्ट) विधि नवाचार :संभावित ऊर्जा ब्लॉक एन्कोडिंग फर्मियन सांख्यिकी का उपयोग करके नवीन डिजाइन V₂ और V₃ विलय कार्यान्वयन का चतुर अनुकूलन चरण ढाल अवस्था पुनः उपयोग संसाधन में कमी प्रायोगिक पूर्णता :कई एल्गोरिदम की व्यवस्थित तुलना वास्तविक भौतिकी पैरामीटर का विस्तृत मूल्यांकन स्पष्ट संसाधन स्केलिंग नियम सत्यापन व्यावहारिक मूल्य :प्रारंभिक दोष-सहिष्णु क्वांटम कंप्यूटर के लिए व्यावहारिक अनुप्रयोग प्रदान करता है संसाधन अनुमान T गेट्स और क्यूबिट्स संख्या तक विशिष्ट है कोड कार्यान्वयन संभावना उच्च (एल्गोरिथम विवरण विस्तृत) लेखन गुणवत्ता :संरचना स्पष्ट, तर्क कठोर तकनीकी विवरण पर्याप्त लेकिन अनावश्यक नहीं चार्ट प्रभावी रूप से तर्क का समर्थन करते हैं श्रेणी सीमा :केवल SU(4) सममिति मामले पर विचार करें स्पिन-आइसोस्पिन ब्रेकिंग इंटरैक्शन शामिल नहीं है त्रि-आयामी स्पेस की धारणा (हालांकि सामान्यीकृत किया जा सकता है) अधूरी लागत अनुमान :प्रारंभिक अवस्था तैयारी का विस्तृत विश्लेषण नहीं अवलोकन योग्य माप लागत पर चर्चा नहीं की गई कुल "घड़ी समय" अनुमान अनुपस्थित है प्रायोगिक सत्यापन अनुपस्थित :छोटे पैमाने पर क्वांटम हार्डवेयर सत्यापन नहीं शास्त्रीय सिमुलेशन के साथ प्रत्यक्ष तुलना नहीं त्रुटि मॉडल का वास्तविक प्रदर्शन परीक्षण नहीं किया गया अनुकूलन संभावना पूरी तरह से अन्वेषण नहीं की गई :सीमित अंतर विधि केवल संक्षेप में उल्लिखित है समानांतर रणनीति विस्तार से चर्चा नहीं की गई विशिष्ट हार्डवेयर के लिए अनुकूलन अपर्याप्त है तुलना की निष्पक्षता :द्वितीय परिमाणीकरण के साथ तुलना करते समय, बाद में आगे अनुकूलन की गुंजाइश हो सकती है विभिन्न विधियों के त्रुटि मॉडल में अंतर हो सकता है क्षेत्र में योगदान :महत्वपूर्ण : पहली बार परमाणु गतिशीलता क्वांटम सिमुलेशन की प्रारंभिक व्यावहारिकता साबित करता हैपरमाणु भौतिकी क्वांटम एल्गोरिथम अनुसंधान के लिए बेंचमार्क प्रदान करता है भौतिकी सिमुलेशन में प्रथम परिमाणीकरण के अनुप्रयोग को आगे बढ़ाता है सैद्धांतिक महत्व :फर्मियन सेमिनॉर्म सामान्यीकरण में सार्वभौमिक मूल्य है ब्लॉक एन्कोडिंग तकनीक अन्य संपर्क इंटरैक्शन प्रणालियों पर लागू की जा सकती है संसाधन स्केलिंग नियम भविष्य के अनुसंधान के लिए मार्गदर्शन प्रदान करते हैं व्यावहारिक मूल्य :उच्च : संसाधन अनुमान से पता चलता है कि प्रारंभिक दोष-सहिष्णु क्वांटम कंप्यूटर व्यावहारिक हैप्रायोगिक भौतिकविदों को क्वांटम कंप्यूटिंग अनुप्रयोग संभावना प्रदान करता है क्वांटम हार्डवेयर विकास निवेश को प्रोत्साहित कर सकता है पुनरुत्पादनीयता :अच्छा : एल्गोरिथम विवरण विस्तृत, परिशिष्ट पूर्ण व्युत्पत्ति प्रदान करता हैपैरामीटर स्पष्ट, सीधे कार्यान्वयन योग्य सुझाव: ओपन सोर्स कोड कार्यान्वयन प्रभाव को बहुत बढ़ाएगा आदर्श अनुप्रयोग :कम ऊर्जा परमाणु बिखरना (E < m_π c²) हल्के परमाणु प्रतिक्रिया (η < 50) बड़ी आयतन आवश्यकता वाली भौतिकी समस्याएं (Ω ≫ η) सीमित परिदृश्य :उच्च घनत्व परमाणु पदार्थ (η/Ω छोटा नहीं) स्पिन-आइसोस्पिन निर्भरता की आवश्यकता वाली प्रक्रियाएं अत्यधिक सटीकता आवश्यकता (ϵ ≪ 10^{-3}) संभावित विस्तार :जाली QCD कम ऊर्जा प्रभावी सिद्धांत ठंडे परमाणु प्रणाली संपर्क इंटरैक्शन अन्य विरल फर्मियन प्रणाली यह कार्य क्वांटम लाभ पथ पर स्थिति:
अल्पकालीन : NISQ के लिए अनुपयुक्त (दोष-सहिष्णुता की आवश्यकता)मध्यकालीन : प्रारंभिक दोष-सहिष्णु क्वांटम कंप्यूटर के आदर्श अनुप्रयोगों में से एकदीर्घकालीन : परमाणु भौतिकी पूर्ण क्वांटम सिमुलेशन के लिए आधार स्थापित करता है39 Roggero et al., PRD 101, 074038 (2020) - पहली बार पायनलेस EFT क्वांटम सिमुलेशन
40 Watson et al., arXiv:2312.05344 - द्वितीय परिमाणीकरण सुधारी विधि
51 Su et al., PRX Quantum 2, 040332 (2021) - रसायन विज्ञान प्रथम परिमाणीकरण
56,57 Low & Chuang, PRL 118, 010501; Quantum 3, 163 - QSP सिद्धांत
62 Su et al., Quantum 5, 495 (2021) - फर्मियन सेमिनॉर्म
63 Childs et al., PRX 11, 011020 (2021) - Trotter त्रुटि सिद्धांत
67 Berry et al., PRA 110, 012612 (2024) - GQSP दक्षता दोगुनी
समग्र मूल्यांकन : यह परमाणु भौतिकी क्वांटम सिमुलेशन क्षेत्र में उच्च गुणवत्ता का सैद्धांतिक क्वांटम एल्गोरिथम पेपर है जो महत्वपूर्ण योगदान देता है। कठोर सैद्धांतिक विश्लेषण और विस्तृत संसाधन अनुमान के माध्यम से, यह पहली बार सरल परमाणु प्रतिक्रिया के क्वांटम सिमुलेशन की प्रारंभिक दोष-सहिष्णु क्वांटम कंप्यूटर पर व्यावहारिकता साबित करता है। विधि नवाचार मजबूत है, विश्लेषण पूर्ण है, और क्षेत्र विकास को आगे बढ़ाने में महत्वपूर्ण भूमिका है। मुख्य कमियां प्रायोगिक सत्यापन और अधूरी अंत-से-अंत लागत विश्लेषण में हैं, लेकिन ये बाद के काम में पूरक हो सकते हैं। उच्च स्तरीय क्वांटम कंप्यूटिंग या परमाणु भौतिकी पत्रिका में प्रकाशन की सिफारिश की जाती है।