شبیه‌سازی مورفولوژی‌های Bulk-Heterojunction با برش‌های الکتریکی دوبعدی

1. مقدمه

Bulk-heterojunctionها (BHJها) ریزساختار اصلی بسیاری از دستگاه‌های الکترونیکی آلی هستند—مشهورترین آن‌ها سلول‌های خورشیدی آلی‌اند، اما همچنین آشکارسازهای نوری، حسگرها، و دیگر سامانه‌های مختلط یونی–الکترونیکی که در آن‌ها مسیرهای تراوا و گلوگاه‌های محلی اهمیت دارند. در یک BHJ، دو فاز در مقیاس 10–100 nm در یکدیگر نفوذ می‌کنند و یک فصل‌مشترک داخلی بزرگ و یک چشم‌انداز انتقال ایجاد می‌کنند که اساساً یک‌بعدی نیست. معاوضه آن، پیچیدگی است: وقتی ریزساختار را resolve می‌کنید، ناگزیر محدودیت‌های پراکلِیشن، بن‌بست‌ها، تنگ‌شدگی‌ها، و تغییرات فضایی شدید در چگالی حامل و بازترکیب را نیز به ارث می‌برید. برای اکثر برازش‌های روزمره دستگاه، از اینجا شروع نمی‌کنید. معمولاً در سطح دستگاه 1D کار می‌کنید—سریع، پایدار، و اغلب کاملاً مناسب—و در صورت لزوم فیزیک‌های اضافی مانند تله‌گذاری و بازترکیب را نیز وارد می‌کنید. برای مثال، آموزش شروع سریع OPV یک گردش‌کار استاندارد 1D (شامل فیزیک تله در صورت نیاز) را نشان می‌دهد که ابزار مناسبی برای اکثر مطالعات روتین و استخراج پارامتر است.

در سوی دیگر نردبان انتزاع، اگر هدف شما پرسش‌های در سطح سامانه مانند مقیاس‌پذیری، تلفات سری/شنت، یا رفتار ماژول/پنل باشد، معمولاً به‌جای resolve کردن مورفولوژی در مقیاس نانو به یک توصیف سطح بزرگ / آگاه از مدار منتقل می‌شوید. آموزش PM6:Y6 سطح بزرگ نمونه‌ای از آن رویکرد مدل‌سازی در سطح بالاتر است.

این آموزش عمداً در جهت مخالف حرکت می‌کند: ریزساختار را (از طریق یک مورفولوژی BHJ صریح) حفظ می‌کند اما با حل مسئله الکتریکی روی برش‌های 2D از نظر محاسباتی عملی باقی می‌ماند. نکته این نیست که هر دستگاهی باید به این شکل برازش شود؛ بلکه داشتن یک گردش‌کار resolveکننده ریزساختار از نظر علمی و فناورانه ارزشمند است. این کار به شما اجازه می‌دهد فرض‌های محیط مؤثر را تحت فشار قرار دهید، ببینید مدل‌های «ساده» کجا از کار می‌افتند، و برای محدودیت‌های ناشی از مورفولوژی بر کارایی—به‌ویژه وقتی پراکلِیشن و تله‌گذاری/بازترکیب موضعی فضایی غالب باشند—شهود بسازید.

2. یک شبیه‌سازی مورفولوژی ایجاد کنید و BHJ را بررسی کنید

از پنجره اصلی، روی New simulation کلیک کنید. این کار کتابخانه نوع دستگاه را که در ?? نشان داده شده است باز می‌کند. روی Morphology دوبار کلیک کنید، سپس پروژه را در جایی روی دیسک محلی خود ذخیره کنید. ذخیره به‌صورت محلی اهمیت دارد زیرا شبیه‌سازی‌های resolveشده از نظر مورفولوژی می‌توانند فایل‌های میانی زیادی تولید کنند (meshها، snapshotها، cacheها)، و ذخیره‌سازی محلی معمولاً بسیار سریع‌تر و قابل‌اعتمادتر از پوشه‌های شبکه‌ای یا همگام‌سازی‌شده با ابر است.

