نیاز به حالت‌های تله در مدل‌های دستگاه آلی

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

1. ساختار فیزیکی و انرژیایی مواد بی‌نظم

نیمه‌رساناهای معدنی سنتی مانند Si بلوری یا GaAs هم بسیار منظم‌اند و هم فوق‌العاده خالص — و اغلب به خلوص “نه-نه” (99.9999999%) می‌رسند. در مقابل، نیمه‌رساناهای آلی به‌ندرت از 99.9% خلوص فراتر می‌روند، که آن‌ها را حدود یک میلیون برابر پُرعیب‌تر از همتایان معدنی‌شان می‌سازد. از نظر ساختاری نیز تفاوت به همان اندازه چشمگیر است: نیمه‌رساناهای معدنی شبکه‌های بلوری منظم تشکیل می‌دهند، مانند تیله‌هایی که مرتب روی یک صفحه سولیتیر چیده شده‌اند (??). خود سیلیکون ساختار مکعبی الماسی اختیار می‌کند، یعنی یک شبکه تقریباً کامل (??). در مقابل، مواد آلی سامانه‌های مولکولی “شل” هستند، با پلیمرهای درهم‌تنیده‌ای که شبیه یک بشقاب اسپاگتی بولونز‌اند — که در آن اسپاگتی‌ها نمایانگر زنجیره‌های پلیمری و سس نمایانگر مولکول‌های کوچک است (??). شبیه‌سازی‌های آمیخته‌های آلی این تصویر را تأیید می‌کنند و نشان می‌دهند بسته‌بندی پلیمرها شبیه اسپاگتی است و مشتقات فولرن در سراسر آن پراکنده‌اند (??).

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

نتیجه روشن است: در حالی که در نیمه‌رساناهای منظم غالباً می‌توان حالت‌های تله را نادیده گرفت، آن‌ها بر فیزیک سامانه‌های بی‌نظم حاکم‌اند. بنابراین هر مدل دستگاه واقع‌بینانه باید یک توصیف دقیق از توزیع‌های تله و سینتیک Shockley–Read–Hall را در بر بگیرد. OghmaNano دقیقاً همین کار را انجام می‌دهد و شبیه‌سازی سامانه‌های بی‌نظمی مانند سلول‌های خورشیدی آلی، OFETها، و پروسکایت‌ها را با دقت فیزیکی ممکن می‌سازد.

برای دیدن اینکه چرا این موضوع مهم است، اکنون به رابطه بین چگالی حامل و تراز فرمی می‌پردازیم تا چند مثال مشخص را بررسی کنیم.

2. چرا حالت‌های تله برای مدل‌سازی دستگاه مهم‌اند (بدون ریاضیات)

یکی از ویژگی‌های مرکزی نیمه‌رساناهای آلی و دیگر نیمه‌رساناهای بی‌نظم این است که چگالی حامل تابعی قوی از هر دو، ولتاژ اعمال‌شده و شدت روشنایی است. با افزایش بایاس یا شدت نور، بار بیشتری به دستگاه تزریق می‌شود یا به‌صورت نوری تولید می‌گردد. چون این مواد تعداد زیادی حالت تله درون گاف نواری دارند، حامل‌ها ابتدا این تله‌ها را پُر می‌کنند و سپس حالت‌های گسترده را اشغال می‌کنند. این پُر شدن تله‌ها به این معناست که حتی تغییرات کوچک در ولتاژ می‌تواند تغییرات بزرگی در چگالی حامل آزاد ایجاد کند.

این موضوع مهم است زیرا بازترکیب در دستگاه مستقیماً به چگالی حامل‌ها وابسته است. شکل کلی نرخ بازترکیب برابر است با

\[ R = k_r \, n(V)\, p(V), \]

که در آن \(k_r\) ثابت بازترکیب است، و \(n(V)\) و \(p(V)\) چگالی‌های الکترون و حفره وابسته به ولتاژ هستند. اگر فرم تابعی \(n(V)\) (و \(p(V)\)) نادرست باشد زیرا حالت‌های تله نادیده گرفته شده‌اند، آنگاه نرخ بازترکیب نیز نادرست خواهد بود. این موضوع مستقیماً به پیش‌بینی‌های نادرست برای ولتاژ مدار باز (\(V_{OC}\)) و دیگر ویژگی‌های کلیدی دستگاه منجر می‌شود.

