הדמיית מורפולוגיות של מבנים הטרוגניים בתלת-ממד באמצעות חתכים חשמליים דו-ממדיים

1. מבוא

מבנים הטרוגניים מרובי-פזות (BHJ) הם המבנה המיקרוסקופי העיקרי המשמש במכשירים אלקטרוניים אורגניים רבים – בעיקר בתאים סולאריים אורגניים, אך גם בגלאי אור, חיישנים ומערכות מעורבות יוניות-אלקטרוניות אחרות, שבהן חשובים מסלולים מקושרים וצווארי בקבוק מקומיים. במבנה BHJ, שתי פזות חודרות זו לזו בסקאלה של 10–100 ננומטר, ויוצרות ממשק פנימי גדול וסביבה של הולכה שאינה חד-ממדית באופן מהותי. המחיר לכך הוא מורכבות: ברגע שמגדירים את המבנה המיקרוסקופי, בלתי נמנע שמתקבלים מגבלות קישוריות, נקודות קצה, חסימות ושונות מרחבית חזקה בצפיפות נושאי המטען ובשילוב מחדש. עבור רוב ההתאמות היומיומיות של המכשיר, לא כדאי להתחיל מכאן. בדרך כלל עובדים ברמת המכשיר החד-ממדי – מהיר, אמין ולעתים קרובות מתאים לחלוטין – ואפשר לכלול גם פיזיקה נוספת, כגון לכידה ושילוב מחדש. לדוגמה, המדריך המהיר להתחלה בתחום תאים סולאריים אורגניים מדגים זרימת עבודה סטנדרטית חד-ממדית (כולל פיזיקת לכידה במידת הצורך), שהיא הכלי המתאים לרוב המחקרים והחילוצים של פרמטרים השגרתיים.

בקצה השני של סולם ההפשטה, אם המטרה שלכם היא לשאול שאלות ברמת המערכת, כגון הרחבת מערך, הפסדים בחיבור טורי/מקבילי, או התנהגות של מודול/פאנל, בדרך כלל תעברו לגישה של... בעל שטח גדול / מודע למעגל חשמלי תיאור במקום פתרון של מבנה בקנה מידה ננומטרי. מדריך מקיף על PM6:Y6, המיועד לשטחים גדולים מהווה דוגמה לגישה המתקדמת יותר בתחום המידול.

מדריך זה בוחר במכוון ללכת בכיוון הפוך: הוא שומר על המבנה המיקרוסקופי (באמצעות מורפולוגיה מוגדרת של תערובת הטרוגנית) אך נשאר מעשי מבחינה חישובית על ידי פתרון הבעיה החשמלית על חתכים דו-ממדיים. המטרה אינה שכל מכשיר יתאים לגישה זו, אלא שיש ערך מדעי וטכנולוגי בשימוש בתהליך עבודה המאפשר ניתוח של המבנה המיקרוסקופי. הוא מאפשר לבדוק את ההנחות של מודלים "אפקטיביים", לראות היכן מודלים "פשוטים" מפסיקים לעבוד, ולפתח הבנה לגבי הגבולות המורפולוגיים המשפיעים על הביצועים – במיוחד כאשר תופעות כמו הולכה ואינטראקציות מקומיות מרחבית דומיננטיות.

2. צרו מודל סימולציה של מורפולוגיה ובדקו את מבנה ה-BHJ.

מהחלון הראשי, לחצו על "הפעלה של סימולציה חדשה". פעולה זו תפתח את ספריית סוגי המכשירים המוצגת בתמונה ??. לחצו פעמיים על "מורפולוגיה", ולאחר מכן שמרו את הפרויקט במקום כלשהו בדיסק המקומי שלכם. שמירה מקומית חשובה מכיוון שסימולציות המבוססות על ניתוח מורפולוגי יכולות ליצור קבצים רבים זמניים (רשתות, צילומי מצב, קבצי מטמון), ואחסון מקומי הוא בדרך כלל מהיר ואמין יותר מאשר תיקיות המסונכרנות ברשת או בענן.