پس از ایجاد پروژه، پنجره اصلی دستگاه باز می‌شود. باید یک شیء قرمز مرکزی را ببینید که نماینده یکی از فازهای bulk-heterojunction (BHJ) است و درون ناحیه محصورکننده لایه فعال قرار گرفته است، همان‌طور که در ?? نشان داده شده است. در حالی که نمای کامل 3D اطلاع‌رسان است، مورفولوژی می‌تواند از نظر بصری متراکم باشد. برای واضح‌تر کردن ساختار داخلی، روی شیء BHJ راست‌کلیک کنید و View → Show cut through Y را انتخاب کنید.

فعال کردن نمای برش‌خورده، رندر نشان‌داده‌شده در ?? را تولید می‌کند که یک برش 2D از میان مورفولوژی 3D را آشکار می‌کند. این کار دیدن اتصال‌پذیری فاز، فصل‌مشترک‌ها، و مسیرهای پراکلِیشن درون BHJ را بسیار آسان‌تر می‌کند.

یک مدل ذهنی مفید در این مرحله چنین است: ناحیه مستطیلی محصورکننده حجم لایه فعال را تعریف می‌کند، در حالی که شیء قرمزِ نهفته یک فاز صریح BHJ (برای مثال، پلیمر) را تعریف می‌کند. فاز مکمل به‌طور ضمنی به‌عنوان باقیمانده حجم لایه فعال که توسط شیء نهفته اشغال نشده است تعریف می‌شود. این جداسازی هنگام بررسی ویژگی‌های شیء در ویرایشگر شیء در بخش‌های بعدی صریح می‌شود.

3. ساختار دستگاه: لایه فعال و مورفولوژی BHJ نهفته

تب Device structure را باز کنید و روی Layer editor کلیک کنید. دستگاه از سه لایه تشکیل شده است: یک تماس بالا، یک لایه فعال منفرد، و یک تماس پایین (شکل ??). این پشته عمداً ساده انتخاب شده است تا درک این‌که مورفولوژی چگونه وارد مسئله الکتریکی می‌شود آسان‌تر شود، بدون آن‌که فصل‌مشترک‌های اضافی یا لایه‌های انتقال معرفی شوند. نکته کلیدی این است که خود لایه فعال یک حجم مستطیلی ساده است که در ویرایشگر لایه تعریف می‌شود. در این مثال این لایه به ماده مولکولی Y6، یک پذیرنده، اختصاص داده شده است و—به‌خودی‌خود—فقط یک جعبه است. تمام پیچیدگی ساختاری مرتبط با bulk-heterojunction به‌صورت جداگانه، با نهفتن یک شیء دوم درون این حجم لایه فعال، وارد می‌شود. آن شیء نهفته یک شبکه CAD پیچیده است که از پایگاه داده Shape وارد شده و نمایانگر فاز پلیمر PM6 است. شبکه پلیمر یک حجم سه‌بعدی بسته را درون لایه فعال تعریف می‌کند. فاز مکمل (ناحیه غنی از پذیرنده) سپس به‌طور ضمنی به‌صورت بخشی از لایه فعال که توسط شبکه پلیمر اشغال نشده است تعریف می‌شود. به بیان دیگر، BHJ با حذف حجمی نمایش داده می‌شود نه با تعریف صریح دو شبکه درهم‌قفل‌شده.

این جداسازی هنگام بررسی اشیا در نمای 3D روشن می‌شود. دو شیء مرتبط وجود دارد: جعبه لایه فعال محصورکننده (که معمولاً به‌صورت یک wireframe آبی نشان داده می‌شود)، و مورفولوژی پلیمر نهفته (که به رنگ قرمز نشان داده می‌شود). برای بررسی ویژگی‌های آن‌ها، روی ناحیه آبی محصورکننده راست‌کلیک کنید و Edit object را انتخاب کنید، همان‌طور که در ?? نشان داده شده است.

در ویرایشگر شیء، تمایز صریح است. برای شیء محصورکننده Active، Object shape یک جعبه ساده است (با غیرفعال بودن meshهای پیچیده). برای شیء Polymer، Object shape به یک ورودی Shape DB اشاره می‌کند، که همان شبکه مثلث‌بندی‌شده تعریف‌کننده مورفولوژی PM6 است. تغییر این مرجع مورفولوژی پلیمر را عوض می‌کند در حالی که لایه فعال پیرامون آن بدون تغییر باقی می‌ماند.

