נוהל ניסויי מבוסס רובוט לפיתוח חישת גז קולורימטרי

Summary

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

Abstract

מאמר זה מציג תוכנית ניסויית מבוססת רובוט שמטרתה לפתח חיישן גז קולורימטרי יעיל ומהיר. התוכנית משתמשת בגישת Design-Build-Test-Learning (DBTL) אוטומטית, המייעלת את תהליך החיפוש באופן איטרטיבי תוך אופטימיזציה של מתכונים מרובים למרווחי ריכוז שונים של הגז. בכל איטרציה, האלגוריתם מייצר אצווה של הצעות מתכונים המבוססות על פונקציות רכישה שונות, ועם הגידול במספר האיטרציות, הערכים של פונקציית המטרה המשוקללת עבור כל מרווח ריכוז משתפרים משמעותית.

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

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

Introduction

היישומים המעשיים של חיישני גז הם נרחבים מאוד ושימשו בתחומים שונים כגון ניטור סביבתי, תעופה וחלל, וטיפול בגז פסולת ^1,2,3. עקרון העבודה של חיישני גז מסתמך בדרך כלל על מספר מנגנונים, כגון אלקטרוכימיה, כרומטוגרפיה של גז ואופטי. מבין מנגנוני זיהוי רבים, מנגנון המבוסס על שינוי צבע התפתח למנגנון חומצה-בסיס הבולט באופן ייחודי. בשל עלותו הנמוכה והיישום הפשוט שלו, הוא נמצא בשימוש נרחב בתכנון חיישני גז ניידים וחד פעמיים רבים, כגון חיישני CO₂ ^1,4,5. חיישן מסוג זה משתמש בשינוי הצבע של כימיקלים מסוימים כדי לזהות ריכוזי גזים. כאשר ריכוז הגז משתנה, הוא גורם לחומר החיישן לחוות תגובות כימיות כגון קומפלקס יוני או שינויי צבע אינדיקטורים, מה שמוביל לשינוי בצבע הצבע הרגיש לגז⁶. על ידי איתור וניתוח השינויים בצבע, ניתן למדוד את ריכוז הגז בעקיפין. בינתיים, למרות היתרונות של עלות נמוכה וניידות, לחיישן מסוג זה עדיין יש כמה חסרונות, כגון מחזור פיתוח ארוך ויעילות נמוכה ^7,8,9. יחד עם זאת, שיטות מסורתיות של תכנון חיישנים מתקשות לעמוד במאפייני חישה מרובים בו זמנית, כגון השגת זמן התגובה, ההפיכות ומגבלת הגילוי הנדרשים. תחת פרדיגמת המחקר והפיתוח המסורתית, קשיים אלו מעכבים מאוד את הייצור והיישום הנרחב של חיישני גז קולורימטריים.

בתגובה לאתגרים שהוזכרו לעיל במחקר ופיתוח לפי דרישה, טכנולוגיית החיישנים הקולורימטריים שפותחה באמצעות תהליך ניסיוני זה יכולה לטפל בכמה מהחסרונות של חישת גז מסורתית. על ידי שימוש בגישת תכנון-בנייה-בדיקה-למידה איטרטיבית (DBTL)^10,11, ניתן לשפר משמעותית את היעילות של פיתוח חיישנים, ובכך להפחית את זמן המחקר והפיתוח ולענות ביעילות על צרכי המחקר והפיתוח ^1,12. במערך פיתוח DBTL טיפוסי, פיתוח חומרים חדשים נלקח כלולאת משוב איטרטיבית. הלולאה מכילה ארבעה שלבים מרכזיים: 1. תכנון פרמטרי האופטימיזציה, היעדים ודגימת מרחב הפרמטרים לניסוי ניסוי; 2. בנה את הדוגמאות של הפרמטרים שנבחרו; 3. בדוק את ערך היעד עבור הדגימות שנבנו; 4. ניתוח למידת מכונה של משוב היעד כדי להנחות את בחירת פרמטרי האצווה הבאה. בתהליך איטרטיבי זה, פלטפורמת הניסוי בעלת התפוקה הגבוהה המאפשרת בנייה ובדיקה מהירה של דגימות, ואלגוריתמי למידת המכונה הם מרכיבי המפתח. פלטפורמת הבדיקה האוטומטית בעלת התפוקה הגבוהה יכולה לבדוק בו זמנית עד 384 יחידות חישה, ולאסוף כמות גדולה של נתוני תגובה באיכות גבוהה. על ידי שימוש באלגוריתמים של למידת מכונה 13,14,15,16,17, כגון אופטימיזציה בייסיאנית רב-אובייקטיבית, ניתן לבצע אופטימיזציה של מדדי חישה מרובים של יחידות החישה (למשל, רגישות, זמן תגובה והפיכות) בו-זמנית ואוטומטית, ובכך לשפר את הביצועים הכוללים של מאפייני חישה שונים. מתכוני יחידות החישה שנוצרו על ידי אלגוריתם האופטימיזציה יכולים להשיג גילוי ריכוז CO₂ כמותי ללא כיול פרטני, ומדד השגיאה הריבועית הממוצעת (RMSE) יכול גם לעמוד באינדיקטורים הנדרשים.

