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

תיאור כללי של משוואת המצב של גז לא אידיאלי
------------------------------------------

התיאור הבסיסי של גז לא אידיאלי הוא משוואת המצב הויריאלית (virial equation of state), הנתונה באופן כללי על ידי:

{\displaystyle {\frac {Pv}{RT}}=1+{\frac {B\left(T\right)}{v}}+{\frac {C\left(T\right)}{v^{2}}}+\ldots }

כאשר v הוא הנפח המולרי, P הלחץ, T הטמפרטורה ו- R קבוע הגזים. הגודל {\displaystyle {\frac {Pv}{RT}}} נקרא מקדם הקומפרסביליות (compressibility factor). {\displaystyle B(T)} נקרא המקדם הויריאלי השני (second virial coefficient) ו-{\displaystyle C(T)} נקרא המקדם הויריאלי השלישי (third virial coefficient). ככל שצפיפות הגז גבוהה יותר, יש צורך במקדמים ויריאליים גבוהים יותר על מנת לתאר את התנהגות הגז במדויק.

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

בגז לא אידיאלי, האנרגיה של מצב מיקרוסקופי s נתונה על ידי:

{\displaystyle E\_{s}=E\_{s,id}+\Phi (r\_{1},r\_{2},...r\_{N})}  
כאשר {\displaystyle E\_{s,id}} היא האנרגיה שהייתה למצב מיקרוסקופי לו הגז היה אידיאלי, ו-\Phi  היא האנרגיה הפוטנציאלית (האנרגיה כתוצאה מאינטראקציה בין החלקיקים).

פונקציית החלוקה הקנונית נתונה על ידי:

{\displaystyle Q\ \left(T,V,N\right)=\sum \_{s}exp\ {\left(-{\frac {E\_{s,id}}{kT}}\right)}{\frac {1}{V^{N}}}\int \_{V}\exp {\left(-{\frac {\Phi }{kT}}\right)d^{3N}r}}

כאשר k הוא קבוע בולצמן. הביטוי {\displaystyle \sum \_{s}exp\ {\left(-{\frac {E\_{s,id}}{kT}}\right)}} הוא פונקציית החלוקה הקנונית עבור גז אידיאלי {\displaystyle Q\_{id}}, והביטוי {\displaystyle \int \_{V}\exp {\left(-{\frac {\Phi }{kT}}\right)d^{3N}r}} מכונה "**אינטגרל קונפיגורציה**", ויסומן ב- {\displaystyle Q\_{N}}, כך שמתקיים: {\displaystyle Q\ \left(T,V,N\right)=Q\_{id}{\frac {Q\_{N}}{V^{N}}}}. ניתן לראות שאינטגרל הקונפיגורציה ה- N מביא לידי ביטוי אינטראקציה בין N חלקיקים בגז.  
המשך הפיתוח יבוצע באנסמבל הגרנד קנוני. פונקציית החלוקה הגרנד קנונית מוגדרת כ:

{\displaystyle \Xi =\sum \_{N=0}^{\infty }\sum \_{s}{\lambda ^{N}\ e^{-\beta E\_{Ns}}}\ =\sum \_{N=0}^{\infty }{\lambda ^{N}\ Q(T,V,N)}}

כאשר \lambda  נקרא אקטיביות אבסולוטית (absolute activity), ושווה ל- {\displaystyle e^{\frac {\mu }{kT}}} . בנוסף, מתקיים:

{\displaystyle PV=kT\ln {\Xi }}

מהקשרים הללו מתקבל:

