Nonlinear dynamic elastic moduli from equilibrium stress fluctuations

Ce papier dérive des expressions de fonctions de corrélation en temps transitoire pour les modules élastiques dynamiques linéaires et non linéaires en utilisant les fluctuations de contrainte à l'équilibre et les équations DOLLS/SLLOD, permettant ainsi le calcul des réponses viscoélastiques anharmoniques à partir de simulations de dynamique moléculaire à l'équilibre sans protocoles de déformation hors équilibre explicites.

L'essentiel

Imaginez que vous avez un trampoline géant et invisible composé de milliards de petits ressorts et de billes (des atomes) qui rebondissent. Vous voulez savoir comment ce trampoline réagit lorsque vous le poussez ou le tirez. Rebondit-il instantanément ? Oscille-t-il ? Devient-il mou ou rigide selon la force de votre poussée ?

Dans le monde de la physique, ces réactions sont appelées les modules élastiques et viscoélastiques. Habituellement, pour les mesurer, les scientifiques doivent étirer ou comprimer physiquement le matériau dans une simulation informatique et observer ce qui se produit. C'est comme essayer de comprendre le fonctionnement d'un moteur de voiture en le percutant contre un mur encore et encore. Cela fonctionne, mais c'est désordonné, coûteux et difficile à contrôler.

Cet article présente une nouvelle astuce ingénieuse pour déterminer ces réactions sans jamais pousser réellement le matériau.

L'astuce du « voyage dans le temps »

Les auteurs (Garbuzov et Beltukov) ont découvert un raccourci mathématique. Ils ont réalisé que si l'on observe simplement le matériau au repos à température ambiante (à l'équilibre), les petits secousses et fluctuations aléatoires des atomes contiennent toutes les informations secrètes dont vous avez besoin.

Pensez-y ainsi : si vous vous tenez dans une pièce bondée et observez les gens se heurter les uns aux autres de manière aléatoire, vous pouvez en fait prédire comment la foule réagirait si quelqu'un commençait soudainement à les bousculer. Vous n'avez pas besoin de commencer à pousser pour connaître la réponse ; les chocs aléatoires contiennent déjà le plan directeur.

Le problème qu'ils ont résolu

Les scientifiques savaient déjà comment utiliser ces « chocs aléatoires » pour prédire :

1. Les réactions statiques : Comment le matériau se comporte lorsque vous le poussez et le maintenez immobile.
2. Les réactions simples et linéaires : Comment il se comporte lorsque vous le poussez doucement et rapidement.

Mais il existait un énorme vide. Personne ne savait comment utiliser les chocs aléatoires pour prédire les réactions complexes et changeantes. Que se passe-t-il si vous poussez le matériau, puis le tirez, puis le poussez plus fort, le tout selon un rythme ? C'est ce qu'on appelle la réponse dynamique non linéaire. C'est comme essayer de prédire le comportement d'un élastique si vous l'étirez, le laissez se détendre, puis l'étirez à nouveau alors qu'il vibre encore. Jusqu'à présent, il n'existait aucune formule pour calculer cela simplement en observant le matériau au repos.

La solution : une nouvelle recette

Les auteurs ont dérivé une nouvelle « recette » (une formule mathématique) qui agit comme un traducteur.

• Les ingrédients : Ils examinent la contrainte (la pression interne) et les termes de Born-Kinetic (une manière élégante de décrire l'énergie combinée des positions des atomes et de leurs vitesses).
• Le processus : Ils calculent comment ces ingrédients sont corrélés entre eux au fil du temps. C'est comme écouter le rythme des chocs aléatoires.
• Le résultat : Ils obtiennent une formule qui vous dit exactement comment le matériau réagira à n'importe quelle poussée ou traction complexe et changeante dans le temps, simplement en analysant les données d'une simulation calme et non perturbée.

Pourquoi c'est important (selon l'article)