این ساختار—جعبه ساده لایه فعال به‌علاوه حجم پلیمر نهفته—در نحوه مدیریت مورفولوژی‌های BHJ در OghmaNano مرکزی است. این کار هندسه را از نظر مفهومی ساده نگه می‌دارد، تضمین می‌کند فاز پلیمر یک ناحیه بسته و به‌خوبی تعریف‌شده تشکیل دهد، و به شما امکان می‌دهد بررسی کنید که چگونه پراکلِیشن، اتصال‌پذیری، و تنگ‌شدگی‌ها در شبکه پلیمر، وقتی حل‌گر الکتریکی اجرا شود، بر انتقال بار اثر می‌گذارند.

4. پارامترهای الکتریکی: تحرک‌ها و باندهای تله SRH

ویرایشگر Electrical parameters را از پنجره اصلی باز کنید. پنل‌های پارامتر جداگانه‌ای برای دو فاز لایه فعال فراهم شده‌اند: ناحیه محصورکننده Active (پذیرنده / ماتریس) و مورفولوژی نهفته Polymer که در ?? و ?? نشان داده شده‌اند. اینجا جایی است که دو فاز درهم‌نفوذکننده از نظر الکتریکی متمایز می‌شوند.

هر دو فاز از یک فرمالیسم پایه انتقال استفاده می‌کنند: drift–diffusion با آمار حامل Maxwell–Boltzmann و tensorهای تحرک متقارن. تفاوت کلیدی، کنتراست تحرک عمدی است. در فاز پلیمر، تحرک حفره بالا است در حالی که تحرک الکترون سرکوب می‌شود؛ در فاز Active مکمل، این وضعیت برعکس است. این کار انتقال انتخابی را اعمال می‌کند: حفره‌ها از طریق شبکه پلیمر پراکله می‌شوند، در حالی که الکترون‌ها از طریق ماتریس پراکله می‌شوند.

این عدم‌تقارن ساده برای ایجاد رفتار فیزیکی معنادار BHJ کافی است. حتی پیش از اجرای شبیه‌سازی، این موضوع تعیین می‌کند که در خروجی‌های فضایی چه چیزی را باید انتظار داشته باشید: مسیرهای پراکلِیشن پیوسته، بن‌بست‌هایی که در آن‌ها حامل‌ها از نظر هندسی به دام می‌افتند، و تجمع بار موضعی در تنگ‌شدگی‌ها یا نواحی با اتصال‌پذیری ضعیف مورفولوژی. علاوه بر انتقال حامل آزاد، تله‌گذاری Shockley–Read–Hall (SRH) غیرتعادلی در هر دو فاز فعال شده است. همان‌طور که در پنل‌های پارامتر نشان داده شده است، چگالی حالات با استفاده از توزیع‌های دُم نمایی برای الکترون‌ها و حفره‌ها نمایش داده می‌شود، با چگالی‌های تله و شیب‌های انرژی یکسان در هر دو فاز. حل‌گر پنج باند تله گسسته را در هر نقطه شبکه ردیابی می‌کند.

هر باند تله از طریق جمله‌های صریح گیراندازی و گسیل (ثابت‌های کوپلینگ آزاد-به-تله که در ویرایشگر دیده می‌شوند) به‌صورت دینامیکی حامل‌ها را با جمعیت‌های آزاد مبادله می‌کند. در نتیجه، سامانه‌ای که حل می‌شود نه‌تنها شامل چگالی‌های آزاد الکترون و حفره است، بلکه شامل چندین جمعیت حامل به‌دام‌افتاده نیز در هر مکان فضایی می‌شود. این کار سختی عددی و هزینه محاسباتی را به‌طور قابل‌توجهی افزایش می‌دهد، اما برای ثبت انتقال پاشنده، بازترکیب یاری‌شده توسط تله، و رفتار گذرای واقع‌گرایانه در نیمه‌رساناهای آلی بی‌نظم ضروری است. مقادیر عددی به‌کاررفته در اینجا به‌عنوان یک مثال پایدار و نمایشی در نظر گرفته شده‌اند، نه یک مدل ماده کالیبره‌شده. وقتی گردش‌کار درک شد، این پارامترها را می‌توان پالایش کرد، در برابر آزمایش اعتبارسنجی کرد، یا در صورت این‌که هدف جداسازی اثرات مورفولوژی و پراکلِیشن باشد و نه فیزیک تله‌گذاری، ساده‌تر کرد.

