3 כחול 1 חום

מכירים את 3Blue1Brown?

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

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

 

 

 

דיריכלה שהופך לנוימן

מהי העתקת דיריכלה לנוימן?
זוהי העתקה שניתן להגדירה על מגוון רחב של אובייקטים, החל מיריעות ועד גרפים.
נציג דוגמה:
בקורס של משוואות דיפרנציאליות חלקיות לומדים על בעיות עם תנאי שפה מסוג דיריכלה וגם על תנאי שפה מסוג נוימן (אם לא למדתם קורס במד”ח, אתם יכולים לדלג על הפסקה הבאה ולקרוא את המשך התאור).
יהי \( \Omega \) תחום ב \( \mathbb{R}^2 \).  ותהי \( f \) פונקציה שמוגדרת על שפת \( \Omega \)  (כלומר, על  \( \partial\Omega \)).
כעת נמצא פונקציה הרמונית \( g \) שמזדהה עם \( f \) על  \( \partial\Omega \), כלומר \( g|_{\partial\Omega}\equiv f \).
נחשב את הנגזרת הנורמלית של   \( g \) על   \( \partial\Omega \) ונסמנה ב   \( \frac{\partial}{\partial \hat{n}}g|_{\partial\Omega} \).  זוהי גם כן פונקציה על \( \partial\Omega \).
ההעתקה ששולחת כל \( f \) שכזו ל \( \frac{\partial}{\partial \hat{n}}g|_{\partial\Omega} \), נקראת העתקת דיריכלה לנוימן (לפעמים גם קוראים לה אופרטור סטקלוב).

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

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

מודלים כמו-מחזוריים

גבישים כמו-מחזוריים (quasi-crystals) התגלו באופן ניסויי על-ידי פרופ’ דן שכטמן עוד ב 1982.
מאז נחקרו ונחקרים שפע של מודלים מתמטיים בהם הכמו-מחזוריות (קווזי-מחזוריות) באה לידי ביטוי.
אחד המודלים שעומד במרכז מחקרנו נקרא אופרטור שרדינגר כמו-מחזורי חד-ממדי…..
ומה זה בדיוק?
אגב, נתאר כאן את המודל המתמטי שבו מתמקד מחקרנו. אם תרצו לקרוא גם על הפיסיקה הרלונטית (שהיא מרתקת, אולם לא נעסוק בה ישירות בפרויקט) אתם מוזמנים לקפוץ לסוף הדף.

נתבונן בפונקציות מהשלמים לממשיים, כלומר ב \(   f\, : \mathbb{Z} \rightarrow \mathbb{R}   \).נגדיר את האופרטור \( H_{\alpha} \) שפועל על פונקציה \( f \) שכזו ומחזיר פונקציה מאותו סוג (שנסמנה \( H_{\alpha}\,f \)) באופן הבא: $$ [H_{\alpha}\,f](n) = f(n-1)+f(n+1)+V(n)f(n) $$ ומיהי \( V(n) \) המצויינת בהגדרת האופרטור?
\(   V\, : \mathbb{Z} \rightarrow \mathbb{R}   \) היא פונקציה מהשלמים לממשיים המכונה פוטנציאל.  בפונקציה זו טמונה הכמו-מחזוריות ובה נמצא כל העניין.דוגמה אחת לפונקציה שכזו היא
$$  1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0  $$ כלומר, רשמנו למעלה 49 ערכים מתוך אינסוף הערכים של \( V \).  הערכים שרשמנו הם בדיוק הערכים \( V(1), V(2), V(3), \ldots , V(49) \).
האם אתם רואים חוקיות?

סדרת הערכים מעלה מבוססת על המספר האי-רציונלי  \( \alpha:=\frac{\sqrt{5}-1}{2} \).
האם המספר הזה מזכיר לכם משהו?  (נסו להוסיף לו אחת ולחפשו ברשת).
נתאר איך ליצור את הפוטנציאל מעלה מתוך המספר \( \alpha:=\frac{\sqrt{5}-1}{2} \).
ראשית, ניצור את סדרת המספרים הבאה   \( z(n):=n\alpha\mod 1 \),
כלומר אנו מחשבים את החלק השברי (שלימין הנקודה העשרונית) של כל אחת מהכפולות השלמות של \( \alpha \).
מתוך הסדרה \( z(n) \), נוכל להגדיר את הפוטנציאל על-ידי:
$$V(n)=   \begin{cases}  1, \quad \text{  if   } \,\, 1-\alpha\leq z(n) < 1 \\ 0,  \quad \text{otherwise} \end{cases}  $$