{\displaystyle {\frac {PV}{kT}}=\ln {\left[\sum \_{N=0}^{\infty }{\lambda ^{N}\ Q(T,V,N)}\right]}=\ln {\left[1+\sum \_{N=1}^{\infty }{\lambda ^{N}\ Q(T,V,N)}\right]}=\ln {\left[1+\sum \_{N=1}^{\infty }{\lambda ^{N}Q\_{id}(T,V,N){\frac {Q\_{N}}{V^{N}}}}\right]}=\ln {\left[1+\sum \_{N=1}^{\infty }{\lambda ^{N}{\frac {Z^{N}}{N!}}{\frac {Q\_{N}}{V^{N}}}}\right]}=\ln {\left[1+\sum \_{N=1}^{\infty }{{\frac {Q\_{N}}{N!}}z^{N}}\right]}}

כאשר השוויון השני נובע מכך שמתקיים {\displaystyle Q={\frac {Z^{N}}{N!}}}, כאשר {\displaystyle Z} היא פונקציית החלוקה הקנונית של חלקיק בודד, לכן האיבר בסכום הנ"ל עבור {\displaystyle N=0} הוא 1, ובשוויון האחרון סומן {\displaystyle z\equiv {\frac {\lambda Z}{V}}}.

באמצעות שימוש בפיתוח טיילור ללוגריתם מתקבל מהביטוי לעיל הקשר:

{\displaystyle {\frac {P}{kT}}={\frac {Q\_{1}}{V}}z+{\frac {1}{2!V}}\ \left(Q\_{2}-Q\_{1}^{2}\right)z^{2}+{\frac {1}{3!V}}\left(Q\_{3}-3Q\_{1}Q\_{2}+2Q\_{1}^{3}\right)z^{3}+\ldots \equiv \sum \_{l=1}^{\infty }{b\_{l}z^{l}}}

כאשר:

{\displaystyle b\_{1}\equiv {\frac {Q\_{1}}{V}}={\frac {1}{V}}\int \_{V}{d^{3}r}=1\ ;b\_{2}\equiv {\frac {Q\_{2}-Q\_{1}^{2}}{2!V}}\ ;b\_{3}\equiv {\frac {Q\_{3}-3Q\_{1}Q\_{2}+2Q\_{1}^{3}}{3!V}}}

המקדם {\displaystyle b\_{l}} מכונה "**אינטגרל מקבץ**" (cluster integral).

על מנת לבטא את המקדמים הויריאליים, יש לבטא את הקומפרסביליות כטור חזקות של {\displaystyle {\frac {1}{v}}}, כלומר יש למצוא קשר בין z ל- v. על מנת לעשות זאת, ניתן להשתמש בקשרים ידועים עבור האנסמבל הגרנד קנוני:{\displaystyle n(z)={\frac {N}{V}}={\frac {\lambda }{V}}\left({\frac {\partial \ln {\Xi }}{\partial \lambda }}\right)\_{T,V}={\frac {z}{V}}\left({\frac {\partial \ln {\Xi }}{\partial z}}\right)\_{T,V}={\frac {z}{kT}}\left({\frac {\partial P}{\partial z}}\right)\_{T,V}=\sum \_{l=1}^{\infty }{{lb}\_{l}z^{l}}}

כאשר במעבר שלפני האחרון נעשה שימוש בקשר {\displaystyle PV=kT\ln {\Xi }} שצוין קודם, ובמעבר האחרון נעשה שימוש בביטוי של {\displaystyle {\frac {P}{kT}}} כטור חזקות של z.

כדי למצוא את הקשר ההפוך - z כפונקציה של n - ניתן לכתוב את z כטור חזקות של n : {\displaystyle z=\sum \_{l=1}^{\infty }{a\_{l}n^{l}}}, ולמצוא את מקדמי הטור על ידי הצבה של {\displaystyle z(n)} הנ"ל בביטוי ל- {\displaystyle n(z)}. כך, ניתן לקבל:

{\displaystyle z=n-2b\_{2}+\left(8b\_{2}^{2}-3b\_{3}\right)n^{3}+\ldots }

מהצבה של קשר זה בהצגה של {\displaystyle {\frac {P}{kT}}} כטור חזקות של z ניתן לקבל את {\displaystyle {\frac {P}{kT}}} כטור חזקות של n:

{\displaystyle {\frac {P}{kT}}=n-b\_{2}n^{2}+\left(4b\_{2}^{2}-2b\_{3}\right)n^{3}+\ldots }

ומכאן ניתן לקבל, תוך שימוש בקשר {\displaystyle n={\frac {N\_{A}}{v}}} :

{\displaystyle {\frac {Pv}{RT}}=1-{\frac {b\_{2}N\_{A}}{v}}+{\frac {\left(4b\_{2}^{2}-2b\_{3}\right)N\_{A}^{2}}{v^{2}}}+\ldots }

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

| {\displaystyle B\left(T\right)=-b\_{2}N\_{A}\ ;C\left(T\right)=\left(4b\_{2}^{2}-2b\_{3}\right)N\_{A}^{2}} |
| --- |

איור 1: האינטראקציות שמבוטאות באמצעות המקדמים הויריאליים השונים - העיגולים מייצגים חלקיקים בגז, והחצים מייצגים את האינטראקציות ביניהם
מביטויים אלו ניתן, באמצעות הצבת הקשרים שנמצאו לעיל בין אינטגרלי מקבץ לאינטגרלי קונפיגורציה, לבטא את המקדמים הויריאליים באמצעות אינטגרלי קונפיגורציה. מאחר שאינטגרל המקבץ ה- i תלוי רק באינטגרלי קונפיגורציה {\displaystyle Q\_{j}} כאשר {\displaystyle j\leq i}, ניתן לראות שהמקדם הויריאלי הראשון לא מביא לידי ביטוי אינטראקציות בין חלקיקים, המקדם הויריאלי השני מביא לידי ביטוי אינטראקציות בין זוגות חלקיקים בלבד, המקדם הויריאלי השלישי מביא לידי ביטוי אינטראקציות בין שלשות וזוגות של חלקיקים וכך הלאה.

הערכת המקדמים הויריאליים
------------------------

על מנת להעריך את המקדמים הויריאליים, יש להעריך את אינטגרלי הקונפיגורציה {\displaystyle Q\_{N}}.

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

{\displaystyle \Phi \ \left(r\_{1},\ r\_{2},\ldots ,\ r\_{N}\right)=\sum \_{i=0}^{N}\sum \_{j=i}^{N}u\_{ij}}

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

{\displaystyle Q\_{N}=\int \_{V}\exp {\left(-{\frac {\Phi }{kT}}\right)d^{3N}r}=\int \_{V}\exp {\left(-{\frac {1}{kT}}\sum \_{i=0}^{N}\sum \_{j=i}^{N}u\_{ij}\right)d^{3N}r}=\int \_{V}\prod \_{i,j=0\ ;i<j}^{N}\exp {\left(-{\frac {u\_{ij}}{kT}}\right)}d^{3N}r}

לשם כך תוגדר **פונקציית המקבץ של מאייר** (Mayer cluster function):

{\displaystyle f\_{ij}=\exp {\left(-{\frac {u\_{ij}}{kT}}\right)}-1}

באמצעות פונקציה זו ניתן לכתוב את אינטגרל הקונפיגורציה כך:

{\displaystyle Q\_{N}=\int \_{V}\prod \_{i,j=0\ ;i<j}^{N}\left(1+f\_{ij}\right)d^{3N}r=\int \_{V}\left[1+\sum \_{i=0}^{N}\sum \_{j=i}^{N}f\_{ij}+\sum \_{i=0}^{N}\sum \_{j=i}^{N}\sum \_{k=0}^{N}\sum \_{l=k}^{N}{f\_{ij}f\_{kl}}+\ldots \right]d^{3N}r}

### המקדם הויריאלי השני

על מנת להעריך את המקדם הויריאלי השני, יש להעריך את האינטגרלים {\displaystyle Q\_{1},\ Q\_{2}}. האינטגרל {\displaystyle Q\_{1}} הוא מיידי:

{\displaystyle Q\_{1}=\int \_{V}{d^{3}r\_{1}=V}}

על מנת להעריך את {\displaystyle Q\_{2}}, נניח, כפי שצוין קודם לכן, שהאינטראקציות בין החלקיקים הן כוחות מרכזיים, כלומר {\displaystyle \Phi \left(r\_{1},r\_{2}\right)=\Phi (r\_{12})}. לכן, לפי ההגדרה של אינטגרל מקבץ, מתקיים:

{\displaystyle Q\_{2}=\iint \_{V}{\exp {\left[-{\frac {\Phi \left(r\_{12}\right)}{kT}}\right]}d^{3}r\_{1}d^{3}r\_{2}}}

{\displaystyle b\_{2}={\frac {Q\_{2}-Q\_{1}^{2}}{2!V}}={\frac {1}{2V}}\iint \_{V}{\left\{\exp {\left[-{\frac {\Phi \left(r\_{12}\right)}{kT}}\right]}-1\right\}d^{3}r\_{1}d^{3}r\_{2}}={\frac {1}{2V}}\iint \_{V}{f\_{12}d^{3}r\_{1}d^{3}r\_{2}}}

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

ניתן להחליף קואורדינטות: {\displaystyle r\_{1},\ r\_{2}\rightarrow \ r\_{12},r\_{cm}}, כאשר {\displaystyle r\_{cm}={\frac {m\_{1}r\_{1}+m\_{2}r\_{2}}{m\_{1}+m\_{2}}}} הוא מיקום מרכז המסה של שני החלקיקים, ו- {\displaystyle r\_{12}\equiv r\_{1}-r\_{2}}.

ניתן לבטא את {\displaystyle b\_{2}} במונחי הקואורדינטות החדשות כ:

{\displaystyle b\_{2}={\frac {1}{2V}}\int \_{V}{d^{3}r\_{cm}}\int \_{V}{f\_{12}d^{3}r\_{12}}={\frac {1}{2}}\int \_{V}{f\_{12}d^{3}r\_{12}}}

ולכן המקדם הויריאלי השני נתון על ידי:

{\displaystyle B\left(T\right)=-b\_{2}N\_{A}=-{\frac {N\_{A}}{2}}\int \_{V}{f\_{12}d^{3}r\_{12}}=-{\frac {N\_{A}}{2}}\int \_{V}{\left\{\exp {\left[-{\frac {\Phi \left(r\_{12}\right)}{kT}}\right]}-1\right\}d^{3}r\_{12}}=2\pi N\_{A}\int \_{0}^{\infty }{\left\{1-\exp {\left[-{\frac {\Phi \left(r\_{12}\right)}{kT}}\right]}\right\}r\_{12}^{2}dr\_{12}}}

כאשר בשוויון האחרון נעשה השימוש בקשר שמתקיים עבור סימטריה ספרית: {\displaystyle d^{3}\ r\_{12}=4\pi r\_{12}^{2}\ dr\_{12}}. לסיכום:

| {\displaystyle B\left(T\right)=2\pi N\_{A}\int \_{0}^{\infty }{\left\{1-\exp {\left[-{\frac {\Phi \left(r\_{12}\right)}{kT}}\right]}\right\}r\_{12}^{2}dr\_{12}}} |
| --- |

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

### המקדם הויריאלי השלישי

בהנחה שהאינטראקציות בין החלקיקים הן כוחות מרכזיים וכן שהפוטנציאל אדטיבי במובן ש:{\displaystyle \Phi \left(r\_{1},r\_{2},r\_{3}\right)=\Phi \left(r\_{12}\right)+\Phi \left(r\_{13}\right)+\Phi \left(r\_{23}\right)}, ניתן להציב את הפוטנציאל הנ"ל בביטוי לאינטגרל הקונפיגורציה {\displaystyle Q\_{3}}, ומהקשר בין המקדם {\displaystyle b\_{3}} לאינטגרלי הקונפיגורציה לקבל:

{\displaystyle b\_{3}={\frac {1}{3!V}}\iiint \_{V}{\left[e^{-{\frac {\Phi \left(r\_{12}\right)}{kT}}}e^{-{\frac {\Phi \left(r\_{13}\right)}{kT}}}e^{-{\frac {\Phi \left(r\_{23}\right)}{kT}}}-e^{-{\frac {\Phi \left(r\_{12}\right)}{kT}}}-e^{-{\frac {\Phi \left(r\_{13}\right)}{kT}}}-e^{-{\frac {\Phi \left(r\_{23}\right)}{kT}}}+2\right]d^{3}r\_{1}d^{3}r\_{2}d^{3}r\_{3}}={\frac {1}{3!V}}\iiint \_{V}{\left[f\_{12}f\_{13}f\_{23}+f\_{12}f\_{13}+f\_{12}f\_{23}+f\_{13}f\_{23}\right]d^{3}r\_{1}d^{3}r\_{2}d^{3}r\_{3}}}

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

על מנת לחשב את המקדם הויריאלי השלישי, יש למצוא ביטוי גם ל- {\displaystyle b\_{2}^{2}} (לפי הביטוי למקדם הויריאלי השלישי שפותח קודם). לשם כך, נתבונן באינטגרל: {\displaystyle \iiint \_{V}{f\_{12}f\_{13}d^{3}r\_{1}d^{3}r\_{2}d^{3}r\_{3}}}. ניתן לבצע מעבר קואורדינטות {\displaystyle r\_{1},\ r\_{2},r\_{3}\rightarrow r\_{cm},\ r\_{12},\ r\_{13}}, ולכתוב את האינטגרל הנ"ל במונחי הקואורדינטות החדשות:

{\displaystyle \iiint \_{V}{f\_{12}f\_{13}d^{3}r\_{1}d^{3}r\_{2}d^{3}r\_{3}}=\int \_{V}{d^{3}r\_{cm}}\int \_{V}{f\_{12}d^{3}r\_{12}}\int \_{V}{f\_{13}d^{3}r\_{13}}=V\int \_{V}{f\_{12}d^{3}r\_{12}}\int \_{V}{f\_{13}d^{3}r\_{13}}=4Vb\_{2}^{2}}

כאשר השוויון האחרון נובע מהקשר בין {\displaystyle b\_{2}} לפונקציית המקבץ של מאייר, שהוצג כאשר חושב המקדם הויריאלי השני.

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

| {\displaystyle C\left(T\right)=\left(4b\_{2}^{2}-2b\_{3}\right)N\_{A}^{2}=-{\frac {N\_{A}^{2}}{3V}}\iiint \_{V}{f\_{12}f\_{13}f\_{23}d^{3}r\_{1}d^{3}r\_{2}d^{3}r\_{3}}} |
| --- |

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

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

משוואת המצב עבור מודלים שונים לאינטראקציות בין החלקיקים - חישוב המקדם הויריאלי השני
-----------------------------------------------------------------------------------

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

איור 2: מודל הכדורים הקשיחים
### מודל הכדורים הקשיחים (hard sphere gas model)

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

{\displaystyle \Phi \left(r\_{12}\right)={\begin{cases}0\ \ \ \ \ \ \ \ r\_{12}>R\\\infty \ \ \ \ \ \ r\_{12}<R\end{cases}}}

לכן המקדם הויריאלי השני הוא:

{\displaystyle B\left(T\right)=2\pi N\_{A}\int \_{0}^{\infty }{\left\{1-\exp {\left[-{\frac {\Phi \left(r\_{12}\right)}{kT}}\right]}\right\}r\_{12}^{2}dr\_{12}}=2\pi N\_{A}\int \_{0}^{R}{r\_{12}^{2}dr\_{12}}={\frac {2\pi R^{3}N\_{A}}{3}}\equiv b\_{0}}

