אנו מציגים מודל אנליטי לחקירת התפשטות גלים מגנטיים-פוטו-תרמואלסטיים מסתובבים במוליכים למחצה מחוזקים בסיבים אניזוטרופיים תחת עירור אופטי. השיטה יכולה לשמש לניתוח טמפרטורה, צפיפות נושאים, תזוזה והתפלגויות מתח בחומרים מתקדמים של מוליכים למחצה.
הרשמה ל-JoVE נדרשת לצפייה בתוכן זה. התחברו או התחילו בגרסת ניסיון בחינם.
מאמר מחקר
אנו מציגים מודל אנליטי לחקירת התפשטות גלים מגנטיים-פוטו-תרמואלסטיים מסתובבים במוליכים למחצה מחוזקים בסיבים אניזוטרופיים תחת עירור אופטי. השיטה יכולה לשמש לניתוח טמפרטורה, צפיפות נושאים, תזוזה והתפלגויות מתח בחומרים מתקדמים של מוליכים למחצה.
מחקר זה בוחן את ההתנהגות המגנטו-פוטו-תרמואלסטית המסתובבת המקושרת של חצי מרחב מוליך למחצה מחוזק בסיבים אניזוטרופיים, הנתון לעירוור אופטי ולשדה מגנטי מיושם. המשוואות השולטות, המשלבות השפעות תרמיות, אלסטיות, צפיפות נשאי, אלקטרומגנטיות וסיבוביות, מנוסחות במסגרת מולטיפיזית מאוחדת. שיטת המצב הנורמלי משמשת כדי להמיר את המשוואות הדיפרנציאליות החלקיות המקושרות למערכת של משוואות דיפרנציאליות רגילות, אשר לאחר מכן נכתבת מחדש כמערכת דיפרנציאלית מטריצית מסדר ראשון. לאחר מכן מפותחת נוסחה מבוססת ערך עצמי לבניית הפתרון האנליטי של הבעיה המקושרת. בעיית הערכים העצמיים המתקבלת, מערכת תנאי הגבול, וכמויות השדות מוערכים נומרית. מוצגות תוצאות נומריות עבור טמפרטורה, צפיפות הנשואים, רכיבי ההעתקה והתפלגויות המאמצים כדי לבחון את השפעות פרמטר הסיבוב ועוצמת השדה המגנטי על השדות הפיזיקליים המקושרים. התוצאות מדגימות את התפקיד החשוב של אינטראקציות סיבובית ואלקטרומגנטיות בשינוי התפשטות גלים, מאפייני החלשה והתפלגויות שדה בתוך התווך הסמי-מוליך האניזוטרופי. בוצעו חישובים נומריים עבור מדיום מוליכים למחצה מבוסס סיליקון.
פוטו-תרמואלסטיות הפכה לתחום מחקר פעיל בזכות יכולתה לתאר את האינטראקציות המקושרות בין שדות תרמיים, מכניים, אופטיים ושדות צפיפות נשאים בחומרים למחצה שנחשפו לקרינת לייזר ועירור פוטותרמי. תופעות משולבות אלו ממלאות תפקיד חשוב במכשירים אופטואלקטרוניים מודרניים, טכנולוגיות מוליכים למחצה, מערכות עיבוד לייזר ויישומים מיקרואלקטרוניים. סעיד ואח' חקרו התפשטות גלים במוצקים למחצה בעלי תכונות תלויות טמפרטורה, בעוד עבאס ואחרים בחנו אינטראקציות פוטו-תרמיות במוליכים למחצה המכילים חללים גליליים ומוליכות תרמית משתנה. בהמשך, איעתישאם אולה ואחרים ניתחו את השפעת מוליכות תרמית משתנה ופולסי לייזר על גלים אלסטיים מוחזרים במדיה מוליכים למחצה. בנוסף, יאדב4 חקר התפשטות גל פלזמה מגנטו-פוטותרמית במוליכים למחצה דיפוזיוניים תחת מסגרת תרמואלסטית רב-פאזית דו-טמפרטורתית. לאחרונה, Lute ואחרים 5 ו-Lute ואחרים.6 פיתחו מודלים מתקדמים של מוליכים למחצה פוטו-תרמואלסטיים הכוללים יצירת חום פוטותרמי, השפעות זיכרון ואינטראקציות לא מקומיות, בעוד שעלאופי ואל-דאלי7 הציעו פורמולציה מולטיפיזית של פוטו-תרמו-היגרואלסטי לחומרים מוליכים למחצה. יתרה מזאת, דווח על מספר מודלים פוטו-תרמואלסטיים כלליים המבוססים על תיאוריות השהיית פאזה ונוסחאות לא מקומיות ששיפרו את תיאור תופעות ההובלה התרמית והנשאית המקושרת במוליכים למחצה 8,9,10,11.
הביקוש הגובר למבנים הנדסיים קלים ובעלי ביצועים גבוהים עורר מחקר נרחב על חומרים אניזוטרופיים וחומרים מחוזקים בסיבים. בזכות מאפייניה המכניים המעולים ועמידותם התרמית המוגברת, קומפוזיטים מחוזקים בסיבים שימשו באופן נרחב ביישומים של תעופה וחלל, אזרחית והנדסה אלקטרונית. חאן ואח' פיתחו מודל תרמואלסטי מיקרופולרי מאוחד לחומרים מחוזקים בסיבים התלויים בטמפרטורה והראו את ההשפעה המשמעותית של פרמטרים מיקרוסטרוקטוריים על מאפייני התפשטות הגלים. Pitarresi ואח' חקרו את התגובה התרמואלסטית של פלסטיקים מחוזקים בסיבים והדגישו את תפקיד ההטרוגניות של החומר בהתפלגות המתח התרמי. Escalante-Solís ואחרים בחנו את השפעת עקמומיות הסיבים על ההתנהגות המיקרומכנית של קומפוזיטים מחוזקים, בעוד אבו-דהאב ואח' חקרו תופעות החזרת גלים במדיומים תרמואלסטיים מחוזקים בסיבים. לאחרונה, פורקייט וקנוריה16 ניתחו מדיום מגנטו-תרמואלסטי מחוזק בסיבים הנתון לגירוי לייזר פולס, בעוד שאקאי ואח' 17 וקווינלן ואח' 18 חקרו תגובות תרמואלסטיות במבני קומפוזיט מחוזק בסיבים בתנאי עומס שונים. מחקרים אלו אישרו כי מאפייני החיזוק והאניזוטרופיה משפיעים באופן משמעותי על ההובלה התרמית, ריכוזי המתח והתנהגות התפשטות הגלים. עם זאת, רוב המחקרים הזמינים הוקדשו לחומרים מרוכבים תרמואלסטיים ולא לקחו בחשבון בו-זמנית את הקישור של מוליכים למחצה פוטו-תרמואלסטיים יחד עם מנגנוני הובלת נשאי.
האינטראקציה בין שדות מגנטיים למדיה תרמואלסטית למחצה זכתה גם היא לתשומת לב רבה בשל חשיבותה במכשירים אלקטרומגנטיים ובחומרים מוליכים. יאדב19 חקר השתקפות גלים במרחב מגנטו-תרמואלסטי אורטוטרופי מסתובב עם אפקטי דיפוזיה, בעוד סלוואמני ואחרים חקרו התפשטות גלים מגנטואלסטיים בננו-קרניים שבריריות לא מקומיות. אבואלרגל ואח' בחנו גירוי משולב של לייזר-תרמומגנטי במבנים נקבוביים ממוגנטים עם השפעות תלויות בזיכרון, בעוד שצ'נדל ואחרים הציעו מסגרת מגנטו-תרמואלסטית לא מקומית המבוססת על הולכת חום מור-גיבסון-תומפסון. בהקשר של מוליכים למחצה, אבואלרגל23 פיתח מודל פוטו-תרמואלסטי חלקי מותאם לחצי מרווח מוליך למחצה מסתובב הנתון לשדה מגנטי, ו-Sur24 חקר אינטראקציות מגנטו-פוטו-תרמואלסטיות במדיה המציגה מאפיינים תורשתיים. לאחרונה, סעידי ואח' וראשיד ואח' 26 ניתחו הפרעות מגנטו-פוטו-תרמואלסטיות במדיה כללית של מוליכים למחצה מסתובבים והראו את ההשפעה המשמעותית של עוצמת מגנטית על התפשטות גלים מקושרים. תצפיות דומות דווחו במודלים מתקדמים של מוליכים למחצה פוטו-תרמו-ויסקואלסטייםושל מור-גיבסון-תומפסון 27,28,29.
השפעות סיבובי הן היבט חשוב נוסף של תרמואלסטיות כללית, שכן הסיבוב מוסיף כוחות קוריוליס וצנטריפטליים נוספים שמשנים משמעותית את מאפייני התפשטות הגלים. חאן ואח' ביצעו ניתוח רגישות של התנהגות גלים במוצקים תרמואלסטיים מסתובבים שנחשפו לעומס תרמי מושרה בלייזר. Yadav31 חקר גלים מגנטותרמואלסטיים במדיה אורתוטרופית מסתובבת עם אפקטי דיפוזיה, בעוד Abouelregal ואחרים 32 בחנו תגובות תרמואלסטיות חלקיות דו-פאזיות במדיה מתוחה מסתובבת המכילות חללים כדוריים. יתרה מזאת, חאן ואח' חקרו התפשטות גלים במרחב תרמואלסטי נקבובי מסתובב תחת מודל שברירי לא מקומי של השהיה תלת-פאזית. השפעות סיבוביות על ננו-מבנים תרמואלסטיים למוליכים למחצה נבדקו גם על ידי אבו-דהאב ואח' 34, שהוכיחו שסיבוב יכול לשנות באופן משמעותי התפלגות שדות תרמית ומכנית. למרות מאמצים אלו, השפעות סיבובי נדירות נחקרו במקביל לאניזוטרופיה, חיזוק סיבים, עירור אופטי, הובלת נשאים ואינטראקציות שדה מגנטי במסגרת אחידה של מוליכים למחצה.
למרות שהושגה התקדמות משמעותית במחקרים שהוזכרו, רוב המחקרים הקיימים התמקדו בשילובים נבחרים של אפקטים פוטותרמואלסטיים, מגנטיים, סיבובי או חיזוק. מודל אנליטי מקיף שמשלב בו-זמנית עירור אופטי, דינמיקת העברת נשאי, אניזוטרופיה, חיזוק סיבים, אינטראקציות שדה מגנטי והשפעות סיבובי בתווך מוליכים למחצה נותר ברובו בלתי זמין. לכן, עדיין יש צורך בניסוח מאוחד שמתאר במדויק את הקישור ההדדי בין המנגנונים הפיזיקליים הללו ומבהיר את השפעתם המשולבת על התפשטות גלים.
בהשראת תצפיות אלו, העבודה הנוכחית מפתחת מודל אנליטי דו-ממדי לחקר התפשטות גל מגנטו-פוטו-תרמואלסטי בחצי מרחב מוליך למחצה אניזוטרופי מסתובב, מחוזק בסיבים, הנתון לעוררות אופטית. המשוואות המנהלות מנוסחות במסגרת מאוחדת המשלבת שדות תרמיים, מכניים, צפיפות נשאים ושדות אלקטרומגנטיים, ונפתרות באמצעות טכניקת מצב נורמלי יחד עם גישת ערכי עצמי. ראוי לציין כי אלאופי ואח' חקרו לאחרונה את התגובה התרמו-פוטואלסטית המקושרת של סיליקון אניזוטרופי מחוזק בסיבים באמצעות נוסחת ערך עצמי. עם זאת, השפעות שדה מגנטיות ותרומות סיבוביות לא נלקחו בחשבון במודל שלהם. החידוש של העבודה הנוכחית טמון בהרחבת הניסוח הזה על ידי שילוב בו-זמנית של אינטראקציות מגנטיות, אפקטים סיבוביים, דינמיקת העברת נשאי, חיזוק אניזוטרופי ועירוור אופטי במסגרת מגנטו-פוטו-תרמואלסטית אחת. כתוצאה מכך, המודל המוצע מספק תיאור מקיף יותר של תופעות התפשטות גלים מקושרים ומאפשר הערכה מפורטת של ההשפעה המשולבת של עוצמת סיבוב מגנטית על השדות הפיזיקליים הרלוונטיים. בניגוד ל-Alaofi ואח' 35, שבחנו את התגובה התרמו-פוטואלסטית של מוליך למחצה מחוזק בסיבים אניזוטרופי ללא שדה מגנטי ואפקטי סיבוב, הנוסחה הנוכחית משלבת אינטראקציה סימולטנית של סיבוב, שדה מגנטי, עירור אופטי, העברת נשאים וחיזוק אניזוטרופי במסגרת מאוחדת. בנוסף, מימוש הערך העצמי מורחב למערכת מקושרת כללית זו, ומאפשר ניתוח של השפעות פיזיקליות מקושרות שלא טופלו במודל הקודם. יתרה מזאת, מאפיין ייחודי של העבודה הנוכחית הוא יישום גישת הערך העצמי. בעוד שמחקרים קשורים רבים, כולל Alaofi ואחרים, מסתמכים על נהלי אלימינציה כדי להפחית את מערכת הממשל לפני קבלת הפתרון, הניסוח הנוכחי נבנה ונפתר באמצעות מסגרת ערכים עצמיים מבוססת מטריצה. מימוש זה מיושם על מערכת מצומדת כללית יותר הכוללת סיבוב, שדה מגנטי, עירור אופטי, העברת נשאים וחיזוק סיבים, ובכך מרחיב את היקף הדגמים שפורסמו קודם.
הגישה מוגבלת. התחברו או התחילו תקופת ניסיון כדי לצפות בתוכן זה.
העבודה הנוכחית מבוססת כולה על מידול אנליטי, חישוב סימבולי וסימולציות נומריות לחקירת אינטראקציות מגנטו-פוטו-תרמואלסטיות מסתובבות במדיות למחצה אניזוטרופיות. לא השתתפו במחקר זה משתתפים אנושיים, ניסויים בבעלי חיים, נתונים קליניים או דגימות ביולוגיות. לכן, לא נדרשו אישור אתי והסכמה מדעת.
ניסוח מתמטי של בעיית המגנטו-פוטו-תרמואלסטי בחצי מרחב סיבוב מחוזק על ידי סיבים למחצה אניזוטרופית
במחקר הנוכחי, נחקר חצי מרחב מוליך למחצה אניזוטרופי דו-ממדי ומסתובב עם סיבים אניזוטרופית תחת עירור אופטי ושדה מגנטי מיושם. המדיום תופס את האזור החצי-אינסופי x ≥ 0, שבו הגבול ב-x = 0 מייצג את המשטח החשוף הנתון לעומס אופטי חיצוני. מערכת הקואורדינטות נבחרת כך שציר ה-x מתארך לתוך התווך, בעוד ציר ה-y נמצא לאורך המשטח, ומייצג את ההתנהגות המישורית של המבנה. השדה המגנטי המופעל וקטור המהירות הזוויתית נלקחים שניהם לאורך ציר ה-z. התווך המוליך למחצה נחשב הומוגני ואלסטי ליניארית, בעוד שאניזוטרופיה מוכנסת דרך סיבים מחזקים מיושרים המוטמעים בכיוון x-. בליעה אופטית בגבול יוצרת נשאי חום מקומיים ונשיאי מטען עודפים, מה שמוביל לאינטראקציות תרמיות, מכניות ונשאיות משולבות בתוך התווך. בנוסף, השדה המגנטי יוצר אפקטים של קישור אלקטרומגנטי, בעוד שהתנועה הסיבובית תורמת להשפעות אינרציאליות המשפיעות משמעותית על התפשטות גלים תרמואלסטיים ועל התגובה הפיזיקלית הכוללת. בהתאם לכך, המצב הפיזי של התווך מיוצג על ידי שדה הטמפרטורה T(x, y, t)(K), צפיפות הנשאים N(x, y, t)(m-3), ורכיבי ההזזה u(x, y,t)(m) ו-v(x,y,t)(m), בהנחה של עיוותים קטנים. איור 1 ממחיש את הגיאומטריה של הבעיה, כולל מערכת הקואורדינטות, עירור אופטי, שדה מגנטי, אפקט סיבובי וכיוון הסיבים. הנוסחה הנוכחית מתאימה למדיה למחצה מחוזקת בסיבים אניזוטרופית הומוגנית הפועלת בתחום העיוות הקטן ומסגרת התרמואלסטיות הליניארית. המודל מניח כיוון סיב קבוע ותכונות חומר קבועות לאורך כל התווך. כתוצאה מכך, התנהגות חומר לא ליניארית, עיוותים גדולים, נזק לחומר ושינויים מרחביים בתכונות החומר אינם נלקחים בחשבון במחקר הנוכחי. לכן, המודל המוצע מיועד לתנאי עומס מתונים שבהם התגובה נשארת בטווח הליניארי. במחקר הנוכחי, ההתעוררות האופטית מדומה באמצעות תנאי גבול שנקבעו לטמפרטורת פני השטח ולצפיפות הנשאים שנוצרת על ידי הפוט. תהליך האינטראקציה המפורט בין לייזר לחומר, כולל בליעה אופטית, עומק חדירה והתפלגות עוצמה, אינו מטופל במפורש. במקום זאת, ההשפעה הכוללת שלו מיוצגת על ידי אמפליטודות הגבול θ0 ו-N0, שמאפיינות את ההתרגשות התרמית והנשאים הנגרמת על ידי השדה האופטי הנכנס.