5. شبکه‌بندی

پیش از اجرای شبیه‌سازی، ارزش دارد شبکه الکتریکی را بررسی کنید تا مطمئن شوید معقول است. در پنجره اصلی، ریبون Electrical را باز کنید (شکل ??) و روی Electrical mesh کلیک کنید. این کار ویرایشگر شبکه نشان‌داده‌شده در ?? را باز می‌کند.

در این مثال، شبکه در صفحه xy برابر با 40 × 40 است، یعنی مسئله الکتریکی در دو بعد حل می‌شود. جهت z عمداً حذف شده است. شبیه‌سازی‌های الکتریکی کاملاً سه‌بعدی در OghmaNano ممکن هستند، اما به یک مورفولوژی 3D واقعی نیاز دارند و با افزایش قابل‌توجه هزینه و پیچیدگی محاسباتی همراه‌اند. ساخت و تفسیر چنین مدل‌هایی فراتر از دامنه این آموزش است.

در عوض، این محاسبه باید به‌عنوان یک مقطع 2D از میان دستگاه در نظر گرفته شود. این برش از مرکز مورفولوژی می‌گذرد و ریزساختار جانبی کامل BHJ را درون آن صفحه حفظ می‌کند. این رویکرد پراکلِیشن، گلوگاه‌ها، و اتصال‌پذیری فاز را ثبت می‌کند، در حالی که سریع، قابل‌تفسیر، و مناسب برای کاوش در مقیاس آموزش باقی می‌ماند.

6. اجرای شبیه‌سازی

برای شروع حل‌گر، روی دکمه آبی Run simulation در پنجره اصلی کلیک کنید. در حالی که اجرا در حال انجام است، ترمینال نقطه بایاس و چگالی جریان را در هر گام چاپ می‌کند (??). در این مثال بایاس اعمال‌شده به تماس بالا از 0.00 V تا تقریباً 0.54 V گام‌برداری می‌شود تا یک جاروب JV تولید شود. در بایاس پایین، چگالی جریان گزارش‌شده تحت روشنایی منفی است، و وقتی بایاس به رژیم مدار باز نزدیک می‌شود، جریان خالص از صفر عبور کرده و تغییر علامت می‌دهد، که با رفتار سلول خورشیدی سازگار است.

وقتی اجرا پایان یافت، تب Output را باز کنید (??) تا فایل‌های تولیدشده را بررسی کنید. خلاصه ترمینال زمان کل اجرا، تعداد معادلات حل‌شده، و آمار پایه I/O دیسک را گزارش می‌کند. در شبیه‌سازی‌های resolveشده از نظر مورفولوژی، پایش خروجی فایل مهم است: نوشتن تعداد زیادی snapshot می‌تواند حتی وقتی خود حل‌گر سریع است، زمان اجرا را غالب کند.

برای یک sanity check هندسی، روی device.csv دوبار کلیک کنید تا نمایش مثلث‌بندی‌شده دستگاهی را که حل‌گر درونی استفاده می‌کند مشاهده کنید (??). این کار تأیید می‌کند که مورفولوژی BHJ همان‌گونه که انتظار می‌رود وارد و گسسته‌سازی شده است.

برای مشاهده کارایی الکتریکی، jv_contact_0.csv یا jv_contact_1.csv را باز کنید تا منحنی JV را رسم کنید (??). برای این شبیه‌سازی 2D، باید از خروجی‌های JV مربوط به هر تماس استفاده شود؛ می‌توان jv.csv را نادیده گرفت.

7.0 بررسی snapshotهای شبیه‌سازی