לאחר יצירת הפרויקט, נפתח חלון המכשיר הראשי. עליך לראות אובייקט מרכזי בצבע אדום, המייצג שלב אחד במבנה ה"בּוּלְק-הֶטֵרוֹ-ג'ַנְקְשֶׁן" (BHJ), הטבוע בתוך אזור השכבה הפעילה המקיפה, כפי שמוצג בתרשים ??. אמנם התצוגה התלת-ממדית המלאה מספקת מידע, אך המורפולוגיה יכולה להיראות צפופה מבחינה ויזואלית. כדי להבהיר את המבנה הפנימי, לחץ לחיצה ימנית על אובייקט ה-BHJ ובחר בתפריט "תצוגה" → "הצג חתך לאורך ציר Y".

הפעלת אפשרות התצוגה החותכת יוצרת את התצוגה המופיעה בתמונה ??, המציגה חתך דו-ממדי של המבנה התלת-ממדי. הדבר מקל מאוד על זיהוי הקשר בין הפאזות, הממשקים ונתיבי ההולכה בתוך תערובת ה-BHJ.

בשלב זה, מודל מנטלי מועיל הוא כדלקמן: האזור המלבני התוחם מגדיר את נפח השכבה הפעילה, בעוד שהאובייקט האדום המשובץ מגדיר את אחת הפאזות המוגדרות של תערובת הטרוגנית בתחום ה-BHJ (למשל, הפולימר). הפאזה המשלימה מוגדרת באופן מרומז כחלק הנותר מנפח השכבה הפעילה, שאינו תפוס על ידי האובייקט המשובץ. הפרדה זו הופכת למפורשת בעת בדיקת תכונות האובייקט בעורך האובייקטים, בסעיפים הבאים.

3. מבנה המכשיר: שכבה פעילה ומורפולוגיה של תערובת הטרו-הצמדה מוטמעת (BHJ).

פתחו את הלשונית "מבנה המכשיר" ולחצו על "עורך שכבות". המכשיר מורכב משלוש שכבות: מגע עליון, שכבה פעילה אחת ומגע תחתון (איור ??). מבנה פשוט זה מקל על הבנת האופן שבו המורפולוגיה משפיעה על הבעיה החשמלית, מבלי להוסיף ממשקים או שכבות הולכה נוספות. הנקודה המרכזית היא שהשכבה הפעילה עצמה היא נפח מלבני פשוט המוגדר בעורך השכבות. בדוגמה זו, היא מוגדרת כחומר מולקולרי Y6, שהוא מקבל, וכשלעצמה – היא פשוט תיבה. כל המורכבות המבנית הקשורה למבנה ה"צומת הטרוגנית" מוצגת בנפרד, על ידי הטמעת אובייקט שני בתוך נפח השכבה הפעילה. האובייקט המושתל הוא רשת CAD מורכבת המיובאת ממסד הנתונים Shape, והיא מייצגת את פאזת הפולימר PM6. רשת הפולימר מגדירה נפח תלת-ממדי סגור בתוך השכבה הפעילה. הפאזה המשלימה (האזור העשיר במקבלים) מוגדרת לאחר מכן באופן מרומז כחלק מהשכבה הפעילה שאינו תפוס על ידי רשת הפולימר. במילים אחרות, מבנה ה-BHJ מיוצג על ידי הימנעות נפח, ולא על ידי הגדרה מפורשת של שתי רשתות משולבות.

ההפרדה הזו הופכת ברורה כאשר בוחנים את האובייקטים בתצוגת התלת-ממד. ישנם שני אובייקטים רלוונטיים: תיבת השכבה הפעילה החיצונית (המוצגת בדרך כלל כקו מתאר כחול), והמורפולוגיה של הפולימר המשובץ (המוצגת באדום). כדי לבחון את התכונות שלהם, לחצו לחיצה ימנית על האזור הכחול החיצוני ובחרו באפשרות "ערוך אובייקט", כפי שמוצג בתמונה ??.

