שרשרת פולימר אידיאלית

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

המודל מציג תיאור פשוט של פולימר כשרשרת של אובייקטים זהים בעלי אורך סופי l המתחברים זה לזה בקצותיהם בקונפורמציות זוויתיות שונות (ראה איור 1).
ממוזער|שמאל|350px|איור 1: תיאור פשוט של פולימר המורכב ממונומרים באורך קבוע b
## עקרונות המודל
N מונומרים זהים מרכיבים את הפולימר שאורכו הוא:L = N\,l 

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

נגדיר וקטור קצה-קצה \vec R של שרשרת פולימרית אידיאלית ווקטורים\vec r_1,\ldots ,\vec r_N
המצביעים המשתייכים למונומרים.נציין כי פיתוח זה מתייחס למצב בו מספר המונומרים N גדול מספיק, כך שמשפט הגבול המרכזי יתקיים.ממוזער|איור המדגים שרטוט של שרשרת פולימרית אידיאלית קצרה.קצוות השרשרת לא מחוברים לזה לזה לכן נקבל שערך התוחלת :

\langle \vec{R} \rangle = \Sigma_{i=1}^N \langle \vec r_i\rangle = \vec 0

היות ש\vec r_1,\ldots ,\vec r_N בלתי תלויים זה בזה, \vec Rמתפלג לפי התפלגות נורמלית (או גאוסיאנית). ולכן בשלשה ממדים נקבל שהשונות תהיה:
: \sigma ^2 = \langle R_x^2 \rangle - \langle R_x \rangle ^2 =\langle R_x^2 \rangle -0 
:ובפרט עבור כל מימד:
:\langle R_x^2\rangle =\langle R_y^2 \rangle = \langle R_z^2\rangle = N\,\frac{l^2}{3}
:ומכאן נובע ש: \langle \vec{R^2} \rangle = N\, l^2 = L \, l~
:ולכן וקטור קצה-קצה יהיה:\sqrt{\langle \vec{R^2} \rangle} = \sqrt N \, l = \sqrt{L \, l}~

כמו כן, וקטור הקצה-קצה מתפלג לפי פונקציית צפיפות ההסתברות הבאה:
:P(\vec R) = \left (\frac{3}{2 \pi N l^2} \right)^{3/2}e^{-\frac{3 \vec R^2}{2 N l^2}}

גודל נוסף, שמחושב פעמים רבות בפיזיקה של פולימרים, הוא רדיוס ההתמדה של הפולימר:
:\mathit{R}_G = \frac{ \sqrt N\, l }{ \sqrt 6\ } 



## רלוונטיות המודל במציאות והכללה של המודל

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

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

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

## דרגות חופש של פולימר

באיור 2 (משמאל) מוצג שרטוט של פולימר המורכב משלושה מונומרים, נקודת המוצא של כל מונומר i מסומנת ב \vec{r_i}.

לפי הגדרת הפולימר כל המונומרים זהים כלומר |\vec{r_1}-\vec{r_2}| = l. ממוזער|איור 2: פרמטרים של פולימר
שתי דרגות החופש של כל מונומר הן:\theta הזווית בין שני מונומריםו-\varphi זווית הסיבוב

## תרמודינמיקה של פולימר
פונקציית החלוקה : z = Tr (e^{-\beta H})האנרגיה החופשית : F = k_B T ln z = U-T S

ממוזער|איור 3: איור נוסף של פולימר, באדום מסומן <math>\vec{R_N}</math>

כדי להסתכל על התרמודינמיקה של הפולימר נגדיר שני משתנים חשובים:

\vec{R_N} הוא המרחק מקצה לקצה של הפולימר ומוגדר על ידי \vec{R_N} = \sum_{i=1}^N \vec{r_i} והזווית בין מונומר i למונומר i+1 תוגדר על ידי (ראה איור 2) cos \theta_{i,i+1} =\frac{ \vec{r_i} \vec{r_{i+1}}}{l^2}נזכיר כי ממוצע תרמי מוגדר ע"פ  = \frac{Tr Ae^{-\beta H}}{Tr e^{-\beta H}} עבור מודל זה: ההמילטוניאן יהיה תלוי בזוויות בלבד : H = H(\theta_i,\varphi_i)_i ובהתאם גם פונקציית החלוקה.

