מודל כללי למגעי חיכוך במערכות קולואידיות

מדוע חלקיקים זעירים והשחיקה ביניהם חשובים

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

איך חלקיקים מחליקים, מתגלגלים ומסובבים

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

שידוך בין חיכוך לרעידות תרמיות

התובנה המרכזית היא שכאשר חיכוך מפחית אנרגיה מהמערכת, כוחות אקראיים מהסביבה התרמית חייבים להזין חזרה אנרגיה, כדי שהחלקיקים יתאזנו לטמפרטורת הנוזל. באמצעות המסגרת המתמטית של משוואות פוקקר–פּלַנֶק, המחברים גוזרים את הצורה המדויקת שחייבים לקבל הכוחות והמערטורים האקראיים כשהחלקיקים חווים חיכוך משיקי במגע. באופן מהותי, הדחיפות האקראיות חייבות להיות מקושרות גם לתרגום וגם לסיבוב באופן מסודר כפי שגם החיכוך עושה. בהתאם לאופן שבו מפרשים חשבונאות סטוכסטית בזמן (איטו, סטרטונוביץ' או סכמות היינגי–קלימונטוביץ'), הרעש מקבל צורות מעט שונות אך מוגדרות במלואן, ויכול להיות פשוט או רעש "מכפילי" מורכב שתלוי במהירות התנועה של החלקיקים.

מה קורה כשחסר רעש החיכוך

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

החיכוך משנה זרימה והפרדת פאזה

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

מסקנות למידול חומרים רכים אמיתיים

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

ציטוט: Hofmann, K., Dormann, KR., Liebchen, B. et al. A general model for frictional contacts in colloidal systems. Commun Phys 9, 139 (2026). https://doi.org/10.1038/s42005-026-02624-5

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