בעורך האובייקטים, ההבחנה ברורה. עבור האובייקט הפעיל הראשי, צורת האובייקט היא קובייה פשוטה (כאשר רשתות מורכבות מושבתות). עבור אובייקט הפולימר, צורת האובייקט מצביעה על רשומה במסד הנתונים של הצורות, שהיא הרשת המשולשת המגדירה את המורפולוגיה של PM6. שינוי הפניה זו מחליף את המורפולוגיה של הפולימר, תוך השארת השכבה הפעילה הסובבת ללא שינוי.

מבנה זה – קופסה פשוטה המכילה שכבה פעילה, בתוספת נפח פולימרי משובץ – הוא מרכזי באופן שבו מטפלים במורפולוגיות BHJ ב-OghmaNano. הוא הופך את הגאומטריה לקלה להבנה, מבטיח שהפאזה הפולימרית יוצרת אזור סגור מוגדר היטב, ומאפשר לחקור כיצד תופעות כמו חדירה, קישוריות והיצרויות ברשת הפולימרית משפיעות על הולכת מטען לאחר הפעלת פותר החשמל.

4. פרמטרים חשמליים: ניידות ואזורי מלכודת SRH

פתחו את עורך הפרמטרים החשמליים מחלון התצוגה הראשי. קיימים פאנלים נפרדים לפרמטרים עבור שתי הפאזות של השכבה הפעילה: אזור הפעילות התוחם (מקבל / מטריצה) והמורפולוגיה הפולימרית המשובצת, כפי שמוצג בתרשימים ?? ו-??. כאן מוגדרות שתי הפאזות החודרות זו לזו כך שהן יהיו מובחנות מבחינה חשמלית.

שתי הפאזות משתמשות באותו פורמליזם בסיסי של הולכה: דיפוזיה עם סטטיסטיקות של נשאי מטען על פי חוק מקסוול-בולצמן וטנזורים סימטריים של ניידות. ההבדל העיקרי הוא בניגודיות מכוונת בניידות. בפאזה הפולימרית, ניידות החורים גבוהה, בעוד שניידות האלקטרונים מופחתת; בפאזה הפעילה המשלימה, המצב הפוך. הדבר מאפשר הולכה סלקטיבית: חורים עוברים דרך רשת הפולימר, בעוד שהאלקטרונים עוברים דרך המטריצה.

א-סימטריה פשוטה זו מספיקה כדי ליצור התנהגות פיזקלית משמעותית במבני BHJ. עוד לפני הפעלת הסימולציה, היא מגדירה מה צפוי לראות בפלט המרחבי: מסלולי הולכה רציפים, נקודות מבוי סתומות שבהן נושאי המטען נלכדים גיאומטרית, והצטברות מטען מקומית בנקודות חניקה או באזורים בעלי קישוריות לקויה במורפולוגיה. בנוסף להולכת נושאי מטען חופשיים, מופעלת לכידת שוקלי-ריד-הול (SRH) שאינה בשיווי משקל בשתי הפאזות. כפי שמוצג בפאנלים של הפרמטרים, צפיפות המצבים מיוצגת באמצעות התפלגויות אקספוננציאליות עבור אלקטרונים וחורים, עם צפיפויות לכידה ושיפועים אנרגטיים זהים בכל פאזה. הפותר עוקב אחר חמישה פסי לכידה בדידים בכל נקודת רשת.

כל רצועת מלכודת מחליפה נושאי מטען באופן דינמי עם אוכלוסיות החופשיות באמצעות מנגנוני לכידה ופליטה מוגדרים (קבועי צימוד בין מצב חופשי למצב כלוא, הניתנים לצפייה בעורך). כתוצאה מכך, המערכת הנפתרת כוללת לא רק את צפיפות האלקטרונים והחורים החופשיים, אלא גם אוכלוסיות מרובות של נושאי מטען כלואים בכל נקודה במרחב. הדבר מגדיל באופן משמעותי את קשיחות החישוב ואת העלות החישובית, אך הוא חיוני לתיאור תופעות של תנועה דיספרסיבית, רקומבינציה בסיוע מלכודות והתנהגות דינמית מציאותית בחומרים אורגניים לא מסודרים. הערכים המספריים המשמשים כאן נועדו לשמש כדוגמה יציבה וממחישה, ולא כמודל חומר מכויל. לאחר שהבנו את תהליך העבודה, ניתן לחדד את הפרמטרים הללו, לאמת אותם באמצעות ניסויים, או לפשט אותם אם המטרה היא לבודד את המורפולוגיה ואת אפקטי הולכה ולא את הפיזיקה של המלכודות.

