DOI: https://doi.org/10.22331/q-2026-04-13-2060
تاريخ النشر: 2026-04-13
المؤلف: Samuele Piccinelli وآخرون
الموضوع الرئيسي: خوارزميات وهندسة الحوسبة الكمومية
نظرة عامة
في هذا البحث، يقدم المؤلفون إطارًا شاملاً لحساب دوال غرين المتأخرة (RGFs) على الحواسيب الكمومية من خلال إعادة صياغة المشكلة كمسألة تمايز الدوائر. تستفيد هذه الطريقة من التطور الزمني الحقيقي ومكونات الدائرة المصممة خصيصًا، والتي تُسمى اضطرابات الدائرة، والتي تمثل قوى اضطرابية خارجية في سياق الاستجابة الخطية. من خلال إنشاء ارتباط مباشر بين مشتقات الدائرة وحسابات RGF، يسمح الإطار بتطبيق تقنيات تمايز متنوعة، بما في ذلك مقدرات عشوائية لا تتطلب اتصالًا إضافيًا بين الكيوبتات لعمليات التطور الزمني. يتم إثبات فعالية هذه الطريقة من خلال التحقق العددي على نماذج الدوران والتفاعل، مما يؤدي إلى علاقات ديناميكية دقيقة حتى في ظل ظروف الضوضاء الواقعية.
يستنتج المؤلفون أن صياغتهم الموحدة تبسط تقييم RGFs عند نقاط زمنية متقطعة إلى تقدير مشتقات ناتج دائرة كمومية مُعلمة بالنسبة لبارامترات الاضطراب القابلة للتحكم. لا يوضح هذا الإطار فقط قابلية تطبيق أدوات تمايز الدوائر المتنوعة على حسابات الاستجابة المتأخرة، بل يوفر أيضًا مسارًا لتطوير عائلات من مقدرات RGF من تمثيل اضطراب مشترك. تتماشى النتائج بشكل جيد مع معايير التشخيص الدقيقة، والإطار جاهز للتكامل مع تقنيات تقدير التدرجات في الأنظمة المقاومة للأخطاء، مما قد يقلل من أعباء العينة. علاوة على ذلك، يفتح الربط الصريح بين دوال الاستجابة ومشتقات الدائرة آفاقًا للتقدم النظري والحاسوبي، بما في ذلك تحسين تصميمات المقدرات واستراتيجيات قياس فعالة لاستخراج العديد من الملاحظات.
مقدمة
تسلط مقدمة ورقة البحث الضوء على أهمية فهم الخصائص الديناميكية للأنظمة الكمومية ذات الجسيمات المتعددة عبر مجالات متنوعة، بما في ذلك الكيمياء الكمومية وفيزياء المادة المكثفة. تؤكد على التحديات التي تطرحها الارتباطات القوية والنمو السريع للتشابك أثناء التطور الزمني، مما يعقد التحقيق في الإثارات ذات الطاقة المنخفضة. لقد ظهرت الحوسبة الكمومية كحل واعد لهذه التحديات، مما يمكّن من حساب دوال الارتباط الزمني، المشار إليها بـ \( C(t – t_0) = \text{Tr}[\rho A(t) B(t_0)] \)، حيث تمثل \( \rho \) حالة متعددة الجسيمات مرجعية و\( A \) و\( B \) هما مشغلان ذوا صلة بعمليات فيزيائية محددة. هذه الدوال حاسمة لربط التنبؤات النظرية بالملاحظات التجريبية، مثل الدوال الطيفية وأطوال التشتت.
تقترح الورقة إطارًا شاملاً لحساب دوال غرين المتأخرة (RGFs) على المعالجات الكمومية، مقدمة خوارزميات كمومية جديدة مناسبة للأجهزة الضوضائية. يوسع المؤلفون التحليلات النظرية السابقة لتأسيس روابط بين دوال الاستجابة الخطية ومشتقات الدوائر الكمومية المعلمة. يسمح هذا الإطار بإعادة بناء دوال الارتباط المعتمدة على الزمن من عائلة دائرة معلمة واحدة، مما يسهل استخراج عوامل الهيكل الديناميكي والكميات الطيفية ذات الصلة. ستفصل الأقسام التالية من المخطوطة الأسس النظرية، وتقنيات تمايز الدوائر، ونتائج المحاكاة، مع الحفاظ على التركيز على القابلية العامة للتطبيق خارج الأنظمة الفيزيائية المحددة.
طرق
في هذا القسم، يقدم المؤلفون منهجية لتعريف فئة من المقدرات لدوال الاستجابة (RGFs) باستخدام مشتقات الدوائر، كما هو موضح في المعادلة (8). يتم توضيح النهج من خلال أمثلة محددة تستخدم دوائر كمومية معلمة للوصول إلى المشتقات الوظيفية اللازمة. يفترض الإطار أنه يمكن تهيئة سجل كمومي في الحالة الأرضية \(|\psi_0\rangle\) وأن كل من المروجات غير المضطربة والمضطربة يمكن تنفيذها بكفاءة، ربما من خلال صيغة المنتج في المحاكاة الكمومية الرقمية.
يستخدم المؤلفون أيضًا تحليل سوزوكي-تروتر لتفكيك تأثيرات الهاميلتونيين \(H_0\) و\(H_s(t’)\)، مما يؤدي إلى مروج إجمالي \(U_s(T, N) \approx \prod_{n=0}^{N-1} e^{-iH_0\tau} e^{-iH_s(n\tau)\tau}\)، حيث \(\tau = T/N\). يسمح ذلك ببناء الوظيفة \(F[s] := U_s^\dagger(T, N) o_R U_s(T, N)\)، والتي تمثل المحاكاة الكمومية الرقمية لقيمة التوقع \(\langle o_R(T) \rangle\) مع دالة قوة محددة \(s(r’, t’)\) عبر المواقع \(r’\) والأوقات \(t’ \in [0, T]\) لوقت تطور إجمالي \(T\).
نتائج
في هذا القسم، تُعرض نتائج المحاكاة التي تركز على الضوضاء الإحصائية، باستخدام وحدات طبيعية حيث يتم تعيين معلمات الاقتران والقفز إلى \( J \equiv 1 \). تشمل التحليلات مساهمات الجسيمين في عامل الهيكل الديناميكي (DSF)، والتي تعتبر حاسمة لفهم طيف الإثارة. تقارن المحاكاة الديناميات الدقيقة المستمدة من تطور تروتر مع النتائج من محاكاة الدوائر الكمومية لكل من طريقة الاضطراب المحلي للدائرة (LCP) وطريقة الاضطراب العشوائي للدائرة (SCP). من الجدير بالذكر أن تباين تقديرات التدرجات مستمد، مما يشير إلى أن تخصيص اللقطات الأمثل يفضل \( S = 1 \)، على الرغم من أن الاعتبارات العملية قد تتطلب \( S \geq 1 \).
تظهر النتائج أن SCP تعيد إنتاج الترددات ذات الصلة لـ DSF بفعالية حتى في ظل ظروف الضوضاء المسببة للانحلال الواقعية، مع الحفاظ على خطأ نسبي أقل من 1%. بالإضافة إلى ذلك، تُظهر طريقة SCP مزايا على LCP من حيث كفاءة العينة، خاصة مع زيادة عدد النقاط الزمنية. تشير النتائج إلى أن SCP يمكن أن تحقق تقارب التقديرات مع أخطاء متبقية تتناسب مع \( 1/\sqrt{S} \)، وأنها تسمح بالحساب المتوازي لجميع دوال الارتباط المطلوبة. بشكل عام، تتماشى النتائج العددية مع الأساليب التحليلية السابقة، مما يؤكد قوة طريقة SCP لمحاكاة الديناميات الكمومية في أنظمة الجسيمات المتعددة.
مناقشة
في هذا القسم، يناقش المؤلفون تقنيتين لتقييم العلاقات الديناميكية في الأنظمة الكمومية: الاضطراب المحلي للدائرة (LCP) والاضطراب المتزامن للدائرة (SCP). تستخدم طريقة LCP دالة دلتا لديراك للتعبير عن المشتق الوظيفي لكمية \( F \) بالنسبة لبارامتر الاضطراب \( s \). من خلال تطبيق اضطراب صغير في موقع وزمن محددين، يستنتج المؤلفون تعبيرًا لدالة غرين \( G \) التي تلتقط استجابة النظام للاضطرابات المحلية. يمكن تعزيز دقة هذه الطريقة باستخدام نهج الفرق المحدود المتماثل، مما يقلل من عدم الاستقرار العددي الذي يتم مواجهته عادةً في الحسابات الكمومية.
من ناحية أخرى، تسمح تقنية SCP بتقدير متزامن للتدرجات عبر اتجاهات اضطراب متعددة، مستفيدة من طرق العينة العشوائية. تعتبر هذه الطريقة مفيدة بشكل خاص لتجميع العينات على مدى طويل لدوال الاستجابة، حيث تتطلب فقط قالب دائرة واحدة ويمكنها التعامل بكفاءة مع اضطرابات متعددة بشكل متوازي. يبرز المؤلفون أن كلا الطريقتين متوافقتان مع الأجهزة الكمومية المعاصرة ويمكن تكييفها لمختلف الأنظمة الكمومية، بما في ذلك النماذج الفيرمونية. يختتم القسم بالتأكيد على إمكانيات هذه التقنيات في حساب كميات مثل عامل الهيكل الديناميكي (DSF)، والذي يعتبر حاسمًا لفهم الخصائص الديناميكية للأنظمة المكثفة.
DOI: https://doi.org/10.22331/q-2026-04-13-2060
Publication Date: 2026-04-13
Author(s): Samuele Piccinelli et al.
Primary Topic: Quantum Computing Algorithms and Architecture
Overview
In this research, the authors introduce a comprehensive framework for calculating Retarded Green’s Functions (RGFs) on quantum computers by reformulating the problem as one of circuit differentiation. This approach leverages real-time evolution and specially designed circuit components, termed circuit perturbations, which represent external perturbative forces in a linear-response context. By establishing a direct correlation between circuit derivatives and RGF computations, the framework allows for the application of various differentiation techniques, including stochastic estimators that do not necessitate additional qubit connectivity for time-evolution operations. The efficacy of this method is demonstrated through numerical validation on interacting spin and fermionic models, yielding accurate dynamical correlations even under realistic noise conditions.
The authors conclude that their unified formulation simplifies the evaluation of RGFs at discrete time points to the estimation of derivatives of a parameterized quantum circuit output concerning controllable perturbation parameters. This framework not only clarifies the applicability of diverse circuit-differentiation tools to retarded response calculations but also provides a pathway to develop families of RGF estimators from a common perturbation representation. The results align well with exact diagonalization benchmarks, and the framework is poised for integration with gradient-estimation techniques in fault-tolerant regimes, potentially minimizing sampling overheads. Furthermore, the explicit mapping between response functions and circuit derivatives opens avenues for theoretical and computational advancements, including improved estimator designs and measurement-efficient strategies for extracting multiple observables.
Introduction
The introduction of the research paper highlights the significance of understanding the dynamical properties of quantum many-body systems across various fields, including quantum chemistry and condensed matter physics. It emphasizes the challenges posed by strong correlations and rapid entanglement growth during time evolution, which complicate the investigation of low-energy excitations. Quantum computing has emerged as a promising solution to these challenges, enabling the computation of time correlation functions, denoted as \( C(t – t_0) = \text{Tr}[\rho A(t) B(t_0)] \), where \( \rho \) represents a reference many-body state and \( A \) and \( B \) are operators relevant to specific physical processes. These functions are crucial for linking theoretical predictions with experimental observables, such as spectral functions and scattering amplitudes.
The paper proposes a comprehensive framework for computing Retarded Green’s Functions (RGFs) on quantum processors, introducing novel quantum algorithms suitable for noisy devices. The authors extend previous theoretical analyses to establish connections between linear response functions and derivatives of parametrized quantum circuits. This framework allows for the reconstruction of time-dependent correlation functions from a single parametrized circuit family, facilitating the extraction of dynamical structure factors and related spectroscopic quantities. The subsequent sections of the manuscript will detail the theoretical underpinnings, circuit-differentiation techniques, and simulation results, while maintaining a focus on general applicability beyond specific physical systems.
Methods
In this section, the authors present a methodology for defining a class of estimators for response functions (RGFs) using circuit derivatives, as outlined in Equation (8). The approach is demonstrated through specific examples that utilize parametrized quantum circuits to access the necessary functional derivatives. The framework assumes that a quantum register can be initialized in the ground state \(|\psi_0\rangle\) and that both the unperturbed and perturbed propagators can be efficiently implemented, potentially through a product formula in digital quantum simulations.
The authors further employ a Suzuki-Trotter decomposition to factorize the actions of the Hamiltonians \(H_0\) and \(H_s(t’)\), leading to a total propagator \(U_s(T, N) \approx \prod_{n=0}^{N-1} e^{-iH_0\tau} e^{-iH_s(n\tau)\tau}\), where \(\tau = T/N\). This allows for the construction of the functional \(F[s] := U_s^\dagger(T, N) o_R U_s(T, N)\), which represents the digital quantum simulation of the expectation value \(\langle o_R(T) \rangle\) with a specified force function \(s(r’, t’)\) over positions \(r’\) and times \(t’ \in [0, T]\) for a total evolution time \(T\).
Results
In this section, the results of simulations focusing on statistical noise are presented, utilizing natural units where the coupling and hopping parameters are set to \( J \equiv 1 \). The analysis includes the two-spinon contributions to the dynamical structure factor (DSF), which are critical for understanding the excitation spectrum. The simulations compare the exact dynamics derived from Trotter evolution with results from quantum circuit simulations for both local circuit perturbation (LCP) and stochastic circuit perturbation (SCP) methods. Notably, the variance of gradient estimates is derived, indicating that optimal shot allocation favors \( S = 1 \), although practical considerations may necessitate \( S \geq 1 \).
The findings demonstrate that SCP effectively reproduces the relevant frequencies of the DSF even under realistic depolarizing noise conditions, maintaining a relative error below 1%. Additionally, the SCP method shows advantages over LCP in terms of sampling efficiency, particularly as the number of time points increases. The results indicate that SCP can achieve convergence of estimates with residual errors scaling as \( 1/\sqrt{S} \), and that it allows for parallel computation of all required correlation functions. Overall, the numerical results align with previous analytic approaches, confirming the robustness of the SCP method for simulating quantum dynamics in many-body systems.
Discussion
In this section, the authors discuss two techniques for evaluating dynamical correlations in quantum systems: Local Circuit Perturbation (LCP) and Simultaneous Circuit Perturbation (SCP). The LCP method utilizes a Dirac delta function to express the functional derivative of a quantity \( F \) with respect to a perturbation parameter \( s \). By applying a small perturbation at a specific position and time, the authors derive an expression for the Green’s function \( G \) that captures the system’s response to local perturbations. The accuracy of this method can be enhanced by using a symmetric finite difference approach, which reduces numerical instabilities commonly encountered in quantum computations.
On the other hand, the SCP technique allows for the simultaneous estimation of gradients across multiple perturbation directions, leveraging stochastic sampling methods. This approach is particularly advantageous for long-time sampling of response functions, as it requires only a single circuit template and can efficiently handle multiple perturbations in parallel. The authors highlight that both methods are compatible with contemporary quantum hardware and can be adapted for various quantum systems, including fermionic models. The section concludes by emphasizing the potential of these techniques for applications in calculating quantities like the Dynamical Structure Factor (DSF), which is crucial for understanding the dynamical properties of condensed matter systems.