אנרגיה חופשית של הלמהולץ

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

==הגדרה==
האנרגיה החופשית של הלמהולץ מוגדרת ע"י:
:<math> F = U-TS </math>
כאשר
*<math>F </math> האנרגיה החופשית של הלמהולץ (לעיתים מסומנת באות <math>A </math>)
*<math>U </math> האנרגיה הפנימית של המערכת
*<math>T </math> הטמפרטורה
*<math>S </math> האנטרופיה

==פיתוח==
האנרגיה החופשית של הלמהולץ מתקבלת על ידי ביצוע התמרת לז'נדר לאנרגיה הפנימית:

:<math>F(T,V, \left\{ N_i \right\})=U(S,V,\left\{ N_i \right\})-\left( {\partial U \over \partial S} \right) _{\left\{ N_i \right\},V} S=U-TS</math>

כאשר <math>N_i</math> מספר החלקיקים מסוג <math>i</math> במערכת ו- <math>\left( {\partial U \over \partial S} \right) _{\left\{ N_i \right\},V}=T</math>.
כלומר, האנרגיה החופשית של הלמהולץ היא פונקציית מצב בה מספר החלקיקים, הנפח והטמפרטורה (במקום האנטרופיה עבור האנרגיה הפנימית) הם המשתנים הבלתי תלויים.

===ביטוי דיפרנציאלי===
באמצעות הביטוי הדיפרנציאלי לאנרגיה הפנימית נקבל:
:<math>dF=dU-d(TS)=TdS-PdV+\sum_{i} \mu_i dN_i - d(TS)=-SdT-PdV+\sum_{i} \mu_i dN_i</math>
כאשר <math>\mu_i</math> הפוטנציאל הכימי.

ממשוואה זו ניתן לקבל ביטויים עבור המשתנים התרמודינמיים הבלתי נשלטים:
{| border="0" cellpadding="4" style="margin: 0 0 1em 1em"
|-
|האנטרופיה –
|align="left"|<math>S = -\left( \frac{\partial F}{\partial T}\right)_{V,\left\{ N_i \right\}} </math>,
|-
|הלחץ –
|align="left"| <math>P =- \left( \frac{\partial F}{\partial V}\right)_{T,\left\{ N_i \right\}} </math>
|-
|והפוטנציאל הכימי –
|align="left"|<math>\mu_i = \left( \frac{\partial F}{\partial N_i}\right)_{T,V,\left\{ N_{j\ne i} \right\}} </math>.
|}

הביטוי שהתקבל עבור הלחץ ניתן לרשום ע"י
:<math>P=-\left( \frac{\partial U}{\partial V}\right)_{T,\left\{ N_i \right\}}+T\left( \frac{\partial S}{\partial V}\right)_{T,\left\{ N_i \right\}}</math>
כאשר הביטוי השמאלי מסמל את "לחץ האנרגיה" והביטוי הימני את "לחץ האנטרופיה".
לחץ האנרגיה <math>-\left( \frac{\partial U}{\partial V}\right)_{T,\left\{ N_i \right\}}</math> דומיננטי ברב החומרים המוצקים ואילו <math>T\left( \frac{\partial S}{\partial V}\right)_{T,\left\{ N_i \right\}}</math> דומיננטי בגזים ובפולימרים אלסטיים (למשל גומי). משמעות הביטוי היא שכאשר הטמפרטורה קבועה, הלחץ כבר לא תלוי רק בביטוי <math>-{\partial U \over \partial V}</math> משום שהאנטרופיה יכולה להשתנות כתוצאה מהשינוי בנפח, אפילו אם האנרגיה הפנימית בלתי תלויה בנפח (למשל גז אידיאלי בטמפרטורה קבועה). 

ניתן לקבל את קשרי מקסוול הבאים :