تحرک نیز به‌طور مشابه تحت تأثیر قرار می‌گیرد. در سامانه‌های بی‌نظم، تحرک حاملِ مؤثر به تعادل بین حامل‌های آزاد و حامل‌های به‌دام‌افتاده وابسته است. یک بیان ساده چنین است

\[ \mu_e(n) = \frac{\mu_e^0 \, n_{\text{free}}}{n_{\text{free}} + n_{\text{trap}}}, \]

که در آن \(\mu_e^0\) تحرک ذاتی الکترون، \(n_{\text{free}}\) چگالی حامل‌های متحرک، و \(n_{\text{trap}}\) چگالی حامل‌های به‌دام‌افتاده است. اگر رابطه چگالی–ولتاژ نادرست باشد، آنگاه وابستگی پیش‌بینی‌شده تحرک–ولتاژ نیز نادرست خواهد بود. در کنار هم، این خطاها به این معنا هستند که منحنی J–V شبیه‌سازی‌شده با آزمایش مطابقت نخواهد داشت، حتی اگر پارامترهای بازترکیب یا تحرک در غیر این صورت معقول باشند.

3. چرا حالت‌های تله برای مدل‌سازی دستگاه مهم‌اند (با ریاضیات)

برای توصیف درست چگالی حامل‌ها، باید چگالی حالت‌های (DoS) زیربنایی را در نظر گرفت. شکل ?? طرحی از DoS را متناظر با ساختارهای نواری منظم و بی‌نظم که پیش‌تر در ?? و ?? نشان داده شدند، ارائه می‌دهد. در یک نیمه‌رسانای منظم، DoS یک لبه نواری تیز دارد (نوار سهموی)، و جمعیت Fermi–Dirac بالای لبه نوار رسانش قرار می‌گیرد. اما در یک نیمه‌رسانای بی‌نظم، DoS یک دنباله از حالت‌های تله موضعی را نشان می‌دهد که تا عمق زیادی در گاف نواری امتداد می‌یابد (که معمولاً با دنباله‌های نمایی یا گاوسی مدل می‌شود). در نتیجه حامل‌ها در این دو حالت، توزیع‌هایی اساساً متفاوت را اشغال می‌کنند.

به‌صورت صوری، چگالی کل الکترون از انتگرال‌گیری DoS وزن‌دهی‌شده با اشغال Fermi–Dirac به دست می‌آید:

\[ n(E_f,T) = \int_{E_{\min}}^{\infty} \rho(E)\, f(E,E_f,T)\, dE, \]

که در آن \(E_f\) تراز شبه‌فرمی، \(\rho(E)\) DoS، و \(f(E,E_f,T)\) توزیع Fermi–Dirac است. در مواد منظم، \(\rho(E)\) تیز است، بنابراین فقط حالت‌های بالای لبه نوار رسانش مشارکت می‌کنند. در مواد بی‌نظم اما \(\rho(E)\) شامل دنباله‌های تله گسترده است، بنابراین حامل‌ها می‌توانند در نزدیکی تراز فرمی ذخیره شوند و چگالی‌های باری یک تا دو مرتبه بزرگی بزرگ‌تر از یک بلور منظم در همان بایاس ایجاد کنند. این موضوع مستقیماً در آزمایش‌های استخراج بار مشاهده می‌شود.

نتیجه‌گیری روشن است: درست‌گرفتن وابستگی چگالی حامل اختیاری نیست. بدون یک مدل حالت تله، هم بازترکیب و هم تحرک بر حسب ولتاژ نادرست خواهند بود و منحنی‌های J–V واقع‌بینانه قابل بازتولید نخواهند بود. OghmaNano این حالت‌های تله را به‌صورت صریح در بر می‌گیرد و مدل‌سازی دقیق دستگاه‌های بی‌نظم مانند آمیخته‌های PM6:Y6 و P3HT:PCBM، و همچنین نیمه‌رساناهای منظم‌تر را هنگام غیرفعال بودن تله‌ها ممکن می‌سازد.