### הפרמטרים האופייניים של המערכת
#### היחס האופייני של פלורי
אידיאליות השרשרת תבוא לידי ביטוי על ידי חוסר קורלציה (מתאם) בין זוויות שונות המופרדות על ידי מרחקים גדולים כלומר lim_{|i-j|\to\infin} =0 **אך שימו לב** שהסכום על כל מיצועי הזוויות לא בהכרח מתאפס אלא שואף לקבועC'_i \equiv \sum_{j = 1}^{N} וסה"כ נקבל את **"היחס האופייני של פלורי"** 
 C_n \equiv \frac{1}{n} \sum_{i = 1}^{n} C'_i 

#### אורך קהן
גודל בעל חשיבות נוסף במערכת הוא אורך קהן (Kuhn length). ניתן להגדירו באמצעות היחס האופייני של פלורי: b = lC_{\infin}.יש הקוראים לו האורך האפקטיבי של המערכת l_{eff} ובספרים שונים הוא יכול להופיע בסימנים שונים.למעשה, באמצעות גודל אופייני זה ניתן לתאר שרשראות מורכבות, בעלות קורלציה (נמוכה!) בין זוויות של מונומרים קרובים, כשרשרת אידיאלית פשוטה עם אורך "מונומר" אפקטיבי b.
#### האורך האופייני של השרשרת
במערכת ישנם שני אורכים חשובים: האורך האופייני של השרשרת R_0 ורדיוס ההתמדה R_G, את הראשון קל לחשב משיקולים תרמודינמיים ואת השני קל למדוד בניסויי פיזור. בין שני הגדלים קיים יחס ליניארי : {R_G}^2 = \frac{1}{6} {R_0}^2 ,כך שע"י הצבת אחד הגדלים ניתן למצוא את השני.

#### רדיוס ההתמדה

באנגלית Radius of gyration, מוגדר ע"פ: R_g^2 \equiv \frac{1}{N} \sum_{i=1}^{N} (\vec{r_i}-\vec{R_{cm}})^2כאשר \vec{r_i} הוא וקטור מיקום של מונומר i ו \vec{R_{cm}} הוא וקטור מיקום מרכז המסה של הפולימר כולו.וקטור מרכז המסה מוגדר ע"פ \vec{R_{cm}} \equiv \frac{1}{N} \sum_{j=1}^{N} \vec{r_j}נכניס את הגדרת מרכז המסה להגדרת רדיוס ההתמדה ונקבל:R_g^2 = \frac{1}{N^2} \sum_{i=1}^{N}\sum_{j=1}^{N} (\vec{r_i}^2-\vec{r_i}\vec{r_j})עבור שרשרת אידיאלית יש סימטריה בין \vec{r_i}, \vec{r_j} ולכן מותר לכתוב: \vec{r_i}^2-\vec{r_i}\vec{r_j} = \vec{r_i}^2-2\vec{r_i}\vec{r_j}+\vec{r_j}^2 = (\vec{r_i}-\vec{r_j})^2.רדיוס ההתמדה הממוצע: = \frac{1}{N^2} \sum_{i=1}^{N}\sum_{j=1}^{N} בגבול הרצף : N \to \infin, l \to 0 ; Nl = const הסכומים הופכים לאינטגרלים  \sum_{i=1}^{N} \to \int_0^N du ,\sum_{j=i}^{N} \to \int_u^N dv

ואז אפשר לחשב את רדיוסי ההתמדה ע"פ הגאומטריה של הפולימר:

### מציאת האורך האופייני של השרשרת <math>R_0</math> משיקולים תרמודינמיים
משיקולי סימטריה ברור כי הממוצע  = 0 . זאת מכיוון שאנו סוכמים על כל הקונפיגורציות האפשריות בבעיה ומאחר ו-\vec{R_N} הוא וקטור, נוכל למצוא עבור כל \vec{r_i} וקטור השווה לו בגודלו והפוך בכיוונו \vec{-r_i}. 

אך השונות לא מתאפסת:  =  =  =\sum_{i,j}^N=\sum_{i,j=1}^N l^2 

ובגבול N>>1 נקבל:  \cong C_{\infin}Nl^2

מאחר שהממוצע על \vec{R_N} מתאפס, השונות : \sigma_{R_N}^2 = -^2 =  \ne 0 נגדיר את האורך האופייני של השרשרת כסטיית התקן : R_0 = \sqrt{}