:<math>\left( {\partial S \over \partial V} \right) _{\left\{ N_i \right\},T}=\left( {\partial P \over \partial T} \right) _{\left\{ N_i \right\},V}</math>
:<math>\left( {\partial S \over \partial N_i} \right) _{T,V,\left\{ N_{l\ne i} \right\}}=-\left( {\partial \mu_i \over \partial T} \right) _{\left\{ N_i \right\},V}</math>
:<math>\left( {\partial P \over \partial N_i} \right) _{T,V,\left\{ N_{l\ne i} \right\}}=-\left( {\partial \mu_i \over \partial V} \right) _{\left\{ N_i \right\},V}</math>
:<math>
\left( {\partial \mu_j \over \partial N_i} \right) _{T,V,\left\{ N_{l\ne i} \right\}}=\left( {\partial \mu_i \over \partial N_j} \right) _{T,V,\left\{ N_{l\ne j} \right\}}</math>

==עיקרון מינימום האנרגיה החופשית ועקרון מקסימום העבודה==
נוכל לתאר מערכת בטמפ' קבועה באמצעות מערכת המצומדת לאמבט חום בנפח קבוע ובעל טמפ' שאינה משתנה כאשר המערכת הכוללת (מערכת + אמבט) סגורה. כלומר, קיים מעבר חום בין המערכת לאמבט והאמבט גדול מאוד ביחס למערכת כך שכל שינוי במערכת לא משפיע על טמפרטורת האמבט. נתאר תהליך כך שהמערכת הכוללת נמצאת בשיווי משקל לפני ואחרי התהליך, לכן המשתנים התרמודינמיים לפני ואחרי התהליך מוגדרים היטב. כיוון שהמערכת הכוללת סגורה מתקיים החוק הראשון של התרמודינמיקה
:<math>\Delta U_{resevoir}+\Delta U+W=0</math>

כאשר <math>\Delta U_{resevoir}</math> השינוי באנרגיית האמבט, <math>\Delta U</math> השינוי באנרגיית המערכת ו- <math>W</math> העבודה שנעשתה על ידי המערכת. האמבט בנפח קבוע ולכן אינו מבצע עבודה. כלומר השינוי באנרגיית האמבט שווה לחום המועבר למערכת
:<math>Q_{resevoir}=\Delta U_{resevoir}=-(\Delta U+W)</math>

טמפרטורת האמבט <math>T</math> קבועה ולכן האמבט נמצא בשיווי משקל ובעל אנטרופיה מוגדרת היטב
:<math>\Delta S_{resevoir}={Q_{resevoir} \over T}=-{\Delta U+W \over T}</math>

המערכת בשיווי משקל בתחילת ובסוף התהליך ולכן בטמפרטורה <math>T</math> והשינוי באנטרופיה הכוללת הוא
:<math>\Delta S_{resevoir}+\Delta S=-{\Delta U+W-T \Delta S \over T}=-{\Delta F +W \over T}
</math>
המערכת הכוללת סגורה ולכן לפי החוק השני של התרמודינמיקה, האנטרופיה הכוללת יכולה רק לגדול. לכן
:<math>W \le -\Delta F</math>

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

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

הוכחה

עבור המערכת שתוארה, נסתכל על האנטרופיה של המערכת הכוללת
:<math>\sigma _{total}=\sigma _R+\sigma _H=\sigma _R (U_{total}-U_H)+\sigma _H</math>

כאשר <math>\sigma _R</math> אנטרופיית האמבט, <math>\sigma _H</math> אנטרופיית המערכת ו- <math>U_H</math> האנרגיה הפנימית של המערכת. האנטרופיה מסומנת כאן ביחידות של אנרגיה: <math>S=k_B\sigma</math> (<math>k_B</math> קבוע בולצמן). אנו מניחים <math>U_{total} \gg U_H</math> ולכן <math>U_R=U_{total}-U_H \approx U_{total}</math> ונוכל לקרב את אנטרופיית האמבט סביב <math>U_{total}</math>
:<math>\sigma_{total} \cong \sigma_R(U_{total})-{\partial \sigma_R \over \partial U_R}U_H+{1 \over 2}{\partial ^2 \sigma_R \over \partial U_R^2 }U_H^2+ \cdots </math>

:<math>{\partial \sigma_R \over \partial U_R}={1 \over \tau}</math>
כאשר <math>\tau =k_B T</math> טמפרטורת האמבט ביחידות של אנרגיה. <math>\tau</math> קבועה ולכן כל הנגזרות מסדר הגבוה מ-1 מתאפסות ומתקבל

