Ogólny model tarcia kontaktowego w układach koloidalnych

Dlaczego drobne ziarnka i ich tarcie mają znaczenie

Wiele codziennych materiałów, od tuszu do drukarek po pastę do zębów i płynne podkłady, składa się z drobnych cząstek zawieszonych w płynie. Cząstki te nieustannie zderzają się i ocierają o siebie. Do niedawna modele komputerowe takich układów w dużej mierze pomijały szczegóły tarcia zachodzącego przy zetknięciu cząstek. Badanie to pokazuje, że to pozornie drobne pominięcie może prowadzić do błędnych przewidywań dotyczących tego, jak takie materiały płyną, zgęstnieją, a nawet spontanicznie rozdzielają się na gęste grudki i rozrzedzone obszary — oraz proponuje ogólny, termodynamicznie spójny sposób naprawy tego problemu.

Jak cząstki ślizgają się, toczą i wirują

Kiedy dwie sferyczne cząstki zderzają się w płynie, nie tylko odbijają się: mogą się między sobą ślizgać i toczyć. Ten ślizg generuje tarcie styczne, które sprzęga ruch translacyjny cząstki z jej obrotem. Tradycyjne modele dla dużych ziaren, jak piasek, już uwzględniają takie tarcie. Jednak w świecie koloidów — cząstek o rozmiarach submikrometrowych do mikrometrowych — losowe potrząsanie przez otaczający płyn (szum termiczny) jest silne i nie można go zignorować. Autorzy postawili sobie za cel zbudowanie modelu, w którym kontakty tarciowe i ta wszechobecna losowa ruchliwość są traktowane razem w sposób respektujący podstawowe zasady termodynamiki.

Połączenie tarcia z termicznym trzęsieniem

Kluczową intuicją jest to, że zawsze gdy tarcie odpływa energię z układu, losowe siły z otoczenia termicznego muszą tę energię z powrotem dostarczać, aby cząstki osiągały właściwą temperaturę otaczającego płynu. Wykorzystując matematyczne ramy równań Fokkera–Plancka, autorzy wyprowadzają precyzyjną postać, jaką muszą mieć te losowe siły i momenty obrotowe, gdy cząstki doświadczają tarcia stycznego przy kontakcie. Istotne jest, że losowe kopniaki muszą być powiązane zarówno z translacją, jak i rotacją w tej samej uporządkowanej formie co samo tarcie. W zależności od interpretacji rachunku stochastycznego w czasie (schematy Itô, Stratonovicha lub Hänggi–Klimontovicha), szum przyjmuje nieco różne, ale w pełni określone formy i może być prosty albo bardziej złożony — „mnożący” się z zależności od prędkości poruszania się cząstek.

Co idzie nie tak, gdy szum tarciowy jest niekompletny

Mając ogólny model, badacze przeprowadzili obliczenia dużej skali, by przetestować jego konsekwencje. Najpierw zbadali pasywne płyny koloidalne z różnymi prawami tarcia i pokazali, że uwzględnienie jedynie deterministycznego tarcia przy jednoczesnym pominięciu odpowiadającej mu części losowej prowadzi do poważnych niespójności. Symulowane cząstki przestają podążać za znanym rozkładem prędkości Maxwella–Boltzmanna, a ich ruch translacyjny i rotacja zdają się mieć różne efektywne temperatury, obie różne od temperatury rozpuszczalnika. Gdy dodano poprawnie skonstruowane losowe siły i momenty, te artefakty znikają: rozkłady prędkości i obrotów zgadzają się z oczekiwaniami teoretycznymi, a temperatura kinetyczna pokrywa się z temperaturą kąpieli.

Tarcie przekształca przepływ i separację faz

Zespół następnie zbadał, jak kontakty tarciowe wpływają na bardziej złożone zachowania. W przepływie napędzanym ciśnieniem przez szczelinowy kanał cząstki toczą się i ślizgają wzdłuż chropowatych ścian. Tarcie sprzęga ścinanie (gradienty prędkości) z obrotem cząstek i wpływa na łatwość poślizgu płynu przy ścianach. Co ciekawe, ogólna lepkość ulega tylko łagodnym zmianom przy umiarkowanych gęstościach, ale długość poślizgu przy powierzchni jest silnie zmniejszona w miarę wzrostu tarcia, podczas gdy ruch toczenia nadal zapobiega całkowitym warunkom bez poślizgu. Przechodząc do materii aktywnej — cząstek samopędnych, które nieustannie konsumują energię — autorzy badali motility-induced phase separation, gdzie aktywne cząstki spontanicznie tworzą gęste klastry. Kontakty tarciowe powiększają zakres warunków, w których to klastrowanie zachodzi. Jednak jeśli związany z tym losowy szum zostanie zaniedbany, przewidywany diagram fazowy zmienia się jakościowo: symulacje mogą pokazywać rozdział faz tam, gdzie model termodynamicznie spójny tego nie przewiduje, lub odwrotnie. To podkreśla, jak wrażliwe są kolektywne zachowania poza równowagą na poprawne odwzorowanie tarcia i szumu.

Co to znaczy dla modelowania realnych miękkich materiałów

Mówiąc prosto, badanie dostarcza brakującego elementu dla wirtualnych laboratoriów, które mają przewidywać zachowanie gęstych zawiesin i układów cząstek aktywnych pod wpływem przepływu lub samopędzenia. Autorzy pokazują, że nie można po prostu włączyć sił tarciowych z modeli granularnych, zostawiając szum termiczny bez zmian; losowe „kopniaki” muszą być starannie dopasowane do tarcia, aby model spełniał podstawowe zasady bilansu energii. Ich ogólna recepta ma zastosowanie do szerokiej klasy praw tarcia i schematów symulacyjnych, i może być użyta w popularnych pakietach dynamiki molekularnej. Otwiera to drzwi do bardziej wiarygodnych symulacji zjawisk takich jak zagęszczanie przy ścinaniu, przepływ wzdłuż teksturowanych powierzchni oraz formowanie wzorów w aktywnych, wirujących lub chiraloidalnych koloidach, przybliżając teorię do złożonego zachowania obserwowanego w rzeczywistych miękkich materiałach.

Cytowanie: 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

Słowa kluczowe: tarcie koloidalne, szum termiczny, zagęszczanie przy ścinaniu, materia aktywna, separacja faz