5. שילוב/חיבור/התאמה

לפני הפעלת הסימולציה, כדאי לבדוק את רשת החשמל כדי לוודא שהיא תקינה. בחלון הראשי, פתחו את הכרטיסייה "חשמל" (איור ??) ולחצו על "רשת חשמל". פעולה זו תפתח את עורך הרשת, כפי שמוצג באיור ??.

בדוגמה זו, הרשת היא בגודל 40 על 40 במישור xy, כלומר הבעיה החשמלית נפתרת בשני ממדים. כיוון z מושמט בכוונה. ניתן לבצע סימולציות חשמליות תלת-ממדיות מלאות ב-OghmaNano, אך הן דורשות מורפולוגיה תלת-ממדית אמיתית, וגורמות לעלייה משמעותית בעלות החישובית ובמורכבות. בנייה וניתוח של מודלים כאלה חורגים ממסגרת הדרכה זו.

במקום זאת, יש להתייחס לחישוב זה כאל חתך דו-ממדי של המכשיר. החתך עובר דרך מרכז המבנה ומכיל את המיקרו-מבנה ה-BHJ המלא במישור זה. גישה זו מאפשרת לזהות תופעות של הולכה, צווארי בקבוק וקישוריות בין הפאזות, תוך שמירה על מהירות, יכולת ניתוח והתאמה טובה לניתוח בקנה מידה לימודי.

6. הפעלת הסימולציה

לחצו על כפתור "הפעלת הדמיה" הכחול בחלון הראשי כדי להתחיל את פתרון הבעיה. במהלך ההרצה, הטרמינל מציג את נקודת ההטיה וצפיפות הזרם בכל שלב (??). בדוגמה זו, ההטיה המופעלת על המגע העליון משתנה בהדרגה מ-0.00 וולט לערך של כ-0.54 וולט, במטרה ליצור סדרת מדידות של זרם-מתח. בערכי הטיה נמוכים, צפיפות הזרם המדווחת היא שלילית תחת תאורה, וכאשר ההטיה מתקרבת למצב של מעגל פתוח, הזרם הכולל חוצה את האפס ומשנה את סימנו, בהתאם להתנהגות של תא סולארי.

לאחר סיום הריצה, פתחו את הלשונית "פלט" (??) כדי לבדוק את הקבצים שנוצרו. סיכום הפלט במסוף מציג את זמן הריצה הכולל, את מספר המשוואות שנפתרו, ונתונים בסיסיים על פעולות קלט/פלט בדיסק. בסימולציות המבוססות על ניתוח מורפולוגי, כדאי לעקוב אחר פלט הקבצים: כתיבה של מספר גדול של תמונות מצב עלולה להאריך את זמן הריצה, גם כאשר הפותר עצמו מהיר.

כדי לבדוק את תקינות הנתונים הגיאומטריים, לחצו פעמיים על device.csv כדי לראות את הייצוג המשולש של הרכיב, כפי שהוא משמש באופן פנימי על ידי הפותר (??). פעולה זו מאשרת שהמורפולוגיה של ה-BHJ יובאה והופרדה ליחידות כצפוי.

כדי לראות את הביצועים החשמליים, פתחו את jv_contact_0.csv או jv_contact_1.csv כדי להציג את גרף ה-JV (??). עבור סימולציה דו-ממדית זו, יש להשתמש בנתוני ה-JV עבור כל מגע; ניתן להתעלם מ-jv.csv.

7.0 בדיקת תמונות מצב מתוך הסימולציה