איור 1: ייצוג סכמטי של חצי מרחב מוליך למחצה אניזוטרופי מסובב, מחוזק בסיבי, הנתון לעוררות אופטית ושדה מגנטי חיצוני. האיור ממחיש את התצורה הפיזיקלית של הבעיה, כולל מערכת הקואורדינטות, עירור אופטי, שדה מגנטי יישום, כיוון סיבים ואפקטים סיבובי שנבחנו בנוסחה הנוכחית. אנא לחצו כאן כדי לצפות בגרסה מוגדלת של הדמות הזו.
הקשר המרכיב של טנזור המתח בתווך חצי-מוליך תרמו-אלסטי אניזוטרופי מחוזק בסיבים ניתן לביטוי בצורה המוכללת כך: 12,15,16.
(1)
כאן, σij מציין את רכיבי טנזור המתח, Cijkl הם מקדמי הקשיחות האלסטיים, ekl מייצג את טנזור העיוות, T הוא עליית הטמפרטורה ביחס לטמפרטורת הייחוס T0, ו-N מציין את צפיפות הנשא העודפת. הטנזורים βij ו-ηij מתאימים למקדמי הקישור התרמואלסטי והנשא, בהתאמה. בהתאם לכך, היחס המרכיב הכולל את ההשפעה המפורשת של חיזוק הסיבים ניתן לכתיבהכ-12,15:
(2)
בנוסחה זו, λ ו-μT הם קבועי האלסטיות של למה, בעוד μL מציין את מודולוס הגזירה האורכי לאורך כיוון הסיב. הפרמטרים α ו-β מתארים את אפקטי החיזוק הקשורים לסיבים המוטמעים. הכמות δij היא סמל דלתא של קרונקר, ו-a i הם הרכיבים של וקטור היחידה המגדיר את כיוון הסיב. במודל הנוכחי, סיבי החיזוק מיושרים לאורך כיוון x כך ש-a = (1,0). המונחים הכוללים אתβij θ ו-ηijN מייצגים השפעות תרמיות וקישור נשאי, בהתאמה. בתצורה הדו-ממדית הנוכחית, רכיבי המאמץ השולטים מצטמצמים לצורות הבאות12,15:
. (3)
. (4)
. (5)
כאן,
ו-
מסמנים את רכיבי ההעתקה לאורך כיווני x - ו-y, בהתאמה, בעוד ש-A ij הם המקדמים האלסטיים האפקטיביים של התווך המחוזק בסיבים אניזוטרופיים. מקדמי הקישור התרמואלסטי והנשא מוגדרים כך
,
,
,
.
ביחסים הנ״ל, αij מייצגים את מקדמי ההתפשטות התרמית, בעוד ξij מסמנים את מקדמי התפשטות הנשא הקשורים לתווך המוליך למחצה. המקדמים האלסטיים האפקטיביים של התווך האניזוטרופי המחוזק בסיבים ניתנים על ידי
,
,
, . 
מקדמים אלה מאפיינים את התגובה האלסטית האניזוטרופית של חומר המוליך למחצה מחוזק ומתארים כיצד כיוון הסיבים משפיע על התנהגות התרמואלסטית המקושרת. כדי להסביר את השפעת האינטראקציות האלקטרומגנטיות בפורמולציה הנוכחית של מגנטו-פוטו-תרמואלסטי, מניחים ששדה מגנטי אחיד מיושם לאורך כיוון z, שהוא נורמלי למישור העיוות x - y. בהתאם לכך, וקטור השדה המגנטי נחשב בצורה23,25 ,
, כאשר H0 מציין את עוצמת השדה המגנטי הקבועה. מכיוון שהניסוח הנוכחי מוגבל לעיוותים דו-ממדיים, שדה ההעתקה של התווך נלקח כ-
, כאשר
ו-
מייצגים את רכיבי ההעתקה לאורך כיווני x ו-y- בהתאמה.
בהנחות של עיוותים קטנים ומדיום מוליך למחצה חשמלי הנע באיטיות, האינטראקציה בין מהירות החלקיק לשדה המגנטי המופעל יוצרת שדה חשמלי מושרה.
בהתבסס על יחסי מקסוול האלקטרומגנטיים לתנועת מדיה מוליך, וקטור השדה החשמלי המושר ניתן לבטאכ-19,23
. (6)
כאשר μ 0 מציין את החדירות המגנטית ו
- הוא וקטור מהירות החלקיקים. החלפת הביטויים של
ו-
ביחס הנ"ל נותנת
. (7)
שמניב
. (8)
בלקחת את הנגזרת בזמן של השדה החשמלי המושרה, נקבל
. (9)
וקטור ההפרעה המגנטית שנוצר עקב עיוות של התווך המוליך למחצה מוגדרכ-23˒25
. (10)
הביטוי לעיל מקיים אוטומטית את תנאי הדיברגנס של מקסוול עבור שדה ההפרעה המגנטית, כלומר
. (11)
כדי לקבוע את צפיפות הזרם החשמלי, מעריכים תחילה את העיקול של וקטור ההפרעה המגנטית. בצורת דטרמיננטה, ניתן לכתוב את אופרטור הסילולכ-23
. (12)
הרחבת הדטרמיננטה מובילה ל-
. (13)
וקטור צפיפות הזרם החשמלי מתקבל אז מהמשוואה האלקטרומגנטית של מקסוול23
. (14)
כאשר ε0 מציין את הפרמיטיביות החשמלית של התווך.
כוח הגוף האלקטרומגנטי הפועל על התווך המוליך למחצה נקבע באמצעות יחס כוח לורנץ23
. (15)
מכפלת
הווקטור ניתנת להערכה בצורה דטרמיננטית כ-
. (16)
בהחלפת הביטויים הקודמים במשוואה (15), רכיבי וקטור כוח הגוף האלקטרומגנטי הופכים ל
. (17)
. (18)
יחסים אלו מצביעים בבירור על כך שהשדה המגנטי המופעל תורם מנגנוני קישור נוספים למשוואות השולטות הן באמצעות מונחים אלקטרומגנטיים דמויי קשיחות והן באמצעות מונחים אינרציאליים מותאמים הפרופורציונליים ל
-. כתוצאה מכך, השדה המגנטי משפיע משמעותית על תכונות ההתפשטות של גלים תרמואלסטיים ועל ההתנהגות הדינמית הכוללת של התווך המוליך למחצה אניזוטרופי מחוזק בסיבים. בנוסף להשפעות האלקטרומגנטיות, השפעת הסיבוב משולבת בניסוח הנוכחי כדי לתאר את התגובה הדינמית של התווך כאשר נצפים ממסגרת התייחסות מסתובבת. המדיום המוליך למחוצה בסיבים נחשב לעבור סיבוב אחיד של גוף קשיח עם וקטור מהירות זוויתי קבוע הניתןעל ידי 33
, כאשר
מסמן את מהירות הזווית הקבועה סביב ציר z. מכיוון שציר הסיבוב הוא נורמלי למישור x - y, כאשר משוואות התנועה מנוסחות במערכת קואורדינטות מסתובבת, נוצרות תאוצות אינרציאליות נוספות בשל אופיה הלא-אינרציאלי של המסגרת המסתובבת. תאוצות אלו מורכבות בעיקר מהאצות קוריוליס וצנטריפוגליות. האצה של קוריוליס משויכת לשדה מהירות החלקיקים ומבוטאתכ-31˒33
. (19)
בהחלפת הביטויים של
ו-
, נקבל
. (20)
לכן, האצה של קוריוליס הופכת ל
. (21)
התאוצה הצנטריפוגלית תלויה ישירות בשדה ההעתקה עצמו ומיוצגת על ידי31.35
. (22)
ראשית, מכפלת
הווקטורים מוערכת כך
. (23)
לאחר מכן, החלפת התוצאה שהושגה ביחס התאוצה הצנטריפוגלית מניבה
. (24)
בהתאם לכך, ניתן לבטא את התרומה הסיבובית הכוללת המופיעה במשוואות השלטון כך:
. (25)
לכן, רכיבי התאוצה הסיבובית בכיווני x ו-y הם
, (26)
. (27)
הביטויים הנ״ל מראים כי תנועת הסיבוב יוצרת קישור נוסף בין רכיבי ההעתקה דרך תאוצת קוריוליס, בנוסף להשפעות האינרציאליות התלויות בתזוזה הנובעות מההאצה הצנטריפוגלית. לכן, הפעולה המשולבת של השדה המגנטי והסיבוב יוצרת שינויים משמעותיים בתגובה הדינמית ובמאפייני הפצת הגלים של התווך המוליך למחצה אניזוטרופי מחוזק על ידי סיבים מסתובבים. כדי לשלב את ההשפעה המשולבת של אינטראקציות אלקטרומגנטיות ותנועת סיבוב, משוואות התנועה של המדיום המוליך למחצה אניזוטרופי מחוזק בסיבים מוכללות לכלול הן את כוח הגוף האלקטרומגנטי והן את ההאצות האינרציאליות הנוספות הנובעות במסגרת התייחסות מסתובבת. בהתאם לכך, משוואת התנועה הכללית עבור רצף מסתובב מתפורר ניתנת לביטוי כך:31˒32
. (28)
כאן, ρ מציין את צפיפות המסה, ו-Fi מייצג את רכיבי כוח הגוף האלקטרומגנטי. יתרה מזאת,
מסמן את וקטור המהירות הזוויתית של המסגרת המסתובבת. ה-EQS. (22) ו-(24) מתאימים לתאוצות קוריוליס וצנטריפוגליות, בהתאמה. משוואות התנועה בכיווני x ו-y ניתנות לכתיבה כך:
. (29)
. (30)
החלפת רכיבי כוח הגוף האלקטרומגנטי שהתקבלו קודם במשוואות לעיל נותנת
. (31)
. (32)
לאחר מכן, החלפת היחסים המרכיבים המתאימים לתווך המוליך למחצה אניזוטרופי מחוזק בסיבים במשוואות הנ״ל מניבה את משוואות התנועה המקושרות, במונחים של רכיבי ההזזה, שדה הטמפרטורה וצפיפות הנשאים32
. (33)
. (34)
לבסוף, באמצעות הביטוי של שדה ההפרעה המגנטית שניתן קודם במשוואה (15), והחלפתו במשוואות הנ״ל, ניתן לכתוב את משוואות התנועה השולטות בצורתן הסופית המקושרת כך
. (35)
. (36)
משוואות אלו חושפות את ההשפעה המשולבת של השפעות סיבובי ושדה מגנטי על התגובה התרמואלסטית של התווך. המונחים הסיבוביים מתחשבים הן בקוריוליס והן בתרומה הצנטריפוגלית, בעוד שהשדה המגנטי מוסיף קישור אלקטרומגנטי נוסף ומשנה את ההתנהגות הדינמית של המערכת. כתוצאה מכך, המשוואות השולטות יוצרות מסגרת מאוחדת לניתוח התפשטות גלים ואינטראקציות מולטיפיזיקליות במוליכים למחצה מחוזקים בסיבים מגנטו-פוטו-תרמואלסטיים מסתובבים. בגירוי אופטי, ההתנהגות התרמית של התווך המוליך למחצה מושפעת משמעותית מהאינטראקציה בין הולכת חום, העברת נשאים ועיוות מכני, מה שמוביל לתהליך תרמו-פוטואלסטי מחובר חזק. בניגוד למודל הולכת החום הקלאסי, התפלגות הטמפרטורה בחומרים למחצה מושפעת לא רק מדיפוזיה תרמית אלא גם משילוב נשאים וקישור תרמואלסטי. כתוצאה מכך, משוואת הולכת החום המוכללת עבור התווך המוליך למחצה מחוזק בסיבים אניזוטרופיים ניתנת לביטוי כך:7˒23.
. (37)
המשוואה לעיל מראה בבירור שהשדה התרמי בתוך המדיום המסתובב של מוליכים למחצה מגנטי-פוטו-תרמואלסטי נשלט על ידי ההשפעה המשולבת של הולכת חום אניזוטרופית, תהליכי רקומבינציית נושאים ואינטראקציות תרמואלסטיות. המונח
מתאר את האנרגיה התרמית הנוצרת כתוצאה מרקומבינציה של נשאים תחת עירור אופטי, בעוד שמונחי הקישור הכוללים
ומציינים
את השפעת העיוות המכני התלוי בזמן על התגובה התרמית של התווך. כתוצאה מכך, שדה הטמפרטורה קשור חזק הן לצפיפות הנשא והן לשדה האלסטי, שמשחק תפקיד חשוב במאפייני ההתפשטות של גלים תרמואלסטיים בחומרים למחצה מחוזקים בסיבים. בניסוח הנוכחי, התפתחות ריכוז הנשא N(x, y, t) בתוך התווך המוליך למחצה נשלטת על ידי ההשפעות המשולבות של דיפוזיית נשאי, תהליכי רקומבינציה והפעלה תרמית הנוצרת על ידי עירור אופטי. בהתאם לכך, משוואת העברת הנשא המתארת את דינמיקת הנשא הלא-שיווי משקל ניתנת לכתיבה כך 4,23
. (38)
כאן, DE מציין את מקדם הדיפוזיה של הנשא, בעוד
שהוא מייצג את אופרטור לפלסיאן הדו-ממדי במישור x-y . המונח
מתייחס לאפקט הרקומבינציה של הנשא הקשור לחיי הנשא τ. יתרה מזאת, κ הוא פרמטר הצימוד התרמו-נשא המוגדר על ידי
, כאשר N0 מציין את ריכוז נשאי שיווי המשקל. מונח הקישור κT מתאר את השפעת שדה הטמפרטורה על יצירת נשאים עודפים בתוך תווך המוליך למחצה. המשוואה שלמעלה מראה כי ריכוז הנשא קשור בחוזקה לשדה התרמי דרך מנגנוני יצירת נשאים המופעלים תרמית. כתוצאה מכך, דינמיקת הנושאים תלויה מאוד הן בהשפעות דיפוזיה תרמית והן בהשפעות רקומבינציה, המשפיעות משמעותית על התגובה הפוטו-תרמואלסטית המקושרת של התווך המוליך למחצה אניזוטרופי מחוזק בסיבים מסתובבים. המשוואות השולטות והניסוח המתמטי של מערכת המוליך למחצה המשולבת מגנטו-פוטו-תרמואלסטי הוקמו כעת לחלוטין. הפרמטרים הפיזיקליים והחומריים של מדיום הסיליקון מסוכמים בטבלה 2, יחד עם הערכים המספריים, היחידות וההפניות המתאימות. פרמטרים אלו משמשים לאחר מכן בחישובים נומריים ובהליך אי-ממדי.
ניסוח לא-ממדי של מודל המוליך למחצה מחוזק סיבים מגנטו-פוטו-תרמואלסטי מסתובב
כדי לפשט את המשוואות השולטות ולקבל ייצוג מתמטי קומפקטי של מערכת המגנטו-פוטו-תרמואלסטית המסתובבת המקושרת, הוכנסו סולמות מאפיינים מתאימים לצורך אי-ממדיזציה של המשתנים הפיזיקליים. הליך אי-ממדיות זה מפחית את מספר הפרמטרים השולטים בחומר ומקל על הטיפול האנליטי והנומרי במשוואות המקושרות. הכמויות האופייניות שנבחרו נבחרות בהתאמה לתכונות התרמואלסטיות, האלקטרומגנטיות, הסיבוביות והנשיאות של התווך המוליך למחצה16,21. בהתאם לכך, מוצגים המשתנים חסרי הממד הבאים:
,
,
,
,
, 



