פרק 9 — משוואת Butler–Volmer ותיאוריית טאפל

9.1 מבוא

בפרק הקודם הכרנו את מושג הקיטוב — הסטייה מפוטנציאל שיווי-משקל — ואת שלושת המנגנונים שיכולים להגביל תגובה אלקטרוכימית. עכשיו נרצה לכמת: אם הפוטנציאל הוסט ב-$\eta$ וולט, מה יהיה הזרם? האם הקשר ליניארי? אקספוננציאלי? ומתי כל אחד מהם תקף?

התשובה היא משוואת Butler–Volmer — ואפשר לטעון שהיא המשוואה המרכזית של כל האלקטרוכימיה הקינטית. היא מקשרת בין קיטוב לבין זרם, מכילה את כל ה”אופי” של תגובה אלקטרוכימית, ומצידה נגזרות שיטות מדידת הקורוזיה הנפוצות ביותר.

---

9.2 משוואת טאפל — ההיגיון האנרגטי

נבחן תגובה אלקטרוכימית כללית:

$$\mathrm{Ox}+n{e}^{-}\leftrightarrow \mathrm{Red}$$

בפוטנציאל שיווי-משקל, שני הכיוונים מתרחשים בקצב שווה. כאשר מסיטים את הפוטנציאל ב-$\eta$ וולט, מה קורה לאנרגיית ה-”קפיצה” שהאלקטרון חייב לעבור?

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

מצב המעבר נוצר במרחק מסוים מפני האלקטרודה. אם עובי שכבת הכפל הוא $r$, ניתן לסמן את קואורדינטת מצב המעבר כ-$\alpha r$, כאשר $\alpha$ הוא גורם הסימטריה (symmetry factor). עבור רוב התהליכים $\alpha =0.5$. הפרמטר $\alpha$ הוא מספר בין 0 ל-1 (ברירת המחדל 0.5). הוא מתאר כיצד מתחלק שינוי הפוטנציאל בין שני הכיוונים: חלק $\alpha$ מהשינוי מקל על מעבר האלקטרון בכיוון החמצון, וחלק $(1-\alpha )$ מקל עליו בכיוון החיזור. זהו מאפיין של אתר ה”מעבר” בין הצורה המחומצנת למחוזרת, ולמעשה מבטא את א-הסימטריה של מחסום האנרגיה.

כאשר מופעל פוטנציאל יתר $\eta$, מחסום האקטיבציה יורד ב:

$\Delta \Delta G^{‡}=-\alpha nF\eta (9.1)$

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

בשיווי משקל קבוע הקצב הוא:

$k_{eq}=k_0{e}^{\frac{E_a}{RT}}=k_0{e}^{\frac{\Delta G^{‡}}{RT}}(9.2)$

תחת פוטנציאל יתר $\eta$:

$k_{\eta }=k_0{e}^{\frac{\Delta G^{‡}-\Delta \Delta G^{‡}}{RT}}=k_0{e}^{\frac{\Delta G^{‡}+\alpha nF\eta }{RT}}=k_{eq}\cdot e^{\frac{-\alpha nF\eta }{RT}}(9.3)$

ובצורה פשוטה יותר:

$k=k_0{e}^{\frac{\alpha nF\eta }{RT}}(9.4)$

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

כעת נתרגם זאת לשפה האלקטרוכימית. מכיוון ש-$\Delta G=-nFE=-nF\eta$, וקבוע הקצב קשור לזרם לפי $k=\frac{i}{AnF}$ (כלומר $rate=kC$), נקבל:

$\frac{i}{AnF}=\frac{i_0}{AnF}{e}^{\frac{-\alpha nF\eta }{RT}}(9.5)$

או בצורה נוחה יותר:

$j=j_0{e}^{\frac{-\alpha nF\eta }{RT}}(9.6)$

זוהי משוואת טאפל המודרנית, כאשר $j_0$ הוא פרמטר הכולל את כל הקבועים לעיל. הוא מכונה צפיפות זרם החלפה (exchange current density) ומאפיין את מהירות הגעת האלקטרודה לשיווי משקל.

