Eppur si muove: Shape of topological defects -- and consequent motion -- in active nematics

1. Problématique et Contexte

Les défauts topologiques sont des caractéristiques fondamentales des phases présentant un ordre de type cristal liquide. Dans les systèmes actifs (où l'énergie est consommée localement, comme dans les tissus biologiques ou les suspensions bactériennes), ces défauts possèdent des propriétés radicalement différentes de leurs homologues à l'équilibre. En particulier, les défauts de charge +1/2 dans les nématiques actifs sont connus pour s'auto-propulser.

Cependant, la littérature théorique précédente reposait souvent sur une hypothèse simplificatrice critique : celle selon laquelle la forme du défaut (la distribution spatiale du paramètre d'ordre) reste inchangée par l'activité et correspond à celle d'un cristal liquide passif à constante de Frank unique (où les déformations de splay et de bend ont le même coût élastique).

Les auteurs soulignent que cette hypothèse conduit à des prédictions incohérentes :

• Elle suggère que dans certains modèles actifs, les défauts +1/2 ne devraient pas s'auto-propulser (si un coefficient d'activité spécifique est nul).
• Elle prédit à tort que des systèmes passifs à deux constantes de Frank pourraient voir leurs défauts s'auto-propulser, ce qui viole la conservation de l'énergie à l'équilibre.
• Elle échoue à expliquer pourquoi les défauts s'auto-propulsent toujours dans les simulations numériques, quelle que soit la nature de l'activité.

L'objectif de cet article est de démontrer que l'activité modifie intrinsèquement la forme des défauts topologiques et que négliger cette modification conduit à des erreurs fondamentales, même pour de très faibles niveaux d'activité.

2. Méthodologie

Les auteurs adoptent une approche combinant théorie analytique et simulations numériques directes :

• Modèle Théorique : Ils utilisent un modèle générique de nématique actif « sec » (sans couplage hydrodynamique explicite avec un fluide, bien que le cas humide soit discuté en annexe). L'équation de mouvement pour le tenseur d'ordre $Q_{ij}$ inclut des termes d'élasticité (constantes de Frank $K_1, K_2$) et des termes non-linéaires actifs ($L_1, L_2$) :
  $$\partial_t Q_{ij} = G_{ij}$$
  où $G_{ij}$ contient les termes d'énergie libre et les termes actifs.
• Développement Analytique :
  • Ils effectuent un développement perturbatif pour déterminer la forme des défauts à la fois dans le cœur (près de $r=0$) et dans le champ lointain ($r \to \infty$).
  • Pour les défauts -1/2, ils obtiennent une expression analytique fermée de la forme du défaut.
  • Pour les défauts +1/2, ils travaillent dans un référentiel se déplaçant avec le défaut (vitesse constante $v_d$) et établissent un lien direct entre la modification du champ angulaire au cœur et la vitesse d'auto-propulsion.
  • Ils utilisent la méthode de Halperin-Mazenko pour calculer la vitesse, mais en utilisant la forme correcte du défaut (modifiée par l'activité) plutôt que la forme passive standard.
• Simulations Numériques :
  • Ils résolvent les équations du modèle (1) dans un domaine circulaire en utilisant une méthode des éléments finis (via MATLAB PDE Toolbox).
  • Ils imposent des conditions aux limites correspondant à la forme d'un défaut $\pm 1/2$ dans un système à constante de Frank unique pour induire un défaut dans le volume.
  • Ils suivent la position du défaut au cours du temps pour mesurer sa vitesse et reconstruisent la forme du paramètre d'ordre pour valider les prédictions analytiques.

3. Résultats Clés

A. Modification de la forme des défauts

L'étude révèle que l'activité modifie significativement la forme des défauts par rapport aux systèmes passifs :

• Défauts -1/2 : La forme dans le champ lointain est décrite par des fonctions $A(\phi)$ et $B(\phi)$ qui dépendent des coefficients actifs et élastiques. Contrairement aux systèmes passifs à constante unique où le champ angulaire est purement $\theta = \pm \phi/2$, l'activité introduit des déviations angulaires et des variations d'amplitude dépendant de l'angle azimutal.
• Défauts +1/2 : La modification de la forme au cœur est directement liée à la vitesse d'auto-propulsion. L'auteur montre que le terme de correction angulaire au cœur, noté $b(\phi)$, satisfait :
  $$b(\phi) = \frac{1}{4}(v_d + a_0 L_2) \sin \phi$$
  (pour une orientation spécifique). Cela prouve que la forme du défaut et sa vitesse sont intrinsèquement couplées.

B. Validation de l'auto-propulsion universelle

Les simulations confirment que les défauts +1/2 s'auto-propulsent dans tous les systèmes actifs étudiés, indépendamment de la présence ou de la valeur des coefficients actifs spécifiques (comme $L_2$).

• La vitesse mesurée $v_d$ correspond parfaitement aux prédictions analytiques obtenues en tenant compte de la forme modifiée (équation 5 du papier), et non à celles basées sur la forme passive (équation 4).
• L'équation (4), basée sur l'hypothèse de forme inchangée, prédit à tort $v_d = 0$ si $L_2=0$, ce qui est contredit par les simulations.

C. Résolution des incohérences théoriques

En utilisant la forme correcte du défaut dans le formalisme de Halperin-Mazenko :

• La vitesse calculée pour un système passif à deux constantes de Frank devient nulle, respectant ainsi la conservation de l'énergie.
• Les prédictions antérieures suggérant que certains défauts actifs ne devraient pas bouger sont invalidées.

D. Implications pour les interactions

Les auteurs soulignent que puisque la forme du défaut détermine les conditions aux limites qu'il impose aux autres défauts, les interactions défaut-défaut dans les nématiques actifs dépendent de l'activité via la modification de la forme. Cela ouvre la porte à des interactions non réciproques, un sujet pour des travaux futurs.

4. Signification et Impact

Cet article apporte une contribution majeure à la physique de la matière active :

1. Correction d'un paradigme erroné : Il démontre que l'hypothèse de « forme de défaut inchangée » est non seulement incorrecte, mais qu'elle mène à des conclusions qualitativement fausses, même à l'ordre le plus bas en activité.
2. Lien forme-mouvement : Il établit une relation analytique directe entre la déformation locale du paramètre d'ordre au cœur du défaut et sa vitesse macroscopique.
3. Universalité : Les résultats suggèrent que l'auto-propulsion des défauts +1/2 est une propriété générique des nématiques actifs, indépendante des détails spécifiques du modèle d'activité, pourvu que la forme du défaut soit correctement prise en compte.
4. Extension potentielle : Bien que l'étude se concentre sur les systèmes « secs », les auteurs indiquent que leurs méthodes s'appliquent également aux systèmes « humides » (avec écoulement fluide) et à d'autres types d'ordres liquides cristallins (hexatiques, smectiques), où la forme des défauts joue un rôle crucial dans l'organisation collective (amas, lignes de défauts).

En résumé, le titre « Eppur si muove » (Et pourtant, elle bouge) fait référence au fait que, contrairement aux prédictions basées sur des modèles simplifiés, les défauts actifs bougent toujours, et que la clé pour comprendre ce mouvement réside dans la compréhension de leur forme déformée par l'activité.