Matière active : quand la turbulence naît sans inertie

Un film rempli de filaments peut se plier, se fragmenter en domaines et produire des écoulements tourbillonnants alors même que l’inertie est négligeable. Appeler cela turbulence ne signifie pas que chaque mécanisme de l’atmosphère ou d’un tuyau rapide se retrouve à l’identique sous le microscope. La question intéressante est inverse : comment une matière qui consomme de l’énergie localement peut-elle fabriquer des mouvements collectifs à grande échelle, puis perdre leur régularité ? Le moteur n’est plus seulement une force appliquée aux frontières ; il est distribué dans les constituants.

Sanchez et ses collègues ont construit en 2012 une matière active à partir de microtubules et de moteurs moléculaires, alimentés par l’ATP. L’assemblage produit des mouvements collectifs et permet de relier une source chimique à des écoulements. Le résultat n’implique pas que toute suspension de filaments ait le même régime : les assemblages, les surfaces et la disponibilité du carburant comptent. ^[1]

Le nombre de Reynolds compare l’ordre de grandeur de l’inertie, ρU²/L, à celui de la viscosité, ηU/L². Son petit nombre ne dit pas que rien ne peut être instable. Il dit que l’accélération inertielle n’est pas le principal terme du bilan des forces. Prenons un exemple numérique entièrement pédagogique : ρ = 1000 kg/m³, U = 10 µm/s, L = 100 µm et η = 10⁻³ Pa·s donnent Re = 10⁻³. Les valeurs ne décrivent pas une préparation expérimentale particulière. Elles montrent seulement comment un écoulement vivant et désordonné peut rester dans le régime où les forces visqueuses dominent.

L’équation de Stokes détermine alors presque instantanément la vitesse correspondant aux contraintes présentes. Pourtant, ces contraintes changent parce que les orientations se déplacent, se tordent et réagissent à l’écoulement. Il existe donc une boucle dynamique même quand le terme d’inertie disparaît : orientation → contrainte → vitesse → nouvelle orientation. Confondre absence d’inertie et absence de mémoire collective ferait manquer le phénomène essentiel. Le temps peut se loger dans la relaxation de la structure plutôt que dans l’accélération d’une masse.

Pour décrire une population de bâtonnets alignés sans choisir une extrémité privilégiée, on utilise un directeur n identifié à −n. Une flèche ordinaire compterait à tort deux états là où la matière en présente un seul. En deux dimensions, le tenseur Q = S(nn − I/2) conserve cette symétrie. S décrit la force de l’alignement ; n décrit sa direction. Quand S devient faible, parler d’une direction précise perd du sens, ce qui est particulièrement important au cœur des défauts. Cette distinction empêche de transformer une carte de directions peu fiables en une belle carte de tourbillons artificiels.

Le cadre des nématiques actifs organise les liens entre orientation, écoulement, défauts et confinement. La revue de Doostmohammadi et ses collègues est utile pour comparer des systèmes dont les détails microscopiques diffèrent. Il faut néanmoins identifier les coefficients et les symétries avant de transférer une conclusion d’un système de moteurs et filaments à un tissu vivant. ^[2]

Faisons le tour d’un défaut en suivant l’orientation des bâtonnets. Si θ augmente de π, le directeur revient à une orientation physiquement équivalente : la charge vaut +1/2. S’il diminue de π, elle vaut −1/2. Les silhouettes dites comète et trèfle viennent de cette géométrie, pas d’un dessin choisi librement. Notre figure construit θ = q arg(x + iy), puis trace des segments sans pointe. Le centre est laissé vide, car le modèle d’orientation seul n’y spécifie pas l’amplitude S. Une couleur de rotation de fluide et une charge de directeur ne mesurent d’ailleurs pas la même chose.

Une paire +1/2 et −1/2 a une charge totale nulle. Sa création ne contredit donc pas la conservation globale imposée par un domaine et ses conditions aux bords. Les deux défauts peuvent se séparer, réorganiser les orientations puis s’annihiler. Mais le nombre de défauts ne suffit pas à déterminer l’intensité du transport : deux films avec une densité comparable peuvent avoir des corrélations temporelles et des interactions avec les parois très différentes. Il faut mesurer simultanément la structure et la dynamique, plutôt que prendre un compteur de singularités pour une description complète.

