DOI: https://doi.org/10.1103/frx3-41bz
تاريخ النشر: 2026-01-05
المؤلف: Lukas Cvitkovich وآخرون
الموضوع الرئيسي: ظواهر النقل الكمي والإلكتروني
نظرة عامة
في هذه الدراسة، تم حساب انقسامات الوادي في سوبرلاتيس Si/SiGe باستخدام نظرية الكثافة الوظيفية من أولى المبادئ (DFT)، والتي توفر إطارًا قويًا لنمذجة الواجهات، والضغوط، والاضطراب الذري بدقة دون الاعتماد على المعلمات التجريبية. تم تقييم نتائج DFT مقابل حسابات الربط الضيق (TB) ونظرية “2k0” ضمن تقريب الكتلة الفعالة (EM). تشير النتائج إلى أنه بينما تؤكد DFT الاستنتاجات الرئيسية لنظرية 2k0، فإنها تكشف أيضًا عن قيود في الطرق شبه التجريبية، لا سيما في تعاملها مع الاضطراب الذري. كشفت حسابات DFT عن تأثيرات خلط الوادي المداري الكبيرة عند انقسامات الوادي الأكبر، والتي لا تلتقطها نظرية 2k0 التقديرية الأصلية ولكن يمكن معالجتها في النماذج غير التقديرية.
أظهر مقارنة انقسامات الوادي DFT مع تلك التي تم الحصول عليها من حسابات TB أن كلا الطريقتين تتماشيان مع الاتجاهات العالمية المتوقعة، على الرغم من أن دقة TB تعتمد على ملف تركيز Ge. بشكل محدد، تؤدي TB بشكل أفضل مع ملفات تركيز أكثر سلاسة مقارنة بتلك التي تحتوي على تباينات سريعة. تنشأ الفجوات من الاختلافات في الأوصاف المجهرية للسبائك، حيث تكون DFT أكثر قدرة على حساب تحويلات الشحنات والتأثيرات بعيدة المدى. على الرغم من هذه القيود، تعتبر TB ونظرية 2k0 مناسبة للتحليلات الإحصائية لانقسامات الوادي في هياكل Si/SiGe المختلفة، لا سيما للتطبيقات في أجهزة الكيوبت الدوراني. تناقش الدراسة أيضًا مقياس بيانات TB مع حجم السوبرسيل وتبرز تأثير الاضطراب الذري على إحصائيات انقسام الوادي، وتخلص إلى أن TB يوفر تقريبًا معقولًا لهذه الظواهر.
مقدمة
تناقش مقدمة الورقة إمكانيات النقاط الكمومية في هياكل Si/SiGe لشبكات الكيوبت الدوراني على نطاق واسع في الحوسبة الكمومية. يبرز المؤلفون جودة الواجهات Si/SiGe المتبلورة مقارنة بواجهات السيليكون/الأكسيد، بالإضافة إلى القدرة على تنقية السيليكون والجرمانيوم نظائريًا للتخفيف من قيود التماسك الناتجة عن الدوران النووي. يشيرون إلى أنه تم تحقيق كيوبتات دورانية متماسكة وعالية الدقة، مع دقة البوابات التي تتجاوز عتبات تصحيح الأخطاء الكمومية. ومع ذلك، لا تزال التحديات قائمة بسبب الحالات المتقاربة تقريبًا بالقرب من حافة نطاق التوصيل، مما يمكن أن يؤدي إلى تباين في انقسامات الوادي التي تؤثر على تماسك الكيوبت ونقل الإلكترونات بين النقاط.
لمعالجة هذه التحديات، يستكشف المؤلفون تصاميم آبار كمومية مبتكرة تهدف إلى تعزيز انقسامات الوادي، مثل نتوءات Ge وتركيزات Ge المتذبذبة. يشيرون إلى نظرية “2k₀”، التي تفترض أن انقسام الوادي يتأثر بمكونات فورييه لإمكانات الاحتجاز. ومع ذلك، فإن هذا النموذج شبه التجريبي له قيود، لا سيما فيما يتعلق باضطراب السبيكة والتأثيرات غير المحلية. يقترح المؤلفون استخدام نظرية الكثافة الوظيفية من أولى المبادئ (DFT) كنموذج أكثر دقة للتحقيق في انقسامات الوادي في سوبرلاتيس Si/Si₀.₇Ge₀.₃، حيث تتضمن DFT تفاعلات معقدة دون الاعتماد على معلمات تجريبية. تتم مقارنة نتائج DFT مع بيانات الربط الضيق، مما يكشف عن رؤى حول قيود النماذج الحالية وأهمية خلط الوادي المداري في الأنظمة ذات انقسامات الوادي الكبيرة. تم هيكلة الورقة لتقديم الهياكل والنموذج أولاً، تليها مناقشة حول انقسامات الوادي في تصاميم آبار كمومية مختلفة.
الطرق
في هذا القسم، يتم تفصيل المنهجية المستخدمة لحساب الانقسام بين ودي ز في الهياكل [001] المسطحة. تركز الدراسة على سوبرلاتيس يتكون من آبار كمومية من السيليكون (Si) بسمك حوالي 10 نانومتر، متداخلة مع حواجز من Si$_{0.7}$Ge$_{0.3}$ بسمك حوالي 8.6 نانومتر. يتماشى تركيز Ge في الحواجز مع القيم التي لوحظت في الدراسات التجريبية الأخيرة. يتم نمذجة الهياكل كخلايا سوبر مع شروط حدود دورية، حيث يتم تعريف المنطقة في الطائرة كـ $S = n_a \times n_a$، مع $a$ تمثل معلمة الشبكة في الطائرة و$n$ عدد صحيح موجب.
يتم تحليل هياكل نموذجية مختلفة، تتضمن توزيعات عشوائية من ذرات Ge وفقًا لملف تركيز محدد. يتم تحديد العدد الإجمالي لذرات Ge في كل طبقة أحادية بواسطة تركيز Ge، مما يضمن الاتساق عبر تجسيدات مختلفة من السوبرسيل. تشمل النماذج تكوينات مع واجهات حادة وسلسة، بالإضافة إلى تركيزات Ge المتذبذبة. يتم تعديل معلمة الشبكة في الطائرة لتتناسب مع تلك الخاصة بوسادة افتراضية من Si$_{0.7}$Ge$_{0.3}$، مما يؤدي إلى إجهاد شد ثنائي المحاور قدره $\epsilon_{\parallel} = 1.15\%$ في طبقة Si النقية. يتم استرخاء المواضع الذرية داخل السوبرسيل لاستيعاب محتوى Ge، مما يؤدي إلى تشوهات محلية في الشبكة في جميع الأبعاد الثلاثة بسبب اضطراب السبيكة. يتم أخذ عينات من تجسيدات عشوائية متعددة لملف تركيز Ge لتقييم تأثير هذا الاضطراب على خصائص النظام.
النتائج
في هذا القسم، تركز النتائج على طاقة انقسام الوادي ($E_{VS}$) في الآبار الملساء والمتعرجة، مقارنة نظرية الكثافة الوظيفية (DFT) مع حسابات الربط الضيق (TB). تستخدم كلا الطريقتين سوبرلاتيس متطابقة، مع الحفاظ على نفس توزيع ومواقع ذرات Ge. من الجدير بالذكر أن فرق نطاق التوصيل (CBO) لنموذج TB قد تم تعديله من 150 مي eV إلى 225 مي eV ليتماشى مع نتائج DFT، مما يبرز الاختلافات الجوهرية بين الطريقتين الحسابيتين.
علاوة على ذلك، تستخدم الدراسة ملفات تعريف مجال كهربائي عمودي متسقة عبر كل من حسابات DFT وTB. يتم دمج الإمكانات الخارجية المستمدة من DFT في هاملتونيان TB، كما هو موضح في الملحق A. يتم تعريف المجالات الكهربائية العمودية ($F_z$) المشار إليها في هذا العمل على أنها تدرج الإمكانات هارتري داخل البئر السيليكوني، مما يضمن مقارنة متماسكة لتأثيرات المجالات الكهربائية على انقسام الوادي في الأنظمة التي تم تحليلها.
المناقشة
في هذا القسم، يناقش المؤلفون المنهجيات والنتائج المتعلقة بانقسامات الوادي في هياكل SiGe، باستخدام نظرية الكثافة الوظيفية (DFT) وحسابات الربط الضيق (TB). تم إجراء حسابات DFT باستخدام كود CP2K، حيث تم تحسين هندسة السوبرلاتيس وحساب هياكلها الإلكترونية، باستخدام الوظيفة التبادلية-التداخلية PBE للاسترخاء الهيكلي والوظيفة الهجينة PBE0 لحسابات انقسام الوادي. يتم تحديد انقسام الوادي، المشار إليه بـ $E_{VS}$، كفرق الطاقة بين أدنى شريحتين فرعيتين عند $k_{\perp} = 0$، دون النظر في اقتران الدوران-المدار. يبرز المؤلفون القيود الحسابية لـ DFT للأنظمة الكبيرة ويقارنون انقسامات الوادي المستمدة من كل من نماذج DFT وTB، مشيرين إلى أنه بينما تنتج كلا الطريقتين اتجاهات مشابهة، تنشأ الفجوات من معالجة اضطراب السبيكة وملفات تركيز Ge المحددة.
يستفيض المؤلفون في توضيح الفروق بين DFT ونظريات الكتلة الفعالة (EM)، لا سيما في سياق انقسامات الوادي. يقدمون إطارًا رياضيًا للتعبير عن انقسام الوادي ضمن نظرية EM ويؤكدون على قيود هذا النهج عند تطبيقه على الإمكانات المتغيرة بسرعة. تشير النتائج إلى أن DFT توفر فهمًا أكثر دقة للبنية الإلكترونية، لا سيما في حساب التأثيرات غير المحلية وتحويلات الشحنات التي تؤثر على مشهد الإمكانات. تكشف التحليلات أنه بينما يمكن لنماذج TB التقاط الاتجاهات الإحصائية في انقسامات الوادي بشكل فعال، إلا أنها قد لا تأخذ في الاعتبار التعقيدات التي يقدمها الاضطراب الذري، خاصة في الأنظمة ذات تدرجات التركيز الحادة. بشكل عام، تؤكد النتائج على أهمية استخدام كل من طرق DFT وTB للحصول على فهم شامل لانقسامات الوادي في هياكل SiGe، لا سيما في سياق تطبيقات النقاط الكمومية.
DOI: https://doi.org/10.1103/frx3-41bz
Publication Date: 2026-01-05
Author(s): Lukas Cvitkovich et al.
Primary Topic: Quantum and electron transport phenomena
Overview
In this study, valley splittings in Si/SiGe superlattices were computed using ab initio density functional theory (DFT), which offers a robust framework for accurately modeling interfaces, strains, and atomistic disorder without relying on empirical parameters. The results from DFT were benchmarked against atomistic tight-binding (TB) and the “2k0” theory within the effective mass (EM) approximation. The findings indicate that while DFT corroborates the primary conclusions of the 2k0 theory, it also exposes limitations in semi-empirical methods, particularly in their treatment of atomistic disorder. DFT calculations revealed significant valley-orbit mixing effects at larger valley splittings, which are not captured by the original perturbative 2k0 theory but can be addressed in non-perturbative models.
The comparison of DFT valley splittings with those obtained from TB calculations showed that both methods align with the expected global trends, although the accuracy of TB is contingent upon the Ge concentration profile. Specifically, TB performs better with smoother concentration profiles compared to those with rapid modulations. The discrepancies arise from the differing microscopic descriptions of alloys, with DFT being more adept at accounting for charge transfers and long-range effects. Despite these limitations, TB and the 2k0 theory are deemed suitable for statistical analyses of valley splittings in various Si/SiGe heterostructures, particularly for applications in spin qubit devices. The study also discusses the scaling of TB data with supercell size and highlights the influence of atomistic disorder on valley splitting statistics, concluding that TB provides a reasonable approximation for these phenomena.
Introduction
The introduction of the paper discusses the potential of quantum dots in Si/SiGe heterostructures for large-scale spin qubit arrays in quantum computing. The authors highlight the superior quality of epitaxial Si/SiGe interfaces compared to silicon/oxide interfaces, as well as the ability to isotopically purify silicon and germanium to mitigate nuclear spin-induced coherence limitations. They note that coherent and high-fidelity spin qubits have been achieved, with gate fidelities surpassing quantum error correction thresholds. However, challenges remain due to nearly degenerate states near the conduction band edge, which can lead to variability in valley splittings that affect qubit coherence and inter-dot electron shuttling.
To address these challenges, the authors explore innovative quantum well designs aimed at enhancing valley splittings, such as Ge spikes and oscillating Ge concentrations. They reference the “2k₀” theory, which posits that valley splitting is influenced by the confinement potential’s Fourier components. However, this semi-empirical model has limitations, particularly regarding alloy disorder and non-local effects. The authors propose using ab-initio density functional theory (DFT) as a more accurate modeling approach to investigate valley splittings in Si/Si₀.₇Ge₀.₃ superlattices, as DFT incorporates complex interactions without relying on empirical parameters. The results from DFT are compared with tight-binding data, revealing insights into the limitations of existing models and the significance of valley-orbit mixing in systems with large valley splittings. The paper is structured to first introduce the heterostructures and modeling framework, followed by a discussion of valley splittings in various quantum well designs.
Methods
In this section, the methodology employed to compute the splitting between Z valleys in planar [001] heterostructures is detailed. The study focuses on a superlattice comprising approximately 10 nm thick silicon (Si) quantum wells, interspersed with barriers of Si$_{0.7}$Ge$_{0.3}$ that are about 8.6 nm thick. The Ge concentration in the barriers aligns with values observed in recent experimental studies. The heterostructures are modeled as supercells with periodic boundary conditions, where the in-plane area is defined as $S = n_a \times n_a$, with $a$ representing the in-plane lattice parameter and $n$ a positive integer.
Various model structures are analyzed, incorporating random distributions of Ge atoms according to a specified concentration profile. The total number of Ge atoms in each monolayer is determined by the Ge concentration, ensuring consistency across different realizations of the supercell. The models include configurations with sharp and smooth interfaces, as well as oscillating Ge concentrations. The in-plane lattice parameter is adjusted to match that of a virtual Si$_{0.7}$Ge$_{0.3}$ buffer, resulting in a tensile, biaxial strain of $\epsilon_{\parallel} = 1.15\%$ in the pure Si layer. The atomic positions within the supercell are relaxed to accommodate the Ge content, leading to local lattice distortions in all three dimensions due to alloy disorder. Multiple random realizations of the Ge concentration profile are sampled to evaluate the impact of this disorder on the system’s properties.
Results
In this section, the results focus on the valley splitting energy ($E_{VS}$) in smoothed and wiggle wells, comparing Density Functional Theory (DFT) with Tight Binding (TB) calculations. Both methods utilize identical superlattices, maintaining the same distribution and atomic positions of Ge atoms. Notably, the conduction band offset (CBO) of the TB model has been adjusted from 150 meV to 225 meV to align with DFT findings, thereby emphasizing the intrinsic differences between the two computational approaches.
Furthermore, the study employs consistent vertical electric field profiles across both DFT and TB calculations. The external potential derived from DFT is integrated into the TB Hamiltonian, as detailed in Appendix A. The vertical electric fields ($F_z$) referenced in this work are defined as the gradient of the Hartree potential within the silicon well, ensuring a coherent comparison of the effects of electric fields on valley splitting in the analyzed systems.
Discussion
In this section, the authors discuss the methodologies and findings related to the valley splittings in SiGe heterostructures, utilizing Density Functional Theory (DFT) and Tight-Binding (TB) calculations. The DFT calculations, performed using the CP2K code, optimize the geometry of superlattices and compute their electronic structures, employing the PBE exchange-correlation functional for structural relaxation and the hybrid PBE0 functional for valley splitting calculations. The valley splitting, denoted as $E_{VS}$, is determined as the energy difference between the lowest two sub-bands at $k_{\perp} = 0$, without considering spin-orbit coupling. The authors highlight the computational limitations of DFT for large systems and compare the valley splittings derived from both DFT and TB models, noting that while both methods yield similar trends, discrepancies arise from the treatment of alloy disorder and the specific Ge concentration profiles.
The authors further elaborate on the differences between DFT and effective mass (EM) theories, particularly in the context of valley splittings. They introduce a mathematical framework to express valley splitting within EM theory and emphasize the limitations of this approach when applied to rapidly varying potentials. The results indicate that DFT provides a more nuanced understanding of the electronic structure, particularly in accounting for non-local effects and charge transfers that influence the potential landscape. The analysis reveals that while TB models can effectively capture statistical trends in valley splittings, they may not fully account for the complexities introduced by atomic disorder, especially in systems with sharp concentration gradients. Overall, the findings underscore the importance of using both DFT and TB methods to gain a comprehensive understanding of valley splittings in SiGe heterostructures, particularly in the context of quantum dot applications.