מחקר עומק
שתי סיבות שונות לא לדעת מה יש בפנים, ורק אחת מהן מתוקנת במכשיר טוב יותר
ספקטרוסקופיית עכבה מודדת מערכת בעשרות תדרים ומחזירה עקומה, ומהעקומה מנסים להסיק מה מתרחש בפנים. זו בעיה הפוכה מהסוג שבו התשובה אינה יחידה: מעגלים חשמליים שונים לחלוטין מייצרים בדיוק אותה עקומה. אבל יש כאן שני סוגים של אי ודאות שקל לבלבל ביניהם, ורק אחד מהם חסין בפני מדידה טובה יותר. ההבחנה הזאת היא גם מה שקובע מתי אלגוריתם אבולוציוני עוזר ומתי הוא רק מייפה את הבעיה.
יש שתי סיבות שונות לגמרי לא לדעת מה קורה בתוך מערכת, והערבוב ביניהן הוא הטעות העמוקה בתחום הזה. הסיבה הראשונה היא שהמדידה אינה טובה מספיק: יש רעש, חלון התדרים צר מדי, שני תהליכים נופלים כמעט באותו מקום. הסיבה השנייה היא שהמידע הדרוש להכרעה אינו קיים בנתונים כלל, בשום איכות מדידה. הראשונה מתוקנת במכשיר טוב יותר, בממוצע על יותר סריקות ובחלון תדרים רחב יותר. השנייה אינה מתוקנת בשום מכשיר, מפני שאין שם מה לתקן.
ספקטרוסקופיית עכבה היא המקום שבו ההבחנה הזאת מתחדדת, והיא עומדת מאחורי בדיקת בריאות של סוללות ליתיום, ניטור קורוזיה, אפיון תאי דלק וחיישנים ביולוגיים. העיקרון פשוט. במקום להפעיל מתח קבוע ולבדוק כמה זרם עובר, מפעילים אות מתנדנד קטן וסורקים תדרים, בדרך כלל מאלפית הרץ ועד מאה קילוהרץ. בכל תדר מודדים שני דברים: פי כמה המערכת מקטינה את הזרם, ובכמה היא מאחרת או מקדימה אותו ביחס למתח. זו הקשבה למערכת בקצבים שונים, כי תהליכים איטיים כמו דיפוזיה מגיבים בתדרים נמוכים, ותהליכים מהירים כמו קיבול בממשק שולטים בתדרים גבוהים.
שני המספרים בכל תדר מתאחדים למספר מרוכב אחד שנקרא עכבה. דוגמה שכדאי להחזיק בראש: נגד של 100 אוהם בטור עם קבל של 10 מיקרופאראד, בתדר 100 הרץ. החלק הממשי הוא 100 אוהם, כלומר מה שהופך לחום. החלק המדומה הוא מינוס 159.2 אוהם, כלומר מה שנאגר ומוחזר. יחד זה גודל של 188.0 אוהם בזווית של מינוס 57.9 מעלות. מי שרוצה להבין מאיפה מגיע החלק המדומה ימצא במהדורה גם דף לימוד על מספרים מרוכבים.
עד כאן הכיוון קדימה, והוא קל: בהינתן המעגל והפרמטרים אפשר לחשב את העקומה בשורה אחת. הקושי הוא בכיוון ההפוך. נתונה העקומה, מהי המערכת. ז'אק אדמר ניסח לפני יותר ממאה שנה שלושה תנאים לבעיה מנוסחת היטב: שקיים פתרון, שהוא יחיד, ושהוא תלוי ברציפות בנתונים. פענוח ספקטרום עכבה מפר את השני ואת השלישי. הקיום אינו הבעיה. הבעיה היא שיש יותר מדי פתרונות, לא מעט מדי.
זה ממקם את הבעיה במשפחה מוכרת, וזו הפרספקטיבה שכדאי להחזיק: שחזור תמונת CT, חידוד תמונה מטושטשת, היפוך גאופיזי, וניתוח זמני הרפיה ב-NMR. כולן בעיות הפוכות לא יציבות, וכולן נפתרות באותה משפחת שיטות. הקרובה ביותר לעכבה היא דווקא האחרונה. חילוץ התפלגות זמני הרפיה מספקטרום עכבה הוא התאום בתחום התדר של ההתמרה ההפוכה שספקטרוסקופיסט NMR מבצע כדי להפריד דעיכות מעריכיות חופפות, והוא קשה מאותה סיבה: מעריכיות בעלות קבועי זמן דומים כמעט אינן ניתנות להבחנה זו מזו. ה-CT הוא דווקא המקרה הקל של המשפחה, וזו הסיבה שמכשיר בבית חולים מחזיר תמונה שאפשר להסתכל עליה בעוד ספקטרום עכבה אינו מחזיר מודל יחיד.
הסוג הראשון של אי הוודאות נקרא בספרות אי זיהוי מבני. הוא תכונה של המודל, לא של המדידה. שרשרת Voigt ומעגל סולם הם שתי צורות קנוניות שונות של אותה פונקציית עכבה, תוצאה סטנדרטית בתורת סינתזת רשתות, ולכן הם מייצרים אותה עקומה בדיוק. לא דומה, זהה. סטיבן פלטשר פרסם ב-1994 ב-Journal of The Electrochemical Society טבלאות שלמות של רשתות מנוונות כאלה. חברת BioLogic מציגה בחומרי ההדרכה שלה מקרה שבו ארבע טופולוגיות אינן ניתנות להבחנה על סמך הספקטרום בלבד. מול ניוון כזה, מכשיר טוב פי אלף אינו מקדם אותך בכלום.
הסוג השני נקרא אי זיהוי מעשי, והוא תכונה של הנתונים שבמקרה נאספו. שני תהליכים עם קבועי זמן קרובים נראים כקשת אחת רחבה; תהליך שקבוע הזמן שלו נופל מחוץ לחלון התדרים אינו נראה כלל; שני פרמטרים יכולים להיות מתואמים כל כך שהנתונים אינם קובעים אף אחד מהם בנפרד. כאן, ורק כאן, מדידה טובה יותר היא בדיוק הפתרון. הורדת רעש, ממוצע על יותר סריקות, הרחבת חלון התדרים ושקלול נכון של השגיאה מצמצמים ממש את קבוצת המודלים שהנתונים מתירים. ההגדרה המקובלת אומרת זאת במפורש: פרמטר הוא מזוהה מעשית אם רווח הסמך שלו סופי, וזה בדיוק מה שנתונים טובים יותר עושים לרווח הסמך.
ברנרד בוקאמפ מנסח את שתי האפשרויות במשפטים עוקבים: מעגלים שקולים שונים יכולים להיות בעלי תגובת עכבה זהה לחלוטין, ומעבר לכך ייתכן שבתוך שגיאת הניסוי שני מעגלים לא קשורים ייתנו גרפים כמעט דומים. זהה לחלוטין מול כמעט דומה. זו כל ההבחנה, והיא קובעת מה כדאי לעשות: מול המקרה השני משפרים את המדידה, ומול המקרה הראשון משנים את הניסוי או מוותרים על ההכרעה.
לכן צריך להפריד שלוש שכבות שהשפה נוטה לערבב. מדידה היא לא התאמה טובה, והתאמה טובה היא לא פירוש. ההחלקה בין השכבות היא המקום שבו מאמרים נשברים. חוקר מודד עקומה, מתאים לה מעגל, מקבל שגיאה קטנה, וכותב שהתהליך הדומיננטי הוא העברת מטען. השגיאה הקטנה אמיתית. המשפט האחרון הוא השערה.
יש כלי שנועד לשמור על השכבה הראשונה, וכדאי לדעת מה הוא לא עושה. יחסי קרמרס-קרוניג קושרים מתמטית את החלק הממשי והמדומה של תגובה של מערכת ליניארית, סיבתית, יציבה וסופית, והמבחן שפרסם בוקאמפ ב-1995 בודק אם הנתונים מקיימים אותם. זו בדיקה של הנתונים, לא של המודל. היא תתפוס אות עירור גדול מדי, ארטיפקט של מכשיר, ובעיקר מערכת ששינתה את עצמה תוך כדי הסריקה. היא לא תכריע בין שני מעגלים. סקירת Nature Reviews Methods Primers מ-2021 אומרת זאת חד: כל ספקטרום שעקבי עם קרמרס-קרוניג ניתן להתאמה על ידי מעגל שקול, והבעיה היא שמודלים כאלה אינם יחידים.
המענה המתמטי לחוסר היציבות הוא רגולריזציה, וכאן נמצא חוט מקשר שקל לפספס. בחילוץ התפלגות זמני הרפיה השיטה הסטנדרטית היא רגולריזציית טיכונוב: במקום לשאול רק איזו התפלגות מתאימה הכי טוב, מוסיפים קנס על חספוס ומבקשים את ההתפלגות החלקה ביותר שעדיין עקבית עם הנתונים. יש לזה מחיר שכדאי להכיר: עוצמת הרגולריזציה משנה את התוצאה, ובחירות שונות נותנות מספר שונה של פסגות. מי שקורא פסגות כתהליכים קורא גם את בחירת הפרמטר.
מכאן הגיוני להגיע ל-Genetic Programming. אלגוריתם גנטי רגיל מקבל מבנה קבוע מראש ומחפש את הערכים הטובים ביותר עבורו. אם החלטתם שהמעגל הוא נגד בטור עם צמד מקבילי, האלגוריתם מחפש שלושה מספרים. Genetic Programming עושה משהו אחר: הפרט באוכלוסייה אינו רשימת מספרים אלא עץ, ולעץ יש גודל ומבנה משתנים. מוטציה יכולה להחליף קבל ברכיב CPE, והכלאה יכולה להעביר ענף שלם ממעגל אחד לאחר. החיפוש רץ בו זמנית על הצורה ועל המספרים, וזה מתאים לבעיה מפני שמה שרוצים לגלות הוא אובייקט סימבולי בעל מבנה משתנה.
אבל חיפוש חופשי על מבנה מייצר בעיה משלו. אם הציון היחיד הוא טיב ההתאמה, האבולוציה תלמד מהר שכל רכיב נוסף מקטין קצת את השגיאה, והעצים יתנפחו. בעכבה זה מסוכן במיוחד, מפני שכל רכיב מזמין פירוש כימי, ומעגל מנופח מייצר סיפור מנופח. לכן ציון רציני מוסיף קנס על מורכבות, וזו אותה תנועה שעושה טיכונוב בהתפלגות זמני ההרפיה: לייבא הנחה מוקדמת שהתשובה צריכה להיות פשוטה. הניסוח המדויק הוא שהמטרה משתנה. לא המעגל שנותן את השגיאה הקטנה ביותר, אלא מעגל מזוהה, סיבתי, מוצדק סטטיסטית וסביר פיזיקלית. עבודת Gene Expression Programming מ-Ghent University, שפורסמה ב-2021 ב-IEEE Transactions on Instrumentation and Measurement, מסירה רכיבים מיותרים במפורש ומחפשת את המעגל הפשוט ביותר שמסביר את המדידות. היא נבדקה בעיקר על סימולציות, ובסופה מקרה בוחן אחד של מדידות ביולוגיות אמיתיות.
צריך לומר בבירור מה זה אומר על ה-AI כאן, כי הניסוח המקובל הפוך. האלגוריתם אינו מגלה מבנים בלי הטיה אנושית. הוא זקוק לדקדוק פיזיקלי שנקבע מראש: אילו רכיבים מותרים, אילו חיבורים חוקיים, אילו גבולות על קבועי זמן, איזו התנהגות חייבת להיות סיבתית. בלי הדקדוק הזה הוא מייצר ביטויים מתמטיים חסרי פשר שמתאימים לנתונים מצוין. מה שהוא משוחרר ממנו הוא ההנחה על הטופולוגיה, לא ההנחות הפיזיקליות. זה הרבה, וזה לא אותו דבר.
ואם המידע להכרעה אינו נמצא בספקטרום, מאיפה הוא כן מגיע. התשובה המקובלת בתחום מנוסחת היטב בסקירת Nature Reviews: מודלים של עכבה אינם יחידים, וספקטרוסקופיית עכבה אינה טכניקה עומדת בפני עצמה. מודל שנבנה כהשערה מכוון את הבחירה בניסויים משלימים שיכולים לחזק אותו או להפריך אותו. במילים אחרות, ההכרעה אינה מגיעה מעיבוד טוב יותר של אותה עקומה, אלא מכך שמודדים שוב תוך שינוי מכוון של תנאי הניסוי ובודקים אם כל רכיב זז לכיוון שהפירוש שלו מחייב.
הבדיקה הראשונה היא טמפרטורה. תהליך אמיתי של העברת מטען מופעל תרמית: קבוע הקצב שלו, ואיתו צפיפות זרם החליפין, מציית לחוק ארניוס עם אנרגיית שפעול, ולכן התנגדות העברת המטען צונחת כשמחממים את התא. קיבול גאומטרי או דיאלקטרי כמעט אינו זז. מדידה של אותו תא בכמה טמפרטורות מפרידה בין שתי הפרשנויות. דיוק שחשוב לשמור עליו: חוק ארניוס חל על צפיפות זרם החליפין, והתנגדות העברת המטען שווה ל-RT חלקי nF כפול אותה צפיפות, כלומר היא נושאת גורם טמפרטורה נוסף. את אנרגיית השפעול קוראים מצפיפות הזרם, לא ישירות מההתנגדות.
הבדיקה השנייה היא גודל וגאומטריה. קיבול השכבה הכפולה הוא תכונה של שטח פנים ואמור לגדול איתו, והתנגדות העברת המטען אמורה לקטון איתו. התנגדות האלקטרוליט מתנהגת אחרת ותלויה בגאומטריית התא: עבור אלקטרודת דיסק היא הולכת כהופכי של הרדיוס ולא של השטח, ואפשר אפילו לחזות אותה מתוך מוליכות ידועה, מה שנותן בדיקה בלתי תלויה על איכות המדידה. הבדיקה השלישית היא מתח: התנגדות העברת מטען תלויה מעריכית במתח היתר דרך קשר באטלר-וולמר ותזוז בסדרי גודל, בעוד קיבול גאומטרי או דיאלקטרי לא. כאן חייבים לדייק במילה: קיבול השכבה הכפולה עצמו כן משתנה עם המתח, ולכן המבחן הזה מפריד קינטיקה מגאומטריה, לא קיבול מהתנגדות. הבדיקה הרביעית היא ריכוז: דילול של אלקטרוליט הרקע מעלה בעיקר את התנגדות התמיסה, ודילול של החומר הפעיל מעלה את מקדם ורבורג וגם את התנגדות העברת המטען. לכל רכיב יש כפתור אחר שמזיז אותו, וזה בדיוק מה שהופך אותם לניתנים להפרדה.
יש בתחום גם דוגמה מצוינת לפירוש שכולם חוזרים עליו ושהנתונים אינם תומכים בו. רכיב הפאזה הקבועה, CPE, נמצא כמעט בכל מעגל מודרני, וההסבר שנלמד בדרך כלל הוא חספוס פני השטח. חומר ההדרכה של Gamry ממליץ להתייחס למעריך כאל קבוע אמפירי בלי בסיס פיזיקלי ממשי, וסקירת Nature Reviews מ-2021 כותבת שמקור ה-CPE אינו ודאי ושההסברים המקובלים כנראה אינם חלים. סקירת הוראה מ-2023 מוסיפה את הממצא הלא נוח: השוואה בין אלקטרודות בדרגות חספוס שונות הראתה שדווקא אלקטרודה מחוספסת יותר התנהגה קרוב יותר לאידאלי, כלומר הפוך מהניבוי של הסבר החספוס. רכיב שמופיע בכל התאמה ושפירושו שנוי במחלוקת הוא תזכורת טובה לכך שהמעגל אינו התא.
הקו הישראלי בסיפור הזה אינו קישוט. קבוצתו של יועד צור בטכניון מריצה מאז תחילת העשור הקודם שיטה בשם ISGP, שבה תכנות אבולוציוני מפענח ספקטרום עכבה. חשוב לדייק במה היא עושה: ISGP אינו מייצר מעגל שקול, ומחבריו מנגידים אותו במפורש לניתוח מעגלים. הוא מחפש את התפלגות זמני ההרפיה, כלומר אילו תהליכים פועלים באילו קצבים, בלי לקבוע מהם. קו העבודה עבר מקתודות של תאי דלק במוצק ב-2012, דרך סופר-קבלים ב-2014 וב-2016, ועד ניתוח קרמיקה מיוננת ב-2017.
מה שקרה השנה מעניין יותר. במאמר שפורסם במרץ 2026 ב-Journal of Solid State Electrochemistry בנו אלאא אגבאריה, שמעיה זברי ויועד צור מחולל של ספקטרומים סינתטיים, כלומר נתונים שבהם התשובה הנכונה ידועה מראש, והריצו עליהם את השיטה של עצמם. ISGP שחזר היטב התנגדות אפקטיבית, קיבול וזמני הרפיה, ושחזר פסגה רחבה גם בנוכחות רעש. פסגה צרה מאוד הוא לא שחזר. היא חזרה כפסגה רחבה ונמוכה יותר, אף שהשטח שמתחתיה נשמר. זהו כישלון שהחוקרים מדווחים על השיטה של עצמם, והוא שווה יותר מהצלחה, מפני שהוא מגדיר את גבול ההפרדה של הכלי במקום להשאיר אותו לניחוש. שימו לב שזו בדיוק הבחנה מהסוג המעשי: היא נקבעת על ידי רעש ועל ידי הייצוג, ולכן היא נושא לשיפור.
סביב הקו הזה מסתדרת חזית שכולה משנת 2024 ואילך, וכדאי לקרוא אותה לפי המעמד של כל פריט. AutoEIS משלב חיפוש אבולוציוני עם הסקה בייסיאנית ובדיקת סבירות סטטיסטית, ותוכנתו פורסמה במאי 2025 ב-JOSS. מחקר קורוזיה מ-2024 ב-npj Materials Degradation השתמש בהסקה בייסיאנית כדי לזהות מתי הספקטרום אינו מכיל מספיק מידע כדי להצדיק מבנה מעגל מסוים, על נתונים מסומלצים. בסוף 2025 הציעו פאן ועמיתיו מערכת דו שלבית שבה אלגוריתם גנטי בוחר את המבנה וריבועים פחותים לא ליניאריים מתאימים את הפרמטרים, עם שגיאה מנורמלת מרבית של 1.07 אחוזים על שתי סוללות של 25 אמפר שעה. ועל הקצה יושבים שניים שהעיתונות נוטה להציג כהישגים גמורים: AutoREC היא preprint שגרסתו האחרונה מ-30 באוגוסט 2026, ו-EIS-GPT התקבל לפרסום ב-Cell Reports Physical Science אך טרם פורסם, וה-preprint שלו מתואר בידי מחבריו כהוכחת היתכנות שנשענת בעיקר על נתונים מסומלצים.
יש גם אילוץ פיזי שאי אפשר לתכנת סביבו. מדידה בתדר נמוך מאוד אורכת זמן, פשוט מפני שמחזור אחד אורך זמן. מחזור בודד באלפית הרץ אורך אלף שניות בדיוק, קרוב לשבע עשרה דקות, וזה גם המינימום שהספרות נוקבת בו למדידה אחת בתדר הזה. הורדת רעש מחייבת לחזור על המדידה, וזה מכפיל את הזמן, ובזמן הזה המערכת הנמדדת ממשיכה לחיות: סוללה נרגעת, משטח ממשיך להיאכל בקורוזיה, חיישן ממשיך לקשור. מדדתם משהו שהשתנה תוך כדי המדידה, וזה בדיוק מה שמבחן קרמרס-קרוניג נועד לתפוס. הממצא הלא אינטואיטיבי במחקר הקורוזיה הוא שאפשר להשמיט חלק ניכר מהנקודות בתדר נמוך בלי לפגוע מהותית באמידת הפרמטרים. אבל הסף תלוי במבנה שמניחים, ומשתנה בסדר גודל בין מערכת למערכת: 0.36 הרץ במקרה אחד לעומת 0.018 הרץ באחר. זו תוצאה תלוית מודל ולא כלל כללי, אבל היא מרמזת שחלק מהזמן שמושקע שם נקנה בלי תמורה.
מה שמצטייר אינו הכרזה על פריצת דרך ואינו קריאה לספקנות. מידול אוטומטי של ספקטרום עכבה עובד היום, במובן שמערכת יכולה לקבל נתונים ולהחזיר מעגל סביר בלי שמומחה יבחר אותו מראש. אוטומציה של הבנה פיזיקלית אינה עובדת, והסיבה מדויקת יותר מאשר חוסר בשלות: חלק מהמידע הדרוש להכרעה אינו נמצא בספקטרום, ולכן שום אלגוריתם לא יחלץ אותו משם. הוא נמצא בגאומטריה של האלקטרודה, בכימיה, בתלות בטמפרטורה ובניסויים משלימים. מה שחסר ב-2026 אינו רכיב, מפני שאבות טיפוס משכנעים קיימים לכל אחד מהם. מה שחסר הוא מערכת שמחברת אותם ומראה, על נתונים ניסויים רחבים ומגוונים, שהיא מייצרת לא התאמה טובה אלא מסקנה פיזיקלית שאפשר לשחזר.