Emergent cohesion via self-caging in maximally entangled rod packings

Résumé Technique : Cohésion Émergente via l'Auto-Encagement dans les Empilements de Bâtonnets Maximalement Enchevêtrés

Énoncé du Problème
Les empilements désordonnés de bâtonnets rigides présentent une forme de rigidité collective distincte du blocage (jamming) observé dans les particules sphériques. Alors que les empilements de sphères bloquées reposent sur des frontières fixes ou la limite thermodynamique pour empêcher l'échappement des grains, les empilements de bâtonnets finis avec des frontières libres peuvent rester cohésifs et auto-équilibrés (par exemple, un nid de brindilles). Cette cohésion provient de l'« enchevêtrement », où les bâtonnets s'obstruction mutuellement le mouvement par la géométrie, le contact et la friction. Cependant, malgré son ubiquité, l'enchevêtrement manque de protocoles systématiques pour sa génération et de liens quantitatifs avec la stabilité mécanique. Le problème central abordé est l'identification du motif géométrique minimal qui sous-tend cette cohésion et la compréhension de la manière dont des interactions purement répulsives et frictionnelles peuvent soutenir une structure autoporteuse sans forces attractives.

Méthodologie
Les auteurs proposent un protocole pour générer des empilements de bâtonnets maximalement enchevêtrés en optimisant une fonctionnelle géométrique soumise à des contraintes de non-pénétration.

1. Définition du Système : Le système consiste en $N$ bâtonnets sphérocylindriques rigides de diamètre $d$ et de longueur $l$, caractérisés par des coordonnées généralisées $q = \{x_i, y_i, z_i, \theta_i, \phi_i\}$. Le rapport d'aspect est défini par $\alpha = l/d$.
2. Quantification de l'Enchevêtrement : L'enchevêtrement est quantifié par le Nombre de Croisements Moyen (ACN, Average Crossing Number), une mesure topologique définie pour des paires de bâtonnets $i$ et $j$ comme :
   $$ACN_{ij} = \frac{1}{4\pi} \int_{-l/2}^{l/2} \int_{-l/2}^{l/2} \frac{|(\mathbf{t}_i \times \mathbf{t}_j) \cdot (\mathbf{r}_i - \mathbf{r}_j)|}{|\mathbf{r}_i - \mathbf{r}_j|^3} ds dt$$
   L'enchevêtrement normalisé total $\tilde{e}(q)$ est la somme des ACN par paires sur le réseau.
3. Protocole d'Optimisation : Partant de configurations d'adsorption séquentielle aléatoire, le système évolue via un flot de gradient (implémenté à l'aide du Fast Inertial Relaxation Engine, FIRE) pour maximiser $\tilde{e}(q)$. Pour garantir le réalisme physique, une potentiel de répulsion harmonique unidirectionnel est imposé pour empêcher le chevauchement des bâtonnets ($d_{ij} < d$).
4. Analyse de Stabilité : La stabilité dynamique des configurations résultantes maximalement enchevêtrées ($q^{**}$) est sondée en soumettant ces dernières à des perturbations aléatoires et en les faisant évoluer sous une dynamique de corps rigide avec des collisions élastiques et une friction de Coulomb. La stabilité est mesurée par le « temps d'échelle de désenchevêtrement » ($t_u$) et la rétention de l'enchevêtrement au cours du temps.

Contributions Clés et Résultats

• Identification de l'Auto-Encagement : Le protocole d'optimisation génère des empilements qui présentent un « auto-encagement » — un état où les contraintes géométriques collectives empêchent l'échappement des bâtonnets, même dans des systèmes finis avec des frontières libres. Contrairement aux sphères bloquées où les grains de surface peuvent toujours s'échapper, les bâtonnets enchevêtrés bloquent mutuellement les mouvements latéraux.
• Régimes Structurels et Crossover : Les empilements résultants affichent deux régimes structurels distincts selon le paramètre $N/(Z\alpha)$, où $Z$ est le nombre de coordination moyen :
  • Régime Hérisson ($N/(Z\alpha) > 1/3$) : Un cœur dense et riche en contacts est entouré d'une couche radialement orientée de type « épineux ».
  • Régime Isotrope ($N/(Z\alpha) < 1/3$) : La structure de la couche disparaît, et l'empilement devient sphérique et compact.
  • Le crossover se produit précisément à $N/(Z\alpha) = 1/3$, ce qui correspond au point où le volume libre dans les espaces de configuration translationnels et rotationnels est minimisé.
• Métriques d'Encagement Géométrique : Les auteurs définissent formellement l'« auto-encagement » en analysant l'espace de configuration accessible ($D_k$) pour chaque bâtonnet. Un bâtonnet est encagé si sa surface de translation libre ($A_k$) et son angle solide de rotation ($\Omega_k$) sont finis. L'étude introduit des rapports de gap normalisés par la longueur, $\tilde{g}_t = \sqrt{A^*}/l$ et $\tilde{g}_r = \sqrt{\Omega^*/2\pi}$, où $A^*$ et $\Omega^*$ sont les volumes libres maximaux dans l'empilement.
• Mécanisme de Cohésion : Il est démontré que la cohésion mécanique provient d'un couplage géométrique-frictionnel :
  • Confinement Géométrique : Dans le régime d'auto-encagement, les bâtonnets sont piégés latéralement, générant des contacts soutenus.
  • Dissipation Frictionnelle : Ces contacts soutenus permettent à la friction de dissiper le mouvement relatif (spécifiquement le glissement axial), arrêtant les degrés de liberté restants.
  • Condition de Stabilité : Les simulations montrent que les empilements avec de faibles $\tilde{g}_t$ et $\tilde{g}_r$ restent cohésifs, tandis que ceux avec de grands volumes libres se désenchevêtrent rapidement. Crucialement, en l'absence de friction, l'enchevêtrement décroît jusqu'à zéro quel que soit le confinement géométrique, démontrant que la géométrie seule est insuffisante ; la friction est nécessaire pour arrêter le mouvement axial permis par la cage.
• Lois d'Échelle : L'étude établit que la transition d'un état cohésif et auto-encagé vers un état désordonné et mobile est gouvernée par l'échelle des volumes libres près du ratio critique $N/(Z\alpha) \approx 1/3$.

Signification
Cet article établit un motif géométrique minimal pour la cohésion induite par l'enchevêtrement dans les empilements de bâtonnets athermiques. Il démontre que la matière granulaire cohésive peut exister sans forces attractives, à condition que le système atteigne une configuration géométrique spécifique (auto-encagement) qui, couplée à la friction, crée une structure autoporteuse mécaniquement stable. Le cadre lie directement les contraintes géométriques locales (volume libre dans l'espace de configuration) aux taux de collision et à la dissipation frictionnelle, offrant une base quantitative pour comprendre la stabilité des systèmes enchevêtrés allant de la matière granulaire aux filaments biologiques. Les auteurs notent que cette perspective s'étend naturellement aux particules non convexes, aux filaments flexibles et aux systèmes actifs où l'enchevêtrement régit la réponse mécanique.