Dimension-dependent continuum limits in tissue mechanics

Résumé technique : Limites de continuum dépendantes de la dimension en mécanique des tissus

Énoncé du problème
Les modèles discrets basés sur les cellules lient efficacement les forces microscopiques à l'évolution à l'échelle du tissu, mais ils sont souvent difficilement exploitables sur le plan computationnel pour un grand nombre de cellules et complexes à paramétrer. Les modèles de continuum offrent une alternative pratique pour l'inférence au niveau de la population et la révélation des liens de transport émergents. Cependant, alors que les limites de continuum en 1D pour la mécanique des tissus épithéliaux sont bien étudiées, les problèmes correspondants en dimensions supérieures restent largement inexplorés. Un fossé critique existe dans la compréhension de la manière dont la dimensionnalité spatiale influence le lien entre les forces à l'échelle cellulaire et le transport à l'échelle tissulaire, particulièrement concernant les diffusivités non linéaires de type milieu poreux, $D(\rho) \propto \rho^\gamma$. La plupart des applications phénoménologiques existantes en mécanique des fluides et en biologie des populations sont limitées à $\gamma > 0$, alors que la mécanique épithéliale motive à la fois $\gamma > 0$ et $\gamma < 0$. Les auteurs étudient si l'exposant de transport macroscopique $\gamma$ est déterminé uniquement par la mécanique microscopique ou s'il est fondamentalement altéré par la dimensionnalité spatiale.

Méthodologie
Les auteurs emploient une stratégie consistant à faire correspondre les taux de relaxation exponentielle à long terme entre les modèles mécaniques discrets et les modèles d'équations aux dérivées partielles (EDP) macroscopiques au voisinage d'un état stationnaire.

1. Modèles discrets :

   • 1D : Le tissu est modélisé comme une chaîne de $N$ liens mécaniques (cellules) avec des limites fixes. Les nœuds évoluent selon un modèle suramorti piloté par une loi de force d'interaction $f(\ell)$. Les auteurs utilisent une famille paramétrée de lois de force, $f(\ell) = k(\ell^\nu - a^\nu)/(\nu a^{\nu-1})$ (avec la limite logarithmique pour $\nu=0$), qui se linéarise en la loi de Hooke près de la longueur de repos $a$.
   • 2D : Les cellules sont modélisées comme un réseau de triangulation où les nœuds interagissent via des forces mécaniques par paires. La dynamique est régie par une équation suramortie similaire sommant les forces sur les nœuds adjacents.
   • Analyse : Les auteurs linéarisent les systèmes discrets autour de leurs états stationnaires uniformément espacés. La relaxation à long terme est gouvernée par la plus grande valeur propre du Jacobien du système (ou de la matrice Laplacienne en 2D), produisant un taux de décroissance exponentielle.

2. Modèles macroscopiques :

   • La densité du tissu $\rho(\mathbf{r}, t)$ est décrite par une équation de diffusion non linéaire : $\partial_t \rho = \nabla \cdot (D(\rho) \nabla \rho)$.
   • Les auteurs supposent une forme générale pour la diffusivité non linéaire : $D(\rho) = a^\beta \rho^\gamma / \tau$, où $a$ est l'échelle de longueur microscopique, $\rho$ la densité, et $\tau$ une échelle de temps caractéristique.
   • En linéarisant l'EDP au voisinage de l'état stationnaire $\bar{\rho}$, les auteurs démontrent que le taux de relaxation dépend de $D(\bar{\rho})$ et de la géométrie du domaine (spécifiquement les valeurs propres du Laplacien).

3. Procédure de grossissement (Coarse-Graining) :

   • Le cœur de la méthode consiste à égaliser les taux de relaxation dérivés des simulations discrètes (dépendants des paramètres microscopiques $k, \eta, a, \nu$ et de la dimension $d$) avec les taux de relaxation prédits par l'EDP macroscopique.
   • En faisant varier la densité à l'état stationnaire $\bar{\rho}$ (via le changement de $N$ ou de la taille du domaine $L$) tout en maintenant les paramètres microscopiques constants, les auteurs infèrent la forme fonctionnelle de $D(\rho)$ et l'exposant $\gamma$.

Résultats clés

• Dépendance dimensionnelle de l'exposant : L'étude établit que l'exposant $\gamma$ dans la diffusivité non linéaire $D(\rho) \propto \rho^\gamma$ est déterminé conjointement par l'exposant de la loi de force microscopique $\nu$ et la dimension spatiale $d$.
• Relation d'échelle : L'analyse dimensionnelle et les données de simulation révèlent une relation d'échelle universelle : $\gamma = \beta - 2/d$, où $\beta$ est un exposant d'échelle lié à la loi de force microscopique.
  • En 1D ($d=1$), pour une loi de force avec $\nu=1$ (loi de Hooke), la diffusivité résultante est $D(\rho) \propto \rho^{-2}$.
  • En 2D ($d=2$), la même loi de force microscopique ($\nu=1$) produit un $D(\rho) \propto \rho^{-1.11}$.
  • Généralement, à mesure que la dimension $d$ augmente, l'amplitude de la non-linéarité $|\gamma|$ diminue.
• Émergence de la diffusion linéaire : Lorsque $d \to \infty$, $\gamma \to 0$, impliquant que la diffusion linéaire devient de plus en plus pertinente dans les dimensions supérieures pour les mêmes mécaniques microscopiques.
• Validation : Des simulations numériques en 1D et 2D confirment que le $D(\rho)$ inféré reproduit fidèlement la dynamique de relaxation discrète. Les résultats sont valables pour diverses conditions initiales et tailles de domaines, confirmant que le transport macroscopique dépend des échelles de longueur locales ($a$ et $\rho$) plutôt que des échelles globales du domaine.

Signification et affirmations
L'article prétend révéler un nouveau mécanisme physique pour les phénomènes de transport macroscopique émergents dans la mécanique des tissus. La signification primaire réside dans la démonstration que la dimensionnalité n'est pas simplement un détail géométrique, mais qu'elle altère fondamentalement la loi macroscopique obtenue après grossissement.

• Limites des modèles 1D : Les enseignements tirés des modèles 1D ne peuvent être supposés valables sans changement dans les dimensions supérieures. Spécifiquement, l'utilisation de modèles linéaires (loi de Hooke) en 1D masque des détails critiques sur la manière dont les échelles de longueur microscopiques se transforment en transport émergent à l'échelle du tissu.
• Nouvelle méthode de grossissement : Les auteurs proposent une méthode pour transformer la dynamique microscopique en fonctions de diffusion non linéaires en comparant les taux de relaxation. Cette approche ne repose pas sur les limites de continuum traditionnelles, ce qui la rend particulièrement adaptée aux systèmes finis.
• Mécanisme physique : Ce travail fournit une explication physique de la raison pour laquelle des diffusivités non linéaires de type milieu poreux apparaissent dans la mécanique des tissus et comment l'exposant spécifique est fixé par l'interaction entre la mécanique microscopique et la dimension spatiale.

Les auteurs concluent que la même loi d'interaction microscopique ne génère pas une limite de continuum indépendante de la dimension ; au contraire, la dimension remodèle la diffusivité non linéaire effective, changeant l'exposant dans $D(\rho) \propto \rho^\gamma$. Ce phénomène est analogue aux comportements dépendants de la dimension dans la théorie des marches aléatoires et les méthodes de groupe de renormalisation.