پس از پایان اجرا، آموزنده‌ترین خروجی‌ها در پوشه snapshots/ یافت می‌شوند. فایل‌های snapshot به‌عنوان متغیرهای حالت داخلی سامانه drift–diffusion–SRH ذخیره می‌شوند، که به‌صورت تابعی از بایاس اعمال‌شده نوشته شده‌اند. آن‌ها حل resolveشده از نظر فضاییِ معادلات انتقال را در هر گام ولتاژ آشکار می‌کنند، نه خلاصه‌های پس‌پردازش‌شده. وقتی یک snapshot باز می‌شود، دو کنترل در دسترس قرار می‌گیرد: یک لغزنده بایاس که در جاروب JV حرکت می‌کند، و یک لغزنده موقعیت برش که برش 2D از میان مورفولوژی را انتخاب می‌کند. در سراسر این آموزش، مسئله الکتریکی روی یک برش 2D حل می‌شود، و تمام داده‌های snapshot متناظر با آن برش هستند.

شکل‌های ??– ?? snapshotهای نماینده 2D را در یک بایاس ثابت نشان می‌دهند. این کمیت‌ها باید در کنار یکدیگر تفسیر شوند. چشم‌انداز باند رسانش (LUMO)، Ec.csv (??)، محیط انرژی برای الکترون‌ها را تعریف می‌کند. در یک شبیه‌سازی resolveکننده مورفولوژی این چشم‌انداز به‌دلیل انرژی‌های وابسته به فاز و پاسخ الکترواستاتیکی به تجمع بار محلی، به‌صورت جانبی تغییر می‌کند. چگالی الکترون آزاد، Q_nfree.csv (??)، نواحی‌ای را برجسته می‌کند که در آن‌ها الکترون‌ها می‌توانند به‌طور مؤثر از طریق شبکه پذیرنده/ماتریس حرکت کنند. تغییرات فضایی شدید نشان‌دهنده مسیرهای پراکلِیشن، تنگ‌شدگی‌ها، و گلوگاه‌های ناشی از مورفولوژی هستند.

چگالی الکترون به‌دام‌افتاده، Q_ntrap.csv (??)، نشان می‌دهد حامل‌ها در کجا در حالت‌های تله تجمع می‌یابند. مقایسه چگالی‌های آزاد و به‌دام‌افتاده اغلب سریع‌ترین راه برای شناسایی نواحی استخراج ضعیف یا تله‌گذاری محلی شدید وابسته به مورفولوژی BHJ است. خروجی مرتبط با تحرک، mu_p_ft.csv (??)، یک تحرک مؤثر میانگین‌گیری‌شده روی حامل‌های آزاد و به‌دام‌افتاده را نشان می‌دهد. تغییرات فضایی هم کنتراست ذاتی فاز و هم تله‌گذاری محلی را بازتاب می‌دهد، و اغلب همبستگی شدیدی با اتصال‌پذیری شبکه زیرین دارد.

در حالی که نماهای 2D اغلب روشن‌ترین راه تحلیل تغییرات فضایی هستند، گاهی تفسیر همین داده‌ها در فضای واقعی آسان‌تر است. با فعال کردن 3D mode در پنجره snapshot، برش 2D دوباره به هندسه دستگاه فرافکنی می‌شود. شکل‌های ??– ?? نمونه‌هایی از این فرافکنی را نشان می‌دهند. تحرک‌های میانگین‌گیری‌شده الکترون و حفره (mu_n_ft.csv، mu_p_ft.csv)، شبه-فرمی آزاد الکترون (Fn.csv)، و چگالی جریان الکترون در جهت انتقال (Jn.csv) برای بصری‌سازی مسیرهای پراکلِیشن و جریان از میان BHJ به‌ویژه مفید هستند.

این خروجی‌های snapshot ارزش علمی اصلی شبیه‌سازی‌های resolveکننده مورفولوژی را فراهم می‌کنند: آن‌ها امکان می‌دهند کمیت‌های قابل‌مشاهده در مقیاس دستگاه به‌طور مستقیم به انرژی‌های موضعی، تله‌گذاری، و اتصال‌پذیری درون bulk heterojunction مرتبط شوند.

8. کارهایی برای امتحان کردن