## מודלים לתיאור שרשרת פולימר אידיאלית
במודלים הבאים נחשב את R_0 תחת מגבלות שונות על הזוויות \theta, \varphi בין המונומרים.
### שרשרת חופשית :Freely Joint Chain
במודל זה נניח כי אין הגבלות על בחירת הזוויות לאורך השרשרת.\theta מתפלגת אחיד (התפלגות אחידה) בין 0 ל-2\pi : \theta_i = \left\{ 0,2\pi \right\} \varphi מתפלגת אחיד (התפלגות אחידה) בין -\pi ל-\pi : \varphi_i = \left\{ -\pi,\pi \right\} 

ולכן : = l^2 = \delta_{ij}כאשר \delta_{ij} = \begin{cases} 0, & \text{if } i \ne j \\ 1, & \text{if } i=j \end{cases} היא הדלתא של קרוניקר.כלומר היחס האופייני של פלורי הוא: C_{\infin} = 1, ואורך קהן הוא: b = lוהאורך קצה-קצה האופייני של השרשרת :

R_0 = \sqrt{}=l\sqrt{N}
### שרשרת מסתובבת חופשית :Freely Rotating Chain
במודל זה נניח כי אחת הזוויות קבועה:\theta_i = \theta = const  \varphi מתפלגת אחיד בין -\pi ל \pi : \varphi_i = \left\{ -\pi,\pi \right\} 

נחשב את :

שכנים קרובים מדרגה ראשונה: = l^2 =l^2שכנים קרובים מדרגה שנייה:= l^2 =l^2מכאן נסיק כי עבור i,j כלליים נקבל : = l^2כלומר : R_0^2 = \sum_{i = 1}^N\sum_{j = 1}^N = \sum_{i = 1}^N+l^2\sum_{i = 1}^N[\sum_{j = 1}^{i-1}(cos\theta)^{i-j}+\sum_{j = i+1}^{N}(cos\theta)^{j-i}]נשים לב כי : \sum_{i = 1}^N=Nl^2  \sum_{j = 1}^{i-1} הוא פשוט סכימה על j\sum_{j = i+1}^{N} הוא פשוט סכימה על j>i 

נחליף למשתנה סכימה קל יותר : k = i-j ונקבל:R_0^2 = Nl^2+l^2\sum_{i = 1}^N[\sum_{k = 1}^{N}(cos\theta)^{k} +\sum_{k=1}^{N-i}(cos\theta)^{k}]תוצאת הסכום עד N כללי: \sum_{k=1}^{N}=\frac{cos \theta}{1-cos\theta} [1-cos^N \theta]ולכן נקבל כי R_0^2 = Nl^2+l^2\sum_{i = 1}^N\frac{cos \theta}{1-cos\theta}[2-cos^{i-1}\theta-cos^{N-i}\theta]

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

 \frac{r_ir_j}{l^2}= (cos\theta)^ = exp[|j-i| ln (cos\theta)] אם נגדיר את אורך הקורלציה : s_p = -\frac{1}{ln(cos \theta)}נקבל כי אורך הקורלציה יורד אקספוננציאלית עם המרחק בין המונומרים.

ובגבול N>>1 וזוויות גדולות נקבל:אורך קורלציה : s_p \cong 1, היחס האופייני של פלורי C_{\infin} = \frac{1+cos\theta}{1-cos\theta} \cong 2 , ואורך קהן b = l\frac{1+cos\theta}{1-cos\theta} \cong 2l 
 והאורך קצה-קצה האופייני של השרשרת : R_0^2 = Nl^2\frac{1+cos\theta}{1-cos\theta}
בגבול הרציפות כאשר N\rightarrow\infty, \theta\rightarrow0:

s_p\thickapprox\frac{2}{\theta^2}\rightarrow\infty

נגדיר את הגודל המקסימלי שהפולימר יכול לקבל: R_{max} = NLcos(\frac{\theta}{2})\thickapprox Nl

לכן נדרוש N\rightarrow\infty, l\rightarrow0, כך ש Nl=const. גבול זה נקרא שרשרת דמוית-תולעת (באנגלית: Worm-like chain).