התוכנית שלנו היא הליך ניסיוני שפותח על בסיס חישת גז קולורימטרי (ראה איור 1 לתרשים הזרימה). עם התפתחות מעבדות מונעות עצמיות, גישת ה-DBTL האוטומטית הראתה סיכויים מצוינים בשל היעילות, המהירות ויכולת החזרה הגבוהה שלה ^5,12. תהליך הפיתוח הידני המסורתי כולל התאמה של משתנה אחד בכל פעםתהליך הפיתוח הידני המסורתי כולל התאמת משתנה אחד בכל פעם, ואחריו שינוי משתנה אחר כדי לייעל את פרמטר היעד ולהשיג את התוצאה הרצויה. החסרונות העיקריים של תהליך זה כוללים יעילות נמוכה בניסויים ידניים, רגישות לטעויות אנוש, קושי בניהול משתנים רב-ממדיים בתרחישים מורכבים רב-ממדיים ונטייה להיתקע באופטימיות מקומית. בהשוואה לתהליך הפיתוח הידני, שיטת DBTL שאומצה בתוכנית ניסיונית זו משתמשת ברובוטיקה בשילוב עם אלגוריתמים מתקדמים של למידה אקטיבית כגון אופטימיזציה בייסיאנית רב-תכליתית. אופטימיזציה בייסיאנית היא גישה הסתברותית לאופטימיזציה של פונקציות אובייקטיביות יקרות להערכה^15,18. הוא בונה מודל חלופי, לעתים קרובות תהליך גאוס, כדי להעריך את פונקציית המטרה ומשתמש בפונקציית רכישה כדי להחליט על הנקודה הבאה לדגימה. פונקציית הרכישה מאזנת בין חקירה (חיפוש אזורים שנדגמו פחות) וניצול (חידוד אזורים ידועים בעלי ביצועים גבוהים) כדי למצוא ביעילות את המקסימום או המינימום הגלובלי. שיטה זו שימושית במיוחד במרחבי חיפוש בעלי מימדים גבוהים ולא קמורים שבהם טכניקות אופטימיזציה מסורתיות נאבקות. לאחר הגדרה גסה של הגבולות העליונים והתחתונים של תוכן הרכיבים, הוא מייעל באופן דינמי ניסויים כדי להשיג את היחס האופטימלי עם הביצועים הטובים ביותר באופן איטרטיבי. שיטה זו משפרת מאוד את היעילות ומפחיתה עלויות ומתפקדת בצורה יעילה יותר במרחב הרב-משתנים לפיתוח המתכון האופטימלי ^5,12.