Dans cette convention, ζ règle l’amplitude de la contrainte active, Γ représente une friction avec le support et fᵉˡ les forces associées à l’élasticité d’orientation. Le signe extensile ou contractile doit toujours être relié à la définition retenue de Q et σᵃ. Une activité constante dans un milieu parfaitement uniforme peut ne produire aucune force volumique : c’est la divergence de la contrainte, donc aussi les gradients d’orientation, qui intervient. Ainsi, une légère courbure peut déclencher un écoulement qui amplifie cette courbure. L’instabilité est collective, même si chaque moteur suit une règle locale.

L’élasticité résiste aux variations d’orientation sur une longueur ℓ, tandis que l’activité tend à les amplifier. Leur comparaison dimensionnelle fait apparaître ℓₐ. Cette relation d’échelle ne fixe ni un coefficient universel, ni la forme précise d’un spectre. Dans notre calcul réduit, multiplier |ζ| par quatre divise ℓₐ par deux ; le graphique ne contient aucune mesure de diamètre réel. Si K ou la concentration changent en même temps que l’ATP, attribuer toute variation de taille au seul ζ serait une erreur d’identification. Un beau collapse de courbes n’annule pas ce problème.

Une autre longueur apparaît quand la friction est présente. Comparer η∇²u et Γu donne ℓₛ = √(η/Γ). Pour des variations beaucoup plus larges que ℓₛ, la friction peut limiter les interactions hydrodynamiques de longue portée. Les deux longueurs ne sont pas interchangeables : ℓₐ compare activité et élasticité ; ℓₛ compare viscosité et friction. Modifier le support peut donc transformer les structures sans modifier directement la chimie des moteurs. En pratique, il faut ajouter la taille du domaine et les conditions de bord à ce bilan.

Le champ de comparaison du graphique provient de ψ = sin x sin y, avec u = (∂ᵧψ, −∂ₓψ). Calculer sa divergence donne immédiatement zéro, car les dérivées croisées se compensent. Sa vorticité n’est pourtant pas nulle : dans cette convention elle vaut 2 sin x sin y. Nous avons donc fabriqué des cellules tourbillonnantes sans résoudre aucune équation de matière active. Cette construction est un contre-exemple utile : la ressemblance d’une image à un film expérimental ne prouve pas que le mécanisme a été identifié. Il faut confronter des prédictions temporelles, des corrélations et des réponses à une perturbation.

Une vitesse estimée par suivi de particules ou corrélation d’images possède une résolution et un bruit. Dériver cette carte pour obtenir la vorticité amplifie les erreurs aux petites échelles. Le protocole devrait annoncer le pas spatial, l’intervalle temporel, le filtrage et la manière de traiter les bords. Un maximum spectral qui se déplace quand le lissage change fortement peut être un artefact. À l’inverse, une longueur robuste entre méthodes et sous perturbations constitue un argument beaucoup plus convaincant qu’une image unique particulièrement spectaculaire.

Le travail de Doostmohammadi et ses collègues sur l’apparition de turbulence mésoscopique en 2017 examine comment des états ordonnés peuvent perdre leur stabilité dans un modèle actif. Il fournit un point d’entrée pour distinguer motifs réguliers, intermittence et désordre développé, sans réduire la transition à un simple équivalent du Reynolds critique d’un tuyau. ^[3]

Hillebrand et Alert ont rapporté en 2025 une transition discontinue dans le modèle étudié, avec coexistence et hystérésis. Leurs résultats rendent le rôle des interactions hydrodynamiques et du confinement particulièrement important. Ce résultat ne doit pas être transformé en loi selon laquelle toutes les matières actives auraient le même seuil ou le même ordre de transition. ^[4]

