اندازهگیریهای CELIV روی سلولهای خورشیدی پروسکایتی: بخش a
1. مرور کلی
در این آموزش با معرفی مفاهیم پایهٔ تکنیک استخراج بار با ولتاژ افزایشی خطی (CELIV) آغاز میکنیم و توضیح میدهیم که چگونه از آن برای بررسی تحرک حامل و چگالی در نیمهرساناهای لایهنازک استفاده میشود. پس از ارائهٔ پسزمینهٔ نظری، بهسرعت به مراحل عملی خواهیم پرداخت و نشان میدهیم چگونه یک شبیهسازی CELIV را در OghmaNano پیکربندی و اجرا کنید. این ترکیب نظریه و شبیهسازی عملی هم یک درک روشن از روش فراهم میکند و هم ابزارهای لازم برای اعمال آن بر دستگاههای خودتان را در اختیار شما میگذارد.
2. پسزمینهٔ CELIV
تکنیک استخراج بار با ولتاژ افزایشی خطی (CELIV) یک روش آزمایشی است که برای مطالعهٔ تحرک و چگالی حامل بار در سلولهای خورشیدی آلی و پروسکایتی استفاده میشود. در یک آرایش معمول، سلول خورشیدی به یک مولد سیگنال متصل میشود که یک بایاس معکوس با افزایش خطی اعمال میکند، در حالی که گذرای جریان حاصل روی یک اسیلوسکوپ ثبت میشود. روشنسازی پیوستهٔ اختیاری یا یک پالس لیزر کوتاه میتواند برای تولید حاملهای بار اضافه شود. شماتیکی از این پیکربندی در ?? نشان داده شده است.
وقتی رمپ ولتاژ اعمال میشود، جریان اندازهگیریشده از دو مؤلفهٔ متمایز تشکیل میشود: یک خط مبنای مستطیلی از ظرفیت خازنی هندسی دستگاه و یک قلهٔ اضافی ناشی از استخراج حاملهای بار متحرک در نیمهرسانا. جملهٔ ظرفیت خازنی هندسی همیشه وجود دارد، در حالی که قلهٔ مربوط به حامل، اطلاعات مستقیمی دربارهٔ تحرک و چگالی حامل فراهم میکند. نمونهای از این گذراهای ولتاژ و جریان در ?? نشان داده شده است.
تحرک حامل را میتوان مستقیماً از گذرای CELIV با اندازهگیری زمانی که جریان استخراج به بیشینهٔ خود میرسد تعیین کرد. در سادهترین حالتِ توزیع یکنواخت حامل و بازترکیب ناچیز، تحرک با معادلهٔ CELIV به صورت زیر داده میشود: \[ \mu = \frac{2 d^{2}}{3 A t_{\text{max}}^{2}} \] که در آن \(d\) ضخامت لایهٔ فعال، \(A\) نرخ رمپ ولتاژ است (\(A = \mathrm{d}V/\mathrm{d}t\))، و \(t_{\text{max}}\) زمانی است که قلهٔ جریان رخ میدهد. این رابطه گذرای مشاهدهشده در آزمایش را به پارامتر بنیادی انتقالِ مورد نظر مرتبط میکند.
3. محدودیتهای CELIV
استخراج معادلهٔ تحرک CELIV بر چند فرض سادهکننده متکی است. نخست، فرض میشود که فقط یک نوع حامل بار بر انتقال غالب است، بهطوریکه تنها یک گونهٔ حامل هنگام اعمال رمپ ولتاژ حرکت میکند. دوم، فرض میشود حاملها بهطور یکنواخت توزیع شدهاند و بهصورت هموار از دستگاه خارج میشوند، مشابه حالت ایدهآلشدهای که در سمت چپ ?? نشان داده شده است. سوم، فرایندهای بازترکیب در طی استخراج ناچیز در نظر گرفته میشوند، بهطوریکه افت قابلتوجهی در حامل رخ نمیدهد. در واقعیت، این شرایط بهندرت برقرار هستند. همانگونه که در سمت راست ?? نشان داده شده است، حاملهای بار اغلب تلههای کمعمق، میانی یا عمیق را اشغال میکنند، و آزادسازی با تأخیر آنها گذرای استخراج پهنشده یا دگرشکلیافته تولید میکند. همانند بسیاری از تکنیکهای آزمایشی، تحرکی که با استفاده از CELIV استخراج میشود باید بنابراین بهعنوان یک تقریب در نظر گرفته شود نه مقدار میکروسکوپی واقعی، و در واقع تحرک ظاهری میتواند در طول خود گذرا نیز تغییر کند (doi:10.1063/1.4818267). این محدودیت بهویژه در نیمهرساناهای آلی نامنظم مهم است، جایی که بینظمی انرژی و تلهگذاری بهشدت بر انتقال اثر میگذارند. در مقابل، مواد پروسکایتی اغلب تحرکهای بالاتر و چگالی تلهٔ کمتری نشان میدهند، که تحلیل استاندارد CELIV را مقاومتر و تفسیر آن را آسانتر میکند.
💡 درک خود را بیازمایید: سعی کنید نظریهٔ CELIV را به یک کاربرد ساده اعمال کنید.
فرض کنید دستگاهی با ضخامت d = 200 nm،
نرخ رمپ ولتاژ A = 2 × 106 V/s،
و قلهٔ جریان مشاهدهشده در tmax = 5 µs دارید.
با استفاده از معادلهٔ تحرک CELIV، تحرک حامل را تخمین بزنید.
نمایش پاسخ
تحرک CELIV به صورت زیر داده میشود
\[
\mu = \frac{2 d^{2}}{3 A t_{\text{max}}^{2}}
\]
با جایگذاری مقادیر:
d = 200 × 10-9 m،
A = 2 × 106 V/s،
tmax = 5 × 10-6 s.
\[ \mu = \frac{2 (200 × 10^{-9})^{2}}{3 (2 × 10^{6}) (5 × 10^{-6})^{2}} \approx 1.1 × 10^{-4} \; \text{cm}^2/\text{Vs} \]
این مثال ساده نشان میدهد که چگونه زمان قله در یک گذرای CELIV میتواند به یک تخمین تحرک تبدیل شود.
4: اجرای یک شبیهسازی CELIV در OghmaNano
روی New simulation کلیک کنید. این کار کتابخانهٔ انواع دستگاههای موجود را باز میکند، که در ?? نشان داده شده است. روی Perovskite cells (که با رنگ قرمز برجسته شده) دوبار کلیک کنید تا پوشهٔ مثالهای پروسکایتی باز شود. فهرستی از شبیهسازیهای ازپیشتنظیمشده، شامل MAPbI₃ device، Perovskite solar cell، و نمایشهای اختصاصی CELIV را خواهید دید، همانطور که در ?? نشان داده شده است. برای این آموزش، Perovskite solar cell – CELIV example را انتخاب کنید. وقتی از شما خواسته شد، شبیهسازی را در پوشهای که دسترسی نوشتن دارید ذخیره کنید.
پس از انتخاب Perovskite solar cell – CELIV example، پنجرهٔ اصلی شبیهسازی باز میشود (نگاه کنید به ??). برای آغاز محاسبه، روی Run simulation (آیکون آبی پخش) کلیک کنید یا F9 را فشار دهید. OghmaNano سپس معادلات drift–diffusion وابسته به زمان را حل کرده و گذرای CELIV را تولید خواهد کرد.
jv.csv (دادههای منحنی JV)، time_j.csv (جریان استخراج برحسب زمان)،
time_v.csv (ولتاژ اعمالشده برحسب زمان)، و optical_output (توزیع میدانها) هستند.
برای مشاهدهٔ نتایج متناظر در ابزارهای ترسیم داخلی، روی هر فایل دوبار کلیک کنید.
گام 5: مشاهدهٔ نتایج CELIV
زبانهٔ Output را باز کنید
(??)
و روی time_v.csv دوبار کلیک کنید. این کار برنامهٔ ولتاژ اعمالشده —
رمپ CELIV — را که در
?? نشان داده شده ترسیم میکند.
سپس، روی time_j.csv دوبار کلیک کنید تا گذرای جریان استخراج نمایش داده شود،
همانطور که در ?? آمده است.
متوجه خواهید شد که گذرا شکل آشنای CELIV را دارد اما **وارونه** به نظر میرسد:
طبق قرارداد علامت در OghmaNano، جریانی که از دستگاه خارج میشود منفی است. وقتی رمپ در
?? اعمال میشود،
بار از دستگاه بیرون کشیده میشود، بنابراین دستگاه جریان منفی تحویل میدهد (قلهٔ استخراج).
قلهٔ بعدی با علامت مخالف (مثبت در
??) هنگامی رخ میدهد که بایاس
بازمیگردد و بار پس از بازگشت گذرا به حالت اولیه، دوباره با شتاب به داخل دستگاه سرازیر میشود.
نمودار در ?? گذرای CELIV را پس از ضرب کردن جریان شبیهسازیشده در 1- نشان میدهد تا با نمایش متعارف مورد استفاده در آزمایشها منطبق شود. در این فرم، خط مبنا متناظر با ظرفیت خازنی هندسی دستگاه است، در حالی که قلهٔ مثبت تیز نشاندهندهٔ استخراج حاملهای بار متحرک در طول رمپ ولتاژ اعمالشده است. اندازه و زمانبندی این قله کمیتهای کلیدی هستند که برای تعیین تحرک و چگالی حامل استفاده میشوند. ارائهٔ گذرا در این جهت متعارف، مقایسهٔ مستقیم شبیهسازیهای OghmaNano با اندازهگیریهای منتشرشدهٔ CELIV را آسانتر میکند. برای استخراج تحرک، زمانی را که جریان به بیشینهٔ خود میرسد (\(t_\text{max}\)) شناسایی کنید و این مقدار را در معادلهٔ تحرک CELIV \(\mu = \tfrac{2d^2}{3At_\text{max}^2}\) جایگذاری کنید، که در آن \(d\) ضخامت لایهٔ فعال و \(A\) نرخ رمپ ولتاژ اعمالشده است.
🧪 مثال حلشده: زمان قلهٔ CELIV را از
?? بخوانید و تحرک را تخمین بزنید.
فرض کنید ضخامت لایهٔ فعال d = 600 nm، نرخ رمپ ولتاژ
A = 5 × 105 V/s، و زمان قلهٔ اندازهگیریشده
tmax = 3 µs باشد. چه تحرکی به دست میآورید؟
نمایش پاسخ
از رابطهٔ CELIV استفاده کنید \[ \mu = \frac{2 d^{2}}{3 A t_{\text{max}}^{2}} \] با \(d = 600\times10^{-9}\,\mathrm{m}\)، \(A = 5\times10^{5}\,\mathrm{V\,s^{-1}}\)، \(t_{\text{max}} = 3\times10^{-6}\,\mathrm{s}\).
\[ \mu = \frac{2(600\times10^{-9})^{2}}{3\,(5\times10^{5})\,(3\times10^{-6})^{2}} \approx 5.3\times10^{-8}\ \mathrm{m^{2}\,V^{-1}\,s^{-1}} \] با تبدیل به \(\mathrm{cm^{2}\,V^{-1}\,s^{-1}}\) (\(1\ \mathrm{m^{2}} = 10^{4}\ \mathrm{cm^{2}}\)): \[ \mu \approx 5.3\times10^{-4}\ \mathrm{cm^{2}\,V^{-1}\,s^{-1}}. \]
این مقدار در بازهٔ معمول گزارششده برای نیمهرساناهای آلی نامنظم قرار دارد. (اگر \({t_\text{max}}\) شما با نرخ رمپ یا روشنسازی تغییر کند، خواهید دید که \(\mu\) نیز متناسب با آن جابهجا میشود.)
👉 گام بعدی: به بخش بعدی شبیهسازی CELIV (C-LIV) بروید — بخش B — جایی که: نرخ رمپ ولتاژ را تنظیم میکنیم، شدت روشنسازی را آزمایش میکنیم، و بررسی میکنیم که تغییر تحرک الکترون/حفره چگونه گذرای CELIV را دگرشکل میکند.