A stochastic simulation of the dislocation-mediated etching of porous GaN distributed Bragg reflectors

1. بيان المشكلة

تعد عملية التآكل الكهروكيميائي لنيتريد الغاليوم (GaN) عملية حاسمة لتصنيع عواكس برغ (DBRs) المستخدمة في الأجهزة الإلكترونية الضوئية (مثل مصابيح LED والليزر). في هذه الهياكل، يتم تكديم طبقات متبادلة من الـ GaN عالي النوع (n-doped) والـ GaN غير المخصب عمداً (NID). عملية التآكل انتقائية للناقلية، مما يعني أنه من المتوقع أن تتآكل طبقات الـ n-doped فقط، بينما تظل طبقات الـ NID سليمة.

ومع ذلك، يوجد تحدٍ جوهري: كيف يصل المادة المآكلة إلى طبقات الـ n-doped المدفونة عندما يكون السطح عبارة عن طبقة NID عازلة؟

• النموذج التقليدي ("نموذج الكباب"): كان يُفترض سابقاً أن الخللول الممتدة (Threading Dislocations - TDs) تعمل كأنابيب رأسية مستمرة (أنابيب نانوية) تسمح للمادة المآكلة بالتدفق مباشرة للأسفل عبر كل طبقة، مما يخلق هيكلاً يشبه "الكباب" من المسامية على طول كل خلل.
• الواقع المرصود: كشفت البيانات التجريبية الحديثة (باستخدام تصوير مقطع FIB العرضي وSTEM) أن "نموذج الكباب" غير كافٍ. بدلاً من ذلك، غالباً ما يتبع التآكل "آلية الشلال" (cascade mechanism): حيث تنتقل المادة المآكلة عبر خلل (TD)، ثم تتحرك عرضياً عبر طبقة مخصبة، وتتقاطع مع خلل (TD) آخر، وتستمر في النزول للأسفل. يؤدي هذا إلى هياكل مسامية معقدة وغير منتظمة حيث يمكن أن تكون الخلول نشطة في بعض الطبقات وغير نشطة في أخرى.
• الفجوة: إن الآليات الفيزيائية التي تتحكم في الانتقال بين مورفولوجيا "الكباب" و"الشلال"، وكيفية تأثير الجهد المطبق على هذه المسارات، ليست مفهومة تماماً. هناك نقص في النماذج التنبؤية التي تربط بيانات الكرونوأمبيرومتري (التيار مقابل الزمن) التجريبية مع البنى المجهرية ثلاثية الأبعاد الناتجة.

2. المنهجية

قام المؤلفون بتطوير محاكاة عشوائية (stochastic simulation) لنمذجة عملية التآكل الكهروكيميائي وقارنوا مخرجاتها بالبيانات التجريبية لـ DBRs من GaN المنماة على ركائز السيليكون.

• بنية المحاكاة:

  • بُنيت باستخدام لغة Python مع مكتبة numba لضمان الكفاءة عبر الترجمة الفورية (just-in-time compilation).
  • نموذج الشبكة: يتم تمثيل هيكل DBR كمصفوفة ثنائية الأبعاد من البكسلات (5 نانومتر × 5 نانومتر).
  • المكونات: يتم تعيين البكسلات كـ: مادة مآكلة، أو GaN غير مخصب (NID)، أو GaN مخصب بالسيليكون (Si-doped)، أو خلل ممتد (TD).
  • نمذجة الخلل (TD): تُنمذج الخلول كخطوط رأسية (بعرض بكسل واحد، وتباعد 500 نانومتر) تخترق التكدس.
  • منطق التآكل:
    • بكسلات NID لديها احتمالية 0% للتآكل.
    • بكسلات الـ doped والـ TD لها احتمالات تآكل محددة: $P_{doped}$ (احتمالية تآكل الطبقة المخصبة) و $P_{disloc}$ (احتمالية تآكل الخلل).
    • آلية الشلال: تسمح الخوارزمية للمادة المآكلة بالتحرك عرضياً. إذا تقاطع مسار عرضي مع خلل (TD) غير نشط، يمكن لهذا الخلل أن يصبح نشطاً، مما يخلق مسار "شلال".
    • القيود: يتم منع اندماج المسام المتجاورة لمحاكاة النطاقات غير المتقاطعة الملاحظة في المجهر.
  • حساب التيار: كل بكسل متآكل يولد وحدة شحنة، مما يسمح للمحاكاة بإنشاء منحنيات الكرونوأمبيرومتري (التيار مقابل التكرار/الزمن).