המטרה הכוללת של מאמר זה היא לבסס הליך ניסיוני המבוסס על שיטת DBTL האוטומטית באמצעות טכנולוגיות מחשב שונות כגון למידת מכונה, אופטימיזציה בייסיאנית מרובת אובייקטים ופלטפורמות בדיקה ניסיוניות, כולל פלטפורמת הטיפול האוטומטית בנוזלים ופלטפורמת בדיקת גז בתפוקה גבוהה. זה יאפשר תכנון ומחקר של חיישני גז קולורימטריים. פלטפורמת הרובוט לטיפול בנוזלים "Opentrons OT-2" המותאמת אישית משמשת להשלמת ניסויים בהתאם להגדרות התוכנית, ומבצעת באופן אוטומטי שלבים כגון סינתזת מתכונים, ערבוב וטבילה. פלטפורמת בדיקת הגז הביתית בעלת תפוקה גבוהה משמשת לבדיקת גז וקריאת חיישנים קולורימטריים בצורה בעלת תפוקה גבוהה, שליטה מדויקת בריכוזי גזי המטרה ורישום שינויי הצבע של יחידות החישה בזמן אמת. בהשוואה למערכות ניסיוניות אחרות שתוכננו על בסיס DBTL, למערכת זו עלות חומרה נמוכה יחסית. במקביל, התייחסנו באופן חלקי להיבטים של המשימה הכרוכים בטעויות אנוש באמצעות גישה חצי אוטומטית. מתן התועלת השולית המקסימלית תוך שמירה על היתרונות של תכנון DBTL.

Protocol

1. ניסוי מקדים (בדיקת היתכנות)

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

1. קבע את טווח הריכוז של גז המטרה וקבע תצורת בדיקת גז.
   הערה: ריכוז גז המטרה בתצורת בדיקת הגז עולה באופן ליניארי או אקספוננציאלי.
2. לפני ואחרי זרימת כל ריכוז של גז המטרה, יש לטהר את מערכת הבדיקה בחנקן. שמור על היחס בין זמן הזרימה של חנקן לגז המטרה על 1:1.
   הערה: אם ריכוז גז המטרה הוא <1 עמודים לדקה, זמן הזרימה הוא ~10 דקות; אם ריכוז גז המטרה הוא ≥1 ppm, זמן הזרימה הוא ~5 דקות.
3. הכינו את פתרונות המקור של משתנים עם ריכוזים מתאימים לתמיסה על סמך גורמים כמו רוויה וצמיגות.
   הערה: בדרך כלל מומלץ לשמור על הריכוז הגבוה ביותר האפשרי של תמיסות המקור.
4. הגדר את הנפח הכולל של התמיסה של ניסוח החיישן הקולורימטרי ל-400 מיקרוליטר; לאחר מכן, הגדר את טווח הנפח עבור כל פתרון מקור של המשתנים בהתאם לספרות.
   הערה: טווחי הנפח עבור תמיסות המקור של צבעים הם בדרך כלל בין 0 μL ל-200 μL, בעוד שטווחי הנפח עבור פתרונות מקור אחרים הם בדרך כלל בין 0 μL ל-100 μL. מרווח הדגימה של תמיסת המקור הוא ~25 μL.
5. צור אצווה של 96 ניסוחים באמצעות פונקציות דגימה אקראיות כדי לאמת את היתכנות גילוי גז המטרה באמצעות קולורימטריה כימית.
6. טען את קובץ ניסוח החיישן הקולורימטרי, פתרון המקור, קצות, צלחת 96 בארות וממברנת PTFE לתוך המטפל בנוזל, וצור ברצף מספרי פרטי זהות בלתי תלויים (איור 2).
7. הגדר את המטפל בנוזל לאופן סימולציה כדי לחקות את הפעולה של סינתזה של נוסחאות חיישנים קולורימטריים, כגון שאיבת נוזלים, חלוקה, ניעור וטפטוף (איור 3).
8. אם אין שגיאות במצב הסימולציה, הגדר את המטפל בנוזלים לניסוי והתחל להפוך את הסינתזה של חיישנים קולורימטריים לאוטומטית.
   הערה: הקוד לאוטומציה של התהליך פותח על בסיס חבילת תוכנת הקוד הפתוח "Opentrons". הסינתזה של חיישנים קולורימטריים צפויה להימשך 3-6 שעות.