4. چرا نباید از بازترکیب Langevin در مدل‌های دستگاه استفاده کنید

نرخ کلاسیک بازترکیب Langevin به صورت زیر تعریف می‌شود

\[ R_{\text{free}} = q \, k_r \, \frac{\mu_e + \mu_h}{2 \, \epsilon_0 \epsilon_r} \, n p , \]

که در آن \(R_{\text{free}}\) نرخ بازترکیب، \(k_r\) ضریب تجربی کاهش Langevin، \(\mu_e\) و \(\mu_h\) تحرک‌های الکترون و حفره، \(n\) و \(p\) چگالی حامل‌ها، و \(\epsilon_0 \epsilon_r\) گذردهی دی‌الکتریک هستند. در نگاه اول این از نظر فیزیکی معقول به نظر می‌رسد: فرض می‌شود بازترکیب هر زمان رخ می‌دهد که یک الکترون و یک حفره تحت حرکت براونی به اندازه کافی به هم نزدیک شوند تا میدان کولنی یکدیگر را حس کنند. این تصویر برای حامل‌های کاملاً آزاد در مایعات ساده یا رساناهای یونی مناسب است. اما در فتوولتائیک آلی (OPV) و دیگر نیمه‌رساناهای بی‌نظم، فرض‌های پشت مدل Langevin از هم می‌پاشند.

چرا این مدل شکست می‌خورد؟ مطالعات تجربی — به‌ویژه از اوایل دهه 2010 — به‌سرعت نشان دادند که بازترکیب Langevin نمی‌تواند هم منحنی‌های J–V در تاریکی و هم در روشنایی را به‌صورت خودسازگار بازتولید کند. این مدل به‌طور نظام‌مند نرخ‌های بازترکیب را بیش‌برآورد می‌کرد، اغلب به اندازه چند مرتبه بزرگی. برای ممکن شدن برازش‌ها، پژوهشگران یک “ضریب کاهش Langevin” \(k_r\) وارد کردند که گاهی تا 10⁻³ کوچک بود. هرچند این کار از نظر عملی مفید بود، اما در واقع اعترافی بود به اینکه خودِ مکانیسم در این سامانه‌ها معتبر نیست.

مشکلات هنگامی روشن می‌شوند که معادله را دقیق‌تر بررسی کنیم:

• وابستگی نادرست به چگالی حامل. نرخ به‌طور سخت‌گیرانه به صورت \(np\)، یعنی مرتبه دوم، مقیاس می‌شود. در آزمایش‌ها، بازترکیب در OPVها اغلب از قوانین توانی غیرصحیح مانند \((np)^{1.5}\) پیروی می‌کند، که بازتاب پویایی محدودشده توسط تله و توزیع حالت‌های موضعی است. فرم Langevin به‌سادگی نمی‌تواند این را ثبت کند.
• برخورد نادرست با تحرک. معادله فرض می‌کند تحرک‌های \(\mu_e\) و \(\mu_h\) ثابت‌اند. در واقع، تحرک در نیمه‌رساناهای بی‌نظم شدیداً به چگالی حامل وابسته است (بحث بالا را ببینید). وارد کردن وابستگی نادرست تحرک، نرخ بازترکیب پیش‌بینی‌شده را بیشتر مخدوش می‌کند.
• تله‌ها و بی‌نظمی را نادیده می‌گیرد. مکانیسم Langevin با همه حامل‌ها به‌عنوان حامل آزاد برخورد می‌کند. در OPVها و مواد مشابه، بیشتر حامل‌ها زمانی را در حالت‌های موضعی به‌دام‌افتاده می‌گذرانند و فقط به‌صورت متناوب در رسانش مشارکت می‌کنند. بنابراین بازترکیب به چشم‌انداز تله وابسته است، نه فقط به حرکت حامل آزاد.

