משפט ליוביל (מכניקה המילטונית)

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

## שימור הנפח במרחב הפאזה
במערכת פיזיקאלית המתוארת על ידי s קואורדינטות מוכללות q_1, q_2, ..., q_s, מרחב הפאזה הוא מרחב בעל 2s ממדים - מימד אחד לכל קואורדינטה מוכללת ומימד אחד לכל תנע הצמוד לקואורדינטה מוכללת: p_i=\frac{\partial \mathcal{L}}{\partial \dot{q_i}}, כאשר \mathcal{L} הלגראנז'יאן של המערכת. כל נקודה במרחב הפאזה מייצגת מצב מוגדר מסוים של המערכת, וכאשר המערכת נעה הנקודה המתאימה לכל מצב מגדירה עקומה הנקראת מסלול הפאזה.

אלמנט נפח אינפיניטסימלי במרחב הפאזה ניתן לרשום כ-d\Gamma=\prod_{i=1}^s{dq_idp_i}, ולכן נפח מסוים במרחב הפאזה, היכול לייצג קונפיגורציות של מערכת מרובת חלקיקים (צבר) הוא \int{d\Gamma}. משפט ליוביל טוען ש**הנפח אינווריאנטי תחת טרנספורמציות קנוניות**. טרנספורמציה קנונית מעבירה את הקוארינטות והתנעים המוכללים לקואורדינטות ותנעים מוכללים חדשים (q_1,...,q_s;p_1,...,p_s;t)\mapsto(Q_1,...,Q_s;P_1,...,P_s;t) באופן שמשמר את משוואות המילטון. בפרט, מכיוון שהתקדמות בזמן - מעבר מהקואורדינטות q(t),p(t) לערכיהן לאחר זמן \tau: q(t+\tau),p(t+\tau) - היא טרנספורמציה קנונית, ניתן לראות ש**הנפח במרחב הפאזה נשמר בזמן**.

### הוכחה
רוצים להראות שהנפח במונחי הקואורדינטות הקודמות \Gamma_{q,p}=\int...\int{dq_1...dq_sdp_1...dp_s} זהה לנפח במונחי הקואורדינטות החדשות \Gamma_{Q,P}=\int...\int{dQ_1...dQ_sdP_1...dP_s}. ניתן להעביר את האינטגרל השני להיות אינטגרל במונחי הקואורדינטות הישנות בעזרת היעקוביאן:

\int{dQ_1...dQ_sdP_1...dP_s}=\int{Jdq_1...dq_sdp_1...dp_s}, ולכן נותר רק להראות שהיעקוביאן הוא 1.

נשתמש בכללים של יעקוביאנים כדי לקבל: J=\frac{\partial(Q_1,...,Q_s,P_1,...,P_s)}{\partial(q_1,...,q_s,p_1,...,p_s)}=\frac{\frac{\partial(Q_1,...,Q_s,P_1,...,P_s)}{\partial(q_1,...,q_s,P_1,...,P_s)}}{\frac{\partial(q_1,...,q_s,p_1,...,p_s)}{\partial(q_1,...,q_s,P_1,...,P_s)}}=\frac{(\frac{\partial(Q_1,...,Q_s)}{\partial(q_1,...,q_s)})_{P=const}}{(\frac{\partial(p_1,...,p_s)}{\partial(P_1,...,P_s)})_{q=const}}.

כל טרנספורמציה קנונית ניתן לרשום במונחי פונקציה יוצרת כלשהי, למשל \Phi(q_i, P_i, t) כך שהקואורדינטות האחרות מתקבלות ממנה: p_i = \frac{\partial\Phi}{\partial q_i}, Q_i=\frac{\partial\Phi}{\partial P_i}. לכן האיבר בדטרמיננטה שבמונה בשורה ה-i ובעמודה ה-k הוא \frac{\partial Q_i}{\partial q_k}=\frac{\partial}{\partial q_k}(\frac{\partial\Phi}{\partial P_i})=\frac{\partial^2\Phi}{\partial q_k P_i}. האיבר ה-i,k בדטרמיננטה שבמכנה הוא \frac{\partial p_i}{\partial P_k}=\frac{\partial}{\partial P_k}(\frac{\partial\Phi}{\partial q_i})=\frac{\partial^2\Phi}{\partial q_i \partial P_k}, כלומר היעקוביאן שבמכנה הוא שחלוף של זה שבמונה, ולכן הדטרמיננטות שלהם זהות. קיבלנו כי J=1 ולכן הנפחים שווים.

## משוואת ליוביל
עבור צבר של מצבים שונים של מערכת פיזיקלית מסוימת, נגדיר את פונקציית צפיפות ההסתברות במרחב הפאזה כך ש- \rho(q_1,...,q_s,p_1,...,p_s,t)dq_1...dq_sdp_1...dp_s היא ההסתברות למצוא את המערכת במצב מסוים שנמצא באלמנט הנפח d\Gamma=dq_1...dq_sdp_1...dp_s. מכיוון שהתקדמות בזמן היא טרנספורמציה קנונית, והנפח במרחב הפאזה קבוע תחת טרנספורמציות קנוניות, נקבל כי הנפח שתופסים מספר מצבים במרחב הפאזה, המתפתחים בזמן באופן קלאסי (על-פי משוואות המילטון) נשמר בזמן, ולכן **צפיפות ההסתברות לאורך מסלול ההתפתחות בזמן, גם היא קבועה בזמן**:

\frac{d\rho}{dt}=\{\rho,\mathcal{H}\}+\frac{\partial\rho}{\partial t}=0

כאשר \mathcal{H} ההמילטוניאן, \{f,g\}=\sum_{i=1}^s{\frac{\partial f}{\partial q_i}\frac{\partial g}{\partial p_i}-\frac{\partial g}{\partial q_i}\frac{\partial f}{\partial p_i}} סוגרי פואסון, והקשר \frac{df}{dt}=\{f,\mathcal{H}\}+\frac{\partial f}{\partial t} ידוע ממכינקה המילטונית.

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

\frac{\partial \rho}{\partial t}=-\{\rho,\mathcal{H}\}

למשוואה זו, הנכונה למכניקה הקלאסית הנשלטת על-ידי משוואות המילטון, יש מקבילה במכניקת הקוונטים, **משוואת פון-נוימן**:

\frac{\partial \rho}{\partial t}=-\frac{1}{i\hbar}[\rho, \mathcal{H}]

כאשר \rho הוא אופרטור הצפיפות.

### במערכות בשיווי משקל
עבור מערכות שנמצאות בשיווי משקל, צפיפות ההסתברות אינה תלויה מפורשות בזמן: \frac{\partial\rho}{\partial t}=0. במקרה זה, משוואת ליוביל הופכת להיות:

\{\rho, \mathcal{H}\}=0

פתרון כללי לכך הוא שצפיפות ההסתברות היא פונקציה של ההמילטוניאן:

\rho(q,p)=\rho(\mathcal{H(q,p)})

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

\rho(q,p)=\frac{\delta(\mathcal{H(p,q)}-U)}{\int{\delta(\mathcal{H(p,q)}-U)d\Gamma}}, וכאשר המצבים בדידים ההסתברות זהה לכל המצבים הזמינים: P_s=\frac{1}{\sum_s{1}}=\frac{1}{g(U)}, זוהי פונקציית ההסתברות בצבר המיקרוקנוני שמראה את ההנחה היסודית של המכניקה הסטטיסטית - בשיווי משקל במערכת סגורה כל המצבים הזמינים הם שווי הסתברות.

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

קטגוריה:מערכות דינמיות
קטגוריה:מכניקה
קטגוריה:משפטים פיזיקליים