איור 3: בור פוטנציאל ריבועי
### בור פוטנציאל ריבועי (square-well potential)

{\displaystyle \Phi \left(r\_{12}\right)={\begin{cases}\infty \ \ \ \ \ \ \ 0<r\_{12}<\sigma \\-\varepsilon \ \ \ \ \ \sigma <r\_{12}<R\\0\ \ \ \ \ \ \ \ \ \ \ \ r\_{12}>R\end{cases}}}

לכן:

{\displaystyle B\left(T\right)=2\pi N\_{A}\int \_{0}^{\infty }{\left\{1-\exp {\left[-{\frac {\Phi \left(r\_{12}\right)}{kT}}\right]}\right\}r\_{12}^{2}dr\_{12}}={\frac {2\pi N\_{A}\sigma ^{3}}{3}}\left[1-\left(R^{3}-1\right)\left(e^{\frac {\varepsilon }{kT}}-1\right)\right]}

איור 4: פוטנציאל לנארד-ג'ונס 6-12
### פוטנציאל לנארד-ג'ונס 6-12 (Lennard-Jones 6-12 potential)

{\displaystyle \Phi \left(r\_{12}\right)=4\varepsilon \left[\left({\frac {\sigma }{r\_{12}}}\right)^{12}-\left({\frac {\sigma }{r\_{12}}}\right)^{6}\right]}

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

עבור מודל זה, המקדם הויריאלי השני הוא (כאשר נסמן {\displaystyle x\equiv {\frac {r\_{12}}{\sigma }}}, {\displaystyle T^{\ast }={\frac {kT}{\varepsilon }}}):

{\displaystyle B\_{2}\left(T\right)={\frac {4b\_{0}}{T^{\ast }}}\int \_{0}^{\infty }{x^{2}\left({\frac {12}{x^{12}}}-{\frac {6}{x^{6}}}\right)exp\left\{-{\frac {4}{T^{\ast }}}\left[\left({\frac {1}{x}}\right)^{12}-\left({\frac {1}{x}}\right)^{6}\right]\right\}dx}}

ואת האינטגרל הזה ניתן לחשב נומרית.

### גז ואן דר ואלס (van der Waals)

למעשה, המודל של גז ואן דר ואלס (van der Waals) הוא מקרה פרטי של הפיתוח הסטטיסטי לגז לא אידיאלי, כאשר הפוטנציאל המתאר את האינטראקציה הוא מהצורה:

{\displaystyle \Phi \left(r\_{12}\right)={\begin{cases}{-u}\_{0}\left({\frac {R}{r\_{12}}}\right)^{-6}\ \ \ \ \ \ r\_{12}>R\\\infty \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ r\_{12}<R\end{cases}}}

עבור פוטנציאל כזה, המקדם הויריאלי השני הוא:

{\displaystyle B\left(T\right)=2\pi N\_{A}\left[\int \_{0}^{R}{r^{2}dr+\int \_{R}^{\infty }{\left\{1-\exp {\left[{\frac {u\_{0}}{kT}}\left({\frac {R}{r}}\right)^{6}\right]}\right\}r^{2}dr}}\right]}

תחת ההנחה ש- {\displaystyle {\frac {u\_{0}}{kT}}\ll 1}, ניתן לקרב את האינטגרנד השני ל- {\displaystyle -{\frac {u\_{0}}{kT}}\left({\frac {R}{r}}\right)^{6}}, כך שמתקבל:

{\displaystyle B\left(T\right)\cong {\frac {2\pi N\_{A}R^{3}}{3}}\left(1-{\frac {u\_{0}}{kT}}\right)}

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

{\displaystyle \left(P+{\frac {2\pi N\_{A}^{2}R^{3}u\_{0}}{3}}{\frac {1}{v^{2}}}\ \right)\left(v-{\frac {2\pi N\_{A}R^{3}}{3}}\right)=kT}