מבחינה היסטורית, משוואה $(9.6)$ התקבלה על ידי טאפל בשנת 1905 בצורתה האמפירית:

$\eta =a+b\log j(9.7)$

9.3 משוואת בטלר-וולמר

תגובת החמצון (האנודית) גדלה אקספוננציאלית עם קיטוב חיובי; תגובת החיזור (הקתודית) גדלה אקספוננציאלית עם קיטוב שלילי. הסכום הכולל נותן את משוואת Butler–Volmer:

$$i=i_0\left[e^{\frac{\alpha nF\eta }{RT}}-e^{-\frac{(1-\alpha )nF\eta }{RT}}\right](9.8)$$

האיבר הראשון הוא הזרם האנודי, השני — הקתודי. בשיווי-משקל ($\eta =0$) שניהם שווים ל-$i_0$, הזרם מתאזן, ו-$i_{net}=0$.

---

9.4 שני גבולות חשובים

משוואת Butler–Volmer מורכבת להתמודד איתה ישירות, אבל בשני גבולות קיצוניים היא מפשטת ל-צורות שימושיות מאוד.

קיטוב קטן — הגבול הליניארי

כאשר $|\eta |$ קטן מאוד (בפועל, פחות מכ-10–15 mV), ניתן לבצע פירוק משוואת בטלר-וולמר לסדרת טיילור מסדר ראשון. התוצאה:

$$i\approx i_0\cdot \frac{nF\eta }{RT}(9.9)$$

כלומר, הזרם פרופורציונלי לקיטוב — קשר ליניארי. זהו הבסיס למדידת LPR (Linear Polarization Resistance): מסיטים את הפוטנציאל בכמה מילי-וולט סביב ה-OCP ומודדים את השיפוע $\Delta i/\Delta \eta$. משיפוע זה ניתן לחלץ את $i_{corr}$.

קיטוב גדול — גבול טאפל

כאשר $|\eta |$ גדול (מעל כ-50–100 mV), אחד מהאיברים האקספוננציאליים שולט לחלוטין, והשני הופך זניח. בקיטוב אנודי גדול:

$$i\approx i_0\cdot e^{\frac{\alpha nF\eta }{RT}}(9.10)$$

נוטלים לוגריתם:

$$\log i=\log i_0+\frac{\alpha nF}{2.303RT}\cdot \eta (9.11)$$

מסדרים כ-$\eta$ בפונקציה של $\log i$:

$$\eta =a+b\cdot \log i(9.12)$$

זוהי משוואת טאפל — יחס ליניארי בין $\eta$ ל-$\log i$, שהיה מוכר ניסויית לפני שהיה לו הסבר תיאורטי. השיפוע $b$ נקרא שיפוע טאפל (Tafel Slope).

---

9.5 שיפוע טאפל ומה הוא מספר לנו

שיפוע טאפל אינו מספר שרירותי — הוא מכיל מידע ממשי על מנגנון התגובה. מנוסחת הגזירה ניתן לראות:

$$b=\frac{2.303RT}{\alpha nF}(9.13)$$

בטמפרטורת החדר, אם $\alpha =0.5$ ו-$n=1$ (תגובה חד-אלקטרונית), מתקבל $b\approx 120mV/decade$. כלומר, כל עלייה של סדר גודל בזרם (פי 10) מחייבת הוספת 120 מילי-וולט של קיטוב.

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

---

9.6 ייצוג לוגריתמי ומדוע הוא נחוץ

כדי לצייר עקומות טאפל, משתמשים תמיד בציר לוגריתמי לזרם. הסיבה מעשית לחלוטין: זרמים אלקטרוכימיים יכולים להשתנות על פני 8–10 סדרי גודל — מ-$10^{-9}A/c{\mathrm{m}}^2$ (ננו-אמפרים, קורוזיה זעירה של מתכות אצילות) ועד $10^{-1}A/c{\mathrm{m}}^2$ ויותר (תהליכי ציפוי, קורוזיה פעילה). על ציר ליניארי, כל הטווח מ-$10^{-9}$ עד $10^{-2}$ ייראה פרקטית אפס ולא ניתן לניתוח.