לאחר סיום הריצה, ניתן למצוא את התוצאות המועילות ביותר בספרייה snapshots/. קבצי הצילום מכילים את המשתנים הפנימיים של מצב המערכת של מודל התפשטות-דיפוזיה-SRH, והם נשמרים כפונקציה של המתח המופעל. הם חושפים את הפתרון המרחבי של משוואות התחבורה בכל שלב מתח, ולא סיכומים מעובדים. כאשר פותחים צילום, זמינים שני פקדים: סרגל מתח, המאפשר מעבר בין ערכי המתח, וסרגל מיקום חתך, המאפשר בחירת חתך דו-ממדי של המורפולוגיה. לאורך מדריך זה, הבעיה החשמלית נפתרת על חתך דו-ממדי, וכל נתוני הצילום מתייחסים לחתך זה.

איורים ??–?? מציגים דוגמאות של תמונות דו-ממדיות בנקודת מתח קבועה. יש לפרש את הנתונים הללו יחד. מפת פס ההולכה (LUMO), Ec.csv (??), מגדירה את הסביבה האנרגטית עבור אלקטרונים. בסימולציה המאפשרת לזהות את המורפולוגיה, מפה זו משתנה בצורה רוחבית בשל תלות אנרגטית בפאזה וכן תגובה אלקטרוסטטית לאגירת מטען מקומית. צפיפות האלקטרונים החופשיים, Q_nfree.csv (??), מדגישה אזורים שבהם אלקטרונים יכולים לנוע ביעילות דרך רשת הקולט/מטריצה. שינויים מרחביים משמעותיים מצביעים על מסלולי הולכה, נקודות חסימה וצווארי בקבוק הנוצרים עקב המורפולוגיה.

צפיפות האלקטרונים הכלואים, Q_ntrap.csv (??), מגלה היכן מצטברים נושאי מטען במצבי מלכוד. השוואה בין צפיפות חופשית לצפיפות כלואה היא לעיתים קרובות הדרך המהירה ביותר לזהות אזורים של חילוץ לקוי או מלכוד מקומי חזק הקשור למורפולוגיה של ה-BHJ. התוצאה הקשורה לניידות, mu_p_ft.csv (??), מציגה ניידות אפקטיבית ממוצעת על פני נושאי מטען חופשיים וכלואים. שינויים מרחביים משקפים הן את הניגודיות הפאזית הטבעית והן את המלכוד המקומי, ולעתים קרובות מתואמים באופן חזק עם הקישוריות של הרשת הבסיסית.

אמנם, תצוגות דו-ממדיות הן לעיתים קרובות הדרך הברורה ביותר לנתח שינויים מרחביים, אך לעיתים קל יותר לפרש את אותם נתונים במרחב תלת-ממדי. על ידי הפעלת מצב תלת-ממדי בחלון הצילום, הפרוסה הדו-ממדית מוקרנת בחזרה אל הגיאומטריה של המכשיר. איורים ??–?? מציגים דוגמאות להקרנה זו. ערכי הממוצע של ניידות האלקטרונים והחורים (mu_n_ft.csv, mu_p_ft.csv), רמת הקוואזי-פרמי של האלקטרונים החופשיים (Fn.csv), וצפיפות זרם האלקטרונים בכיוון התנועה (Jn.csv) שימושיים במיוחד לצורך המחשת מסלולי הולכה וזרימת זרם דרך שכבת ה-BHJ.

תצלומי המסך הללו מספקים את הערך המדעי העיקרי של סימולציות המאפשרות ניתוח מורפולוגי: הם מאפשרים לקשור ישירות תכונות הניתנות למדידה בקנה מידה של רכיב לאנרגטיקה מקומית, לתהליכי לכידה ולקישוריות בתוך המבנה ההטרוגני.

8. משימות שכדאי לנסות

הדרך המהירה ביותר לפתח אינטואיציה היא לשנות את המודל ולראות אילו שינויים מתרחשים. המשימות המפורטות להלן נועדו להיות קצרות: החליפו את רשת המורפולוגיה, הפעילו שוב את המודל, והשוו את עקומת ה-JV (ואם יש לכם זמן, צלמו כמה תמונות מצב באותו מתח).