אם נסמן {\displaystyle a\equiv {\frac {2\pi N\_{A}R^{3}u\_{0}}{3}}}, {\displaystyle b\equiv {\frac {2\pi N\_{A}R^{3}}{3}}}, נקבל את משוואת המצב של גז ואן דר ואלס.

### טבלת סיכום למודלים השונים

|  | מודל הכדורים הקשיחים | בור פוטנציאל ריבועי | פוטנציאל לנארד-ג'ונס 6-12 | גז ואן דר ואלס |
| --- | --- | --- | --- | --- |
| הפוטנציאל {\displaystyle \Phi \left(r\_{12}\right)} | 

Φ

(

r

12

)

=

{

0
 
 
 
 
 
 
 
 

r

12

>
R

∞
 
 
 
 
 
 

r

12

<
R



{\displaystyle \Phi \left(r\_{12}\right)={\begin{cases}0\ \ \ \ \ \ \ \ r\_{12}>R\\\infty \ \ \ \ \ \ r\_{12}<R\end{cases}}}

{\displaystyle \Phi \left(r_{12}\right)={\begin{cases}0\ \ \ \ \ \ \ \ r_{12}>R\\\infty \ \ \ \ \ \ r_{12}<R\end{cases}}} | 

Φ

(

r

12

)

=

{

∞
 
 
 
 
 
 
 
0
<

r

12

<
σ

−
ε
 
 
 
 
 
σ
<

r

12

<
R

0
 
 
 
 
 
 
 
 
 
 
 
 

r

12

>
R



{\displaystyle \Phi \left(r\_{12}\right)={\begin{cases}\infty \ \ \ \ \ \ \ 0<r\_{12}<\sigma \\-\varepsilon \ \ \ \ \ \sigma <r\_{12}<R\\0\ \ \ \ \ \ \ \ \ \ \ \ r\_{12}>R\end{cases}}}

{\displaystyle \Phi \left(r_{12}\right)={\begin{cases}\infty \ \ \ \ \ \ \ 0<r_{12}<\sigma \\-\varepsilon \ \ \ \ \ \sigma <r_{12}<R\\0\ \ \ \ \ \ \ \ \ \ \ \ r_{12}>R\end{cases}}} | 

Φ

(

r

12

)

=
4
ε

[

(

σ

r

12

)

12

−

(

σ

r

12

)

6

]

{\displaystyle \Phi \left(r\_{12}\right)=4\varepsilon \left[\left({\frac {\sigma }{r\_{12}}}\right)^{12}-\left({\frac {\sigma }{r\_{12}}}\right)^{6}\right]}

{\displaystyle \Phi \left(r_{12}\right)=4\varepsilon \left[\left({\frac {\sigma }{r_{12}}}\right)^{12}-\left({\frac {\sigma }{r_{12}}}\right)^{6}\right]} | 

Φ

(

r

12

)

=

{


−
u

0

(

R

r

12

)

−
6

 
 
 
 
 
 

r

12

>
R

∞
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

r

12

<
R



{\displaystyle \Phi \left(r\_{12}\right)={\begin{cases}{-u}\_{0}\left({\frac {R}{r\_{12}}}\right)^{-6}\ \ \ \ \ \ r\_{12}>R\\\infty \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ r\_{12}<R\end{cases}}}

{\displaystyle \Phi \left(r_{12}\right)={\begin{cases}{-u}_{0}\left({\frac {R}{r_{12}}}\right)^{-6}\ \ \ \ \ \ r_{12}>R\\\infty \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ r_{12}<R\end{cases}}} |
| המקדם הויריאלי השני {\displaystyle B\left(T\right)} | 



2
π

R

3

N

A

3

{\displaystyle {\frac {2\pi R^{3}N\_{A}}{3}}}

{\displaystyle {\frac {2\pi R^{3}N_{A}}{3}}} | 



2
π

N

A

σ

3

3

[

1
−

(

R

3

−
1

)

(

e

ε

k
T

−
1

)

]

{\displaystyle {\frac {2\pi N\_{A}\sigma ^{3}}{3}}\left[1-\left(R^{3}-1\right)\left(e^{\frac {\varepsilon }{kT}}-1\right)\right]}

{\displaystyle {\frac {2\pi N_{A}\sigma ^{3}}{3}}\left[1-\left(R^{3}-1\right)\left(e^{\frac {\varepsilon }{kT}}-1\right)\right]} | 



4

b

0

T

∗


∫

0

∞

x

2

(

12

x

12

−

6

x

6


)

e
x
p

{

−

4

T

∗


[

(

1
x

)

12

−

(

1
x

)

6

]

}

d
x

{\displaystyle {\frac {4b\_{0}}{T^{\ast }}}\int \_{0}^{\infty }{x^{2}\left({\frac {12}{x^{12}}}-{\frac {6}{x^{6}}}\right)exp\left\{-{\frac {4}{T^{\ast }}}\left[\left({\frac {1}{x}}\right)^{12}-\left({\frac {1}{x}}\right)^{6}\right]\right\}dx}}

{\displaystyle {\frac {4b_{0}}{T^{\ast }}}\int _{0}^{\infty }{x^{2}\left({\frac {12}{x^{12}}}-{\frac {6}{x^{6}}}\right)exp\left\{-{\frac {4}{T^{\ast }}}\left[\left({\frac {1}{x}}\right)^{12}-\left({\frac {1}{x}}\right)^{6}\right]\right\}dx}} | 



2
π

N

A

R

3

3

(

1
−

u

0

k
T

)

{\displaystyle {\frac {2\pi N\_{A}R^{3}}{3}}\left(1-{\frac {u\_{0}}{kT}}\right)}

{\displaystyle {\frac {2\pi N_{A}R^{3}}{3}}\left(1-{\frac {u_{0}}{kT}}\right)} |

מאפיינים של גז לא אידיאלי
-------------------------

לכל מאפיין של גז לא אידיאלי {\displaystyle X(T,V)} (למשל פונקציות המצב התרמודינמיות), ניתן להגדיר את התיקון למאפיין זה ביחס לגז אידיאלי {\displaystyle \Delta X(T,V)}, כך ש:  
{\displaystyle \Delta X\left(T,V\right)\equiv X\left(T,V\right)-X\_{id}(T,V)}

כאשר {\displaystyle X\_{id}(T,V)} הוא אותו המאפיין, עבור מודל של גז אידיאלי. מאחר שבנפח אינסופי ניתן להתייחס לכל גז כאל גז אידיאלי, ניתן לכתוב:

{\displaystyle \Delta X\left(T,V\right)\equiv \left[X\left(T,V\right)-X\left(T,\infty \right)\right]-\left[X\_{id}(T,V)-X\_{id}\left(T,\infty \right)\right]=\int \_{\infty }^{V}\left[\left({\frac {\partial X}{\partial V}}\right)\_{T}-\left({\frac {\partial X\_{id}}{\partial V}}\right)\_{T}\right]dX}

ראו גם
------

* גז אידיאלי
* גז ואן דר ואלס

ביבליוגרפיה
-----------

* Laurendeau, N. M. **Statistical Thermodynamics Fundamentals and Applications**. New York: Cambridge University Press, 2005
* Pathria, R. K. Beale, Paul D. **Statistical Mechanics**. 3rd ed. Amsterdam: Elsevier, 2011
* Reichl, L. E. **A Modern Course in Statistical Physics**. Weinheim, Germany: Wiley‐VCH Verlag GmbH & Co. KGaA, 2016
* Goodwin, A. R. H. Sengers, J. V. Peters, Cor J. **Applied Thermodynamics of Fluids**. Cambridge: RSC Pub. 2010