9. מכניסים את החיישנים הקולורימטריים לתנור בחום של 40 מעלות ומחממים למשך 50 דקות.
10. הניחו את החיישנים הקולורימטריים המיובשים בתא הגזים ובדקו את אחידות התאורה בסביבת הבדיקה ואת אטימות תא הבדיקה. לאחר אישור שהוא ללא בעיות, השתמש בבקרי זרימת המסה (MFC) כדי לשלוט באופן אוטומטי על קצבי הזרימה של גז האנליט (עם ריכוז c_a) והחנקן (ריכוז c_n), תוך ביצוע תצורת בדיקת הגז. נניח שקצב זרימת הגז הכולל הוא S בנפח/דקה, וריכוז האנליט היעד הוא c. קצבי הזרימה של MFC גז האנליט ו-MFC חנקן בנפח/דקה הם
    (1)
    (2)
11. במהלך תהליך הזרימה, מקם מצלמה מעל תא הגזים כדי לצלם תמונה כל 5 שניות כדי לתעד את שינויי הצבע של החיישנים הקולורימטריים (איור 4).
    הערה: בדיקת הזרימה צפויה להסתיים בעוד כשעתיים.
12. שימו לב שהמחשב טוען את התמונות שצולמו בסדר כרונולוגי, מחלץ את ערכי הצבע RGB מכל חיישן קולורימטרי בכל תמונה, ולאחר מכן מחשב את ההבדלים ביחס לצבע הבסיס שנמדד לפני החשיפה לגז המטרה, ובכך מתווה נתון וריאציה צבעונית של החיישנים הקולורימטריים לאורך זמן הזרימה. הנוסחה לחישוב הפרש הצבעים ΔE מוצגת כ-Eq (1):
    (3)
13. שימו לב אם ישנם חיישנים קולורימטריים המציגים שינויי צבע משמעותיים בגז המטרה והאם ערכי שינוי הצבע עולים בשילוב עם עלייה בריכוז גז המטרה. אם כן, היתכנות השימוש בחיישנים קולורימטריים כימיים לגילוי גז המטרה מאומתת.
14. בהתבסס על התוצאות שלפני הניסוי, התאם את הגבולות העליונים והתחתונים של פתרונות המקור של המשתנים והסר פתרונות מקור בעלי השפעה לא מהותית.

2. השתמש בפלטפורמה ניסויית רובוטית לביצוע תהליך האופטימיזציה האיטרטיבי של Design-Build-Test-Learn (DBTL)

1. עיצוב: הגדר פונקציות רב-אובייקטיביות וצור ניסוחים של חיישנים קולורימטריים.
   1. הגדר פונקציה רב-אובייקטיבית לחישוב ניקוד ההערכה המשוקלל עבור מספרי כשרון מרובים של החיישנים הקולורימטריים (איור 5).
      הערה: ניתן למצוא את הפירוט של ציון ההערכה המשוקלל בהפניה¹.
   2. אם אופטימיזציה של DBTL נמצאת בסיבוב הראשוני, צור באופן אקראי 96 ניסוחים ראשוניים של חיישנים קולורימטריים וצור מזהה קמפיין עבור משימת אופטימיזציה זו.
   3. אם ה-DBTL נמצא בסבב השני או מאוחר יותר, תכנן את האצווה הבאה של 96 ניסוחי חיישנים קולורימטריים תוך שימוש בפונקציות רכישה שונות (כגון גבול ביטחון עליון, הסתברות לשיפור ושיפור צפוי). בנוסף, כוונן את ההיפרפרמטרים של פונקציית הרכישה בכל סיבוב. בדרך כלל, κ כהיפרפרמטר עבור UCB, עם ערך שאינו עולה על 5 ומשמש כהיפרפרמטר עבור EI ו-POI, שיורד בהדרגה לקרוב ל-0 ככל שאיטרציות האופטימיזציה גדלות. שלוש המשוואות הבאות (4), (5) ו-(6) בהתאמה ממחישות את פונקציות הרכישה עבור UCB, EI ו-POI:
      (4)
      (5)
      (6)
      כאשר μ(x) הוא הממוצע; σ(x) הוא השונות; κ, הם פרמטרים המאזנים בין חקירה וניצול, הנשלטים על ידי סיבוב האיטרציה וקצב הדעיכה; f(x) הוא ציון ההערכה המשוקלל; f_currentMax הוא ציון ההערכה המשוקלל הגבוה ביותר הנוכחי.
   4. לאחר יצירת ניסוחים של חיישנים קולורימטריים בכל סיבוב שאינו ראשוני, שים לב אם ניסוחים אלה של חיישנים קולורימטריים מציגים ציפיות גבוהות או אי ודאות גבוהה. אם לא, חזור על הפעולה של שלב 2.1.3. אם כן, בצע את שלב הבנייה.