• التحقق التجريبي:

  • العينات: تم إنماء ثلاثة عينات عبر MOVPE:
    • العينة A: DBR قياسي بطبقات متساوية السمك (63 نانومتر مخصب / 51 نانومتر NID).
    • العينة B و C: DBRs ذات سماكات طبقات مخصبة متفاوتة (تزداد أو تنقص عبر التكدس) لاختبار متانة النموذج.
  • التآكل: تم إجراء التآكل في حمض الأوكساليك بجهود تتراوح من 6 فولت إلى 12 فولت دون وجود أخاديد ليتوغرافية.
  • التوصيف: تم تسجيل الكرونوأمبيرومتري أثناء التآكل. وتم تحليل البنى المجهرية بعد التآكل عبر SEM وFIB tomography (بالإشارة إلى بيانات سابقة لـ Thornley et. al).

3. المساهمات الرئيسية

1. تطوير نموذج شلال عشوائي: نجح المؤلفون في إنشاء أول محاكاة قادرة على إعادة إنتاج مورفولوجيا "الشلال" حيث تنتقل مسارات التآكل بين الخلول (TDs) عرضياً، متجاوزين بذلك "نموذج الكباب" البسيط.
2. الربط بين الجهد والاحتمالية: تربط الدراسة بشكل مباشر بين الجهد المطبق تجريبياً ومعاملات المحاكاة. حيث يقابل زيادة الجهد المطبق في التجارب زيادة في احتمالية تآكل الطبقات المخصبة ($P_{doped}$) في المحاكاة.
3. تكميم المورفولوجيا: تنجح المحاكاة في التنبؤ بنسبة هياكل "K5" (الخلول نشطة في جميع الطبقات الخمس) إلى هياكل "C" (الشلالات) كدالة لظروف التآكل، وهو ما يطابق الاتجاهات الملاحظة في تصوير FIB.
4. التوسع في الهياكل غير المنتظمة: تم التحقق من صحة النموذج باستخدام عينات ذات سماكات طبقات غير منتظمة (العينتان B و C)، مما يثبت قابليته للتطبيق بما يتجاوز DBRs الدورية القياسية.
5. رؤية حول تآكل الأخاديد: تم تكييف المحاكاة لنمذجة التآكل من الأخاديد المحددة ليتوغرافياً، مما يشير إلى أنه حتى في النهج القائم على الأخاديد، فإن الشلالات الناتجة عن الخلل ستسيطر في النهاية على مورفولوجيا المسام بعيداً عن جدار الخندق.

4. النتائج الرئيسية

• مطابقة الكرونوأمبيرومتري:
  • تعيد المحاكاة إنتاج ذروات التيار المتذبذبة المميزة التي تُرى في التجارب، حيث تقابل كل ذروة تقريباً تآكل طبقة مخصبة.
  • تأثير الجهد: تؤدي الجهود المنخفضة (low $P_{doped}$) إلى تذبذبات مخمدة وتيارات منخفضة، بينما تنتج الجهود العالية (high $P_{doped}$) ذروات واضحة وعالية السعة.
  • التباين: لا تلتقط المحاكاة طفرة التيار الأولية العالية المشاهدة في التجارب (والتي تُعزى إلى التآكل الجانبي في العينات الحقيقية، والذي تستبعده المحاكاة).