, , . 
כאן, C T מציין את מהירות הגל האלסטי האופיינית, בעוד ש-t* מייצג את זמן ההרפיה התרמית האופייני לתהליך התרמואלסטי המקושר. יתרה מזאת, הפרמטר מגדיר את הפרמטר
הסיבובי הלא-ממדי שמאפיין את השפעת תנועת הסיבוב על ההתנהגות הדינמית של התווך. החלפת הכמויות חסרות הממד לעיל במשוואות השלטון שנגזרו קודם לכן ממירה את המערכת המקושרת לצורה מנורמלת. טרנספורמציה זו מפשטת משמעותית את המבנה המתמטי של המשוואות ומספקת מסגרת מתאימה לחקירת ההשפעות המשולבות של שדה מגנטי, עירור אופטי, סיבוב, אניזוטרופיה ואינטראקציות נשא והעברה. לשם הפשטות, הסימון הראשוני הקשור לכמויות חסרות ממד מושמט בניתוח הבא. בהתאם לכך, ניתן לכתוב את המשוואות השולטות של מערכת מגנטו-פוטו-תרמואלסטית המסתובבת המקושרת בצורה הלא-ממדיתהבאה 16,20:
, (39)
, (40)
, (41)
. (42)
רכיבי המתח הלא-ממדיים המתאימים של התווך הסיבי האניזוטרופי המסתובב מתקבלים כך:
, (43)
, (44)
. (45)
המקדמים הלא-ממדיים a i (i = 1,2,...,18) מייצגים שילובים של הפרמטרים הפיזיקליים, התרמיים, האלקטרומגנטיים, הנשאים, והסיבובים של התווך המוליך למחצה המקושר. מקדמים אלו מאפיינים את השפעת האניזוטרופיה, חיזוק הסיבים, השדה המגנטי, הקישור התרמואלסטי, העברת נשאים ותנועת סיבוב על ההתנהגות הכוללת של המערכת. כתוצאה מכך, משוואות השלטון הלא-ממדיות שהתקבלו מספקות מודל מתמטי קומפקטי ויעיל לניתוח תופעות התפשטות גלים מקושרים ואינטראקציות מולטיפיזיות במדיות מוליכים למחצה מחוזקות בסיבים מגושים-פוטו-תרמואלסטיים מסתובבים. הפרמטרים הלא-ממדיים ai, γi ו-δi הוצגו כדי לייצג שילובים קומפקטיים של תכונות פיזיקליות וחומריות השולטות בהתנהגות המגנטו-פוטו-תרמואלסטית המסתובבת המקושרת של המדיום המוליך למחצה אניזוטרופי מחוזק בסיבים. כל מקדם משקף מנגנון אינטראקציה ספציפי בתוך המערכת המקושרת ומספק תובנות לגבי ההשפעה היחסית של התהליכים הפיזיקליים הבסיסיים. לשם בהירות, הפרמטרים חסרי הממד וההגדרות המתאימות שלהם מסוכמים בטבלה 1.
טבלה 1: הגדרות ופרשנויות פיזיקליות של הפרמטרים חסרי הממד המשמשים בנוסח הנוכחי. הטבלה מסכמת את הפרמטרים חסרי הממד המופיעים במשוואות השולטות יחד עם משמעותם הפיזיקלית ותפקידם בתיאור האינטראקציות התרמואלסטיות, האלקטרומגנטיות, צפיפות הנשאים והסיבוביות המקושרות. אנא לחצו כאן להורדת הטבלה הזו.
פתרון אנליטי באמצעות טכניקת המצב הנורמלי
כדי לגזור את הפתרון האנליטי של מערכת המגנטו-פוטו-תרמואלסטית המסתובבת המקושרת, נעשה שימוש בטכניקת המצב הנורמלי. שיטה זו משמשת רבות בתיאוריות תרמואלסטיות ומוליכים למחצה מוכללות בשל יעילותה בהמרת משוואות דיפרנציאליות חלקיות מקושרות למערכת מצומצמת של משוואות דיפרנציאליות רגילות. גישה כזו שימושית במיוחד בניתוח התפשטות גלים, החלשה ואינטראקציות מולטיפיזיקליות במדיות מוליכים למחצה אניזוטרופיות. לאחר ניתוח המצב הנורמלי, מניחים שכל גדלי השדה הפיזיקליים משתנים הרמונית ביחס לזמן ולקואורדינטה המרחבית הרוחבית . בהתאם לכך, שדה הטמפרטורה, צפיפות הנשאים, רכיבי ההזזה וכמויות המאמץ מיוצגים בצורה האקספוננציאלית24,27
. (46)
כאן, ω מציין את פרמטר התדר המרוכב השולט בשוני הזמן של השדות הפיזיקליים. החלק הממשי של ω קשור להחלשה (או גדילה זמנית) של משרעת הגל, בעוד שהחלק המדומה מייצג את ההתנהגות התנודתית של מצב ההעברה. פרשנויות אלו תואמות את ניתוח המצב הנורמלי הקונבנציונלי שאומץ במחקר הנוכחי, בעוד ש-a מייצג את מספר הגל הקשור לשונות המרחבית לאורך כיוון y. הכמויות
, ו
- מתאימות למשרעת השדה התלויות רק בקואורדינטת המרחבית x. בהחלפת ייצוגי המצב הנורמלי לעיל במשוואות השלטון הלא-ממדיות שהושגו קודם לכן ופישוט הביטויים המתקבלים, מערכת הדיפרנציאלית החלקית המקורית המקושרת מומרת לקבוצת משוואות דיפרנציאליות רגילות ביחס לקואורדינטה המרחבית x. כתוצאה מכך, המשוואות השולטות בתחום הטרנספורמטיבי מקבלות את הצורה הבאה:
, (47)
, (48)
, (49)
. (50)
יתרה מזאת, רכיבי המאמץ המוטרנספורמציה המתאימים מתקבלים כ-
, (51)
, (52)
. (53)
כאן,
, מציין את האופרטור הדיפרנציאלי ביחס לקואורדינטה המרחבית . המערכת המומרת שמתקבלת מהווה את הבסיס המתמטי לבניית משוואת האופייניות ולהפקת הפתרון האנליטי המלא של בעיית המגנטו-פוטו-תרמואלסטי המסתובבת המקושרת. המקדמים המופיעים במשוואות המומרת מוגדרים כך:
,
,
,
,
,
, 
,
. 