:<math>\sigma_{total} \cong \sigma_R(U_{total})-{1 \over \tau}U_H=\sigma_R(U_{total})-{U_H-\tau \sigma_H \over \tau}=\sigma_R(U_{total})-{F_H \over \tau}</math>

המערכת הכוללת סגורה ולכן קונפיגורציית שיווי המשקל שלה (כלומר המסתברת ביותר) מתקבלת עבור <math>\sigma_{total}</math> מקסימלי ולכן <math>{F_H \over \tau}</math> מינימלי. מינימום זה מושג על ידי שינוי כל הפרמטריים הפנימיים של המערכת החופשיים להשתנות וגם על ידי שינויים של האנרגיה הפנימית <math>U_H</math>. אם מוצאים תחילה את המינימום של <math>{F_H \over \tau}</math> כאשר <math>U_H</math> קבוע, אז את המעבר האחרון של המינימיזציה ניתן לכתיבה ע"י

:<math>{1 \over \tau}F(\tau ,V)= \min_U{[{U \over \tau}-\sigma_H(U,V)]}</math>

<math>\sigma_H(U,V)</math> בנוסחא זו היא האנטרופיה המקסימלית כאשר <math>U,V</math> נתונים.

==במכניקה סטטיסטית==
האנרגיה החופשית של הלמהולץ שימושית בעבודה עם צבר קאנוני, בה הטמפ', הנפח ומספר החלקיקים קבוע, והיא פונקציית המצב האופיינית לצבר זה. בעזרת הביטויים של האנטרופיה והאנרגיה הפנימית ניתן להביע את האנרגיה החופשית בעזרת פונקציית החלוקה <math>Z </math>:
:<math>F=-Tk_B \ln Z </math>.

====הוכחת הקשר בין האנרגיה החופשית של הלמהולץ לפונקציית החלוקה====
במקרה של צבר קאנוני, הסיכוי למצוא את המערכת במצב מיקרוסקופי <math>r</math> עם אנרגיה <math>\epsilon_r</math> נתון ע"י:
:<math>P(\epsilon_r)={e^{-\beta \epsilon_r} \over Z}</math>
כאשר <math>\beta={1 \over k_B T}</math> (<math>k_B</math> קבוע בולצמן) ו- <math>Z</math> פונקציית החלוקה:
:<math>Z=\sum_r{e^{-\beta \epsilon_r}}</math>
והסכימה היא על כל המצבים המיקרוסקופיים הזמינים במערכת עם אנרגיה <math>\epsilon_r</math>.
האנרגיה הפנימית הממוצעת היא:
:<math>U \equiv \langle E \rangle
= \sum_r P_r E_r
= \sum_r \frac{e^{-\beta E_r} E_r}{Z}
= \sum_r \frac{-\frac{\partial}{\partial \beta} e^{-\beta E_r}}{Z}
= \frac{-\frac{\partial}{\partial \beta} \sum_r e^{-\beta E_r}}{Z}
= -\frac{\partial \log Z}{\partial \beta}</math>

בנוסף
:<math>S=-\left( {\partial F \over \partial T} \right) _{V,\left\{ N_i \right\}}
 = k\beta ^2 \left( {\partial F \over \partial \beta} \right) _{V,\left\{ N_i \right\}}</math>
ולכן:
:<math>F=U-TS
 = -{\partial log{Z} \over \partial \beta} - \beta \left( {\partial F \over \partial \beta} \right) _{V,\left\{ N_i \right\}}</math>
:<math>F+\beta \left( {\partial F \over \partial \beta} \right) _{V,\left\{ N_i \right\}} = -{\partial \log{Z} \over \partial \beta}</math>
:<math>\left( {\partial (F\beta) \over \partial \beta} \right) _{V,\left\{ N_i\right\}} = -{\partial \log{Z} \over \partial \beta}</math>
:<math>\beta F = -\log{Z} + \psi(N_i,V)</math>

