נקע (הנדסה)

ממוזער|נקעים בסיליקון, במישור 100.
[ממוזער|תצלום ב[מיקרוסקופ אלקטרונים חודר (TEM) של נקעים בפלדה. הנקעים הם הקווים הכהים על הרקע הבהיר במרכז התמונה.]]
**נקע** (אנגלית: **Dislocation**) בהנדסת חומרים הוא סוג של פגם במבנה הגבישי של חומר. קיום הנקעים הוא פתרון לבעיות שלא ניתן להן מענה על ידי המבנה הגבישי המושלם בחומר. כך למשל מוסבר תהליך הדפורמציה הפלסטית בעזרת תנועת נקעים בתוך החומר, וכן הולכה חשמלית. תכונות אלו תלויות גם בצפיפות הנקעים בגביש הנעה בין 10^2\thicksim10^{12}cm^{-2}. את השפעת הנקעים על אלסטיות החומר פיתח בשנת 1905 ויטו וולטרה, ואת תנועת הנקעים בסדר גודל אטומי הסביר סר ג'פרי טיילור בשנת 1934.

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

### אלסטיות
על פי התורה האלסטית, חומר הנתון תחת כוח משני קצותיו יעבור עיוות, אשר נתון ביחס ישר לכוח המופעל באופן הבא:

Y{\operatorname{\delta}\!L\over\operatorname{}\!L}=\frac{F}{A}

כאשר A ו-L הם שטח הפאה ואורכו של החומר, F הכוח המופעל ו-Y מודול יאנג. בעזרת נוסחה זו ניתן להסיק על הכוח הדרוש לעיוות מסוים, אולם מדידות על חומרים הניבו ערכים קטנים מאלו התאורטיים עד כי בחלק מן החומרים השוני מגיע לסדר גודל של פי 10^4. השוואה בין מדידות נבחרות לערכן התאורטי ניתן לראות בטבלה 1. במציאות הכוח שיש להפעיל כדי לצור עיוות בחומר קטן משמעותית מן הכוח התאורטי – את העובדה הזאת אפשר להסביר על ידי קיומם של נקעים, המאפשרים עיוות של החומר על ידי תנועתם.

### אי-סדר גבישי בשני ממדים
סדר ארוך טווח של המבנה הגבישי אינו מתקיים בחומרים דו־ממדיים, כיוון שהוא מופר על ידי תנודות (פלקטואציות) תרמיות של האטומים. נגדיר את סטיית המיקום של כל נקודה בגביש על ידי השדה הווקטורי u(\vec{r}) כך שווקטור המיקום יינתן כך: \vec{r}+\vec u(\vec{r}).

מכאן שהאנרגיה הפוטנציאלית של גביש דו־ממדי היא

U=\int d^2rC\underset{\alpha \beta}{\sum}{\partial u_\alpha\over\partial r_\beta}{\partial u_\alpha\over\partial r_\beta}

C הוא קבוע בעל ממדים של אנרגיה ליחידת שטח.

על ידי התמרת פורייה של השדה הווקטורי u(\vec{r}), האנרגיה הפוטנציאלית תינתן כך:

U=\int d^2r\frac{C}{2}\underset{\alpha \beta \vec{k}\vec{ k'}}{\sum}k_\beta k'_\beta e^{i(\vec{k'}-\vec{ k})\vec{r}} u_\alpha(\vec{ k}) u^*_\alpha(\vec{ k'})=
\frac{VC}{2}\underset{\alpha \vec{k}}{\sum}k^2\left \vert u_\alpha(\vec{k}) \right \vert ^2

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

בעזרת האנרגיה הפוטנציאלית ניתן לחשב את השונות בסטיית המיקום, כלומר השונות של השדה הווקטורי \vec u(\vec{r}):

\langle u^2\rangle=\underset{\beta \vec{k}}{\sum}\langle\left \vert u_\beta(\vec{k})\right \vert ^2 \rangle

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

\langle\left \vert u_\beta(\vec{k})\right \vert ^2 \rangle = \frac{K_BT}{VCk^2}

בעזרת ביטוי זה נוכל לחשב את השונות

\langle u^2\rangle=\underset{\alpha \vec{k}}{\sum}\frac{K_BT}{VCk^2}=2\int \frac{d^2k}{(2\pi)^2}\frac{K_BT}{Ck^2}

במעבר האחרון נעשה שימוש בפונקציית צפיפות המצבים במרחב התנע עבור גוף דו־ממדי.

השדה הווקטורי \vec u(\vec{r}) מתאר תנודות הגדולות ממרחק כלשהו D אשר גדול מן המרחק הבין-אטומי. לכן מתקיים התנאי \vec u(\vec{k})=0 עבור k>\frac{1}{D}.

מכאן שהשונות תסתכם כך:

\langle u^2\rangle=2\int\limits_{0}^{\frac{1}{D}} \frac{dk}{2\pi k}\frac{K_BT}{C}\rightarrow \infty 

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

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

## המבנה הגאומטרי של נקעים
ממוזער|שני סוגי הנקעים העיקריים: נקע קצה, בשמאל, ונקע בורג, בימין. באיור משורטטים גם וקטור בורגרס, באדום, וקו הנקע, בירוק, עבור כל נקע. למעלה במרכז מוצג גביש (דו־ממדי) ללא נקע.

כדי להבין את הגאומטריה של נקע, נשתמש במודל התלת-ממדי של המבנה המסודר של גביש, כאשר האטומים מיוצגים בעזרת כדורים. עבור נקע המצוי בסריג הגבישי, מוגדר וקטור בורגרס – \vec{b} – המייצג את הנקע, ומציין את הכיוון ואת הגודל של התקדמותו בחומר. כמו כן מוגדר וקטור קו הנקע \vec{t}, המקביל לציר הנקע ומפריד בין
ממוזער|'''איור 1''' איור סכימטי של נקע קצה במישורי הגביש. וקטור קו הנקע מסומן בכחול וכיוונו לפנים האיור, ווקטור בורגרס האנכי לו מסומן בשחור.
האזור של הגביש שבו הסריג מושלם ובין האזור שהנקע משפיע עליו. כדי למצוא את וקטור בורגרס, יש להתבסס על כלל יד ימין כאשר וקטור קו הנקע משמש ציר, ולסגור סביבו לולאה בעלת מספר קבוע של אתרי סריג בכל אחת מצלעותיה (איור 3). על הלולאה להשתרע רחוק מספיק מווקטור קו הנקע כך שהאטומים שהיא מחברת נתונים בדפוס הסריג המושלם. עבור סריג אידיאלי ללא נקעים נקודות ההתחלה והסיום של הלולאה יהיו באותו אתר סריג. אך אם הלולאה מקיפה נקע, יהיה הפרש בין נקודות ההתחלה והסיום של הלולאה. הפרש זה וכיוונו מהווים את וקטור בורגרס.

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

### נקע קצה
נקע קצה (באנגלית Edge dislocation) ניתן לתיאור על ידי הוספת חצי מישור עודף של אטומים לתוך מבנה הגביש, כאשר הוא גורם לעיוות מקומי של המישורים האטומיים מסביבו (איור 1). מכיוון שבאזור הנקע יש עודף מישורי גביש, מתקיים לחץ רב באזור זה. במקרה של הפעלת לחץ חיצוני המישור העודף יתקדם דרך מישורי הגביש, כאשר הוא יוצר ומנתק קשרים עם מישורים אלו, עד שיגיע לגבול החיצוני של החומר. עבור נקע קצה וקטור בורגרס מאונך לווקטור קו הנקע (\vec{b}\perp\vec{t}) כפי שמוצג באיורים 1 ו-3.

### נקע בורג
נקע בורג (Screw dislocation) מוגדר כאשר באזור שמעבר לווקטור קו הנקע חלה החלקה של אתרי הסריג בצידו האחד של מישור המקביל לווקטור קו הנקע, ביחס לאתרי הסריג שנמצאים מעבר למישור זה. צורה נוספת לחשוב על כך היא הוספת מישור של אתרי סריג מעל לחלק מן האתרים הקיימים (איור 2). באופן זה הלולאה הדמיונית מתחילה ומסתיימת בשני מישורים שונים ומכאן "צורת הבורג". במקרה זה וקטור בורגרס מקביל לווקטור קו הנקע (\vec{b}\parallel\vec{t}) כפי שניתן לראות באיור 3.
ממוזער|'''איור 3''' בחלק העליון סריג ובו נקע קצה אל מול סריג מושלם. וקטור בורגרס מסומן, וכן וקטור קו הנקע בצורת T, וכיוונו אל פנים מישור התמונה. בחלק התחתון סריג ובו נקע בורג אל מול סריג מושלם. וקטור בורגרס ווקטור קו הנקע מקבילים זה לזה.

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

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

כדי למצוא כוח זה, נניח מודל של סריג חד־ממדי – שרשרת של מסות זהות המחוברות בקפיצים בעלי קבוע קפיץ k והנתונות בפוטנציאל פרבולי מחזורי בעל קבוע קפיץ \Kappa (איור 4). על השרשרת מופעל כוח חיצוני f כך שהפוטנציאל הפרבולי המחזורי הוא

U(x)=\frac {1}{2}\Kappa\left [ x-a\cdot int(x/a) \right ]^2-fx

a הוא המרחק בין המסות בעת שיווי משקל. הכוח הפועל על המסה n הוא

f_n=k\left [ x_{n+1}-x_n-a \right ]+k\left [ x_{n-1}-x_n+a \right ]-{\operatorname{d}\!U\over\operatorname{d}\!x}

כעת בהינתן נקע בחומר, אחת המסות חסרה. ללא הגבלת הכלליות ניתן להניח כי בור הפוטנציאל במיקום n=0 ריק, כך שהמסה 0 נתונה בבור -a והמסה 1 נתונה בבור a. מכאן שהכוח על מסה n הוא

f_n=k\left [ x_{n+1}-x_n-a \right ]+k\left [ x_{n-1}-x_n+a \right ]+f-\Kappa\left [ x_{n}-a(n-1) \right ] ; n\leq0

f_n=k\left [ x_{n+1}-x_n-a \right ]+k\left [ x_{n-1}-x_n+a \right ]+f-\Kappa\left [ x_{n}-na \right ] ; n>0

בשיווי משקל הכוח השקול מתאפס כך שהפתרונות לשתי המשוואות הם
x_n=\frac{f}{\Kappa}+a(n-1)+A_le^{qn} ; n\leq 0

x_n=\frac{f}{\Kappa}+an+A_re^{-qn} ; n>0

על ידי הצבת הפתרונות במשוואות הכוחות, וכן השוואת שני הפתרונות עבור המיקומים n=0,1 נקבל
ממוזער|486px|'''איור 4''' מערך חד־ממדי של מסות המחוברות בקפיצים, הנתונות תחת פוטנציאל מחזורי חיצוני המשורטט מתחתיהן. המסה במיקום 0 חסרה, כך שקיים נקע במערך.
\frac{\Kappa}{k}=e^q+e^{-q}-2

A_l=\frac {a}{e^q+1}

A_r=\frac {-a}{e^{-q}+1}

כאשר הכוח החיצוני המופעל גדול מספיק, מסה 0 הנתונה בשיווי משקל במיקום -a תתגבר על בור הפוטנציאל ותנוע מעבר למיקום -\frac{a}{2} כך שהנקע, כלומר בור הפוטנציאל הריק, ינוע בכיוון זה גם הוא.

כוח קריטי זה יחושב על ידי הפתרון עבור n\leq0

x_0=-\frac{a}{2}=\frac{f_c}{\Kappa}-a+A_l 

f_c=\frac{a\Kappa}{2}tanh(\frac{q}{2}) 

הכוח המקסימלי הדרוש להזזת המסות ממקומן אל בור הפוטנציאל הסמוך הוא \frac{a\Kappa}{2}  בין שיש נקעים ובין שאין. קיום הנקעים מוביל לפונקציית tanh(\frac{q}{2}) , וכך מקטין את הכוח הדרוש לעיוות הגביש.

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

## אנרגיה של נקעים ותנאים היווצרם

### האנרגיה הפוטנציאלית הדרושה ליצירת נקעים
כדי למצוא את האנרגיה הדרושה ליצירת נקע נתייחס לגביש תלת־ממדי ובו נקע בורג שנקבע כראשית הצירים. וקטור קו הנקע מקביל לציר \widehat{z} ולכן סטיית המיקום גם היא בכיוון זה \bigl(\vec u(\vec{r})=u_z(x,y)\bigr). נתייחס אל תורת האלסטיות הליניארית עבורה סטיית המיקום רציפה בכל המרחב מלבד בסביבת הנקע.

אם כך האנרגיה הפוטנציאלית היא:

U=\frac{a\mu}{2}\int d^2r(\nabla u)^2

a הוא קבוע בעל יחידת אורך שנבחר להיות אורך וקטור בורגרס של הנקע (שאורכו כקבוע הסריג) ו-μ קבוע לאמה. בשיווי משקל סטיית המיקום מקיימת את משוואת לפלס \nabla ^2 u=0, ולכן הפתרון שייבחר הוא פונקציה אנליטית מרוכבת, בעלת אי־רציפות בעת הקפה סביב מרכז הנקע. אי־רציפות זו מייצגת את וקטור בורגרס.

u_z(x,y)=\frac {a}{2\pi}Im\{ln(x+iy)\}

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

U=\frac{a^3\mu}{8\pi ^2}\int\limits_{a}^{R} dr2\pi r\frac{1}{r^2}=\frac{a^3\mu}{4\pi}ln\bigl(\frac{R}{a}\bigr)+w

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

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

u_z(x,y)=\frac {a}{2\pi}Im\{ln(x+iy)-ln(x-x_0+iy)\}

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

U=\frac{a^3\mu}{2\pi}ln\bigl(\frac{x_0}{a}\bigr)+2w

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

### טמפרטורה קריטית ליצירת נקעים
ניתן לראות כי האנרגיה של נקע יחיד מתבדרת עבור גביש גדול מאוד \bigl(R\gg1\bigr), אך עם זאת קיומם של נקעים מתאפשר מתוך שיקולי אנטרופיה והאנרגיה החופשית של המערכת.

בהנחה כי בריבוע באורך צלע R וקבוע סריג a יש \frac {R^2}{a^2} אתרים למיקום מרכז הנקע, האנטרופיה היא

S=2k_Bln\bigl(\frac {R}{a}\bigr)

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

F=U-TS=\bigl(\frac{a^3\mu}{4\pi}-2k_BT\bigr)ln\bigl(\frac{R}{a}\bigr)

בטמפרטורה הקריטית האנרגיה החופשית מתאפסת. מכאן שהטמפרטורה הקריטית היא

T_c=\frac{a^3\mu}{8\pi k_B}

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

## קישורים חיצוניים

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

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