2. Build: סינתזה אוטומטית של חיישנים קולורימטריים. ראה פרטים בשלבים 1.7-1.10.
3. בדיקה: תן לבדיקה האוטומטית של חיישנים קולורימטריים להמשיך. ראה פרטים בשלבים 1.11-1.13.
4. למד: התאם מודל פונדקאית
   1. השתמש במשתני ניסוח החיישן הקולורימטרי ובציוני ההערכה המשוקללים שלהם כקלט ופלט של המודל החלופי, בהתאמה.
   2. השתמש ברגרסיה של תהליך גאוס עם הפונקציה הממוצעת עבור k(x_i, x_j) ופונקציית השונות (גרעין) k(x_i, x_j) כדי להתאים למודל החלופי.
      הערה: איור 6 מציג את המודל החלופי למשתנים חד-ממדיים ודו-ממדיים. באלגוריתם התאמת התהליך הגאוסי שפותח, נבחר גרעין מטרן אנזוטרופי בשילוב עם גרעין רעש לבן (ראה שתי משוואות הבאות, Eq (7) ו-Eq (8)) כדי להבטיח את הכלליות של פונקציית הקרנל.
      (7)
      (8)
5. קריטריוני סיום אופטימיזציה: שים לב שאופטימיזציה של DBTL נעצרת כאשר מספר האיטרציות מגיע לערך מוגדר מראש, או כאשר אין שיפור משמעותי בציונים המשוקללים של החיישנים הקולורימטריים שנוצרו.

3. בנייה ואפיון של מערך החיישנים הקולורימטריים האופטימלי

1. עבור ששת מרווחי הריכוז של גז המטרה, בצע קמפיינים של אופטימיזציה של n DBTL כדי לגלות שישה ניסוחים של חיישנים קולורימטריים אופטימליים או כמעט-אופטימליים. בתוך כל מרווח ריכוז, ניסוח החיישן הקולורימטרי האופטימלי מציג את ציון ההערכה המשוקלל הגבוה ביותר בבדיקת גז היעד.
2. השתמש במטפל בנוזלים כדי לבנות בערך את מערכי החיישנים הקולורימטריים המורכבים משש ניסוחי החיישנים הקולורימטריים המותאמים. שלבים 1.7-1.10 מפרטים את הפעולה הספציפית באמצעות המטפל בנוזלים.
3. בדיקת חיי מדף:
   1. בנה 14 מערכי חיישנים קולורימטריים וחלק אותם לשתי קבוצות. אחסן קבוצה אחת במצב פתוח בטמפרטורה של 25 מעלות צלזיוס והשנייה בוואקום.
   2. שמור על תנאי בדיקה עקביים ובצע בדיקות תגובה יומיות לאורך תקופה של 7 ימים כדי להעריך את ההשפעה של שני תנאי האחסון על ביצועי מערך החיישנים הקולורימטרי, ובכך להסיק את חיי המדף בשתי אסטרטגיות האחסון.

4. כיול מערך החיישנים הקולורימטריים

1. לבניית מערך חיישנים קולורימטרי, ראה פרטים בשלבים 3.1-3.2.
2. דגימת נתונים:
   1. בחר 5-10 ערכי ריכוז במרווחים שווים עבור כל טווח ריכוזים של גז המטרה, עם סך של לא פחות מ-20 ריכוזים שונים.
   2. חשוף את מערך החיישנים הקולורימטריים לגז המטרה לפי ערך הריכוז בסדר עולה או יורד ורשום את ערכי התגובה של מערך החיישנים הקולורימטרי.
      הערה: נדרשו חמישה עד 10 מחזורי CO₂/N₂ עבור כל ריכוז CO₂ .