L'article affirme qu'il s'agit d'une amélioration majeure car :

1. C'est plus sûr et moins cher : Vous n'avez pas besoin d'exécuter des simulations de « déformation » coûteuses et difficiles où vous étirez physiquement le matériau. Vous exécutez simplement une simulation standard du matériau au repos.
2. C'est plus précis : Lorsque vous essayez d'étirer des matériaux très légèrement dans une simulation, le signal est souvent faible et bruyant (comme essayer d'entendre un murmure dans une tempête). En utilisant la méthode des « chocs aléatoires », vous obtenez une image plus claire sans le bruit.
3. Cela unifie tout : Leur formule est une « clé universelle ». Si vous tournez les boutons vers une fréquence nulle, elle devient l'ancienne formule statique. Si vous désactivez les parties complexes, elle devient l'ancienne formule linéaire. Mais elle déverrouille également la porte vers le monde complexe et non linéaire qui était auparavant fermé.

L'essentiel

Cet article offre aux scientifiques un nouvel outil pour prédire le comportement des matériaux sous des forces complexes et changeantes. Au lieu de « briser » le matériau dans un ordinateur pour voir comment il réagit, ils peuvent maintenant simplement « écouter » les vibrations naturelles et aléatoires du matériau pour prédire son comportement futur. Cela transforme une pièce chaotique et bruyante d'atomes rebondissants en un manuel d'instructions clair sur la façon dont le matériau répondra au monde.

Résumé technique

Résumé technique : Modules élastiques dynamiques non linéaires à partir des fluctuations de contrainte à l'équilibre

Énoncé du problème
Alors que les formules de fluctuation pour le calcul des modules élastiques et viscoélastiques à partir de simulations de dynamique moléculaire (DM) à l'équilibre sont bien établies pour les cas linéaires et non linéaires quasi-statiques (fréquence nulle), ainsi que pour les modules dynamiques linéaires (dépendants du temps), un écart théorique majeur subsiste. Plus précisément, aucune formule de fluctuation n'existe pour les modules dynamiques non linéaires. Ces modules sont nécessaires pour décrire la réponse viscoélastique anharmonique des matériaux soumis à des déformations dépendantes du temps, petites mais finies. Les méthodes existantes pour obtenir ces propriétés reposent généralement sur des protocoles de déformation hors équilibre explicites, qui peuvent être coûteux en calcul ou difficiles à contrôler. Ce travail vise à dériver les expressions de fluctuation manquantes pour ces modules dynamiques non linéaires.

Méthodologie
Les auteurs dérivent les expressions requises en partant des équations microscopiques du mouvement pour un système de particules en interaction.

1. Équations du mouvement : La dérivation utilise les équations du mouvement de DOLLS/SLLOD. Les auteurs restreignent leur analyse aux déformations irrotationnelles (où le tenseur gradient de vitesse $L$ est symétrique), condition sous laquelle les formulations DOLLS et SLLOD coïncident. Cela permet d'utiliser un cadre hamiltonien augmenté d'un terme de couplage entre les moments des particules et le champ de vitesse d'écoulement.
2. Variables « ramenées en arrière » (Pulled-Back) : Pour gérer la dépendance temporelle de la déformation, les auteurs emploient une représentation « ramenée en arrière » des variables de l'espace des phases (coordonnées $\vec{R}_a$ et moments $\vec{P}_a$). Cette transformation mappe l'état déformé actuel vers une configuration de référence, garantissant que les variables restent continues même lorsqu'une déformation est appliquée soudainement.
3. Distribution hors équilibre : Les auteurs construisent la fonction de distribution de probabilité hors équilibre en faisant évoluer la distribution d'équilibre canonique initiale. Ils expriment la moyenne hors équilibre du tenseur des contraintes comme une moyenne sur l'ensemble d'équilibre, pondérée par l'exponentielle du changement d'énergie dû à l'évolution du point de phase ($\Delta H_\Gamma$).
4. Développement perturbatif : Le tenseur des contraintes et le changement d'énergie sont développés en série de puissances par rapport au tenseur de déformation finie de Lagrange $E$. Les auteurs séparent les contributions en termes linéaires et non linéaires, identifiant des tenseurs « Born–cinétiques » ($C^{BK}$ et $N^{BK}$) qui représentent des constantes élastiques instantanées dérivées des termes d'énergie potentielle et cinétique.
5. Dérivation de la fonction de corrélation : En substituant ces développements dans l'expression de la moyenne des contraintes hors équilibre et en conservant les termes jusqu'au second ordre en déformation, les auteurs dérivent des formules explicites pour les modules dynamiques en tant que fonctions de corrélation dépendantes du temps impliquant le tenseur des contraintes et les termes Born–cinétiques.

Contributions et résultats clés
Le résultat principal de l'article est la dérivation d'expressions de fluctuation sous forme fermée pour les modules élastiques dynamiques linéaires et non linéaires :

• Modules dynamiques linéaires ($C_{ijkl}$) : Les auteurs retrouvent la formule de fluctuation viscoélastique linéaire connue (Éq. 45), qui exprime le module comme une combinaison de la moyenne d'équilibre du tenseur Born–cinétique et des fonctions de corrélation temporelle du tenseur des contraintes.
• Modules dynamiques non linéaires ($N_{ijklmn}$) : La contribution centrale est la dérivation de la formule de fluctuation pour le tenseur de module dynamique non linéaire d'ordre six (Éq. 46). Cette expression implique :
  • La moyenne d'équilibre du tenseur Born–cinétique non linéaire ($N^{BK}$).
  • Des corrélations croisées entre le tenseur des contraintes et le tenseur Born–cinétique linéaire.
  • Des corrélations d'ordre supérieur impliquant des produits de tenseurs des contraintes et la différence entre le tenseur Born–cinétique à différents instants.
  • Des termes proportionnels à $\beta^2$ (l'inverse du carré de la température) impliquant des corrélations triples de contraintes.

Les tenseurs dérivés satisfont des propriétés de symétrie spécifiques : les tenseurs Born–cinétiques possèdent à la fois des symétries mineures et majeures, tandis que les modules dynamiques dépendants du temps complets ne satisfont que les symétries mineures (et une symétrie de permutation spécifique pour le tenseur non linéaire).

Signification et affirmations
L'article affirme que ces formules dérivées constituent une extension naturelle des travaux antérieurs, unifiant des résultats séparés pour les modules linéaires quasi-statiques, non linéaires quasi-statiques et dynamiques linéaires en un cadre cohérent unique.

Les auteurs soulignent l'utilité pratique de ces résultats :

• Calcul uniquement à l'équilibre : Les formules permettent le calcul de la réponse anharmonique dépendante du temps complète des matériaux entièrement à partir de simulations DM à l'équilibre. Cela élimine le besoin de simulations de déformation hors équilibre explicites.
• Applicabilité : Cette approche est particulièrement avantageuse pour les systèmes où les protocoles hors équilibre sont coûteux, difficiles à contrôler, ou souffrent de rapports signal/bruit faibles à petites déformations.
• Accessibilité : Les fonctions de corrélation requises n'impliquent que des moyennes d'équilibre des contraintes, des termes Born–cinétiques et de leurs corrélations temporelles, tous facilement accessibles dans les simulations d'équilibre standard.

Les auteurs notent que, bien que le travail actuel soit restreint aux écoulements irrotationnels, le formalisme fournit une base pour de futures extensions aux écoulements rotationnels et aux modules d'ordre supérieur, ainsi que pour des applications aux solides amorphes tels que les polymères et les matériaux biologiques. Ils identifient également la mise en œuvre numérique et la validation par rapport à des simulations hors équilibre directes comme la prochaine étape nécessaire.