מקדמים אלו מכילים את התרומה המשולבת של אלסטיות אניזוטרופית, אינטראקציה עם שדה מגנטי, קישור תרמי, העברת נשאים ואפקטים סיבובית. לכן, המערכת המומרת מספקת ייצוג קומפקטי המתאים לקבלת שורשי האופייניים ולחקירת התנהגות התפשטות הגל המקושר בתוך התווך הסמי-מוליך המחוזק הסיבים המסתובב.
קובץ משלים 1: פתרון אנליטי באמצעות צורת מטריצה. קובץ זה מכיל את ניסוח המטריצות המפורט, את תהליך פתרון הערך העצמי, את גזירת משוואות האופייניות, ואת השלבים האנליטיים הביניים המשמשים להשגת הפתרון הכללי של מודל מגנטו-פוטו-תרמואלסטי המקושר. אנא לחצו כאן להורדת הקובץ הזה.
ניסוח המטריצות המפורט, הליך פתרון הערך העצמי והגזרת משוואות האופייניות מסופקים בקובץ המשלים 1.
תנאי גבול וקביעת הקבועים הלא ידועים
להשלמת הניסוח האנליטי, הפתרונות הכלליים שהושגו הוחלפו בתנאי הגבול שנקבעו על פני השטח x = 0. החלפה זו יצרה מערכת אלגברית מקושרת הכוללת את קבועי
המשרעת הלא ידועים . כל דרישת גבול הקשורה לטמפרטורה, צפיפות נושאים, הזזה מכנית ואילוצי מאמץ הובאה במונחים של מצבי העצמי המותרים, מה שיצר קבוצה ליניארית של משוואות הקשורה בין המקדמים
. לנוחות, המערכת האלגברית שהתקבלה נכתבה מחדש בצורת מטריצה קומפקטית כ-BC = D, כאשר B מציין את מטריצת המקדמים שנבנתה מרכיבי וקטור העצמי שהוערכו על פני
השפל, מייצג את וקטור הקבועים הלא ידועים, ו-D מתאים לוקטור שנוצר מתנאי הגבול שהוטלו, כולל פרמטר העומס התרמי θ0, תקופת הגירוי של הנשא N0, ותנאי ההעתקה שנקבעו. לאחר הערכת הקבועים הללו, הם הוחזרו לביטויים הכלליים של משתני השדה כדי לקבל את הפתרונות האנליטיים המלאים. ביטויים אלו שימשו לאחר מכן בחישובים מספריים ובהדמיה גרפית של שדות התרמואלסטיים, צפיפות הנשאים וההזזה בתוך המדיום המסתובב של המוליך למחצה אניזוטרופי מחוזק על ידי סיבים. תנאי הגבול שאומצו מייצגים משטח מוליך למחצה מואר אופטית הנתון לגירורים תרמיים ונשיאים בו-זמנית. תנאי הטמפרטורה שנקבע מדמים את העומס התרמי שנוצר על ידי השדה האופטי הנכנס, בעוד שתנאי צפיפות הנשא מתחשב בנשאים העודפים שנוצרים על ידי האור האור האור. בנוסף, מגבלת ההעתקה הרוחבית מייצגת כליאה מכנית של פני השטח בכיוון -, בעוד שתנאי מתח הגזירה המתאפס מתייחס לגבול חופשי ממשיכית ללא אחיזה. כתוצאה מכך, תנאי הגבול התרמיים, האלקטרוניים והמכניים המעורבים שנבחרו מספקים ייצוג פיזיקלי עקבי של אינטראקציות פוטו-תרמואלסטיות מקושרות על פני השטח של המוליך למחצה ומספקים את המגבלות הנדרשות לקביעת הקבועים הלא ידועים של התמיסה. תנאי הגבול שהוטלו הם כדלקמן:
הגבלת טמפרטורה:
. (78)
תנאי זה מייצג טמפרטורת פני שטח משתנה הרמונית הנגרמת על ידי חימום אופטי תקופתי. הוא פועל כגירוי תרמי ראשוני שמניע את תהליכי ההובלה התרמואלסטית והנשיאה המקושרים בתוך התווך. המשרעת θ0 מאפיינת את עוצמת העומס התרמי המופעל.
מגבלת צפיפות הנשאים:
. (79)
תנאי גבול זה מתאר את צפיפות הנשאים שנוצרת על ידי הפוט הנובעת מהאורה האופטית. הוא משקף את ההתרגשות האלקטרונית הנובעת מספיעת פוטונים ואת המודולציה ההרמונית שלה התואמת את השדה האופטי הנכנס.
מגבלת הזזה:
. (80)
תנאי זה מציין שהגבול מוגבל מכנית בכיוון הרוחבי. לכן, לא מתרחשת הזזה לאורך כיוון ה-על פני השטח.
מגבלת מאמץ גזירה:
. (81)
תנאי זה מתאים לגבול חופשי מאחיזה ביחס למאמץ גזירה. הוא מבטיח שאין כוחות משיקים לפעול על פני השטח, וזה תואם לגבול מכני חופשי בכיוון המשיק. בנוסף לתנאי הגבול ב-x = 0, הדרישה הפיזיקלית באינסוף הוטלה כך:
הבטחת פתרונות פיזיקליים חסומים בתוך התחום החצי-אינסופי. כדי לספק סקירה ברורה של ההליך האנליטי והחישובי שאומץ במחקר הנוכחי, השלבים העיקריים של מתודולוגיית הפתרון מסוכמים באיור 2.

