פרק 11 — ספקטרוסקופיית אימפדנס אלקטרוכימית (EIS)

11.1 מבוא

השיטות שעסקנו בהן עד כה — Tafel Analysis, LPR, עקומות קיטוב — פועלות על ידי הפעלת שינוי משמעותי בפוטנציאל ומדידת תגובת הזרם. הן מספקות מידע מצוין על פרמטרים מרכזיים כגון $i_{corr}$ ו-$E_{corr}$, אבל יש להן מגבלה מובנית: הן “רואות” את הממשק כמכלול — כסכום של כל התהליכים ביחד, ללא הפרדה ביניהם.

ספקטרוסקופיית אימפדנס אלקטרוכימית (EIS) פועלת אחרת לחלוטין. במקום לדחוף את המערכת בכוח, מפעילים עליה תנודה קטנה מאוד של פוטנציאל AC — כמה מילי-וולט בלבד — ומודדים את תגובת הזרם בתדירויות שונות, לרוב על פני 5–6 סדרי גודל (מ-$10^{-3}$ הרץ ועד $10^5$ הרץ). העובדה שתהליכים שונים “מגיבים” בתדירויות שונות מאפשרת להפריד ביניהם ולקבל תמונה הרבה יותר עשירה של מה שקורה על פני האלקטרודה.

---

11.2 מהו אימפדנס?

בזרם ישר (DC), ההתנגדות של מעגל היא $R=U/i$ — ערך פשוט שאינו תלוי בזמן. בזרם חלופי (AC), המצב מורכב יותר: הקבלים והמשרנים מגיבים בצורה שונה לפי תדירות, ושינוי הפוטנציאל גורם לזרם שאינו בהכרח באותה פאזה. לכן משתמשים באימפדנס $Z$ — הכללה של ההתנגדות לאות AC:

$$Z(\omega )=\frac{\tilde{U}(\omega )}{\tilde{I}(\omega )}$$

האימפדנס הוא גודל מרוכב:

$$Z=Z^{\prime }+j{Z}^{\prime \prime }$$

כאשר $Z^{\prime }$ הוא החלק הממשי (קשור לאיבוד אנרגיה — כמו נגד) ו-$Z^{\prime \prime }$ הוא החלק המדומה (קשור לאגירת אנרגיה — כמו קבל). הזווית בין $Z$ לציר הממשי היא זווית הפאזה $\varphi$, שאומרת כמה “מפגר” הזרם אחרי הפוטנציאל. הערה: ב-EIS משתמשים ב-$j$ לסימון $\sqrt{-1}$ (ולא ב-$i$ כנהוג במתמטיקה) — זוהי מוסכמה היסטורית באלקטרוכימיה ובהנדסת חשמל.

---

11.3 מדוע מערכות אלקטרוכימיות מתנהגות כמו מעגלים

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

ההתנגדות האוהמית של האלקטרוליט ($R_s$) מתנהגת בדיוק כנגד: אין תלות בתדירות, $Z=R_s$.

השכבה הכפולה החשמלית (Double Layer), שנוצרת מהפרדת מטענים על פני המתכת, מתנהגת בקירוב כמו קבל — $Z_C=1/(j\omega C_{dl})$. ככל שהתדירות עולה, האימפדנס של הקבל יורד.

מעבר המטען דרך הגבול — האלקטרון “קופץ” מהמתכת ליון — מתנהג כנגד $R_{ct}$ שמחובר במקביל לקבל השכבה הכפולה.

דיפוזיה של מגיבים ותוצרים לפני השטח ומהם מתוארת על ידי אימפדנס Warburg — גודל שיש לו חלק ממשי וחלק מדומה שווים, ומתנהג כקו בזווית 45° בדיאגרמת Nyquist.

---

11.4 מעגל Randles — המודל הפשוט ביותר

מעגל Randles הוא נקודת המוצא לניתוח EIS של תא פשוט. הוא מורכב מ-$R_s$ בטור עם חיבור מקביל של $C_{dl}$ ו-$R_{ct}$, ולעיתים עם אימפדנס Warburg בטור עם $R_{ct}$:

$$R_s\text{ — }(C_{dl}\parallel R_{ct})\text{ — }[W]$$

בניתוח מעגל זה, בתדירות גבוהה מאוד הקבל “מקצר” את $R_{ct}$ ורואים רק את $R_s$. בתדירות בינונית שולט מעגל ה-RC ומתקבל הקשת (semicircle) שקוטרה שווה ל-$R_{ct}$. בתדירות נמוכה מאוד שולטת הדיפוזיה ומופיע ה”זנב” של וורבורג.

11.5 דיאגרמות Nyquist

דיאגרמת Nyquist מציגה את $-Z^{\prime \prime }$ (הרכיב המדומה, בסימן הפוך) כפונקציה של $Z^{\prime }$ (הרכיב הממשי), כאשר כל נקודה בגרף מייצגת תדירות אחת. התדירות עצמה אינה מצוינת על הגרף — היא “נעלמת” — ולכן לעיתים מסמנים כמה נקודות ייחוס.

עבור מעגל Randles פשוט (ללא דיפוזיה), דיאגרמת Nyquist מראה חצי עיגול: ה-X intercept הימני (בתדירות נמוכה) הוא $R_s+R_{ct}$, ה-X intercept השמאלי (בתדירות גבוהה) הוא $R_s$, וקוטר החצי עיגול הוא $R_{ct}$. הנקודה בראש החצי עיגול (שיא ה-$-Z^{\prime \prime }$) מתרחשת בתדירות $\omega =1/(R_{ct}{C}_{dl})$.

ה-$R_{ct}$ שמתקבל הוא פרמטר מרכזי: $R_{ct}$ גדול אומר שמעבר המטען קשה — כלומר, הקורוזיה איטית. מהקשר לחוק Butler–Volmer ניתן לקשר ישירות:

$$R_{ct}=\frac{RT}{nF\cdot i_0}$$

כלומר, $R_{ct}$ גדול = $i_0$ קטן = קינטיקה איטית.

הערה: נוסחה זו נגזרת מלינאריזציה של משוואת Butler–Volmer בסביבת פוטנציאל שיווי המשקל. במדידות אימפדנס, משתמשים באות AC בעלת אמפליטודה קטנה מאוד, כך שהסטייה מפוטנציאל שיווי המשקל (פוטנציאל היתר $\eta$) קטנה מאוד. בגבול $\eta \to 0$ ניתן לפתח את שתי האקספוננטות בטור טיילור ולשמור רק את האיבר הלינארי:

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

ומכאן:

$R_{ct}={\frac{d\eta }{di}|}_{\eta =0}=\frac{RT}{nF{i}_0}$

---

11.6 דיאגרמות Bode

דיאגרמת Bode מציגה שני גרפים זה לצד זה כפונקציה של $\log \omega$: את $|Z|$ (מודול האימפדנס) ואת $\varphi$ (זווית הפאזה). בניגוד ל-Nyquist, כאן התדירות מופיעה ישירות על ציר ה-X, ולכן קל יותר לקשר בין אזורים בגרף לתהליכים פיזיקליים בסקאלות זמן מוגדרות.

ב-Bode של Randles פשוט: בתדירות גבוהה מאוד, $|Z|\to R_s$ (קבוע, הפאזה ≈ 0). בתדירות ביניים, $|Z|$ עולה והפאזה מגיעה לשיא שלילי (ה-RC מניב העברת פאזה של עד −90°). בתדירות נמוכה מאוד, $|Z|\to R_s+R_{ct}$ (קבוע שוב, הפאזה ≈ 0). כאשר מופיע Warburg, בתדירויות הנמוכות $|Z|$ ממשיך לעלות עם שיפוע של −½ ב-log-log, ובמקביל הפאזה שואפת ל-−45°.

---

11.7 אימפדנס Warburg ותהליכי דיפוזיה

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

$$Z_W=\frac{\sigma }{\sqrt{\omega }}(1-j)$$