כמו כן, על ציר לוגריתמי אזורי הקיטוב הגדול (גבול טאפל) נראים כקווים ישרים — ומזיהויים המפורש ניתן לחלץ את $b$ ואת $i_0$.

---

9.7 צפיפות זרם ולא זרם מוחלט

במדידות אלקטרוכימיות משתמשים תמיד בצפיפות זרם $j=i/A$ (אמפר למ”ר או לסמ”ר), ולא בזרם המוחלט $i$. אם נמדוד שני דגמים שונים בגודל — אחד של 1 סמ”ר ואחד של 10 סמ”ר — ונחבר אותם לאותו פוטנציוסטט, הזרם המוחלט יהיה שונה פי עשר, אבל צפיפות הזרם — שמייצגת את מה שקורה על פני יחידת שטח — תהיה זהה. צפיפות הזרם היא תכונה אינטרינזית של הממשק המתכת-תמיסה, לא של גודל הדגם.

---

9.8 דיאגרמות Evans — פוטנציאל הקורוזיה והזרם מהצטלבות עקומות

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

נקודת החיתוך של שתי העקומות — כאשר $i_{anodic}=i_{cathodic}$ — נותנת ישירות:

• $E_{corr}$: פוטנציאל הקורוזיה (המיקום על ציר ה-$E$),
• $i_{corr}$: זרם הקורוזיה (המיקום על ציר ה-$\log i$).

הכוח של הדיאגרמה הוא שהיא מאפשרת לחשוב ויזואלית על האפקטים של שינויים במערכת. למשל: מה יקרה אם נגביר את ריכוז החמצן? העקומה הקתודית “תעלה” (יותר תגובת חיזור בכל פוטנציאל) — ונקודת החיתוך תעבור לפוטנציאל אנודי יותר ולזרם גבוה יותר. תוצאה: יותר חמצן → יותר קורוזיה. זהו דוגמה לכך שתוצאה אינטואיטיבית-לכאורה (“חמצן שורף ברזל”) קיבלה הסבר כמותי דרך הדיאגרמה.

להפך: מה יקרה אם נוצרת פסיבציה? העקומה האנודית “תשקע” — $i_{corr}$ ירד בסדרי גודל. זהו עקרון ההגנה בפסיבציה מודרנית.

---

9.9 זרם קורוזיה וחוק פאראדיי

ה-$i_{corr}$ אינו מספר אבסטרקטי — ניתן להמיר אותו ישירות לקצב אובדן מסה באמצעות חוק פאראדיי:

$$\frac{dm}{dt}=\frac{i_{corr}\cdot M}{nF}(9.14)$$

כאשר $M$ הוא המשקל המולקולרי של המתכת, $n$ מספר האלקטרונים בתגובה, ו-$F$ קבוע פאראדיי. מכאן ניתן לחשב ישירות את קצב הנסיגה של פני השטח (mm/year), גודל סטנדרטי ונפוץ בהנדסת קורוזיה. כלומר, מדידה חשמלית בניסוי של שעה יכולה לתת תחזית לאורך חיים של עשרות שנים — בתנאי שהמודל מתאים למציאות.

---

9.10 מגבלות המודל

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

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

---

9.11 סיכום

משוואת Butler–Volmer מקשרת בין קיטוב לזרם, ומכילה שני פרמטרים מרכזיים: $i_0$ (פעילות האלקטרודה בשיווי-משקל) ו-$\alpha$ (א-הסימטריה של מחסום האנרגיה). מגבולותיה נגזרים שני יחסים שימושיים: קשר ליניארי בקיטוב קטן (בסיס ל-LPR) ומשוואת טאפל בקיטוב גדול. דיאגרמות Evans מאפשרות להציג את המערכת גרפית ולהבין בצורה ישירה כיצד שינויים בקינטיקה האנודית או הקתודית משפיעים על $E_{corr}$ ו-$i_{corr}$.

בפרק הבא נשתמש בכל אלה ונעבור מהתיאוריה למדידה מעשית: Tafel Analysis, LPR, ושיטות נוספות למדידת קצב קורוזיה ממשי.