כאשר <math>\psi(N_i,V)</math> פונקציה כלשהי הנובעת מהאינטגרציה. עבור <math>T\to 0</math>, נשתמש בהגדרת האנטרופיה במכניקה סטטיסטית: <math>S=k_B \log{g(n,U)}</math> כאשר <math>g(n,U)</math> פונקציית הריבוי ונקבל
:<math>\beta F(T=0)=\beta U(T=0)-{1 \over k_B} S = -\log{g_0}</math>
כאשר <math>g_0</math> פונקציית הריבוי של אנרגיית מצב היסוד. לכן:
:<math>-\log{Z} + \psi(N_i, V) = -\log{g_0}</math>
אבל <math>Z(T \to 0)=g_0</math> מהגדרתו ולכן <math>\psi (N_i, V) = 0</math> וקיבלנו:
:<math>F=-k_BT\log{Z}</math>

==הכללת האנרגיה החופשית עבור גופים==
באופן יותר כללי, הביטוי <math>PdV</math> תלוי במאמץ ובמעוות הפועל על הגוף ולכן נוכל לרשום את האנרגיה החופשית:
:<math>dF=-SdT+\sigma_ {i j}d\epsilon_{i j} + \mu_i dN_i </math>
כאשר <math>\sigma_{i j}</math> טנזור המאמץ, <math>\epsilon_{i j}</math> טנזור המעוות ואנו משתמשים בהסכם הסכימה של איינשטיין.
כעת, מהמשוואה הדיפרנציאלית ניתן למצוא ביטוי עבור טנזור המאמץ:
:<math>\sigma_{ij}=\left( \frac{\partial F}{\partial \epsilon_{ij}}\right)_{T,\left\{ N_l \right\}}</math>

חוק הוק

עבור גוף איזוטרופי הנמצא בטמפרטורה קבועה ועובר מעוות, נרצה למצוא את האנרגיה החופשית כתלות בטנזור המעוות. נניח כי המעוות קטן ולכן נוכל לרשום את האנרגיה החופשית כחזקות של טנזור המעוות. <math>F</math> סקלר ולכן גם הביטויים התלויים בטנזור המעוות חייבים להיות סקלרים. <math>\sigma_{ij}=\left( \frac{\partial F}{\partial \epsilon_{ij}}\right)_{T,\left\{ N_i \right\}}</math> ולכן לא ייתכנו ביטויים ליניאריים של טנזור המעוות. טנזור המעוות סימטרי ולכן ניתן לקבל באמצעותו 2 סקלרים בלתי תלויים: סכום ריבועי האלכסון הראשי <math>(\epsilon_{ii}^2)</math> וסכום הריבועים של כל הרכיבים <math>(\epsilon_{ij}^2)</math> 
:<math>F=F_0 + {1 \over 2}\lambda \epsilon_{ii}^2 + \mu \epsilon_{ij}^2</math>
כאשר <math>F_0</math> האנרגיה החופשית של הגוף בטמפרטורה הנתונה ללא מעוות ו- <math>\lambda, \mu</math> הם קבועי לאמה.
הסכום <math>(\epsilon_{ii}^2)</math> מייצג את המעוות הנפחי ולכן כאשר <math>(\epsilon_{ii}^2)=0</math> פועל מעוות גזירה בלבד. נוח לרשום את האנרגיה החופשית באמצעות ביטוי שקול:

:<math>F=\mu \left( \epsilon_{ij} - {1 \over 3} \delta_{ij} \epsilon_{ii}\right)^2 + {1 \over 2}K\epsilon_{ll}^2</math>
כאשר <math>K=\lambda + {2 \over 3}\mu</math> מודול הנפח ו- <math>\mu</math> מודול הגזירה.
נשים לב שכעת הביטוי הימני מייצג מעוות נפחי והביטוי השמאלי הוא טנזור בעל עקבה <math>0</math> ולכן מייצג מעוות גזירה.

בשיווי משקל האנרגיה החופשית מינימלית. אם לא פועלים כוחות חיצוניים על הגוף, ל- <math>F</math> כפונקציה של <math>\epsilon_{ij}</math> חיים להיות מינימום ב- <math>\epsilon_{ij} = 0</math>. לכן הביטוי של <math>F</math> כפונקציה של <math>\epsilon_{ij}</math> הוא פונקציה ריבועית חיובית. אם <math>\epsilon_{ll}=0</math>, רק הביטוי השמאלי בנוסחא נותר. באותו אופן, אם <math>\epsilon_{ij} = const * \delta_{ij}</math> רק הביטוי הימני נותר. לכן נקבל
:<math>K>0 \quad,\quad \mu>0</math>