איור 2: תרשים זרימה של הליך הפתרון האנליטי שאומץ במחקר הנוכחי, כולל ניסוח המשוואות השולטות, ניתוח מצבים נורמליים, פתרון ערכים עצמיים, יישום תנאי שפה והערכת משתני השדה הפיזיקליים. הדיאגרמה מסכמת את השלבים החישוביים העיקריים לקבלת הפתרון האנליטי ולאחר מכן להערכת התגובה המקושרת בין מגנטו-פוטו-תרמואלסטי באופן נומרי. אנא לחצו כאן כדי לצפות בגרסה מוגדלת של הדמות הזו.
הגישה מוגבלת. התחברו או התחילו תקופת ניסיון כדי לצפות בתוכן זה.
כדי לקבל את התוצאות המספריות שהוצגו במחקר זה, נבחר הסיליקון (Si) כחומר המוליך למחצה המייצג בשל יישומיו הרחבים ביישומים של מוליכים למחצה והנדסת תרמואלסטיות. התכונות הפיזיקליות, התרמיות, האלסטיות, האלקטרומגנטיות והחומרים הקשורים לנשאים שנעשה בהם שימוש בחישובים מסוכמות בטבלה 2. קבועים חומריים אלו הוחלפו במשוואות השולטות חסרות הממדים כדי לייצר את התוצאות המספריות שנדונות בחלקים הבאים. בנוסף, החישובים המספריים בוצעו על פני התחום
הגישה מוגבלת. התחברו או התחילו תקופת ניסיון כדי לצפות בתוכן זה.
ניתוח קבילות, יציבות וגבולות ערכים עצמיים
קבילות הערכים העצמיים שהתקבלו נשלטת על ידי הדרישה הפיזיקלית שכל משתני השדה יישארו חסומים בתוך התחום החצי-אינסופי x ≥ 0. הפתרון העצמי הכללי מבוטא במונחים של מצבים אקספוננציאליים מהצורה
. כתוצאה מכך, ההתנהגות האסימפטוטית של הפתרון תלויה ישירות בחלק הממשי של הערכים
הגישה מוגבלת. התחברו או התחילו תקופת ניסיון כדי לצפות בתוכן זה.
המחברים מצהירים שאין להם אינטרסים מתחרים.
המחברים מביעים את הערכתם לדיקן המחקר והלימודים המתקדמים באוניברסיטת המלך ח'אלד על מימון העבודה הזו באמצעות תוכנית קבוצות מחקר גדולות תחת מענק מספר RGP2/230/47.
הגישה מוגבלת. התחברו או התחילו תקופת ניסיון כדי לצפות בתוכן זה.
| שם | חברה | מספר קטלוג | הערות |
|---|---|---|---|
| תוכנת חישוב (ניתוח סימבולי ומספרי) | Wolfram Research | נעשה שימוש ב-Wolfram Mathematica (גרסה 12.0) | נועד לחישוב ערכים עצמיים, יישום פתרונות אנליטיים והערכה מספרית |
| מרחב פרמטרים של חומרים (מאפייני מוליכים למחצה מסוג סיליקון) | מקורות ספרותיים שונים | לא חל | קבועים פיזיקליים (אלסטיים, תרמיים, קשורים לנשאים) המשמשים בחישובים (טבלה 1) |
| מחשב אישי/תחנת עבודה | HP | לא חל | החישובים בוצעו במחשב אישי סטנדרטי שהפעיל מערכת הפעלה Windows עם זיכרון מספק להדמיה מספרית |
| עורך משוואות | Microsoft Word ו-MathType | לא חל | נועד לעיצוב והצגת ביטויים מתמטיים בכתב היד |
| תוכנת ניהול הפניות | Elsevier | לא חל | נועד לניהול הפניות ועיצוב ציטוטים (בסגנון Vancouver) |
הגישה מוגבלת. התחברו או התחילו תקופת ניסיון כדי לצפות בתוכן זה.