משימה 1 – החלף את רשת המורפולוגיה

בתצוגה הראשית של המכשיר, לחצו לחיצה ימנית על האובייקט האדום BHJ (רשת הפולימר המשובצת) ובחרו באפשרות "עורך רשת". בעורך, ודאו שהאפשרות הנוכחית היא "צור צורה מתוך מסד הנתונים". שנו את הצורה שנבחרה ל-morphology/1, morphology/2 ו-morphology/3 (איורים ??, ??, ??). הפעילו מחדש את הסימולציה לאחר כל שינוי והשוו את עקומת ה-JV.

מה שאתם מחפשים הוא כיצד קישוריות וצווארי בקבוק משפיעים על תהליכי חילוץ ושילוב מחדש: חלק מהצורות הפיזיות יוצרות מסלולי זרימה נקיים יותר (זרם גבוה יותר), בעוד שאחרות מכוונות את הנשאים דרך צווארים צרים ונקודות קצה (הצטברות רבה יותר, לכידה רבה יותר, ובדרך כלל ביצועים גרועים יותר של התא).

כחלק מניסוי מכוון שנועד להיות "שגוי אך מלמד", תמצאו גם צורה בשם "קומקום" במסד הנתונים. החליפו את רשת ה-BHJ בצורה של קומקום, הפעילו שוב, וצרו עקומת JV עבור התקן עם תוספת בצורת קומקום. זה אינו BHJ בעל משמעות פיזיקלית, אך הוא מדגים שהתהליך מונע באופן מהותי על ידי הגיאומטריה: ניתן להכניס כל רשת CAD סגורה ולנתח אותה מבחינה חשמלית.

משימה 2 – שינוי עוצמת לכידת האלקטרונים (רצועות לכידה מסוג SRH).

בפרמטרים החשמליים, התאימו את הגדרות מלכודות SRH שאינן במצב שיווי משקל, ולאחר מכן הפעילו שוב את הסימולציה. נסו לשנות פרמטר אחד בכל פעם:

• הגדילו או הקטינו את צפיפות המלכודות בסדר גודל של 1-2, והשוו בין Q_nfree ל-Q_ntrap.

מטרת המשימה הזו היא לבחון כיצד תופעת הפיזור/הלכידה משפיעה על תנועת החלקיקים ומתחרה בה, וכיצד היא משפיעה על המבנה הפיזי: לכידת חלקיקים נוטה להגביר את השפעתם של אזורים "בעייתיים" (נקודות קצה, אזורים מצומצמים), מכיוון שהמטען החשמלי אינו יכול להתפזר במהירות.

משימה 3 – שנה את קבועי קצב הרקומבינציה

עדיין תחת הפרמטרים החשמליים, שנה את ערכי הרקומבינציה והפעל מחדש. שתי הפרעות פשוטות:

• שנו את קבוע קצב הרקומבינציה בין מצבים אנרגטיים זהים (n_free ל-p_free) כלפי מעלה או כלפי מטה בשיעור של 1–2 סדרי גודל. שימו לב כיצד עקומת ה-JV והתפלגות המרחבית של המטען מגיבים לשינוי.
• אם הפרויקט שלכם חושף פרמטרים של לכידה/צימוד של SRH (קבועי צימוד בין מצב חופשי למצב לכודה), הגדילו או הקטינו את הפרמטרים הללו יחד, והשוו את היחס בין Q_nfree ל-Q_ntrap באותו מתח הטיה.

הרגל מועיל לביצוע משימות 2–3 הוא לשמור על המורפולוגיה קבועה, לשנות פרמטר אחד, ואז לבדוק בכל פעם את שני הדברים הבאים: (i) כיצד עקומת ה-JV משתנה, ו-(ii) האם מסלולי הזרם (למשל, Jn ושדות רמת פרמי) הופכים לחלקים יותר או "משובשים" יותר בשל המיקרו-מבנה.