Efficient Simulation of Gate Leakage Currents in Ultra-scaled Transistors from First Principles
تقدم هذه الورقة طريقة حدود النقل الكمي (QTBM) الموسعة التي تحسب بكفاءة تيارات تسرب البوابة من المبادئ الأولية في ترانزستورات التأثير المجالي من نوع WS2-hBN فائقة القياس، وذلك من خلال التغلب على التحديات في الشروط الحدية المفتوحة والأنظمة الخطية المتفرقة غير المنظمة، متفوقةً بذلك على التقنيات الحديثة للأجهزة التي تتجاوز أطوال بواباتها 25 نانومتر.
المؤلفون الأصليون:Mauro Dossena, Nicolas Vetsch, Alexander Maeder, Anders Winka, Vincent Maillou, Denghui Lu, Alexandros Nikolaos Ziogas, Jiang Cao, Mathieu Luisier
في السعي الدؤوب لجعل الأجهزة الإلكترونية أصغر وأسرع، يواجه المهندسون عقبة فيزيائية جوهرية: البوابة. في الترانزستور الحديث، تعمل البوابة كمفتاح، يتحكم في تدفق الكهرباء عبر قناة. ولجعل المفتاح أكثر حساسية واستجابة، يجب جعل الطبقة العازلة تحت هذه البوابة رقيقة للغاية. ومع ذلك، عندما تصبح هذه الطبقة رقيقة جداً، تبدأ الإلكترونات في عملية "النفق". فبدلاً من الالتزام بقواعد المفتاح والبقاء ضمن مسارها المحدد، فإنها تنفق كمياً عبر العازل وتتسرب إلى البوابة. هذا التسرب يهدر الطاقة ويولد حرارة، مما قد يجعل الأجيال القادمة من الحواسيب فائقة الصغر غير فعالة أو حتى عديمة الفائدة. ولتصميم أجهزة أفضل، يجب أن يكون العلماء قادرين على التنبؤ بدقة بمقدار التيار الذي سيتسرب، لكن حساب ذلك للمواد الصغيرة والمعقدة المستخدمة في الرقائق من الجيل التالي كان لفترة طويلة كابوساً حاسوبياً.
لقد طور فريق من الباحثين في المعهد الفيدرالي السويسري للتكنولوجيا في زيورخ (ETH Zurich) الآن طريقة جديدة لحل هذه المشكلة، مما يسمح لهم بمحاكاة تسرب البوابة في ترانزستورات أكبر بكثير وأكثر واقعية مما كان ممكناً في السابق. يركز عملهم على نوع محدد من الترانزستورات المستقبلية المبنية من "ساندوتش" من المواد: قناة مصنوعة من ثنائي كبريتيد التنغستن، وطبقة عازلة من نيتريد البورون السداسي، وبوابة معدنية مصنوعة من البلاتين. في هذه الأجهزة، تقع البوابة مباشرة فوق القناة، مما يخلق مساحة شاسعة يمكن للإلكترونات أن تهرب منها. التحدي الذي يواجه علماء الحاسوب هو أن المعادلات الرياضية المطلوبة لتتبع هذه الإلكترونات الهاربة تصبح معقدة وكثيرة لدرجة أن الحواسيب الفائقة القياسية تنفد منها الموارد من الوقت والذاكرة قبل أن تتمكن من إنهاء الحساب، خاصة عندما يكون طول البوابة أكثر من بضعة نانومترات.
عالج الباحثون هذه العقبة عن طريق تغيير كيفية مقاربتهم لرياضيات المشكلة، حيث قاموا بتفكيكها إلى جزأين يمكن إدارتهما. أولاً، أدركوا أن التعامل مع البوابة بأكملها ككتلة واحدة ضخمة وصلبة هو أمر غير ضروري وغير فعال. وبدلاً من ذلك، قاموا بتقسيم البوابة إلى وحدات أصغر متكررة. ومن خلال حساب سلوك وحدة صغيرة واحدة ثم دمج هذه النتائج رياضياً، استطاعوا نمذجة البوابة بأكملها دون الحاجة إلى حل المعادلات المستحيلة للكل دفعة واحدة. يشبه هذا النهج كيف يمكن للمرء تقدير الوزن الإجمالي لسلسلة طويلة عن طريق وزن حلقة واحدة وضربها، بدلاً من رفع السلسلة بأكملها.
ثانياً، عالجوا مسألة البيانات غير المنظمة والفوضوية التي تنشأ عندما تتفاعل الإلكترونات مع البوابة. فالطرق القياسية لحل هذه المعادلات تخلق تأثير "الملء"، حيث تمتلئ ذاكرة الحاسوب ببيانات غير ضرورية، مما يؤدي إلى إبطاء كل شيء. قدم الفريق تقنية تبسط وصف تأثير البوابة، حيث تقوم بتجريد المعلومات الزائدة مع الحفاظ على الفيزياء الأساسية سليمة. سمح هذا لهم باستخدام خوارزمية أكثر كفاءة لحل المعادلات المتبقية، محولين مشكلة كانت ستتطلب ساعات أو أياماً من وقت الحوسبة إلى مشكلة يمكن حلها في دقائق.
نتائج هذه الطريقة الجديدة مذهلة. اختبر الباحثون نهجهم على ترانزستور بطول بوابة يبلغ 6.4 نانومتر، وقارنوا حساباتهم المبسطة مقابل محاكاة كاملة ودقيقة. طابقت النتائج المحاكاة بشكل شبه مثالي، مما أثبت أن اختصاراتهم لم تضح بالدقة. وبدافع من هذا التحقق، دفعوا بالمحاكاة إلى طول بوابة يبلغ 25.6 نانومتر، وهو حجم يمثل الأجهزة الواقعية ولكنه كان كبيراً جداً في السابق لمحاكاته بهذا المستوى من التفصيل. لقد نجحوا في حساب كيفية نفق الإلكترونات من المصدر والمصب إلى البوابة، وحددوا مستويات طاقة معينة يكون فيها التسرب أكثر احتمالاً للحدوث.
يُظهر هذا العمل أنه أصبح من الممكن الآن نمذجة تسرب البوابة في الترانزستورات شديدة التصغير بدقة عالية، حتى عندما تمتد البوابة على جزء كبير من القناة. ومن خلال الجمع بين هذين الابتكارين الرياضيين، أزال الفريق الحاجز الحسابي الذي كان يحصر مثل هذه الدراسات في هياكل صغيرة واصطناعية. تفتح هذه القدرة الباب أمام المهندسين لاختبار تركيبات مختلفة من المواد وتصاميم البوابات على الحاسوب قبل بنائها في المختبر، مما يضمن إمكانية جعل الترانزستورات المستقبلية أصغر وأكثر كفاءة دون الوقوع في مشكلة التيار المتسرب.
ملخص تقني: المحاكاة الفعالة لتيارات تسرب البوابة في الترانزستورات فائقة التصغير من المبادئ الأولية
بيان المشكلة مع استمرار تصغير أبعاد الترانزستور وسماكة الأكسيد المكافئة (EOT)، بدأت تيارات تسرب البوابة تبرز كعامل حدي حرج للتقنيات المستقبلية. وتعد هذه المشكلة حادة بشكل خاص في المفاتيح المنطقية التي تستخدم ثنائيات أشباه الموصلات ثنائية الأبعاد من مركبات المعادن الانتقالية (TMDCs)، مثل WS2، والتي تفتقر إلى الأكاسيد الذاتية. وبينما تم استخدام أدوات النقل الكمي (QT) المستمدة من المبادئ الأولى (ab initio) لدراسة الواجهات بين أشباه الموصلات والعوازل، إلا أنه لم يتم تطبيقها بفعالية على تسرب البوابة من قبل. عادةً ما تقوم أدوات النقل الكمي التقليدية بنمذجة حقن الحوامل عند المصدر والمصب فقط. إن إدخال مسار حقن ثالث عند البوابة — وهو منطقة فضاء كبيرة تمتد عبر القناة بأكملها — يخلق عقبتين حسابيتين هامتين:
الشروط الحدية المفتوحة (OBCs): يتطلب حساب طاقة الارتداد الذاتي ومصفوفات الحقن لتلامس البوابة حل مسائل القيم الذاتية غير الخطية (NEVPs). وتتضاعف التكلفة الحسابية لهذه المسائل بشكل تكعيبي (O(N3)) مع عدد المدارات، مما يجعل الحسابات لأقطاب بوابة طويلة باهظة التكلفة بشكل لا يُحتمل.
حل النظام الخطي: يؤدي وجود تلامس البوابة إلى تغيير نمط التشتت (sparsity pattern) لمصفوفة النظام (ES−H−ΣB)، مما يجعلها غير منتظمة الهيكل. ويتسبب هذا النمط غير المنتظم في حدوث "امتلاء" (fill-in) مفرط أثناء عملية التفكيك المباشر (LU factorization)، ويمنع التطبيق الفعال لتقنيات حذف غاوس الكتلي القياسية مثل خوارزمية الدوال الخضراء المتكررة (RGF).
المنهجية لمعالجة هذه التحديات، قام المؤلفون بتوسيع طريقة الحدود الناقلة الكمية (QTBM) باستخدام تقريبين أساسيين:
توسع بلوخ للشروط الحدية المفتوحة: للتخفيف من التوسع التكعيبي لمسألة القيم الذاتية غير الخطية (NEVP)، تم تقسيم تلامس البوابة الكامل (الخلية الكاملة، FC) إلى خلايا وحدة (UC) أصغر وقابلة للتكرار. ومن خلال افتراض الدورية في الاتجاهات المستعرضة، طبق المؤلفون توسع بلوخ. حيث قاموا بحساب طاقة الارتداد الذاتي للحد (σGB) ومصفوفة الحقن (injG) لخلية وحدة واحدة عبر شبكة من ناقلات الموجات المساعدة (ka,kb). ثم تمت إعادة بناء هذه المصفوفات لتشكيل مصفوفات التلامس الكاملة (ΣGB و InjG) باستخدام نواتج كرونيكر (Kronecker products). ومع أن هذا يعد تقريبًا دقيقًا للخلايا الدورية اللانهائية، إلا أنه يعمل كتقريب تم التحقق من صحته للخلايا ذات أطوال البوابات المحدودة.
التفكيك منخفض الرتبة لمصفوفة النظام: للتعامل مع التشتت غير المنتظم، تم تقريب مصفوفة طاقة الارتداد الذاتي للبوابة ΣGB باستخدام تفكيك منخفض الرتبة (ΣGB≈RGΛG−1ϕˉG)، حيث تكون الرتبة nG أصغر بكثير من أبعاد المصفوفة NG. ومن خلال الجمع بين هذا التفكيك ومتطابقة وودبري (Woodbury identity)، أعاد المؤلفون صياغة معادلة QTBM. يعمل هذا التعديل على نقل حدود طاقة الارتداد الذاتي من مصفوفة النظام إلى متجه الحقن في الطرف الأيمن. وبذلك يستعيد مصفوفة النظام نمط التشتت المنظم، ويصبح هرميًا (Hermitian)، مما يجعله قابلاً للتطبيق الفعال لحذف غاوس الكتلي (وتحديدًا خوارزمية توماس الكتلي).
المساهمات الرئيسية
توسيع إطار عمل QTBM: يقدم البحث توسعة مبتكرة لإطار عمل QTBM قادرة على التعامل مع الأجهزة متعددة الأطراف ذات تلامسات البوابة الكبيرة على مستوى المبادئ الأولى (ab initio).
الكفاءة الخوارزمية: يؤدي الجمع بين توسع بلوخ والتقريب منخفض الرتبة إلى تقليل التعقيد الحسابي لحسابات OBC، ويسمح بحل الأنظمة الخطية الكبيرة غير المنتظمة التي كانت مستعصية سابقًا.
التطبيق من المبادئ الأولى: تم تطبيق الطريقة على ترانزستور ذو تأثير مجالي (FET) من نوع WS2 أحادي الطبقة مع عازل hBN بسمك 1.8 نانومتر وبوابة من البلاتين (Pt)، باستخدام مصفوفات نظرية الكثافة الوظيفية (DFT) المولدة بواسطة CP2K ومعالجتها بواسطة الحل مفتوح المصدر quatrex.
النتائج تم التحقق من الطريقة المقترحة واختبار كفاءتها على هيكل ترانزستور WS2-hBN-Pt:
التحقق من الصحة: أظهرت المقارنات بين خط الأساس كامل الرتبة، والتقريب منخفض الرتبة، والنهج المشترك لتوسع بلوخ/الرتبة المنخفضة لطول بوابة (Lg) يبلغ 6.4 نانومتر، توافقًا ممتازًا في دالات النقل من المصدر إلى البوابة ومن المصب إلى البوابة، مما يؤكد صحة التقريبات.
القابلية للتوسع: نجحت الطريقة في محاكاة أجهزة بطول بوابة يصل إلى Lg=25.6 نانومتر. وفي المقابل، فإن معاملة البوابة ككتلة كاملة واحدة حدت من عمليات المحاكاة لتكون Lg≤12.8 نانومتر بسبب قيود الذاكرة والوقت.
مكاسب الأداء:
حساب OBC: وفرت طريقة توسع بلوخ تسريعًا يصل إلى 5.0× وتوفيرًا في الذاكرة يصل إلى 5.0× لطول بوابة Lg=12.8 نانومتر مقارنة بنهج الخلية الكاملة.
حل النظام: حقق التقريب منخفض الرتبة، الذي تم حله عبر خوارزمية توماس الكتلي المخصصة، تسريعًا في وقت التشغيل يصل إلى 7.8× وتوفيرًا في الذاكرة بنسبة 2.8× عند Lg=22.4 نانومتر مقارلة بالمحلل كامل الرتبة (cuDSS).
الرؤى الفيزيائية: كشفت المحاكاة عن قمم متميزة في دالات النقل من المصدر إلى البوابة ومن المصب إلى البوابة حول 7.55 إلكترون فولت، وهي تتوافق مع حالات حزمة التوصيل في قناة WS2. كما لوحظ أن النقل من المصب إلى البوابة أكبر من النقل من المصدر إلى البوابة، وهو ما يتفق مع هبوط الجهد الأكثر حدة في جانب المصب.
الأهمية يزعم المؤلفون أن هذا العمل يتيح الحساب من المبادئ الأولى (ab initio) لدالات نقل تسرب البوابة للأجهزة الإلكترونية النانوية الواقعية ذات القنوات ثنائية الأبعاد وأطوال بوابات تتجاوز 25 نانومتر. ومن خلال التغلب على العوائق الحسابية المرتبطة بتلامسات البوابة الكبيرة، يمهد هذا النهج الطريق لاستقصاء تأثير طبقات العزل المختلفة على خصائص التسرب في الترانزستورات ذات التأثير المجالي فائقة التصغير، وهي قدرة لم تكن متاحة سابقًا في حلول QT القياسية.