תורת שדה ממוצע

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

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

הרעיון הופיע לראשונה בפיזיקה בעבודתו של פייר קירי[1] ופייר ווייס כדי לתאר מעברי פאזה[2]. גישות נוספות שאימצו את השיטה תרמו למחקר ויישומים במודלים אפידמיולוגיים[3], תורת התורים[4], ביצועי רשתות מחשבים גדולות ותורת המשחקים[5].

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

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

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

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

ניסוח פורמלי עריכה

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

 

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

 

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

 

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

 

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

 

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

 

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

 

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

 

כאשר השדה הממוצע נתון על ידי

 

דוגמאות עריכה

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

מודל איזינג עריכה

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

 

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

 

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

 

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

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

 

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

 

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

 

כאשר   הוא מספר האתרים בסריג ושוב השתמשנו בסימון  .

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

גביש נוזלי עריכה

גביש נוזלי הוא תיאור לאוסף של מצבי צבירה בין נוזל למוצק שמתרחשים בצבירים של מולקולות מוארכות. ניתן לתאר את מעבר הפאזה מנוזל איזוטרופי - בו כל המולקולות מצביעות בכיוונים אקראיים, לגביש נוזלי נמטי - בו המולקולות מצביעות לכיוון אחיד, באמצעות תורת לנדאו ובחירת פרמטר הסדר כשדה ממוצע. לשם הדוגמה נציג מודל של אלפרד זָאוּפֶּה ווילהלם מָאיֶר[7] לתיאור מעבר הפאזה איזוטרופי-נמטי, בו מניחים אינטראקציות בין המולקולות שפרופורציונליות לחתך הפעולה של שתי המולקולות, ומזניחים את אפקט הדחיה הסטרית התלוי בטמפרטורה.

 
אילוסטרציה של הסתדרות המולקולות בפאזה הנמטית

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

 
דרגות חופש במוקולה בגביש נוזלי

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

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

נניח כמובן כי   מנורמלת במשתנה האקראי  .

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

בחירה טבעית לפרמטר סדר כזה היא פונקציית הקומולנט השני:

 

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

 

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

 

כלומר הפרש האנטרופיות:

 

הפוטנציאל שעלינו לעשות לו מינימיזציה הוא האנתלפיה:

 

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

 

כאשר   הוא פרמטר נרמול של ההתפלגות.

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

הרחבה לשדות ממוצעים תלויים בזמן עריכה

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

ראו גם עריכה

הערות שוליים עריכה

  1. ^ Kadanoff, L. P. (2009).
  2. ^ Weiss, Pierre (1907).
  3. ^ Boudec, J. Y. L.; McDonald, D.; Mundinger, J. (2007).
  4. ^ Baccelli, F.; Karpelevich, F. I.; Kelbert, M. Y.; Puhalskii, A. A.; Rybko, A. N.; Suhov, Y. M. (1992).
  5. ^ Lasry, J. M.; Lions, P. L. (2007).
  6. ^ HE Stanley (1971).
  7. ^ Maier W. and Saupe A., "Eine einfache molekulare theorie des nematischen kristallinflussigen zustandes", Z. Naturforsch. A
  8. ^ Lars Onsager, "THE EFFECTS OF SHAPE ON THE INTERACTION OF COLLOIDAL PARTICLES", Annals of the New York Academy of Sciences