שיטת DFT-1/2 מדויקת למצבים פגומים רדודים: חישוב יעיל של אנרגיות קשירת תורם בסיליקון

מדוע התאמות זעירות בסיליקון חשובות

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

אטומים התורמים אלקטרונים נוספים

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

מדוע חישובים סטנדרטיים לא מספיקים

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

תיקון פשוט עם תועלת גדולה

המחברים בונים על טכניקה הנקראת DFT-1/2, שמוסיפה תיקון אנרגיית עצמי מקורב ישירות אל תוך DFT סטנדרטי. באופן פרקטי, הם משנים קמעה את הפוטנציאל היעיל של אטומים נבחרים על ידי הסרת חצי אלקטרון בצורה מושגית מאורביטלות אטומיות מסוימות. תחילה הם מיישמים את התיקון על סיליקון במצב מוצק כך שפתח הרצועה המחושב יתאים לניסוי הרבה יותר טוב ויספק רפרנס מהימן לרצועת ההולכה. לאחר מכן הם בוחנים את אופיו האלקטרוני של מצב התורם ומגלים שעבור כל דופנטי קבוצת V הוא נשלט על ידי אורביטל s של אטום הדופנט. הם מיישמים אז תיקון מותאם של חצי אלקטרון לאורביטל הזה ומכוונים רדיוס חיתוך יחיד כדי למקסם את ההפרדה בין רמת התורם לבין מצב ההולכה הריק הקרוב ביותר. חשוב שהתיקון המאופטם נשאר תקף כשמגדילים את תיבת הסימולציה, ולכן ניתן להשתמש בו מחדש בסופראלאים שמכילים עד אלפי אטומים.

עד כמה השיטה עובדת טוב

עם תיקון דו‑שלבי זה — קודם למארח הסיליקון ואחר כך לתורם עצמו — השיטה מניבה אנרגיות קשירת תורם שתואמות במידה רבה לערכים הנמדדים. עבור ארסן בסיליקון, האנרגיה החוזה נבדלת מהניסוי רק ב-0.3 מיליאלקטרון וולט, התאמה כמעט מושלמת ושווה ערך לשיטות היברידיות יקרות בהרבה. עבור אנטימון וזרחן השגיאות הן בערך 5 ו-8 מיליאלקטרון וולט בהתאמה, שיפור משמעותי על פני DFT לא מתוקן. עבור ביסמוט, דופנט כבד מאוד, המחברים כוללים גם קשירה ספין-מסלולית, אפקט יחסותי שמעט מעצב מחדש את רמות האנרגיה. זה מקטין את אנרגיית הקשירה המחושבת עד לגבול של כ-5 מיליאלקטרון וולט מניסוי ומדגיש פיזיקה ששיטות מוקדמות ויקרות יותר הזניחו. כדי להראות שהגישה אינה מוגבלת לסיליקון, הם מיישמים בהצלחה את אותו תהליך על תורם מימן בחמצן-אבץ (ZnO), שוב משחזרים אנרגיות קשירה שנמדדו בתוך כמה מיליאלקטרון וולט.

כלי מעשי לעיצוב שבבים של העתיד

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

ציטוט: Claes, J., Partoens, B., Lamoen, D. et al. An accurate DFT-1/2 approach for shallow defect states: efficient calculation of donor binding energies in silicon. npj Comput Mater 12, 153 (2026). https://doi.org/10.1038/s41524-026-02003-2

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