ניזכר כי: R_0^2 = \sum_{i = 1}^N\sum_{j = 1}^N=\sum_{i = 1}^N\sum_{j = 1}^N(cos(\theta))^{\left\vert i-j \right\vert}=l^2\sum_{i = 1}^N\sum_{j = 1}^N
exp(-\frac{\left\vert i-j \right\vert}{s_p})

נקבל לאחר מעבר לגבול הרצף: R_0^2=l^2\textstyle \int_{0}^{R_{max}} \displaystyle \frac{du}{l}\textstyle \int_{0}^{R_{max}} \displaystyle \frac{dv}{l}
exp(-\frac{\left\vert u-v \right\vert}{ls_p})

ולבסוף נקבל: R_0^2=2ls_pR_{max}-2(ls_p)^2[1-exp(-\frac{R_{max}}{ls_p})]

נקבל שני מקרים: (בכל מקרה R_{max}\gg l)

# המקרה הקשיח: ls_p\gg R_{max}\gg l נפתח את האקספוננט לטור ונקבל R_0\thickapprox R_{max}=lN
# המקרה הגמיש: R_{max}\gg ls_p\gg l נקבל שהאקספוננט שואף לאפס, ולכן R_0^2\thickapprox 2ls_pR_{max}-2(ls_p)^2\thickapprox 2ls_pR_{max} לכן מקבלים R_0=\sqrt{lN2ls_p}= \sqrt{Ll_{eff}}, כלומר שרשרת אידיאלית עם l_{eff} = 2ls_p

### שרשרת מסתובבת מוגבלת: Hindered Rotation
במודל זה נניח כי אחת הזווית קבועה:
\theta_i = \theta = const  
\varphi מתפלגת ע"פ פוטנציאל U(\varphi_i) בין -\pi ל \pi, וההסתברות לקבלת זווית מסוימת מתפלגת ע"פ התפלגות בולצמן 
 = \frac{\int_0^{2\pi} (cos \varphi) exp(-U(\varphi_i)/k_B T)d\varphi }{\int_0^{2\pi}exp(-U(\varphi_i)/k_B T)d\varphi }

ובגבול N>>1 נקבל פשוטאורך קורלציה : s_p \cong \frac{1+}{1-}, היחס האופייני של פלורי C_{\infin} = \frac{1+cos\theta}{1-cos\theta} \frac{1+}{1-} , ואורך קהן b = l\frac{1+cos\theta}{1-cos\theta} \frac{1+}{1-}  
 והאורך קצה-קצה האופייני של השרשרת : R_0^2 = Nl^2\frac{1+cos\theta}{1-cos\theta}\frac{1+}{1-}

### שרשרת מסתובבת איזומרית Rotational Isomeric State
איזומר היא מולקולה שבה סידור האטומים יכול להתקיים בצורות שונות. מודל השרשרת המסתובבת האיזומרית מתייחס לקבוצת איזומרים מוגדרת:
# בעלי שלושה מצבים
# בין המצבים ניתן לעבור באמצעות סיבוב סביב קשר קוולנטי בודד בין שני אטומים.
# נניח מחסום פוטנציאל גבוה, כלומר \Delta E>>k_BT. כאשר \Delta E הוא פער האנרגיה בין המצבים.
ההנחות הזוויתיות של המודל מסתכמות ל:\theta_i = \theta = const  \varphi מתפלג אחיד בין מספר מצומצם של זוויות בדידות. 

מספר המצבים של המערכת נקבע ע"פ מספר הקשרים הראשיים בשרשרת, כלומר עבור n קשרים ראשיים יהיו n-2 זוויות סיבוב ולכן סה"כ יהיו 3^{n-2} מצבי סיבוב איזומריים.

לדוגמה ל-n-pentane יש 4 קשרים ראשיים ולכן 3^{4-2}=9 מצבי סיבוב.

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

### סיכום והשוואה בין מודלי שרשרת הפולימר האידיאלית

## לקריאה נוספת

* Rubinstein & Colby, **Polymer Physics**, OXFORD university press, 26 June 2003
* P.J.Flory, **Statistical Mechanics of Chain Molecules**, HANSER PUBLISHERS,
* P.J.Flory, ** Principles of Polymer Chemistry**, Cornell university press
* M.DOI & S.F.EDWARDS, **Polymer Physics**, OXFORD CLARENDON press

2016 06
קטגוריה:פולימרים