הצגת המקור של הדף "פולימרים ואינטגרלים מסלוליים"

[[קובץ:Syndiotactic_polypropene.png|ממוזער|250 פיקסלים|פולימר מורכב משרשראות ענק מולקולריות הבנויות ממולקולות קטנות - מונומרים]]
[[קובץ:Single_Polymer_Chains_AFM.jpg|ממוזער|250 פיקסלים|שרשראות פולימר ליניאריות אמיתיות איך שנצפו ע" [[מיקרוסקופ כוח אטומי]]]]
[[פולימרים]] הם [[מולקולה|מולקולות]] ענק המורכבות מיחידות חוזרות הקשורות ביניהן, הנפוצים הן בעולם החי כלו למשל ב-[[DNA]] והן בתעשייה כמו למשל ב[[פלסטיק]]. פולימרים שמקורם מהטבע מכונים [[ביופולימרים|פולימרים טבעיים]]. בנוסף להם קיימת קבוצה גדולה נוספת של פולימרים מלאכותיים או חצי מלאכותיים שנקראים פולימרים סינתטיים או חצי-סינתטיים בהתאמה. המבנה המיוחד והארוך של פולימרים מעניק להם תכונות פיזיקליות מיוחדות, כמו [[קשיות]], קשיחות, ויסקו-אלסטיות, ונטייה ליצור מבנים זכוכיתיים וקריסטליים למחצה. הפולימרים נחקרים בכל מיני ענפים, והחקר הפיזיקלי שלהם ב[[פיזיקה של חומר מעובה]] מתעניין בפלוקטואציות התרמיות שלהם, בתכונות המכניות, ובקיניטיקה של התגובות הכוללות תהליכי פירוק הפולימרים ל[[מונומר]]ים או [[פילמור]] המונומרים אל פולימרים.
משום שהפולימרים הם מולקולות גדולות, תכונותיהם הפיזיקליות בדרך כלל הן מסובכות בשיטות ניתוח דטרמיניסטיות. על כן, הגישה הסטטיסטית בדרך כלל היא היעילה במקרה שלהם, במיוחד אם מדובר בסקלה גדולה. הסיבה העקרונית להצלחה של תיאורם בגישה זו היא גודלם הגדול, תכונה המאפשרת את תקפות [[הגבול התרמודינמי]], אף על פי שמדובר במולקולות סופיות.
הפלוקטואציות התרמיות משפיעות באופן רציף על הצורה של הפולימרים בתמיסות. כדי למדל את התופעה משתמשים בעקרונות של [[פיזיקה סטטיסטית|מכניקה סטטיסטית]] ודינמיקה. גישת [[אינטגרלי מסלול|האינטגרל המסלולי]] היא גישה יעילה לטיפול בבעיות כאלו ומציעה דרך נוחה לחשב ממוצעים סטטיסטיים. האינטגרל המסלולי בעולם הפולימרים מתחשב בכל התצורות המרחביות האפשריות לפולימר, ממשקל אותן סטטיסטית, גם בנוכחות שדה חיצוני (ובצימוד לאמבט חום תרמי שקובע את הטמפרטורה של המערכת). הרבה בעיות שלא נפתרו קודם לגישה זו נפתרו כבר באמצעותה בצורה מוצלחת, כמו הנפח האסור, שזירה, חיבורים (links) וקשרים (knots), וכו'. מדענים בולטים תרמו לפיתוח תחום זה כוללים [[פייר ז'יל דה-ז'ן]], מאסאו דוי, פרידריך וויגל ו[[קלינרט האגן]].

==הקדמה לגישת האינטגרל המסלולי==
[[קובץ:Three_paths_from_A_to_B.png|ממוזער|250 פיקסלים|האינטגרל המסלולי בתורת הקוונטים: מתוארים רק שלושה מסלולים אפשריים לחלקיק לעבור מנקודה A בזמן מסוים ומקובע אל הנקודה B בזמן מרובע אחר. כל מסלול מאלה יתרום לאמפליטודת ההסתברות למעבר כזה]]
שיטת [[אינטגרלי מסלול|אינטגרלי המסלול]] (לא לבלבל עם [[אינטגרל קווי|אינטגרלים קווים]]) היא תיאור של [[מכניקת הקוונטים]] שעל פיה, למשל, ניתן להציג את [[הסתברות]] המעבר של חלקיק מנקודה אחת לאחרת כסכום ([[אינטגרל]]) על כל המסלולים האפשריים בהם החלקיק יכול לעבור בין שתי הנקודות. שיטה זו פותחה על ידי הפיזיקאי חתן [[פרס נובל]] [[ריצ'רד פיינמן]] ב-1948. גישה זו מצטרפת לאלו של [[ארווין שרדינגר|שרדינגר]] ו[[ורנר הייזנברג|הייזנברג]] כתיאור נוסף של מכניקת הקוונטים. היא מדגישה את השוני בין המכניקה הקוונטית ו[[מכניקה קלאסית|המכניקה הקלאסית]], כיוון שבזו האחרונה, מעבר חלקיק בין שתי נקודות נתונות מתבצע לאורך מסלול אחד ויחיד. אינטגרלי מסלול הוכיחו את עצמם כיעילים מאוד לצורך פיתוח התמונה הפיזיקלית, ופתרון של בעיות רבות בתורת הקוונטים, וההכללה שלהם ל[[תורת השדות הקוונטית|תורת שדות]] היא כיום השפה המקובלת של ה[[פיזיקה]] התאורטית, במיוחד כאשר רוצים לדבר על [[תורת ההפרעות (מכניקת הקוונטים)|חישובים הפרעתיים]].
כאמור, אינטגרלים מסוג זה הומצאו מתוך ניסיון להבנה עמוקה יותר של [[מכניקת הקוונטים]] והם שונים בתכלית מ[[אינטגרל קווי|אינטגרלי המסלול]] ב[[מתמטיקה]]. במתמטיקה אינטגרלי מסלול (או אינטגרלים מסילתיים) מתייחסים לאינטגרלים של [[פונקציה|פונקציות]] ([[סקלר (פיזיקה)|סקלריות]] או [[וקטור (פיזיקה)|וקטוריות]]) לאורך מסילות במרחב. לעומת זאת כאן האינטגרציה היא של [[פונקציונל]]ים והמסילות במרחב הן למעשה משתני האינטגרציה. כשם ש[[חשבון וריאציות]] הוא הכללה של פעולת ה[[נגזרת]] עבור פונקציונלים, אינטגרלי מסלול, בהקשר בו נדון כאן, מהווים הכללה של פעולת ה[[אינטגרציה (מתמטיקה)|אינטגרציה]].
אז בפורמליזם זה דנים באינטגרציה פונקציונלית, כאשר מספר המשתנים שמבצעים עליהם אינטגרציה הוא [[אינסוף]] רציף. אנחנו נראה את זה פה גם- בחקר ה[[פולימר]]ים באמצעות פורמליזם זה. במקרה שלנו, הפונקציה שמתבצע עליה האינטגרל הפונקציונלי (שהיא האינסוף הרציף שעליו מבצעים אינטגרציה) תהיה הצורה (או יותר מדויק ה[[קונפורמציה (כימיה)|קונפרמציה]]) של הפולימר.
==מודל הפולימר עם הקשר המעופף אקראית==
לפולימרים הגמישים ישנו מספר עצום של תצורות (או קונפיגורציות) שנובעות מדרגות החופש של כיווני [[קשר כימי|הקשרים הכימיים]]. צורת הפולימרים אפוא ניתנת לתיאור ביעילות על ידי [[פיזיקה סטטיסטית]] וממוצעים [[תרמודינמיקה|תרמודינמיים]]. נתאר בחלק זה את החקר של תכונות סטטיסטיות של פולימר אחד במצב [[שיווי משקל תרמודינמי]].
===השרשרת האידיאלית===
[[קובץ:Ideal_chain_random_walk.svg|ממוזער|250 פיקסלים|איור המדגים שרטוט של שרשרת פולימרית אידיאלית קצרה.]]
[[קובץ:linear_poly.png|ממוזער|שמאל|250 פיקסלים| תיאור פשוט של פולימר המורכב ממונומרים באורך קבוע b]]
[[קובץ:RN_for_polymer.png|ממוזער|250 פיקסלים|איור נוסף של פולימר, באדום מסומן <math>\mathbf{R}</math>]]
כדי לחקור פולימרים גמישים, מתחילים עם מודל פשטני שבו אנחנו מניחים <math>N</math> קשרים, כל אחד באורך קבוע <math>l</math>, וזווית מרחבית מקרית ב[[התפלגות אחידה רציפה|התפלגות אחידה]], וכל הקשרים [[אי תלות (סטטיסטיקה)|בלתי תלויים (הסתברותית)]]. מודל זה נקראה [[שרשרת פולימר אידיאלית]].
ה[[קונפורמציה (כימיה)|קונפורמציה]] של פולימר כזה, על כן, מתוארת על ידי <math>N+1</math> וקטורי מקום <math>{\left\{{\mathbf{R}}_n\right\}}^N_{n=0}</math>, או, באופן שקול, על ידי וקטורי הקשר <math>{\left\{{\mathbf{r}}_n\right\}}^N_{n=1}</math>, כאשר <math>{\mathbf{r}}_n={\mathbf{R}}_n-{\mathbf{R}}_{n-1}</math> (ו-<math>{\mathbf{R}}_0</math> אינו חשוב). נסמן ב-<math>\psi \left({\mathbf{r}}_n\right)</math> את [[פונקציית צפיפות|פונקציית צפיפות ההסתברות]] של וקטור הקשר <math>{\mathbf{r}}_n</math>, ובמקרה שלנו הפונקציה הזו היא
:<math>\psi \left(\mathbf{r}\right)=\frac{1}{4\pi l^2}\delta \left(\left|\mathbf{r}\right|-l\right)</math>
(מתקיים הנרמול <math>\int{d^3r\ \psi \left(\mathbf{r}\right)}=1</math>). לכן, צפיפות ההסתברות של קונפורמציה מסוימת היא (מאי-תלות נוכל לכפול)
:<math>\mathrm{\Psi }\left({\left\{{\mathbf{r}}_n\right\}}^N_{n=1}\right)=\prod^N_{n=1}{\psi \left({\mathbf{r}}_n\right)}</math>
כדי לאפיין את גודל הפולימר נגדיר את וקטור הקצה-קצה של הפולימר להיות <math>\mathbf{R}={\mathbf{R}}_N-{\mathbf{R}}_0=\sum^N_{n=1}{{\mathbf{r}}_n}</math>, ואז ברור כי <math>\left\langle \mathbf{R}\right\rangle =0</math>, והאורך הטיפוסי והאופייני של השרשרת מוגדר על ידי <math>\overline{R}=\sqrt{\left\langle {\mathbf{R}}^2\right\rangle }</math>, שבמקרה שלנו הוא:
:<math>\overline{R}=\sqrt{\left\langle {\mathbf{R}}^2\right\rangle }=\sqrt{\sum^N_{n=1,m=1}{\left\langle {\mathbf{r}}_n\cdot {\mathbf{r}}_m\right\rangle }}=\sqrt{\sum^N_{n=1,m=1}{l^2{\delta }_{nm}}}=\sqrt{N}\ l</math>
===השרשרת המסתובבת חופשית===
[[קובץ: A-sketch-of-a-freely-rotating-chain-The-bond-angle-th-is-fixed-whereas-the-torsional.png|ממוזער|250 פיקסלים|סקיצה לשרשרת עם קשר חופשי אזימוטלית. הקשר הוא בזווית קבועה <math>\theta </math> עם הקשר הקודם, בעוד שיש חופש בזווית האזימוטלית]]
נראה פה מודל אחר שבו גם נקבל <math>\overline{R}\propto \sqrt{N}</math>. במודל זה, הקשר ה-<math>n</math> מהווה זווית <math>\theta </math> קבועה עם הקשר הקודם (הקשר ה-<math>n-1</math>), והזווית האזימוטלית מתפלגת אחיד (ועדיין האורך הוא קבוע-<math>l</math>). מודל זה נקרא [[שרשרת פולימר אידיאלית| שרשרת מסתובבת חופשית]].
במקרה כזה, <math>{\mathbf{r}}_n</math> הוא [[שרשרת מרקוב|מרקובי]], ובהינתן <math>{\mathbf{r}}_{n-1}</math>, התוחלת של <math>{\mathbf{r}}_n</math> היא:
:<math>\left\langle {\mathbf{r}}_n|{\mathbf{r}}_{n-1}\right\rangle ={cos \theta \ }\ {\mathbf{r}}_{n-1}</math>
ולכן, מיד מבינים (מ[[משפט ההחלקה]]) כי
:<math>\left\langle {\mathbf{r}}_n\cdot {\mathbf{r}}_m\right\rangle ={\left({cos \theta \ }\right)}^{\left|n-m\right|}l^2</math>
ואז עבור <math>N\gg 1</math>:
:<math>\left\langle {\mathbf{R}}^2\right\rangle \cong Nl^2\frac{1+{cos \theta \ }}{1-{cos \theta \ }}</math>
כלומר יצא לנו:
:<math>\overline{R}=\sqrt{N}\ l\sqrt{\frac{1+{cos \theta \ }}{1-{cos \theta \ }}}\equiv \sqrt{N}\ l_{eff}</math>
כאשר <math>l_{eff}\equiv \frac{\overline{R}}{\sqrt{N}}=l\sqrt{\frac{1+{cos \theta \ }}{1-{cos \theta \ }}}</math>.
והנה פירוט לחישוב:
:<math>\left\langle {\mathbf{R}}^2\right\rangle =\sum^N_{n=1,m=1}{\left\langle {\mathbf{r}}_n\cdot {\mathbf{r}}_m\right\rangle }=\sum^N_{n=1,m=1}{{\left({cos \theta \ }\right)}^{\left|n-m\right|}l^2}=</math>
:<math>=l^2\left(2\sum^N_{ \begin{array}{c}
n=2, \\
m<n \end{array}
}{{\left({cos \theta \ }\right)}^{n-m}}+\sum^N_{n=1}{1}\right)=l^2\left(2\sum^N_{m=2}{\sum^{n-1}_{m=1}{{\left({cos \theta \ }\right)}^{n-m}}}+N\right)=</math>
:<math>=l^2\left(2\sum^N_{m=2}{{\left({cos \theta \ }\right)}^n\sum^{n-1}_{m=1}{{\left({cos \theta \ }\right)}^{-m}}}+N\right)=l^2\left(2\sum^N_{m=2}{{\left({cos \theta \ }\right)}^n{\left({cos \theta \ }\right)}^{-1}\frac{1-{\left({cos \theta \ }\right)}^{-(n-1)}}{1-{\left({cos \theta \ }\right)}^{-1}}}+N\right)=</math>
:<math>=l^2\left(2\sum^N_{m=2}{\frac{{\left({cos \theta \ }\right)}^{n-1}-1}{1-{\left({cos \theta \ }\right)}^{-1}}}+N\right)=l^2\left(\frac{2}{1-{\left({cos \theta \ }\right)}^{-1}}\left[-N+1+\sum^{N-1}_{m=1}{{\left({cos \theta \ }\right)}^n}\right]+N\right)=</math>
:<math>=l^2\left(\frac{2}{1-{\left({cos \theta \ }\right)}^{-1}}\left[-N+1+{cos \theta \ }\frac{1-{\left({cos \theta \ }\right)}^{(N-1)}}{1-{cos \theta \ }}\right]+N\right)=</math>
:<math>= l^2\frac{1}{1-{\left({cos \theta \ }\right)}^{-1}}\left(-N-N{\left({cos \theta \ }\right)}^{-1}+2+\frac{2{cos \theta \ }}{1-{\left({cos \theta \ }\right)}^{-1}}\frac{1-{\left({cos \theta \ }\right)}^{(N-1)}}{1-{cos \theta \ }}\right)\cong </math>
:<math>\cong l^2\frac{1}{1-{\left({cos \theta \ }\right)}^{-1}}\left(-N-N{\left({cos \theta \ }\right)}^{-1}\right)=-Nl^2\frac{1+{\left({cos \theta \ }\right)}^{-1}}{1-{\left({cos \theta \ }\right)}^{-1}}=Nl^2\frac{1+{cos \theta \ }}{1-{cos \theta \ }}</math>
כאשר בשלב מסוים לקחנו <math>N\gg 1</math>.
===שרשרת עם קשר מעופף אקראית===
באופן כללי, אם מדובר בתעופה אקראית מרקובית, אז פונקציית הדפורמציה ניתנת לרישום כך:
<math>\mathit{\Psi}\left({\left\{{\mathbf{r}}_n\right\}}^N_{n=1}\right)={\psi }_0\left({\mathbf{r}}_1\right)\prod^N_{n=2}{\psi \left({\mathbf{r}}_n|{\mathbf{r}}_{n-1}\right)}</math>
וניתן לראות שעבור <math>N\gg 1</math> גם מקבלים <math>\overline{R}\propto \sqrt{N}</math>. אפשר גם בקלות להראות שזה עדיין יתקיים עם תנאי יותר חלש: סמי מרקובית- בהינתן התפלגות עד מספר שכנים קבוע מלפני, כבר השאר לא חשוב (אי תלות הסתברותית). כאשר מסמנים <math>\overline{R}\cong \sqrt{N}l_{eff}</math>, קבוע הפרופורציה <math>l_{eff}</math> נקרא האורך האפקטיבי של הקשר. הקבוע <math>C_{\infty }=\frac{l^2_{eff}}{l^2}</math> מייצג את הקשיחות של השרשרת.
===הפילוג של וקטור הקצה-קצה===
בחזרה לשרשרת האידיאלית, פונקציית צפיפות הפילוג של וקטור הקצה-קצה ניתנת ע"י
:<math>\mathrm{\Phi }\left(\mathbf{R},N\right)=\int{d^3r_1\int{d^3r_2}\dots \int{d^3r_N}\delta \left(\mathbf{R}-\sum^N_{n=1}{{\mathbf{r}}_n}\right)\mathit{\Psi}\left({\left\{{\mathbf{r}}_n\right\}}^N_{n=1}\right)}=\frac{1}{{\left(2\pi \right)}^3}\int{d^3r_1\int{d^3r_2}\dots \int{d^3r_N}\int{d^3k}e^{i\mathbf{k}\cdot \left(\mathbf{R}-\sum^N_{n=1}{{\mathbf{r}}_n}\right)}\prod^N_{n=1}{\psi \left({\mathbf{r}}_n\right)}}=</math>
:<math>=\frac{1}{{\left(2\pi \right)}^3}\int{d^3ke^{i\mathbf{k}\cdot \mathbf{R}}}\prod^N_{n=1}{\int{d^3r_ne^{-i\mathbf{k}\cdot {\mathbf{r}}_n}\psi \left({\mathbf{r}}_n\right)}}=\frac{1}{{\left(2\pi \right)}^3}\int{d^3ke^{i\mathbf{k}\cdot \mathbf{R}}}{\left(\int{d^3re^{-i\mathbf{k}\cdot \mathbf{r}}\psi \left(\mathbf{r}\right)}\right)}^N=\frac{1}{{\left(2\pi \right)}^3}\int{d^3ke^{i\mathbf{k}\cdot \mathbf{R}}}{\left(\frac{{sin kl\ }}{kl}\right)}^N\cong </math>
:<math>\cong \frac{1}{{\left(2\pi \right)}^3}\int{d^3ke^{i\mathbf{k}\cdot \mathbf{R}}}e^{-\frac{Nl^2k^2}{6}}={\left(\frac{3}{2\pi Nl^2}\right)}^{\frac{3}{2}}{exp \left(-\frac{3{\mathbf{R}}^2}{2Nl^2}\right)\ }</math>
כאשר השתמשנו בהצגת פורייה של הדלתא של דיראק, וגם בכך שעבור <math>N\gg 1</math>, <math>{\left(\frac{{sin kl\ }}{kl}\right)}^N</math> נהיה מאוד קטן אלא אם כן <math>kl</math> הוא מאוד קטן, ואז תקף הקירוב <math>{\left(\frac{{sin kl\ }}{kl}\right)}^N\cong \ {\left(1-\frac{{\left(kl\right)}^2}{6}\right)}^N\cong e^{-\frac{Nl^2k^2}{6}}</math>, וגם השתמשנו ב[[אינטגרל גאוסיאני|אינטגרל הגאוסי]] בסוף. קיבלנו אפוא שהפילוג הוא [[התפלגות נורמלית|גאוסי]]. פילוג זה תמיד יתקבל עבור שרשרת עם קשר מעופף רנדומית אם
<math>N</math> הוא מספיק גדול, כאשר השונות תהיה פרופורציונלית ל-
<math>N</math>.
==השרשרת הגאוסיאנית==
===לפני גבול הרצף===
[[קובץ: The-equivalent-picture-for-a-Gaussian-chain-is-a-number-of-effective-monomers-connected.png|ממוזער|250 פיקסלים|השרשרת הגאוסיאנית]]
ראינו שבפילוג של וקטור הקצה-קצה, המבנה הלוקאלי של הקשר מתבטא רק דרך האורך האפקטיבי של הקשר <math>l_{eff}</math>, וזה נכון באופן כללי: המבנה הלוקאלי משנה רק את האורך האפקטיבי, ולא מתבטא בצורות אחרות. לכן, כדי לחקור תכונות גלובאליות של פולימרים, נוכל להתחיל רק מהמודל הפשוט שקיים לנו.
נחשוב על שרשרת שוקטור הקשר הבסיסי שלה מתפלג [[התפלגות נורמלית|גאוסית]] (בתלת מימד):
:<math>\psi \left(\mathbf{r}\right)={\left(\frac{3}{2\pi l^2}\right)}^{\frac{3}{2}}{exp \left(-\frac{3{\mathbf{r}}^2}{2l^2}\right)\ }</math>
כאשר <math>\left\langle {\mathbf{r}}^2\right\rangle =l^2</math>. פונקציית צפיפות ההסתברות של הקונפורמציה עכשיו היא
:<math>\mathit{\Psi}\left({\left\{{\mathbf{r}}_n\right\}}^N_{n=1}\right)=\prod^N_{n=1}{{\left(\frac{3}{2\pi l^2}\right)}^{\frac{3}{2}}{exp \left(-\frac{3{\mathbf{r}}^2_n}{2l^2}\right)\ }}={\left(\frac{3}{2\pi l^2}\right)}^{\frac{3}{2}N}{exp \left(-\sum^N_{n=1}{\frac{3{\left({\mathbf{R}}_n-{\mathbf{R}}_{n-1}\right)}^2}{2l^2}}\right)\ }</math>
וזה נקרא [[שרשרת פולימרית גאוסיאנית|שרשרת גאוסית]], והיא מתארת נכונה לא את המבנה הלוקאלי של הפולימר, אלא את התכונות שלו בסקלה גדולה. היתרון בשרשרת הגאוסית הוא שהיא מודל שהוא הרבה יותר קל לטיפול מתמטי מאשר המודלים הקודמים.
תכונה חשובה למשל לשרשרת הגאוסית היא שהפילוג של <math>{\mathbf{R}}_n-{\mathbf{R}}_m</math> הוא גם כן גאוסי (תלת־ממדי), אך עם שונות <math>\left\langle {\left({\mathbf{R}}_n-{\mathbf{R}}_m\right)}^2\right\rangle =\left|n-m\right|l^2</math>, כלומר
:<math>\mathrm{\Phi }\left({\mathbf{R}}_n-{\mathbf{R}}_m,n-m\right)={\left(\frac{3}{2\pi \left|n-m\right|l^2}\right)}^{\frac{3}{2}}{exp \left(-\frac{3{\left({\mathbf{R}}_n-{\mathbf{R}}_m\right)}^2}{2\left|n-m\right|l^2}\right)\ }</math>
===גבול הרצף===
כאשר לוקחים את [[הגבול התרמודינמי|גבול הרצף]] מקבלים שפונקציית צפיפות ההסתברות של הקונפורמציה היא
:<math>\mathrm{\Psi }\left({\mathbf{R}}_n\right)=C{\mathrm{exp} \left[-\frac{3}{2l^2}\int^N_0{dn{\left(\frac{\partial {\mathbf{R}}_n}{\partial n}\right)}^2}\right]\ }</math>
עכשיו <math>n</math> (שהיה מספר המונומרים או מספר הצעדים) הוא [[תיאור פרמטרי של עקום|משתנה פרמטריזציה]] לאורך הפולימר.
== קונפורמציית השרשרת תחת שדה חיצוני ואינטגרלים מסלוליים==
===פונקציית המשקל הסטטיסטי או פונקציית גרין===
[[קובץ:Three_paths_from_A_to_B.png|ממוזער|250 פיקסלים|שלושה מסלולים אפשריים לפולימר שמתחיל בנקודה A ומסתיים בנקודה B]]
כדי לחקור את התכונות הסטטיסטיות של פולימר תחת שדה חיצוני, צריך לדעת מהו המשקל הסטטיסטי של כל [[קונפורמציה (כימיה)|קונפורמציה]] או צורה של הפולימר, ואז לחשב את [[פונקציית חלוקה (פיזיקה)|פונקציית החלוקה]], את [[אנרגיה חופשית של הלמהולץ|פונקציית האנרגיה החופשית של הלמהולץ]], ואז לחשב כל מיני גדלים תרמודינמיים.
אם ישנו פוטנציאל חיצוני <math>U\left(\mathbf{r}\right)</math> שפועל על כל מקטע, הפילוג ב[[שיווי משקל תרמודינמי]] מתוקן על ידי [[התפלגות בולצמן|פקטור בולצמן]]
:<math>\mathit{\Psi}\left({\mathbf{R}}_n\right)=C{exp \left[-\frac{3}{2l^2}\int^N_0{dn{\left(\frac{\partial {\mathbf{R}}_n}{\partial n}\right)}^2}-\frac{1}{k_BT}\int^N_0{dn\ U\left({\mathbf{R}}_n\right)}\right]\ }</math>
פונקציית המשקל הסטטיסטי, אפוא, היא
:<math>P_N\left(\mathbf{R},{\mathbf{R}}^'\right)=\frac{\int^{{\mathbf{R}}_N=\mathbf{R}}_{{\mathbf{R}}_0={\mathbf{R}}^'}{\mathcal{D}{\mathbf{R}}_n\ e^{-\frac{3}{2l^2}\int^N_0{dn{\left(\frac{\partial {\mathbf{R}}_n}{\partial n}\right)}^2}-\frac{1}{k_BT}\int^N_0{dn\ U\left({\mathbf{R}}_n\right)}}}}{\int{d^3R\ \int^{{\mathbf{R}}_N=\mathbf{R}}_{{\mathbf{R}}_0={\mathbf{R}}^'}{\mathcal{D}{\mathbf{R}}_n\ e^{-\frac{3}{2l^2}\int^N_0{dn{\left(\frac{\partial {\mathbf{R}}_n}{\partial n}\right)}^2}}}}}</math>
כאשר זה משמעותו שאנחנו בוחנים (או לוקחים בחשבון) את כל הקונפורמציות שמתחילות ב- <math>{\mathbf{R}}_0={\mathbf{R}}^'</math>. בנוכחות פוטנציאל חיצוני, <math>P_N\left(\mathbf{R},{\mathbf{R}}^'\right)</math> מייצגת את המשקל הסטטיסטי של השרשרת שמתחילה ב- <math>{\mathbf{R}}_0={\mathbf{R}}^'</math> ומסתיימת ב- <math>{\mathbf{R}}_N=\mathbf{R}</math>. נשים לב שהמכנה אינו תלוי לא ב- <math>\mathbf{R}</math> ולא ב-<math>{\mathbf{R}}^'</math>.
מסיבות שיתבררו בחלק הבא, נקרא לפונקציה זו גם [[פונקציית גרין]], ונסמן אותה ב-<math>G\left(\mathbf{R},{\mathbf{R}}^';N\right)</math>:
:<math>G\left(\mathbf{R},{\mathbf{R}}^';N\right)\equiv P_N\left(\mathbf{R},{\mathbf{R}}^'\right)</math>
בהיעדר כוח חיצוני (<math>U\left(\mathbf{r}\right)=0</math>), <math>G\left(\mathbf{R},{\mathbf{R}}^';N\right)</math> תהיה פשוט הפילוג הגאוסי:
:<math>G_0\left(\mathbf{R},{\mathbf{R}}^';N\right)=G_0\left(\mathbf{R}-{\mathbf{R}}^';N\right)={\left(\frac{3}{2\pi Nl^2}\right)}^{\frac{3}{2}}{exp \left(-\frac{3{\left(\mathbf{R}-{\mathbf{R}}^'\right)}^2}{2Nl^2}\right)\ }</math>
פונקציית החלוקה (המתחשבת בכל הקונפורמציות האפשריות) היא פשוט:
:<math>Z=\int{d^3R\ \int{d^3R^'\ G\left(\mathbf{R},{\mathbf{R}}^';N\right)}}</math>
===המשוואה הדיפרנציאלית שפונקציית גרין מקיימת===
אף על פי שלפונקציית המשקל הסטטיסטי הקודמת <math>P_N\left(\mathbf{R},{\mathbf{R}}^'\right)\equiv G\left(\mathbf{R},{\mathbf{R}}^';N\right)</math> יש משמעות פיזיקלית, נגדיר אותה גם עבור <math>N</math> שלילי להיות אפס, ואז, נסיק כי היא מקיימת את המשוואה הדיפרנציאלית הבאה:
:<math>\left[\frac{\partial }{\partial N}-\frac{l^2}{6}\frac{{\partial }^2}{\partial {\mathbf{R}}^2}+\frac{1}{k_BT}U\left(\mathbf{R}\right)\right]G\left(\mathbf{R},{\mathbf{R}}^';N\right)=\delta \left(\mathbf{R}-{\mathbf{R}}^'\right)\delta \left(N\right)</math>
ולכן קראנו לה [[פונקציית גרין]]. עכשיו נראה את זה.
קודם כל, בהיעזר ב[[התמרת פורייה]] של <math>G_0\left(\mathbf{R},{\mathbf{R}}^';N\right)</math> (כאשר זו היא [[פונקציית דלתא של דיראק|פונקציית הדלתא של דיראק]]),
:<math>G_0\left(\mathbf{R},{\mathbf{R}}^';N\right)=\frac{1}{{\left(2\pi \right)}^3}\int{d^3k\ e^{-i\mathbf{k}\cdot \left(\mathbf{R}-{\mathbf{R}}^'\right)}e^{-\frac{Nl^2k^2}{6}}}\mathrm{\Theta }\left(N\right)</math>
נוכל בקלות לוודא כי
:<math>\frac{\partial }{\partial N}G_0\left(\mathbf{R},{\mathbf{R}}^';N\right)=\frac{1}{{\left(2\pi \right)}^3}\int{d^3k\ e^{-i\mathbf{k}\cdot \left(\mathbf{R}-{\mathbf{R}}^'\right)}\left[-\frac{l^2k^2}{6}e^{-\frac{Nl^2k^2}{6}}\mathrm{\Theta }\left(N\right)+\delta \left(N\right)\right]}=\frac{l^2}{6}\frac{{\partial }^2}{\partial {\mathbf{R}}^2}G_0\left(\mathbf{R},{\mathbf{R}}^';N\right)+\delta \left(\mathbf{R}-{\mathbf{R}}^'\right)\delta \left(N\right)</math>
וכדי לגזור את המשוואה בנוכחות פוטנציאל חיצוני <math>U</math>, נשים לב שמתוך הגדרת <math>G\left(\mathbf{R},{\mathbf{R}}^';N\right)</math>, מתקיים:
[[קובץ:r_tag_tag.png|ממוזער|250 פיקסלים|שרשרת שמתחילה ב-<math>{\mathbf{R}}^'</math>, עוברת דרך <math>{\mathbf{R}}^{''}</math> אחרי <math>n</math> "צעדים" (או יותר מדויק: אחרי ריצת <math>n</math> של [[תיאור פרמטרי של עקום|משתנה הפרמטריזציה]] אחרי [[הגבול התרמודינמי|גבול הרצף]]), ומסתיימת ב-<math>\mathbf{R}</math> אחרי <math>N</math> "צעדים" בסה"כ.]]
:<math>G\left(\mathbf{R},{\mathbf{R}}^';N\right)=\int{d^3R^{''}\ G\left(\mathbf{R},{\mathbf{R}}^{''};N-n\right)G\left({\mathbf{R}}^{''},{\mathbf{R}}^';n\right)}</math>
לכל <math>0<n<N</math>. פיזיקלית, הפקטור <math>G\left(\mathbf{R},{\mathbf{R}}^{''};N-n\right)G\left({\mathbf{R}}^{''},{\mathbf{R}}^';n\right)</math> מייצג את המשקל הסטטיסטי לשרשרת שמתחילה ב-<math>{\mathbf{R}}^'</math>, עוברת דרך <math>{\mathbf{R}}^{''}</math> אחרי <math>n</math> "צעדים" (או יותר מדויק: אחרי ריצת <math>n</math> של [[תיאור פרמטרי של עקום|משתנה הפרמטריזציה]] [[הגבול התרמודינמי|גבול הרצף]]), ומסתיימת ב-<math>\mathbf{R}</math> אחרי <math>N</math> "צעדים" בסה"כ. סכימה (או אינטגרציה במקרה זה) על כל ה-<math>{\mathbf{R}}^{''}</math> תיתן את המשקל הסטטיסטי לשרשרת שמתחילה ב-<math>{\mathbf{R}}^'</math> ומסתיימת ב-<math>\mathbf{R}</math> אחרי <math>N</math> "צעדים".
נשתמש בתכונה זו עבור <math>\mathrm{\Delta }N</math> קטן:
:<math>G\left(\mathbf{R},{\mathbf{R}}^';N+\mathrm{\Delta }N\right)=\int{d^3R^{''}\ G\left(\mathbf{R},{\mathbf{R}}^{''};\mathrm{\Delta }N\right)G\left({\mathbf{R}}^{''},{\mathbf{R}}^';N\right)}</math>
בטווח ה-<math>n</math>-ים שבין <math>N</math> ו-<math>N+\mathrm{\Delta }N</math>, בסדר הראשון ב-<math>\mathrm{\Delta }N</math> נקבל:
:<math>\int^{N+\mathrm{\Delta }N}_N{dn\ U\left({\mathbf{R}}_n\right)}\cong \mathrm{\Delta }N\ U\left(\mathbf{R}\right)</math>
ולכן:
:<math>G\left(\mathbf{R},{\mathbf{R}}^{''};\mathrm{\Delta }N\right)\cong e^{-\frac{\mathrm{\Delta }N}{k_BT}U\left(\mathbf{R}\right)}G_0\left(\mathbf{R},{\mathbf{R}}^{''};\mathrm{\Delta }N\right)</math>
ואז:
:<math>G\left(\mathbf{R},{\mathbf{R}}^';N+\mathrm{\Delta }N\right)\cong e^{-\frac{\mathrm{\Delta }N}{k_BT}U\left(\mathbf{R}\right)}\int{d^3R^{''}\ G_0\left(\mathbf{R}-{\mathbf{R}}^{''};\mathrm{\Delta }N\right)G\left({\mathbf{R}}^{''},{\mathbf{R}}^';N\right)}</math>
ועבור <math>\mathrm{\Delta }N</math> קטן, ל-<math>G_0\left(\mathbf{R}-{\mathbf{R}}^{''};\mathrm{\Delta }N\right)</math> ישנו פיק חד ב-<math>\mathbf{R}={\mathbf{R}}^'</math>, ולכן, נוכל לפתח סביב <math>\mathbf{r}\mathrm{=}\mathbf{R}-{\mathbf{R}}^{''}</math>:
:<math>G\left(\mathbf{R},{\mathbf{R}}^';N+\mathrm{\Delta }N\right)\cong e^{-\frac{\mathrm{\Delta }N}{k_BT}U\left(\mathbf{R}\right)}\int{d^3R^{''}\ G_0\left(\mathbf{R}-{\mathbf{R}}^{''};\mathrm{\Delta }N\right)G\left({\mathbf{R}}^{''},{\mathbf{R}}^';N\right)}=</math>
:<math>=e^{-\frac{\mathrm{\Delta }N}{k_BT}U\left(\mathbf{R}\right)}\int{d^3r\ G_0\left(\mathbf{r};\mathrm{\Delta }N\right)G\left(\mathbf{R}-\mathbf{r},{\mathbf{R}}^';N\right)}\cong e^{-\frac{\mathrm{\Delta }N}{k_BT}U\left(\mathbf{R}\right)}\int{d^3r\ G_0\left(\mathbf{r};\mathrm{\Delta }N\right)\left(1-r_{\alpha }\frac{\partial }{\partial R_{\alpha }}+\frac{1}{2}r_{\alpha }r_{\beta }\frac{\partial }{\partial R_{\alpha }}\frac{\partial }{\partial R_{\beta }}\right)G\left(\mathbf{R},{\mathbf{R}}^';N\right)}=</math>
:<math>=e^{-\frac{\mathrm{\Delta }N}{k_BT}U\left(\mathbf{R}\right)}\left(1-0+\frac{1}{2}\frac{1}{3}\mathrm{\Delta }Nl^2{\delta }_{\alpha \beta }\frac{\partial }{\partial R_{\alpha }}\frac{\partial }{\partial R_{\beta }}\right)G\left(\mathbf{R},{\mathbf{R}}^';N\right)=e^{-\frac{\mathrm{\Delta }N}{k_BT}U\left(\mathbf{R}\right)}\left(1+\frac{1}{6}\mathrm{\Delta }Nl^2\frac{{\partial }^2}{\partial {\mathbf{R}}^2}\right)G\left(\mathbf{R},{\mathbf{R}}^';N\right)</math>
ולכן, אם נרשום את הכל עד הסדר הראשון ב-<math>\mathrm{\Delta }N</math>, נקבל:
:<math>\left(1+\mathrm{\Delta }N\frac{\partial }{\partial N}\right)G\left(\mathbf{R},{\mathbf{R}}^';N\right)=\left(1-\frac{\mathrm{\Delta }N}{k_BT}U\left(\mathbf{R}\right)\right)\left(1+\frac{1}{6}\mathrm{\Delta }Nl^2\frac{{\partial }^2}{\partial {\mathbf{R}}^2}\right)G\left(\mathbf{R},{\mathbf{R}}^';N\right)=\left(1-\frac{\mathrm{\Delta }N}{k_BT}U\left(\mathbf{R}\right)+\frac{1}{6}\mathrm{\Delta }Nl^2\frac{{\partial }^2}{\partial {\mathbf{R}}^2}\right)G\left(\mathbf{R},{\mathbf{R}}^';N\right)</math>
כלומר
:<math>\frac{\partial }{\partial N}G\left(\mathbf{R},{\mathbf{R}}^';N\right)=\left(-\frac{\mathrm{1}}{k_BT}U\left(\mathbf{R}\right)+\frac{l^2}{6}\frac{{\partial }^2}{\partial {\mathbf{R}}^2}\right)G\left(\mathbf{R},{\mathbf{R}}^';N\right)</math>
וזה רוב מה שרצינו להוכיח. זה היה נכון עבור כל <math>N>0</math>. כדי לראות מה קורה באפס, נסתכל על <math>N</math>-ים קטנים:
:<math>G\left(\mathbf{R},{\mathbf{R}}^';N\right)\cong e^{-\frac{N}{k_BT}U\left(\mathbf{R}\right)}G_0\left(\mathbf{R},{\mathbf{R}}^';N\right)</math>
(הסברנו את זה לפני). ומכאן מגיע ה-<math>\delta \left(\mathbf{R}-{\mathbf{R}}^'\right)\delta \left(N\right)</math> ולסיכום, הוכחנו את מה שרצינו להוכיח:
:<math>\left[\frac{\partial }{\partial N}-\frac{l^2}{6}\frac{{\partial }^2}{\partial {\mathbf{R}}^2}+\frac{1}{k_BT}U\left(\mathbf{R}\right)\right]G\left(\mathbf{R},{\mathbf{R}}^';N\right)=\delta \left(\mathbf{R}-{\mathbf{R}}^'\right)\delta \left(N\right)</math>
וכך ברור למה קראנו לפונקציית המשקל הסטטיסטי בשם [[פונקציית גרין]].
===אנלוגיה עם תורת הקוונטים===
נשים לב כי עבור <math>N>0</math>, המשוואה על <math>P_N\left(\mathbf{R},{\mathbf{R}}^'\right)</math> היא:
:<math>-\frac{\partial P_N\left(\mathbf{R},{\mathbf{R}}^'\right)}{\partial N}=\left[-\frac{l^2}{6}\frac{{\partial }^2}{\partial {\mathbf{R}}^2}+\frac{1}{k_BT}U\left(\mathbf{R}\right)\right]P_N\left(\mathbf{R},{\mathbf{R}}^'\right)</math>
וזה מאוד דומה ל[[משוואת שרדינגר]]:
:<math>i\hslash \frac{\partial {\psi }_t\left(\mathbf{R},{\mathbf{R}}^'\right)}{\partial t}=\left[-\frac{{\hslash }^2}{2m}\frac{{\partial }^2}{\partial {\mathbf{R}}^2}+V\left(\mathbf{R}\right)\right]{\psi }_t\left(\mathbf{R},{\mathbf{R}}^'\right)</math>
אז <math>\frac{i\hslash }{t}</math> הפך ל-<math>\mathrm{-}N</math>, <math>\frac{{\hslash }^2}{2m}</math> הפך ל-<math>\frac{l^2}{6}</math>, ו-<math>V\left(\mathbf{R}\right)</math> הפך ל-<math>\frac{1}{k_BT}U\left(\mathbf{R}\right)</math>. לכן, בפולימרים נקבל [[דעיכה מעריכית|פתרונות דועכים]] ולא [[גל]]יים (כי הזמן [[מספר מדומה|מדומה]] כבר), ומשום שזו עדיין [[תורת שטורם-ליוביל|בעיית שטורם ליוביל]], עדיין מדובר בפתרונות עם [[הפרדת משתנים]] (בין "זמן" ומרחב, בדיוק איך שפותרים את משוואת שרדינגר כאשר ה[[המילטוניאן]] ב"ת בזמן), כאשר סט הפתרונות המרחביים הוא פתרון [[פונקציה עצמית|בעיית הערכים העצמיים]] של ההמילטוניאן (שהוא גם פה [[אופרטור הרמיטי]]):
:<math>\left[-\frac{l^2}{6}\frac{{\partial }^2}{\partial {\mathbf{R}}^2}+\frac{1}{k_BT}U\left(\mathbf{R}\right)\right]{\phi }_k=E_k{\phi }_k</math>
:<math>\int{d^3r}\ {\phi }^*_k{\phi }_{k^'}={\delta }_{k,k^'}</math>
:<math>\sum^{\infty }_{k=1}{{\phi }^*_k\left(\mathbf{r}\right){\phi }_k\left({\mathbf{r}}^'\right)}=\delta \left(\mathbf{r}-{\mathbf{r}}^'\right)</math>
:<math>P_N\left(\mathbf{R},{\mathbf{R}}^'\right)=\sum^{\infty }_{k=1}{{\phi }^*_k\left(\mathbf{R}\right){\phi }_k\left({\mathbf{R}}^'\right)e^{-NE_k}}</math>
===דוגמה: שרשרת בקופסה===
ניקח פולימר בקופסה בעלת ממדים <math>L_x,L_y,L_z</math>. הפוטנציאל הוא אינסופי מחוץ לקופסה, ואפס בתוכה. לכן, וכמו בבעיית [[בור פוטנציאל אינסופי|חלקיק בקופסה]], זה מתבטא בתנאי השפה ש-<math>G\left(\mathbf{R},{\mathbf{R}}^';N\right)=0</math> על השפה (ומחוצה לה גם אפס), והיא (פונקציית המשקל הסטטיסטי או [[פונקציית גרין]] איך שקראנו לה גם) מקיימת בתוך הקופסה:
:<math>-\frac{\partial G\left(\mathbf{R},{\mathbf{R}}^';N\right)}{\partial N}=-\frac{l^2}{6}\frac{{\partial }^2}{\partial {\mathbf{R}}^2}G\left(\mathbf{R},{\mathbf{R}}^';N\right)</math>
והפתרון הוא הפתרון המוכר לבעיית [[בור פוטנציאל אינסופי|החלקיק בקופסה]] שמגיע מ[[הפרדת משתנים]] (ועם החלפת <math>\frac{i\hslash }{t}</math> ב-<math>\mathrm{-}N</math> ו-<math>\frac{{\hslash }^2}{2m}</math> ב-<math>\frac{l^2}{6}</math>):
:<math>P_N\left(\mathbf{R},{\mathbf{R}}^'\right)=G\left(\mathbf{R},{\mathbf{R}}^';N\right)=\prod_{i\in \left\{x,y,z\right\}}{g_i\left(R_i,R^'_i;N\right)}</math>
כאשר
:<math>g_i\left(R_i,R^'_i;N\right)=\frac{2}{L_i}\sum^{\infty }_{k=1}{{\mathrm{sin} \left(\frac{k\pi R_i}{L_i}\right)\ }{\mathrm{sin} \left(\frac{k\pi R^'_i}{L_i}\right)\ }e^{-\frac{k^2{\pi }^2Nl^2}{6L^2_i}}}</math>
ואז [[פונקציית חלוקה (פיזיקה)|פונקציית החלוקה]] היא:
:<math>Z=\int{d^3R\ \int{d^3R^'\ \prod_{i\in \left\{x,y,z\right\}}{g_i\left(R_i,R^'_i;N\right)}}}=\prod_{i\in \left\{x,y,z\right\}}{Z_i}</math>
כאשר
:<math>Z_i=\int{dR_i\int{dR^'_i\ g_i\left(R_i,R^'_i;N\right)}}=\frac{8}{{\pi }^2}L_i\sum^{\infty }_{k=0}{\frac{1}{{\left(2k+1\right)}^2}e^{-\frac{{\pi }^2Nl^2}{6}\ \frac{{\left(2k+1\right)}^2}{L^2_i}}}</math>
ו[[אנרגיה חופשית של הלמהולץ|האנרגיה החופשית של הלמהולץ]] היא
:<math>F=-k_BT{ln Z\ }</math>
כך שה[[לחץ]] על הדופן המאונכת לציר <math>i</math> הוא:
:<math>P_i=-\frac{L_i}{L_xL_yL_z}\frac{\partial F}{\partial L_i}=\frac{L_ik_BT}{L_xL_yL_z}\frac{1}{Z_i}\frac{\partial Z_i}{\partial L_i}</math>
נסתכל על זה בשני הגבולות הבאים:
[[קובץ:confined_polymer_a.png|ממוזער|250 פיקסלים|פולימר בקופסה, כאשר הפולימר הרבה יותר קטן מהקופסה (כלומר <math>\sqrt{N}\ l\ll L_x,L_y,L_z</math>)]]
'''(1) כאשר הפולימר הרבה יותר קטן מהקופסה (כלומר <math>\sqrt{N}\ l\ll L_x,L_y,L_z</math>):'''
<math>Z_i</math> הוא בקירוב
:<math>Z_i\cong \frac{8}{{\pi }^2}L_i\sum^{\infty }_{k=0}{\frac{1}{{\left(2k+1\right)}^2}}=L_i</math>
ואז
:<math>P_i\cong \frac{L_ik_BT}{L_xL_yL_z}\frac{1}{L_i}=\frac{k_BT}{V}</math>
וזאת משוואת המצב של [[גז אידיאלי]], כצפוי.
[[קובץ:confined_polymer_b.png|ממוזער|250 פיקסלים|פולימר בקופסה, כאשר הפולימר הרבה יותר גדול מהקופסה (כלומר <math>\sqrt{N}\ l\gg L_x,L_y,L_z</math>)]]
'''(2) כאשר הפולימר הרבה יותר גדול מהקופסה (כלומר <math>\sqrt{N}\ l\gg L_x,L_y,L_z</math>):'''
<math>Z_i</math> הוא בקירוב
:<math>Z_i\cong \frac{8}{{\pi }^2}L_ie^{-\frac{{\pi }^2Nl^2}{6L^2_i}}</math>
ואז
:<math>P_i\cong \frac{L_ik_BT}{L_xL_yL_z}\frac{e^{-\frac{{\pi }^2Nl^2}{6L^2_i}}+L_ie^{-\frac{{\pi }^2Nl^2}{6L^2_i}}\frac{{\pi }^2Nl^2}{6}\frac{2}{L^3_i}}{L_ie^{-\frac{{\pi }^2Nl^2}{6L^2_i}}}=\frac{k_BT}{L_xL_yL_z}\left(1+\frac{{\pi }^2Nl^2}{3L^2_i}\right)\cong \frac{{\pi }^2Nl^2}{3L^2_i}\frac{k_BT}{V}</math>
ועקרונית פה, הלחץ על הדופן בציר מסוים אינו שווה ללחץ על דופן בציר אחר (אלא אם כן ממדי הקופסה שווים בצירים האלה). [[אנאיזוטרופיות]] כזו אחראית על כל מיני תכונות מיוחדות לפולימרים.
===עוד יישומים===
השיטה הזו של [[פונקציית גרין]] מיושמת בהרבה בעיות שקשורות לפולימרים במצבים לא [[הומוגניות (פיזיקה)|הומוגניים]] כגון ספיחה של פולימרים, פולימרים בקרבת [[מעבר פאזה|מעברי פאזה]], פולימרים בתווך [[נקבוביות|נקבובי]], ועוד.
==אפקט הנפח האסור==
===הקדמה===
במודלים של הפולימרים שהזכרנו עד כה ה[[אינטראקציה|אינטראקציות]] בין קטעי הפולימר היו רק אינטראקציות בטווח קצר, ולא התחשבנו כלל באינטגרציות לטווח ארוך (קצר וארוך פה מדברות על המרחק לאורך השרשרת לפי [[תיאור פרמטרי של עקום|משתנה הפרמטריזציה]] <math>n</math>, ולא על המרחק המרחבי). בפועל, יש אינטראקציה כזו כאשר חלקים רחוקים על הפולימר הם קרובים מרחבית. אינטראקציה ברורה היא [[אפקט סטרי|האפקט הסטרי]], הרי משום שכל מקטע בפולימר תופס [[נפח]] במרחב, מקטעים אחרים לא יכולים לקחת אותו. אינטראקציה זו מנפחת את הפולימרים (ביחס לזה שאם חושבים עליהם בלי האינטראקציה), ואז האורך הטיפוסי של הפולימרים יהיה יותר גדול. גם אם ישנם כוחות משיכה (כמו [[קשרי ואן דר ואלס| ואן דר ואלס]]), כל עוד שהכוח הדוחה הזה שולט, הפולימרים מתנפחים. אפקט זה נקרא אפקט הנפח האסור.
[[קובץ:excluded_volume.png|ממוזער|250 פיקסלים|אינטראקציית הנפח האסור]]
אפקט זה נחקר בהתחלה על ידי קהן, והייתה לו התקדמות מודרנית שהתחילה עם פלורי. התקבלה התובנה אז שהאינטראקציה ארוכת הטווח משנה את התכונות הסטטיסטיות של הפולימר בצורה משמעותית. למשל, האורך הטיפוסי של הפולימר לא הולך כמו <math>\sqrt{N}</math>, אלא כמו <math>N^{\nu }</math>, כאשר <math>\nu </math> הוא <math>\frac{3}{5}</math> בערך, כך שאפקט הנפח האסור הוא מאוד חשוב עבור שרשראות ארוכות.
כאשר כוללים את האינטראקציה ארוכת הטווח בחשבון, חישוב מדויק נהיה בלתי-אפשרי. ולכן, משתמשים בכל מיני גישות כדי לעסוק בבעיה זו, ואנחנו נדבר פה על הגישה של [[פונקציית גרין]] שהגדרנו, כאשר ההגדרה שלה כללה [[אינטגרלי מסלול]] כזכור- בעצם זו הגישה כאמור.
===מודל השרשרת עם הנפח האסור===
בפועל, האינטראקציה ארוכת הטווח היא מסובכת- יש בה [[אפקט סטרי|אפקטים סטריים]], [[קשרי ואן דר ואלס|משיכות ואן דר ואלס]], וכל מיני אינטראקציות אחרות המועברות על ידי [[תווך]] ה[[ממס]] (ו[[מולקולה|מולקולותיו]]). בכל אופן, כל עוד מדובר בחקר תכונות בסקלה גדולה, הפרטים של האינטראקציה לא כל כך משנים, משום שאפקט הנפח האסור נשלט על ידי האינטראקציה בין חלקים רחוקים של הפולימר.
ניתן אז למדל את האינטראקציה הזו בין המקטע (הנקודתי) ה-<math>n</math> לבין המקטע ה-<math>m</math> על ידי <math>k_BT\ \tilde{v}\left({\mathbf{R}}_n-{\mathbf{R}}_m\right)</math>, כאשר <math>\tilde{v}</math> היא פונקציית [[אינטראקציה]] קצרת טווח שבדרך כלל מקורבת ל[[פונקציית דלתא של דיראק|פונקציית דלתא]] <math>\tilde{v}\left({\mathbf{R}}_n-{\mathbf{R}}_m\right)=v\ \delta \left({\mathbf{R}}_n-{\mathbf{R}}_m\right)</math> (<math>v</math> הוא הנפח האסור). פונקציית הדלתא לא תגרום לבעיות, כי גדלים מדידים מערבים במקרה הזה אינטגרלים שמתנהגים כראוי. ה[[אנטרופיה]] אכן מתבדרת, אך הפרש האנטרופיה לא, וזה הדבר המעניין פיזיקלית ובעל משמעות. דבר כזה יש לו הרבה מקרים מקבילים או דומים ב[[תורת השדות הקוונטית|תורת השדות הקוונטיים]].
אנרגיית האינטראקציה הכוללת אם כן ניתנת לתיאור ע"י
:<math>U_1=\frac{1}{2}vk_BT\int^N_0{dn\int^N_0{dm\ }\delta \left({\mathbf{R}}_n-{\mathbf{R}}_m\right)}</math>
ואם נסמן את [[ריכוז (כימיה)|הריכוז]] הלוקאלי של המקטעים ב-<math>c\left(\mathbf{r}\right)</math> (כלומר <math>c\left(\mathbf{r}\right)=\int^N_0{dn\ \delta \left(\mathbf{r}-{\mathbf{R}}_n\right)}</math>), נרשום את אנרגיית האינטראקציה בצורה הזאת:
:<math>U_1=\int{d^3r\ \frac{1}{2}vk_BT\ c^2\left(\mathbf{r}\right)}</math>
כאשר צורה זו של כתיבה מראה שהביטוי שרשמנו לאנרגיה הוא הראשון ב[[משוואת מצב|פיתוח הויריאלי]] ביחס לריכוז הלוקאלי <math>c\left(\mathbf{r}\right)</math>. לכן, הפרמטר <math>v</math> נחשב למקדם הויריאלי בין המקטעים.
עקרונית, [[משוואת מצב|הפיתוח הויריאלי]] יכול להמשיך ולכלול עוד חזקות יותר גבוהות של <math>c\left(\mathbf{r}\right)</math>:
:<math>U_1=\int{d^3r\ \left[\frac{1}{2}vk_BT\ c^2\left(\mathbf{r}\right)+\frac{1}{6}wk_BT\ c^3\left(\mathbf{r}\right)+\ .\ .\ .\right]}</math>
אבל למרות זאת, האיברים הבאים זניחים, כי צפיפות המקטעים קטנה בפולימר, והיא מוערכת להיות:
:<math>\overline{c}\cong \frac{N}{{\overline{R}}^3}\propto N^{1-3\nu }=N^{-\frac{4}{5}}</math>
(כאשר לקחנו
<math>\nu =0\mathrm{.6}</math>, שמתקבל מחישובים בהסתמך רק על האיבר הראשון), וזה מאוד קטן כאשר <math>N</math> מאוד גדול. לכן, התכונות העיקריות של אפקט הנפח האסור יכולות להתקבל בצורה מספקת מתוך החישובים שמתחשבים בצורת האנרגיה שרשמנו בהתחלה כדי לתאר את תרומת האפקט הזה:
:<math>U_1=\frac{1}{2}vk_BT\int^N_0{dn\int^N_0{dm\ }\delta \left({\mathbf{R}}_n-{\mathbf{R}}_m\right)}</math>
ואז פונקציית הפילוג של [[קונפורמציה (כימיה)|הקונפורמציה]] <math>{\mathbf{R}}_n</math> (בהתחשב ב[[אינטראקציה]] זו, ובהיעדר פוטנציאל חיצוני) היא:
:<math>\mathit{\Psi}\left({\mathbf{R}}_n\right)=C{exp \left[-\frac{3}{2l^2}\int^N_0{dn{\left(\frac{\partial {\mathbf{R}}_n}{\partial n}\right)}^2}-\frac{1}{2}v\int^N_0{dn\int^N_0{dm\ }\delta \left({\mathbf{R}}_n-{\mathbf{R}}_m\right)}\right]\ }</math>
ולמודל הזה יש רק שני פרמטרים: <math>l</math> שמייצג את האינטראקציה קצרת הטווח, ו-<math>v</math> שמייצג את האינטראקציה ארוכת הטווח. ההנחה הבסיסית למודל זה היא ההפרדה החדה בין הטווחים של אינטראקציות אלו, ואף על פי שאפשר לדון בתקפות של הנחה זו, בדרך כלל מאמינים שהיא נקודת מוצא נכונה לחקר הבעיה של הנפח האסור.
===כיוונים תאורטיים===
נדבר פה בקצרה על כמה גישות תאורטיות לטיפול בבעיה (מבחינה פיזיקלית-סטטיסטית), ונסביר יותר על החלק הקשור ב[[אינטגרלי מסלול|אינטגרל המסלולי]]. נגביל את הדיון ל-<math>v>0</math> (אינטראקציה דוחה). מקרים אחרים גם ניתנים לטיפול.
====תורה פשוטה של פלורי ושדרוג לתורת שדה ממוצע====
לפלורי היה רעיון פשוט כדי לחשב את הנפח של פולימר: להסתכל על האיזון בין שני האפקטים: האפקט הדוחה של הנפח האסור שגורם לפולימר להתנפח, והאפקט [[אלסטיות|האלסטי]] שנובע מההתקשרות בין [[מונומר|המונומרים]] לאורך השרשרת שרוצה מהפולימר להתכווץ. אפשר להסתכל על מודל פשוט אז של [[אנרגיה חופשית של הלמהולץ|האנרגיה החופשית של הלמהולץ]] ולהגיע לצורה הזו:
:<math>F\cong k_BT\left(\frac{3R^2}{2Nl^2}+v\frac{N^2}{R^3}\right)</math>
ואז הנפח הממוצע של הפולימר מוערך להיגזר מ-<math>R</math> שמביא את <math>F</math> למינימום:
:<math>\overline{R}=\sqrt{N}l{\left(\frac{v\sqrt{N}}{l^3}\right)}^{\frac{1}{5}}\propto N^{\frac{3}{5}}</math>
וזה נותן את האקספוננט <math>\nu </math> שדיברנו עליו. ואף על פי ש-<math>\nu </math> שיוצא הוא קרוב לערך האמפירי, אלא שעוד תחזיות של גדלים אחרים מתורה זו לא מספקות. למשל, הביטוי שרשמנו לאנרגיה החופשית אינו עקבי עם צפיפות ההסתברות לווקטור הקצה-קצה שמתקבלת מסימולציות במחשב. גם יש התנהגות מוזרה של ה[[אנטרופיה]]. תורה זו בכל מקרה נחשבת למין [[תורת שדה ממוצע]], אבל לא ממש, כי לא מחשבים בה שום ממוצע לשדה. תורת השדה הממוצע תתואר בהמשך.
====חישוב הפרעתי והיעזרות בפונקציית גרין שמוגדרת על ידי אינטגרלים מסלוליים====
בהנחה שהנפח האסור <math>v</math> הוא קטן, פונקציית פילוג הקונפורמציה ניתנת לפיתוח כ[[טור חזקות]] של <math>v</math>, ואז נוכל לבצע [[תורת ההפרעות (מכניקת הקוונטים)|חישוב הפרעתי]] של הגודל <math>\left\langle {\mathbf{R}}^2\right\rangle </math>. נראה כעת חישוב מפורט לדוגמה של איך מחשבים את התיקון הראשון.
יהי <math>{\left\langle \ \cdot \ \right\rangle }_0</math> מיצוע לפי הפילוג של [[שרשרת פולימר אידיאלית|השרשרת האידיאלית]]. אז עבור השרשרת עם הנפח האסור:
:<math>\left\langle {\mathbf{R}}^2\right\rangle =\frac{{\left\langle {\mathrm{exp} \left(-U_1\right)\ }{\left({\mathbf{R}}_N-{\mathbf{R}}_0\right)}^2\right\rangle }_0}{{\left\langle {\mathrm{exp} \left(-U_1\right)\ }\right\rangle }_0}</math>
כאשר פה (ניקח את <math>k_BT</math> להיות אחד):
:<math>U_1=v\int^N_0{dm\int^m_0{dn\ }\delta \left({\mathbf{R}}_n-{\mathbf{R}}_m\right)}</math>
ומפיתוח עד הסדר הראשון ב-<math>U_1</math>:
:<math>\left\langle {\mathbf{R}}^2\right\rangle \cong \frac{{\left\langle \left(1-U_1\right){\left({\mathbf{R}}_N-{\mathbf{R}}_0\right)}^2\right\rangle }_0}{{\left\langle 1-U_1\right\rangle }_0}\cong {\left\langle {\left({\mathbf{R}}_N-{\mathbf{R}}_0\right)}^2\right\rangle }_0\left(1+{\left\langle U_1\right\rangle }_0\right)-{\left\langle U_1{\left({\mathbf{R}}_N-{\mathbf{R}}_0\right)}^2\right\rangle }_0=</math>
:<math>=Nl^2\left(1+v\int^N_0{dm\int^m_0{dn\ {\left\langle \delta \left({\mathbf{R}}_n-{\mathbf{R}}_m\right)\right\rangle }_0}}\right)-v\int^N_0{dm\int^m_0{dn\ {\left\langle {\left({\mathbf{R}}_N-{\mathbf{R}}_0\right)}^2\ \delta \left({\mathbf{R}}_n-{\mathbf{R}}_m\right)\right\rangle }_0}}</math>
ועכשיו נחשב איברים בנפרד:
:<math>{\left\langle \delta \left({\mathbf{R}}_n-{\mathbf{R}}_m\right)\right\rangle }_0=\int{d^3R_N\ d^3R_m\ d^3R_n\ G_0\left({\mathbf{R}}_N-{\mathbf{R}}_m;N-m\right) G_0\left({\mathbf{R}}_m-{\mathbf{R}}_n;m-n\right)G_0\left({\mathbf{R}}_n-{\mathbf{R}}_0;n\right)\ \delta \left({\mathbf{R}}_n-{\mathbf{R}}_m\right)}=</math>
:<math>=\int{d^3R_N\ d^3R_n\ G_0\left({\mathbf{R}}_N-{\mathbf{R}}_n;N-m\right)G_0\left(0;m-n\right)G_0\left({\mathbf{R}}_n-{\mathbf{R}}_0;n\right)}=G_0\left(0;m-n\right)</math>
:<math>{\left\langle {\left({\mathbf{R}}_N-{\mathbf{R}}_0\right)}^2\ \delta \left({\mathbf{R}}_n-{\mathbf{R}}_m\right)\right\rangle }_0=</math>
:<math>=\int{d^3R_N\ d^3R_m\ d^3R_n\ G_0\left({\mathbf{R}}_N-{\mathbf{R}}_m;N-m\right) G_0\left({\mathbf{R}}_m-{\mathbf{R}}_n;m-n\right)G_0\left({\mathbf{R}}_n-{\mathbf{R}}_0;n\right){\left({\mathbf{R}}_N-{\mathbf{R}}_0\right)}^2\ \delta \left({\mathbf{R}}_n-{\mathbf{R}}_m\right)}=</math>
:<math>=\int{d^3R_N\ d^3R_n\ G_0\left({\mathbf{R}}_N-{\mathbf{R}}_n;N-m\right) G_0\left(0;m-n\right)G_0\left({\mathbf{R}}_n-{\mathbf{R}}_0;n\right){\left({\mathbf{R}}_N-{\mathbf{R}}_0\right)}^2}=</math>
:<math>=\int{d^3R_N\ d^3R_n\ G_0\left({\mathbf{R}}_N-{\mathbf{R}}_n;N-m\right) G_0\left({\mathbf{R}}_n-{\mathbf{R}}_0;n\right){\left[\left({\mathbf{R}}_N-{\mathbf{R}}_n\right)+\left({\mathbf{R}}_n-{\mathbf{R}}_0\right)\right]}^2}=</math>
:<math>=G_0\left(0;m-n\right)\int{d^3R_n\ G_0\left({\mathbf{R}}_n-{\mathbf{R}}_0;n\right)\left[\left(N-m\right)l^2+{\left({\mathbf{R}}_n-{\mathbf{R}}_0\right)}^2\right]}=\left(N-m+n\right)l^2\ G_0\left(0;m-n\right)</math>
ואז עד הסדר הראשון ב-<math>v</math> נקבל:
:<math>\left\langle {\mathbf{R}}^2\right\rangle \cong Nl^2\left(1+v\int^N_0{dm\int^m_0{dn\ G_0\left(0;m-n\right)}}\right)-vl^2\int^N_0{dm\int^m_0{dn\ \left(N-m+n\right)G_0\left(0;m-n\right)}}=</math>
:<math>=Nl^2+vl^2\int^N_0{dm\int^m_0{dn\ \left(m-n\right)G_0\left(0;m-n\right)}}=Nl^2+vl^2\int^N_0{dm\int^m_0{dn\ \left(m-n\right){\left(\frac{3}{2\pi \left(m-n\right)l^2}\right)}^{\frac{3}{2}}}}=</math>
:<math>=Nl^2\left[1+\frac{v}{l^3}{\left(\frac{3}{2\pi }\right)}^{\frac{3}{2}}\int^N_0{dm\int^m_0{dn\ {\left(m-n\right)}^{-\frac{1}{2}}}}\right]=Nl^2\left[1+\frac{1}{N}\frac{v}{l^3}{\left(\frac{3}{2\pi }\right)}^{\frac{3}{2}}\int^N_0{dm\ 2\sqrt{m}}\right]=</math>
:<math>=Nl^2\left[1+\frac{1}{N}\frac{v}{l^3}{\left(\frac{3}{2\pi }\right)}^{\frac{3}{2}}\frac{4}{3}N^{\frac{3}{2}}\right]=Nl^2\left[1+\frac{4}{3}{\left(\frac{3}{2\pi }\right)}^{\frac{3}{2}}\frac{v\sqrt{N}}{l^3}\right]\equiv Nl^2\left(1+\frac{4}{3}z\right)</math>
כאשר מגדירים <math>z\equiv {\left(\frac{3}{2\pi }\right)}^{\frac{3}{2}}\frac{v\sqrt{N}}{l^3}</math>.
מפיתוח דומה וארוך ניתן למצוא את הסדרים הבאים, ומתקבל בסוף:
:<math>\frac{\left\langle {\mathbf{R}}^2\right\rangle }{Nl^2}=1+\frac{4}{3}z-2.075z^2+6.297z^3-25.057z^4+116.135z^5-594.717z^6+\ .\ .\ .</math>
נשים לב שבפועל פרמטר הפיתוח הוא <math>z</math> (ולא <math>v</math>), ו-<math>z</math> הוא פרופורציונלי ל-<math>\sqrt{N}</math>, כך ש[[תורת ההפרעות (מכניקת הקוונטים)|הפיתוח ההפרעתי]] נהיה חסר משמעות כאשר <math>N</math> גדול. המצב הזה שונה ממש מ[[משוואת מצב|הפיתוח הויריאלי]] של גז לא מושלם, כאשר שם פרמטר הפיתוח הוא הריכוז הממוצע <math>\overline{c}</math> שפרופורציוני ל-<math>N^{1-3\nu }</math> שהוא קטן בפולימרים. הסיבה להבדל זה הוא שבפולימרים ההתנגשות בין שני מקטעים משפיעה על וקטור הקצה-קצה בצורה משמעותית (אפקט מסדר גודל של 1), כאשר בגזים לא מושלמים האפקט הוא בסדר גודל של <math>N^{-1}</math>. [[תורת ההפרעות (מכניקת הקוונטים)|הפיתוח ההפרעתי]] במקרה של בעיית הפולימרים מנומק רק כאשר המקטעים של הפולימר אינם סבירים להתנגש זה בזה. ולכן מה שצריך להיות קטן בעצם הוא המספר הממוצע של ההתנגשויות שקורות לאורך כל הפולימר:
:<math>Nv\overline{c}\cong Nv\frac{N}{{\left(\sqrt{N}l\right)}^3}=\frac{v\sqrt{N}}{l^3}</math>
וזה בעצם <math>z</math> (עד כדי קבוע).
אפשר להראות ש[[טור חזקות|הטור]] ([[תורת ההפרעות (מכניקת הקוונטים)|ההפרעתי]]) שרשמנו עבור <math>\frac{\left\langle {\mathbf{R}}^2\right\rangle }{Nl^2}</math> הוא אסימפטוטי, ושיש עלייה חדה במקדמיו שהולכת כמו <math>n^n</math>. בכל מקרה, על ידי טכניקות סכימה מסוימות ניתן להפיק אינפורמציה מספקת עבור ההתנהגות האסימפטוטית. בפרט, האקספוננט <math>\nu </math> מוערך להיות
<math>\nu =0.588\pm 0.001</math>
====מודל הפיתוח האחיד====
הקושי בחישובים עם [[תורת ההפרעות (מכניקת הקוונטים)|הפיתוח ההפרעתי]] ניתן לעקיפה יעילה מספיק על ידי תוכנית אחרת של חישובים. לפי שיטה זו, מניחים כי הפיתוח של השרשרת מיוצג על ידי הפיתוח של אורך קשר, כלומר שצפיפות ההסתברות של הקונפורמציה <math>{\mathbf{R}}_n</math> היא:
:<math>{\mathit{\Psi}}^'\left({\mathbf{R}}_n\right)=C{exp \left[-\frac{3}{2{l^'}^2}\int^N_0{dn{\left(\frac{\partial {\mathbf{R}}_n}{\partial n}\right)}^2}\right]\ }</math>
כאשר <math>l^'</math> הוא אורך קשר שיפותח. הממוצע <math>\left\langle {\mathbf{R}}^2\right\rangle </math> של הקונפורמציה
:<math>\mathit{\Psi}\left({\mathbf{R}}_n\right)=C{exp \left[-\frac{3}{2l^2}\int^N_0{dn{\left(\frac{\partial {\mathbf{R}}_n}{\partial n}\right)}^2}-\frac{1}{2}v\int^N_0{dn\int^N_0{dm\ }\delta \left({\mathbf{R}}_n-{\mathbf{R}}_m\right)}\right]\ }</math>
יכול להיכתב גם ככה:
:<math>\left\langle {\mathbf{R}}^2\right\rangle =\frac{{\left\langle {\mathrm{exp} \left(-B\left[{\mathbf{R}}_n\right]\right)\ }{\left({\mathbf{R}}_N-{\mathbf{R}}_0\right)}^2\right\rangle }^'}{{\left\langle {\mathrm{exp} \left(-B\left[{\mathbf{R}}_n\right]\right)\ }\right\rangle }^'}</math>
כאשר <math>B\left[{\mathbf{R}}_n\right]</math> הוא הפונקציונל:
:<math>B\left[{\mathbf{R}}_n\right]=\frac{3}{2}\left(\frac{1}{l^2}-\frac{1}{{l^'}^2}\right)\int^N_0{dn{\left(\frac{\partial {\mathbf{R}}_n}{\partial n}\right)}^2}+\frac{1}{2}v\int^N_0{dn\int^N_0{dm\ \delta \left({\mathbf{R}}_n-{\mathbf{R}}_m\right)}}</math>
ואם ההנחה שלנו היא קירוב טוב, אז <math>B</math> הוא קטן, ולכן:
:<math>\left\langle {\mathbf{R}}^2\right\rangle \cong \frac{{\left\langle \left(1-B\left[{\mathbf{R}}_n\right]\right){\left({\mathbf{R}}_N-{\mathbf{R}}_0\right)}^2\right\rangle }^'}{{\left\langle 1-B\left[{\mathbf{R}}_n\right]\right\rangle }^'}</math>
ומחישוב דומה למה שעשינו בחלק הקודם, נקבל:
:<math>\left\langle {\mathbf{R}}^2\right\rangle =N{l^'}^2+N{l^'}^2\left(\frac{{l^'}^2}{l^2}-1-\frac{4}{3}\frac{l^3}{{l^'}^3}z\right)</math>
ואם ההנחה שרשמנו היא הטובה ביותר (<math>l^'</math> הוא הטוב ביותר), אז התיקון מסדר ראשון חייב להתאפס, כלומר:
:<math>\frac{{l^'}^2}{l^2}-1-\frac{4}{3}\frac{l^3}{{l^'}^3}z=0</math>

ואז אם נגדיר את מקדם הניפוח <math>{\alpha }^2\equiv \frac{\left\langle {\mathbf{R}}^2\right\rangle }{Nl^2}</math>, נקבל <math>\alpha =\frac{l^'}{l}</math>, ומתקיים:

:<math>{\alpha }^5-{\alpha }^3=\frac{4}{3}z</math>
כלומר, לסיכום קיבלנו כי:
:<math>\left\langle {\mathbf{R}}^2\right\rangle =N{\alpha }^2l^2</math>
כאשר <math>\alpha </math> מקיים:
:<math>{\alpha }^5-{\alpha }^3=\frac{4}{3}z</math>
ועבור <math>z</math> קטן התאוריה הזו נותנת תשובה כמו [[תורת ההפרעות (מכניקת הקוונטים)|הפיתוח ההפרעתי]] שעשינו. מצד שני, עבור <math>z</math> גדול התורה הזו נותנת <math>\overline{R}\propto N^{\frac{3}{5}}</math>- מה שמקבלים מתורת פלורי.
תאוריה זו נותנת תחזיות לא טובות עבור גדלים מסוימים- בכלל (למשל) היא מניחה ש-<math>\mathbf{R}</math> הוא גאוסי. בדרך כלל, הבחירה הטובה ביותר עבור <math>l^'</math> תלויה בגודל אותו חוקרים בשיטה דומה (התאפסות של התיקון הראשון).
====תורת השדה הממוצע====
אפשר לשדרג את תורת פלורי אל [[תורת שדה ממוצע]] על ידי בחינת אנסמבל השרשראות שמתחילות ומסתיימות בנקודות מקובעות במרחב, ואז הצפיפות הממוצעת של מקטע תהיה:
:<math>\overline{c}\left(\mathbf{r}\right)=\int^N_0{dn\ \left\langle \delta \left(\mathbf{r}-{\mathbf{R}}_n\right)\right\rangle }</math>
ואז המקטע ב-<math>\mathbf{r}</math> מרגיש פוטנציאל של <math>\overline{c}\left(\mathbf{r}\right)vk_BT</math> בממוצע, ופילוג של הקונפורמציה נתון ע"י
:<math>{\mathit{\Psi}}_{MF}\left({\mathbf{R}}_n\right)=C{exp \left[-\frac{3}{2l^2}\int^N_0{dn{\left(\frac{\partial {\mathbf{R}}_n}{\partial n}\right)}^2}-v\int^N_0{dn\ \overline{c}\left({\mathbf{R}}_n\right)}\right]\ }</math>
ו[[פונקציית גרין]] אז מקיימת
:<math>\left[\frac{\partial }{\partial n}-\frac{l^2}{6}\frac{{\partial }^2}{\partial {\mathbf{R}}^2}+v\overline{c}\left(\mathbf{R}\right)\right]G\left(\mathbf{R},0;n\right)=\delta \left(\mathbf{R}\right)\delta \left(N\right)</math>
ובהינתן <math>G\left(\mathbf{R},0;n\right)</math>, הצפיפות הממוצעת של המקטעים מחושבת כך:
:<math>\overline{c}\left(\mathbf{r}\right)=\frac{1}{G\left(\mathbf{R},0;N\right)}\int^N_0{dn\ G\left(\mathbf{R},\mathbf{r};N-n\right)G\left(\mathbf{r},0;n\right)}</math>
ושתי המשוואות האחרונות ביחד נותנות משוואה סגורה על <math>G\left(\mathbf{R},0;N\right)</math>.
תוצאות שמתקבלות מחישובים ע"פ תורה זו נותנות תחזיות דיי מעניינות לשרשראות עם נפח אסור. למשל, מקדם המבנה <math>g\left(\mathbf{k}\right)</math> באזור <math>k</math>-ים גדולים (<math>kR_g\gg 1</math>, כאשר <math>R_g</math> הוא אורך הג'ירציה או רדיוס ההתמדה) מתנהג כמו <math>k^{-\frac{5}{3}}</math>, ותוצאה זו גם מתקבלת מהמודל של הפיתוח האחיד.
למודל זה נוספו פלוקטואציות על ידי קוסמס ופריד, אבל התכונות המהותיות שאפשר לקבל אותן בסופו של דבר לא שונות עקרונית ממה שמקבלים ממודל זה.
====חבורת הרנורמליזציה====
[[תורת שדה ממוצע|תורת השדה הממוצע]] שהצגנו אינה מתחשבת בעובדה שההתנגשויות בין מקטעי הפולימר הן בעצם קורלטיביות: הרי הסגמנטים קשורים זה לזה, וכשיש התנגשות, זה מעודד התנגשות אחרת. כדי להתחשב בקורלציות כאלו מיישמים את טכניקת [[רנורמליזציה|חבורת הרנורמליזציה]] על פולימרים.
תוצאות ברורות ומעניינות ניתנו על ידי וילסון ופישר בהסתמך על ההבחנה שבמרחב עם מימד יותר גבוה, האפקט של הקורלציות נהיה יותר חלש, ואז [[תורת ההפרעות (מכניקת הקוונטים)|הפיתוח ההפרעתי]] הפשוט נהיה יעיל בתוצאותיו. למשל, בבעיית הנפח האסור, אם הפולימר הוא במרחב <math>D</math>-ממדי (אוקלידי), פרמטר הפיתוח <math>z</math> יהיה פרופורציונלי ל-<math>z\sim \frac{N^2v}{{\left(\sqrt{N}l\right)}^D}\propto N^{\frac{4-D}{2}}</math>, שנהיה קטן אם <math>D>4</math> (המימד 4 נקרא המימד הקריטי). ואז אפשר לתת פיתוח סביב <math>\epsilon =4-D</math> עבור <math>\nu </math>, ופיתוח כזה תוצאתו היא:
:<math>\nu =\frac{1}{2}\left(1+\frac{1}{8}\epsilon +\frac{15}{256}{\epsilon }^2+\ .\ .\ .\right)</math>
ועבור <math>D=3</math> זה ייתן <math>\nu =0.592</math>, שהוא מספיק קרוב לערך [[תורת ההפרעות (מכניקת הקוונטים)|ההפרעתי]] המדויק שהראינו קודם (<math>\nu =0.588\pm 0.001</math>). חישובים [[רנורמליזציה|בתורת חבורת הרנורמליזציה]] לפולימרים (ובכלל) הם מסובכים יחסית, וניתן למצוא אותם בספרות.
==ראו גם==
* [[פולימר]]
* [[שרשרת פולימר אידיאלית]]
* [[שרשרת פולימרית גאוסיאנית]]
* [[אינטגרלי מסלול]]
==לקריאה נוספת==
* M. Rubinstein & H. Colby, ‘’’Polymer Physics’’’, New York: 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, ‘’’The Theory of Polymer Dynamics’’’, New York: Oxford University Press; 1986
* A.Y. Grosberg & A.R. Khokholv, ‘’’Statistical physics of macromolecules’’’, American Institute of Physics; 1994
* Wax, N., ‘’’Noise and Stochastic Processes’’’. Dover Publishing Co., New York A954) represents a collection of classic papers.
* Lax, M., Rev. Mod. Phys. 32, 25 A960; 38, 541 A966.
* Kubo, R., Rep. Prog. Phys. 29, 255 A966.
* Batchelor, G. K., ‘’’An Introduction to Fluid Dynamics’’’, Chap. 4. Cambridge Univ. Press A970.
* Doi, M., Chem. Phys. 79, 5080 A983.
* van Kampen, N. G. ‘’’Stochastic Processes in Physics and Chemistry’’’. North-Holland, Amsterdam A981.
* Risken, H., ‘’’The Fokker-Planck Equation’’’. Springer, Berlin A984.
* Zwanzig, R., Adv. Chem. Phys. 15, 325 A969.
* Fixman, M., Chem. Phys. 69, 1527 A978.
* Rainville, E. D., ‘’’Special Functions’’’. Macmillan, New York A960.
* Happel, J., and Brenner, H., ‘’’Low Reynolds Number Hydrodynamics’’’. Prentice Hall, Englewood Cliffs, N. J. A965.
* Batchelor, G. K., ‘’’Fluid. Mech. ‘’’ 74, 1 A976.
* Felderhof, B. U., Physica 89A, 373 A977); /. Phys. All, 929 A978.
* Mazur, P., and van Saarloos, W., Physica 110A, 147 A982.
* Kirkwood, J. G., Rec. Trav. Chim. 68, 649 A949.
* Kirkwood, J. G., ‘’’Macromolecules, in Documents on Modem Physics’’’. Gordon & Breach, New York A967.
* Burgers, J. M., ‘’’Second Report on Viscosity and Plasticity’’’, Chap. 3. Amsterdam Academy of Sciences, North-Holland, Nordermann A938,
* Kuhn, W., and Kuhn, H., Helv. Chim. Acta 37, 97 A944; 38, 1533 A945; 39, 71 A946.
* Kramers, H. A., Chem. Phys. 14, 415 A946.
* Batchelor, G. K., ‘’’Fluid Mech. ’’’ 41, 545 A970.
* Peterson, J. M., and Fixman, M., Chem. Phys. 39, 2516 A963.
* Meyer, R. E., (ed.) ‘’’Theory of Dispersed Multiphase Flow’’’. Academic Press, New York, A983.
* Goldstein, H., ‘’’Classical Mechanics’’’. Addison-Wesley, Reading, Mass A959.
* Erpenbeck, J. J., and Kirkwood, J. G., Chem. Phys. 29, 909 A958; 38, 1023 A963.
* Fixman, M., and Kovac, J., Chem. Phys., 61, 4939 A974; ibid. 61, 4950 A974; ibid. 63, 935 A975).
* Hassager, O., Chem. Phys. 60, 2111, 4001 A974.
* Nakajima, H., and Wada, Y., ‘’’Biopolymers ’’’ 16, 875 A977; 17, 2291 A978.
* Doi, M., Nakajima, H., and Wada, Y., ‘’’Colloid Polym. Sci.’’’ 253, 905 A975; 254, 559 A976.
* Iwata, K., Chem. Phys. 71, 931 A979.
* Yamakawa, H., Ann. Rev. Phys. Chem. 35, 23 A984.
 
[[קטגוריה:פיזיקה של חומר מעובה]]
[[קטגוריה:פולימרים]]