• التطور المورفولوجي:
  • الجهد المنخفض (Low $P_{doped}$): تُظهر المحاكاة مزيجاً من قطع "الكباب" القصيرة وشلالات متكررة. جبهة التآكل غير منتظمة، مما يؤدي إلى قمم تيار غير منتظمة.
  • الجهد العالي (High $P_{doped}$): تتنبأ المحاكاة بنسبة أعلى من هياكل "K5" (كباب كامل) وهياكل "C" (شلالات) أقل. وهذا يطابق بيانات FIB التي تُظهر أن الجهود الأعلى تفضل التآكل الرأسي الأكثر استمرارية.
  • طبقة تلو الأخرى مقابل الشلال: عند الاحتمالات العالية، يبدو التآكل أكثر شبهاً بـ "طبقة تلو الأخرى". أما عند الاحتمالات المنخفضة، فإن التآكل العرضي يتنافس مع الانتشار الرأسي، مما يجعل جبهة التآكل غير منتظمة ويؤدي إلى تلاشي التذبذبات.
• عينات السماكة المتغيرة:
  • تنبأت المحاكاة بنجاح باتجاهات الكرونوأمبيرومتري للعينتين B و C.
  • حددت بشكل صحيح أن العينات ذات الطبقات الأكثر سمكاً بالقرب من السطح (العينة B) تُظهر انخفاضاً في سعة الذروة بمرور الوقت، بينما تُظهر العينات ذات الطبقات السطحية الأرق (العينة C) ذروات أكثر استقراراً.
  • محدودية: تبالغ المحاكاة قليلاً في تقدير تآكل الطبقات العميقة مقارنة بالتجارب، مما يشير إلى أن الوصول للمادة المآكلة للطبقات العميقة في الواقع أصعب مما يفترضه النموذج.
• مورفولوجيا المسام:
  • تعيد المحاكاة إنتاج هياكل المسام "الدقيقة والمتعرجة" عند الجهود المنخفضة.
  • محدودية: تفشل المحاكة في إعادة إنتاج "توسع المسام المخروطي" الملاحظ في طبقات NID عند الجهود العالية (حيث تقوم المادة المآكلة بتوسيع طبقة NID حول قلب الخلل)، حيث يمنع النموذج بصرامة تآكل مناطق NID بعيداً عن قلب الخلل.

5. الأهمية

• الرؤية الميكانيكية: توفر هذه الدراسة دليلاً قوياً على أن التآكل الكهروكيميائي لـ GaN هو عملية عشوائية تحكمها احتمالات محلية بدلاً من كونها عملية حتمية ومنتظمة. وهي تثبت صحة فرضية "الشلال" كآلية مهيمنة في التآكل بوساطة الخلل (TD).
• تحسين العمليات: من خلال ربط الجهد المطبق باحتمالية التآكل ($P_{doped}$)، تقدم المحاكاة أداة للمهندسين للتنبؤ والتحكم في المسامية والانعكاسية لـ DBRs دون الحاجة إلى إجراء تجارب تآكل مكثفة قائمة على التجربة والخطأ.
• القابلية للتطبيق الواسع: تشير قدرة النموذج على التعامل مع سماكات الطبقات غير المنتظمة وهياكل التآكل البديلة (الأخاديد) إلى إمكانية استخدامه كأداة متعددة الاستخدامات لتصميم أجهزة فوتونية نتريدية معقدة.
• الفيزياء المستقبلية: يسلط العمل الضوء على أنه بينما يمكن نمذجة السلوك العياني (macroscopic) عشوائياً، فإن الفيزياء المجهرية (مثل نفق Zener، أو النفق بمساعدة المصائد عند نوى الخلل في طبقات NID) تتطلب مزيداً من البحث لفهم التباينات في الاحتمالية مع الجهد بشكل كامل.