נמצא את הדיפרנציאל של <math>F</math> (בטמפרטורה קבועה)
:<math>dF=K\epsilon_{ll} d\epsilon_{ll} +2\mu \left( \epsilon_{ij} - {1 \over 3} \epsilon_{ll} \delta_{ij} \right) 
d\left( \epsilon_{ij} - {1 \over 3} \epsilon_{ll} \delta_{ij} \right) </math>
:<math>dF= \left[ \epsilon_{ll} \delta_{ij} + 2\mu \left( \epsilon_{ij} - {1 \over 3} \epsilon_{ll} \delta_{ij} \right) \right] d\epsilon_{ij}</math>
מכאן שטנזור המאמץ הוא
:<math>\sigma_{ij} = \epsilon_{ll} \delta_{ij} + 2\mu \left( \epsilon_{ij} - {1 \over 3} \epsilon_{ll} \delta_{ij} \right)</math>
נשים לב ש <math>\epsilon_{ii}=\sigma_{ii} / 3K</math> ונקבל שטנזור המעוות כתלות בטנזור המאמץ הוא
:<math>\epsilon_{ij}= {\delta_{ij} \sigma_{ll} \over 9K} + {\sigma_{ij} - {1 \over 3} \delta_{ij} \sigma_{ll} \over 2\mu}</math>
מהביטוי עבור טנזור המאמץ נסיק שהמאמץ הנפחי תלוי רק בסכום <math>\sigma_{ii}</math> והקשר בין <math>\sigma_{ii}</math> ל- <math>\epsilon_{ii}</math> תלוי רק במודול הנפח <math>K</math>. בדחיסה הידרוסטטית <math>\sigma_{ij} = -P\delta_{ij}</math> ולכן:
:<math>\epsilon_{ii} = -{P \over K}</math>
כיוון שהמעוות קטן, כך גם <math>\epsilon_{ii},P</math> קטנים ונוכל לרשום באופן דיפרנציאלי:
:<math>{1 \over K} = -{1 \over V} \left( {\partial V \over \partial P} \right) _T</math>
כאשר הביטוי <math>1 / K</math> נקרא מקדם הדחיסות.

קיבלנו שטנזור המעוות תלוי ליניארית בטנזור המאמץ. זהו חוק הוק, התקף עבור מעוותים קטנים.

כיוון שהאנרגיה החופשית נתונה כפונקציה ריבועית, מתקיים:
:<math>\epsilon_{ij} {\partial F \over \partial \epsilon_{ij} } = 2F</math>
וכיוון ש - <math>\sigma_{ij} = {\partial F \over \partial u_{ij} }</math> נקבל
:<math>F={1 \over 2} \sigma_{ij} \epsilon_{ij}</math>
על ידי הצבה של <math>\sigma_{ij}</math> בנוסחא מהביטוי של <math>\epsilon_{ij} (\sigma_{ij})</math>, נקבל שגם
<math>F (\sigma_{ij})</math> פונקציה ריבועית ולכן באופן דומה <math>\sigma_{ij} {\partial F \over \partial \sigma_{ij} } = 2F</math> וע"י השוואה לנוסחא האחרונה נקבל:
:<math>\epsilon_{ij}=\left( \frac{\partial F}{\partial \sigma_{ij}}\right)_{T,\left\{ N_i \right\}}</math>
חשוב לציין שנוסחא זו נכונה רק כאשר חוק הוק תקף. בשונה מהביטוי 
:<math>\sigma_{ij}=\left( \frac{\partial F}{\partial \epsilon_{ij}}\right)_{T,\left\{ N_i \right\}}</math>
שנכון באופן כללי.

==לקריאה נוספת==
* Charles Kittel and Herbert Kroemer, Thermal Physics, San Francisco: W.H. Freeman and Company
* L.E. Reichl, A Modern Course in Statistical Physics, Wiley-Interscience Publication
* ד. ברגמן, פיזיקה תרמית, דיונון-אוניברסיטת תל אביב
* M. W. Zemansky and chard H. Dittman, Heat and Thermodynamics, McGraw-Hill

==קישורים חיצוניים==
* 
* 

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

קטגוריה:פוטנציאלים תרמודינמיים