שימו לב שהיה מעין סדר בהגדרה של סדרת הערכים שמהווה את \( z(n) \) אולם סדרת הערכים של הפוטנציאל \( V(n) \) אינה נראית כה מסודרת וודאי אינה מחזורית.
יחד עם זאת, היא בעלת תכונות יפות (שיעזרו לנו במחקר) ומכילה סדר ותבניות חוזרות שאינן גלויות ממבט ראשון ובשל כך היא מכונה כמו-מחזורית.  בהתאמה, האופרטור  \( H_{\alpha} \) שהגדרנו מעלה נקרא אופרטור כמו-מחזורי (קווזי-מחזורי).

אז מה חוקרים?
הבעיות המרתקות שעומדות בפנינו כרוכות בנסיון להבין את הספטקרום של האופרטור \( H_{\alpha} \).
לשם פשטות, ניתן לחשוב על הספקטרום של האופרטור בתור אוסף הערכים העצמיים (המוכללים) של האופרטור.
מסתבר שכאשר לוקחים את הערך האי-רציונלי \( \alpha:=\frac{\sqrt{5}-1}{2} \) מקבלים שהספקטרום של האופרטור \( H_{\alpha} \) הוא קבוצת נקודות פרקטלית, מעין קבוצת קנטור.
לעומת זאת, אם נבחר \( \alpha:=\frac{p}{q} \) רציונלי (כלומר, \( p \) ו \( q \) טבעיים), ונציב את \( \alpha \) מעלה בנוסחה של הפטנציאל \( V(n) \) נקבל שהפוטנציאל יהיה באמת מחזורי (ולא רק “כמו-מחזורי”). נסו לחשוב מה יהיה המחזור במקרה כזה.
אם נבחן את האופרטור \( H_{\alpha} \) עבור פוטנציאל מחזורי שכזה, נקבל שהפעם הספקטרום שלו הוא פשוט הרבה יותר – הספקטרום יהיה אוסף סופי של קטעים סגורים (למעשה יהיו בדיוק \( q \) קטעים).

כעת, נוכל לצייר את הספקטרומים עבור כל ערכי \( \alpha \) השונים בין \( 0  \)  ל \( 1 \) ולקבל את האיור הבא:

תמונה זו לקוחה מתוך מאמר של Ostlund ו Kim משנת 1985… וכמו שאתם רואים, יש פה מקום לשיפור 😉

איך קוראים את האיור?
הציר האנכי מסמן את ערכי \( \alpha \) השונים בין \( 0  \) ל \( 1 \) .
לכל ערך של \( \alpha \) מסומן הספקטרום של האופרטור \( H_{\alpha} \) המתאים.  תוכלו לנסות לשים לב שעבור ערכים של \( \alpha=\frac{p}{q} \) מקבלים שהספקטרום מורכב מ \( q \) קטעים שונים…
ומה עם ערכי \( \alpha \) אי-רציונליים?  שם נצפה לראות אוסף פרקטלי (שאינו בן מניה) של נקודות.  רק שקצת קשה לראותו באיור הזה 😉
כדי להבין מה קורה בערכים אי-רציונליים, נוכל להתקרב אליהם באמצעות הערכים הרציונליים ולמעשה נצטרך להעזר במושג שנקרא שברים משולבים.

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

ואם הכתוב מעלה מעניין אתכם ואתם רוצים ללמוד עוד על הנושא, אבל אין לכם כוח או זמן לעבודת המחקר, אני ממליץ על הספר הפופולרי למחצה הבא:
Butterfly in the Quantum World: The story of the most fascinating fractal

ומה לגבי הפיסיקה שמאחורי המודל?

האופרטור \( H_{\alpha} \) משמש לתאור התנהגותו של חלקיק קוונטי שמוגבל לנוע בתווך חד-ממדי מסוים. אופרטור זה נקרא גם המילטוניאן (וזה המקור לסימונו באות H) ויש לו משמעות של אנרגיה. שני המחוברים הראשונים בהגדרה של \( H_{\alpha} \), כלומר \( f(n-1)+f(n+1) \) מהווים את האנרגיה הקינטית והמחובר האחרון \( V(n)f(n) \) הינו האנרגיה הפוטנציאלית. בתור דוגמה, תוכלו לקרוא על הפיסיקה שמאחורי המודל במאמרים של פרופ’ אריק אקרמנס בפקולטה לפיסיקה בטכניון.

 

פונקציות עצמיות של הלפלסיאן על יריעות וגרפים

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

Quantum Graph Symbol
פונקציה עצמית של גרף קוונטי

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

הפרויקטים יכללו אנליזה מתמטית משולבת בתכנות נומרי.

פונקציה עצמית של רבוע
פונקציה עצמית של רבוע

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

Neumann Domains
פונקציה עצמית של טורוס דו-ממדי