כאשר $\sigma$ הוא מקדם Warburg (תלוי בריכוז ומקדם דיפוזיה). שני הדברים האופייניים לוורבורג: קו ב-45° ב-Nyquist (החלק הממשי שווה למדומה), ושיפוע −½ ב-Bode של $|Z|$. הופעת וורבורג בספקטרום EIS היא “חתימה” ברורה לכך שמעבר מסה (ולא מעבר מטען) הוא הצוואר הבקבוק.

---

11.8 CPE — כאשר הקבל “לא מושלם”

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

במקרים אלה מחליפים את $C_{dl}$ ב-Constant Phase Element (CPE):

$$Z_{CPE}=\frac{1}{Q(j\omega {)}^n}$$

כאשר $n$ הוא מספר בין 0 ל-1 (קבל אידאלי: $n=1$; נגד: $n=0$). ה-$n$ מכמת את מידת הא-אידאליות — ולעיתים קרובות עצמו הוא מידע פיזיקלי חשוב על אחידות פני השטח.

---

11.9 בחירת מעגל תמורה ו-Fitting

הניתוח המעשי של ספקטרום EIS מורכב משלושה שלבים: ראשית, בחירה של מעגל תמורה (Equivalent Circuit) שמתאר פיזיקלית את המערכת. שנית, fitting — התאמת פרמטרי המעגל לנתונים הניסויים, בדרך כלל בשיטת ריבועים פחותים. שלישית, פרשנות הפרמטרים שהתקבלו — $R_s$, $R_{ct}$, $C_{dl}$, מקדמי וורבורג וכו’.

הסכנה הגדולה ביותר ב-EIS היא Overfitting: הוספת רכיבים רבים יותר למעגל השקול תמיד תשפר את ה-${\chi }^2$ של ה-fitting — אבל מעגל מורכב מדי ייתן פרמטרים חסרי משמעות פיזיקלית. כלל האצבע: מספר הרכיבים במעגל צריך להיות מוצדק על ידי הבנה פיזיקלית של המערכת, לא רק על ידי שיפור ה-fit. ספקטרום EIS אינו חידה מתמטית — הוא מדידה של מערכת פיזיקלית.

---

11.10 EIS לניטור ציפויים

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

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

---

11.11 EIS לאורך זמן — Time-Domain EIS

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

11.12 זיהוי דפוסים — ה”ECG” של קורוזיה

אחד הדברים שמייחדים מומחי EIS מנוסים הוא יכולת זיהוי דפוסים בספקטרום עוד לפני fitting מלא — בדיוק כמו קריאת ECG ברפואה: הרופא מזהה פתולוגיה בצורת הגל לפני שמחשב פרמטרים. בEIS יש “חתימות” אופייניות: לולאה אינדוקטיבית בתדירות נמוכה — לרוב סימן לפיטינג פעיל. שני קשתות ברורות — מתכת מצופה עם שני ממשקים. חצי עיגול “שקוע” — פיזור קבועי זמן. זנב Warburg — דיפוזיה שולטת. עם ניסיון, הספקטרום “מדבר” לפני שמגיעים ל-fitting.

---

11.13 יתרונות וחסרונות — השוואה

שיטה	יתרון מרכזי	חסרון מרכזי
Tafel Analysis	$i_{corr}$, שיפועי טאפל ישירות	פולשנית; משנה פני שטח
LPR	מהירה, לא הורסת; לניטור רציף	פחות מידע; דורשת $B$ מוקדם
EIS	עשירה מאוד במידע; מפרידה תהליכים	מדידה איטית; פרשנות מורכבת

שלוש השיטות אינן מתחרות — הן משלימות. LPR לניטור יומיומי, Tafel לאפיון מהיר, EIS לחקירה מעמיקה.

---

11.14 סיכום

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

הכלים המרכזיים: דיאגרמת Nyquist (מציגה אחת מול האחרת — ממשי מול מדומה), דיאגרמת Bode (מציגה כנגד תדירות — מאפשרת לזהות אזורים), ומעגל Randles כנקודת מוצא. פרמטרים כגון $R_s$, $R_{ct}$, $C_{dl}$ ומקדמי Warburg מתקבלים מ-fitting ומאפשרים קישור ישיר לתהליכים פיזיקליים — לרבות $i_{corr}$ דרך קשר Butler–Volmer.

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