Cette équation est notre modèle jouet d’hystérésis ; elle ne reproduit pas les paramètres du travail cité. L’état A = 0 est stable pour μ < 0. Une branche d’amplitude non nulle peut déjà exister pour −1/4 < μ < 0. Sur cet intervalle, deux états stables coexistent, séparés par une branche instable. En augmentant μ depuis l’état calme puis en le diminuant depuis l’état agité, on peut suivre deux chemins différents. Le calcul donne un sens concret à « l’histoire du système » : ce sont les bassins d’attraction et les perturbations, pas une mémoire mystérieuse.

Saw et ses collègues ont relié des défauts nématiques à des contraintes et à l’extrusion de cellules dans des épithéliums étudiés. C’est un lien mécanobiologique important, mais son transfert à un autre tissu exige de considérer adhésion, division, mort et régulation biochimique. La forme d’un défaut ne suffit pas à déterminer le destin d’une cellule. ^[5]

Une explication physique est utile lorsqu’elle précise à quel niveau elle s’applique. Le tenseur d’orientation peut résumer des milliers de détails et prédire certaines contraintes ; il ne remplace pas pour autant une description des voies de signalisation. Réciproquement, connaître une protéine ne donne pas automatiquement la dynamique collective d’un tissu. L’enjeu scientifique est de construire des passerelles mesurables entre ces niveaux : modifier une activité, observer la contrainte, puis vérifier la réponse cellulaire avec des contrôles qui évitent de confondre corrélation spatiale et causalité.

Un protocole convaincant ferait varier séparément le carburant, la friction du support et la taille du domaine. Il mesurerait orientation, vitesse, densité de défauts et corrélations avec les mêmes règles d’analyse. Il répéterait la montée et la descente du contrôle, en laissant un temps d’attente annoncé. Il comparerait plusieurs conditions initiales, puis des domaines de tailles différentes. Une boucle d’hystérésis qui disparaît lorsqu’on attend davantage peut signaler un transitoire lent plutôt qu’une coexistence stationnaire. Voilà pourquoi un diagramme de phase nécessite un protocole temporel aussi soigneux que le contrôle des paramètres.

Le gain conceptuel est considérable : un fluide peut devenir désordonné sans emprunter la voie inertielle habituelle. Les tourbillons visibles sont la conséquence d’un couplage entre consommation d’énergie, orientation et forces. La spirale n’est donc ni un moteur universel ni une explication autonome. Elle devient intéressante quand ses courbures, ses singularités et sa dynamique permettent de tester ce couplage. C’est à cette condition qu’un motif séduisant devient un outil pour comprendre le transport, l’auto-organisation et certaines mécaniques du vivant.

Multiplions le bilan de Stokes par u, puis intégrons sur un domaine incompressible périodique. La pression ne fournit pas de puissance nette ; l’intégration par parties transforme le terme visqueux en une dissipation positive lorsqu’on le place de l’autre côté du bilan. Dans un modèle réduit où le travail élastique est négligé, la puissance active injectée doit équilibrer dissipation visqueuse et friction. Pour un modèle complet avec orientation évolutive, il faut aussi suivre le travail élastique et sa relaxation : supprimer ce terme ferait passer une variation d’énergie stockée pour de la puissance active.

Pour un mode variant sur une longueur ℓ, le rapport entre dissipation visqueuse et friction est de l’ordre de η/(Γℓ²) = (ℓₛ/ℓ)². Si ℓ vaut dix fois ℓₛ, ce rapport vaut environ un centième dans cette estimation ; si ℓ vaut un dixième de ℓₛ, il vaut environ cent. Ce changement de dominance explique pourquoi le même matériau sur deux supports peut produire des corrélations différentes. Ce sont des rapports d’échelle, pas des mesures. Leur utilité est de désigner les termes à garder avant d’ajuster une simulation à un film.

Pour f(A) = μA + A³ − A⁵, une petite perturbation δA autour d’un état A* évolue au premier ordre comme dδA/dt = f′(A*)δA. Le signe de f′ détermine donc croissance ou décroissance. À A* = 0, f′ = μ. Pour un état non nul, poser y = A*² et utiliser μ = y² − y donne f′ = 2y(1 − 2y). La branche y > 1/2 est stable, la branche 0 < y < 1/2 instable. Cette vérification indépendante justifie les traits pleins et pointillés du graphique, sans prendre son apparence pour une preuve.