در کنار هم، این مسائل به این معنا هستند که بازترکیب Langevin در بهترین حالت یک تقریب خام و در بدترین حالت یک تقریب گمراه‌کننده است. حتی با یک ضریب کاهش برازش‌شده \(k_r\)، این مدل در ثبت فیزیک درستِ بازترکیب کمک‌گرفته از تله و وابستگی تحرک به ولتاژ ناکام می‌ماند. بنابراین استفاده از بازترکیب Langevin در مدل‌های دستگاه مانند جا دادن یک میخ مربعی در سوراخ گرد است: شاید به شما یک عدد بدهد، اما نتیجه‌ای فیزیکی معنادار نخواهد داد.

5. چگونه می‌توان بازترکیب Langevin را در مدل‌های دستگاه “کارا” کرد

مشکلات کلیدی بازترکیب کلاسیک Langevin وابستگی نادرست آن به چگالی حامل و نیاز به یک ضریب کاهش دلخواه است. یک راهی که پژوهشگران برای “کارا” کردن بازترکیب Langevin آزموده‌اند این است که خودِ تحرک‌ها را وابسته به چگالی حامل وارد کنند، به صورت:

\[ R_{\text{free}} = q \, k_r \, \frac{\alpha \mu_e(n) + \beta \mu_h(n)}{2 \, \epsilon_0 \epsilon_r} \, n_{\text{tot}} p_{\text{tot}} , \]

در اینجا، یک لبه تحرک تعریف می‌شود: حامل‌های بالای لبه تحرک در رسانش مشارکت می‌کنند، در حالی که حامل‌های پایین‌تر به‌عنوان حامل‌های به‌دام‌افتاده در نظر گرفته می‌شوند. در این صورت، تحرک‌های میانگین را می‌توان به شکل زیر نوشت

\[ \mu_e(n) = \frac{\mu_e^0 \, n_{\text{free}}}{n_{\text{free}} + n_{\text{trap}}}, \qquad \mu_h(p) = \frac{\mu_h^0 \, p_{\text{free}}}{p_{\text{free}} + p_{\text{trap}}}. \]

اگر چگالی حامل‌های آزاد بسیار کمتر از چگالی حامل‌های به‌دام‌افتاده باشد، این موضوع به یک نرخ بازترکیب مؤثر به صورت زیر منجر می‌شود

\[ R(n,p) = q \, k_r \, \frac{\alpha \mu_e^0 \, n_{\text{free}} p_{\text{trap}} + \beta \mu_h^0 \, p_{\text{free}} n_{\text{trap}}} {2 \, \epsilon_0 \epsilon_r}. \]

به این ترتیب، بازترکیب Langevin عملاً بر حسب برهم‌کنش‌های بین حامل‌های آزاد و به‌دام‌افتاده (\(n_{\text{free}}p_{\text{trap}}\) و \(p_{\text{free}}n_{\text{trap}}\)) بازتفسیر می‌شود. این، در اصل، معادل تصویر Shockley–Read–Hall (SRH) از بازترکیب است: بازترکیب حامل‌های آزاد با حامل‌های به‌دام‌افتاده.

هرچند این رویکرد در حالت ماندگار به‌طور معقولی عمل می‌کند، اما بر یک فرض قوی تکیه دارد: اینکه همه حامل‌ها در یک مکان معین یک تراز شبه‌فرمی یکتا دارند، یعنی آن‌ها در تعادل موضعی با سرعت گرمایی‌سازی بی‌نهایت هستند. این فرض ممکن است در شرایط حالت ماندگار، زمانی که حامل‌ها برای رسیدن به تعادل فرصت دارند، قابل قبول باشد، اما در حوزه زمانی از بین می‌رود. در واقع، توزیع چگال حالت‌های تله در نیمه‌رساناهای آلی این احتمال را کم می‌کند که حامل‌ها بتوانند مانند یک گاز یگانه و متعادل رفتار کنند. در مقابل، صورت‌بندی SRH از این فرض اجتناب می‌کند، و بنابراین توصیفی از نظر فیزیکی معتبرتر از بازترکیب و تله‌گذاری در مواد بی‌نظم ارائه می‌دهد.