3. הקלטת נתוני חישה: לאחר חשיפה לגז המטרה בריכוזים שונים למשך 5 או 10 דקות, חלץ ערכי ערוץ R, G ו-B של n חיישנים קולורימטריים ממערך החיישנים הקולורימטריים כתכונות להזנה למודל הכיול. וקטור התכונה מכיל 3 × n ממדים.
4. חלוקת מערך נתונים: בהתבסס על מספר הריכוזים של גז המטרה, חלקו את נתוני התגובה בריכוזים שונים לערכת האימון, ערכת האימות וערכת הבדיקה ביחס של 7:1.5:1.5.
5. אימון מודלים: ראשית, השתמש בחבילות Python, כגון scikit-learn ו-torch, כדי לבנות ארבעה מודלים של רגרסיה של למידת מכונה לכיול: Ridge, Random Forest, Xgboost ו-Deep Neural Network. לאחר מכן, הגדר את שורש השגיאה הריבועית הממוצעת כפונקציית ההפסד. לבסוף, טען את מערך הנתונים והתחל לאמן את מודל הכיול.
6. בחירת מודל: אמת את הביצועים של ארבעת מודלי הרגרסיה של למידת מכונה באמצעות ערכת האימות ובחר את המודל עם השגיאה הריבועית הממוצעת הקטנה ביותר כמודל הכיול הסופי.
7. בדיקת דגם: השתמש בערכת הבדיקה כדי לבצע בדיקת ביצועים במודל הכיול הסופי, להעריך את הביצועים של מערך החיישנים הקולורימטריים ומודל הכיול לניתוח כמותי של גז המטרה.

תוצאות

דוגמה טיפוסית למערך ניסיוני זה היא "מערך חיישני CO 2 קולורימטרי רחב טווח בעל רגישות גבוהה"¹². ראשית, הניסוי מייצר תרשים המבוסס על השינוי ב-ΔE לאורך זמן בריכוז_CO2 קבוע לאחר אופטימיזציה רב-אובייקטיבית בייסיאנית (איור 7A). בהתבסס על זמן הת?...

Discussion

מאמר זה מציע תכנון ניסיוני שיכול לפתח חיישני גז קולורימטריים במהירות ובדייקנות רבה יותר. ניתן להשתמש בתהליך ניסיוני זה לפיתוח חיישנים קולורימטריים לגזים שונים, כגון לחות, CO₂ ואמוניה ^1,4,5. באמצעות השיטה של פלטפורמה ...

Materials

Name	Company	Catalog Number	Comments
96-Well Deep Well Plate	NEST	NEST 2 mL 96-Well Deep Well Plate, V Bottom
96-Well PCR Plate	NEST	NEST 0.1 mL 96-Well PCR Plate
cresol red	sigma aldrich	1.05225	Dyes for colorimetric reagents
Ethyl cellulose	sigma aldrich	200689	Dyes for colorimetric reagents
Ethyl cellulose	Aladdin	E110670-100g	Additive
Industrial Camera	HKVision	MV-CS060-10UM/C-PRO	used for recording color changes
Liquid handler	Opentrons	OT2	liquid handler
Mass Flow Controller	ASERT	AST10-DLCMX-500C-042-A2B2-48VY	used in controlling analytes gas mixtures
m-cresol purple	sigma aldrich	1.05228	Dyes for colorimetric reagents
Opentrons OT-2 Tips	Opentrons	OT-2 Tips, 300µL
Opentrons OT-2 Tips	Opentrons	OT-2 Tips, 20µL
phenol red	sigma aldrich	1.07241	Dyes for colorimetric reagents
polyethylene glycol	sigma aldrich	P1458	Dyes for colorimetric reagents
PTFE film	Interstate Specialty Products	PM15M	PTFE mambrane
Tetrabutylammonium hydroxide	sigma aldrich	86854	Base for colorimetric reagents
thymol blue	sigma aldrich	1.08176	Dyes for colorimetric reagents