سریع‌ترین راه برای ساخت شهود، اغتشاش دادن مدل و دیدن آن‌چه تغییر می‌کند است. کارهای زیر طوری طراحی شده‌اند که سریع باشند: شبکه مورفولوژی را عوض کنید، دوباره اجرا کنید، و منحنی JV را مقایسه کنید (و، اگر وقت داشتید، چند snapshot را در همان بایاس نیز بررسی کنید).

Task 1 — تعویض شبکه مورفولوژی

در نمای اصلی دستگاه، روی شیء قرمز BHJ (شبکه پلیمر نهفته) راست‌کلیک کنید و Mesh editor را انتخاب کنید. در ویرایشگر، گزینه فعلی باید Shape from database باشد. شکل انتخاب‌شده را به morphology/1، morphology/2، و morphology/3 تغییر دهید (شکل‌های ??، ??، ??). پس از هر تغییر شبیه‌سازی را دوباره اجرا کنید و منحنی JV را مقایسه کنید.

چیزی که باید به آن نگاه کنید این است که اتصال‌پذیری و گلوگاه‌ها چگونه استخراج و بازترکیب را تغییر می‌دهند: بعضی مورفولوژی‌ها مسیرهای پراکلِیشن تمیزتری تولید می‌کنند (جریان بیشتر)، در حالی که بعضی دیگر حامل‌ها را مجبور می‌کنند از گردن‌های باریک و بن‌بست‌ها عبور کنند (تجمع بیشتر، تله‌گذاری بیشتر، و معمولاً عملکرد JV ضعیف‌تر).

به‌عنوان یک آزمایش عمداً «نادرست اما آموزنده»، در پایگاه داده شکلی به نام teapot نیز خواهید یافت. شبکه BHJ را به teapot تغییر دهید، دوباره اجرا کنید، و برای دستگاهی با inclusion به شکل قوری یک منحنی JV تولید کنید. این یک BHJ از نظر فیزیکی معقول نیست، اما نشان می‌دهد که این گردش‌کار واقعاً هندسه‌محور است: هر شبکه CAD بسته‌ای را می‌توان وارد کرد و از نظر الکتریکی تحلیل کرد.

Task 2 — تغییر شدت تله‌گذاری (باندهای تله SRH)

در Electrical parameters، تنظیمات تله‌های SRH غیرتعادلی را تغییر دهید و دوباره اجرا کنید. هر بار یک تغییر را امتحان کنید:

• چگالی تله را 1-2 مرتبه بزرگی افزایش/کاهش دهید و Q_nfree را با Q_ntrap مقایسه کنید.

هدف این کار دیدن این است که انتقال پاشنده/محدودشده به‌وسیله تله چگونه با مورفولوژی رقابت می‌کند: تله‌گذاری معمولاً نواحی «بد» (بن‌بست‌ها، تنگ‌شدگی‌ها) را تشدید می‌کند زیرا بار نمی‌تواند به همان سرعت تخلیه شود.

Task 3 — تغییر ثابت‌های نرخ بازترکیب

همچنان در Electrical parameters، بازترکیب را تغییر دهید و دوباره اجرا کنید. دو اغتشاش ساده:

• ثابت نرخ بازترکیب آزاد-به-آزاد (n_free به p_free) را 1–2 مرتبه بزرگی بالا/پایین ببرید. ببینید منحنی JV و توزیع فضایی بار چگونه پاسخ می‌دهند.
• اگر پروژه شما پارامترهای گیراندازی/کوپلینگ SRH (ثابت‌های کوپلینگ آزاد↔تله) را نشان می‌دهد، آن‌ها را به‌طور همزمان بالا/پایین مقیاس دهید و توازن Q_nfree و Q_ntrap را در همان بایاس مقایسه کنید.

یک عادت مفید برای Taskهای 2–3 این است که مورفولوژی را ثابت نگه دارید، یک پیچ تنظیم را تغییر دهید، و سپس هر بار همان دو چیز را بررسی کنید: (i) منحنی JV چگونه جابه‌جا می‌شود، و (ii) آیا مسیرهای جریان (برای مثال Jn و میدان‌های شبه-فرمی) هموارتر می‌شوند یا بیشتر توسط ریزساختار «شکسته» می‌شوند.