Nonlinear dynamic elastic moduli from equilibrium stress fluctuations

本文利用平衡态应力涨落以及DOLLS/SLLOD方程，导出了线性和非线性动态弹性模量的瞬态时间关联函数表达式，从而使得能够从平衡态分子动力学模拟中计算非谐粘弹性响应，而无需显式的非平衡形变协议。

要点

想象你有一个由数十亿个微小弹簧和球体（原子）弹跳组成的巨大、无形的蹦床。你想知道当你推或拉这个蹦床时，它会如何反应。它是瞬间弹回？还是会摇晃？或者根据你施加力度的不同，它会变得柔软还是坚硬？

在物理学中，这些反应被称为弹性模量和粘弹性模量。通常，为了测量它们，科学家必须在计算机模拟中实际拉伸或挤压材料，并观察会发生什么。这就像试图通过反复将汽车撞向墙壁来了解发动机的工作原理。这种方法可行，但混乱、昂贵且难以控制。

本文介绍了一种巧妙的新技术，可以在不实际推动材料的情况下确定这些反应。

“时间旅行”技巧

作者（Garbuzov 和 Beltukov）发现了一个数学捷径。他们意识到，如果你只是观察材料在室温下静止（处于平衡态）时的状态，原子那微小、随机的抖动和涨落就包含了你所需的所有秘密信息。

可以这样理解：如果你站在一个拥挤的房间里，观察人们随机地互相碰撞，你实际上可以预测，如果突然有人开始推搡他们，人群会如何反应。你不需要开始推搡就能知道答案；随机的碰撞已经包含了蓝图。

他们解决的问题

科学家们已经知道如何利用这些“随机碰撞”来预测：

1. 静态反应：当你推动材料并保持静止时，材料的感觉如何。
2. 简单、线性的反应：当你轻柔且快速地推动材料时，它的感觉如何。

但存在一个巨大的空白。没有人知道如何利用随机碰撞来预测复杂、变化的反应。如果你先推材料，再拉它，然后更用力地推它，所有这些动作都按照某种节奏进行，会发生什么？这被称为非线性动态响应。这就像试图预测一根橡皮筋的行为：如果你拉伸它，让它弹回，然后在它仍在振动时再次拉伸它。直到现在，还没有一种公式可以通过仅观察材料静止不动来计算这种情况。

解决方案：新配方

作者推导出了一个新的“配方”（数学公式），它充当翻译的角色。

• 原料：他们观察应力（内部压力）和Born-动能项（一种描述原子位置及其速度组合能量的复杂方式）。
• 过程：他们计算这些原料随时间相互关联的方式。这就像聆听随机碰撞的节奏。
• 结果：他们得到了一个公式，可以精确告诉你材料将如何对任何复杂、随时间变化的推或拉做出反应，而这一切只需通过分析来自平静、未受干扰的模拟数据即可实现。

为什么这很重要（根据论文所述）

该论文声称这是一个重大升级，因为：

1. 更安全、更经济：你不需要运行昂贵且困难的“变形”模拟，即实际拉伸材料。你只需运行材料静止状态的标准模拟。
2. 更准确：当你尝试在模拟中非常轻微地拉伸材料时，信号往往微弱且充满噪声（就像试图在风暴中听到耳语）。通过使用“随机碰撞”方法，你可以在没有噪声的情况下获得更清晰的图像。
3. 统一一切：他们的公式是一把“万能钥匙”。如果你将旋钮调至零频率，它就变成了旧的静态公式。如果你关闭复杂部分，它就变成了旧的线性公式。但它也解锁了此前被锁住的复杂非线性世界的大门。

核心结论

这篇论文为科学家提供了一种新工具，用于预测材料在复杂、变化力作用下的行为。以前，他们需要在计算机中“破坏”材料以观察其反应；而现在，他们只需“聆听”材料自然、随机的振动，即可预测其未来的行为。它将一个充满弹跳原子的混乱、嘈杂的房间，转化为一份关于材料将如何响应世界的清晰说明书。

技术摘要

技术摘要：基于平衡态应力涨落的非线性动态弹性模量

问题陈述
虽然利用平衡态分子动力学（MD）模拟计算弹性和粘弹性模量的涨落公式，在准静态（零频）线性和非线性情形以及线性动态（时间相关）模量方面已确立，但仍存在显著的理论空白。具体而言，目前尚无针对非线性时间相关模量的涨落公式。这些模量对于描述材料在微小但有限的时间相关应变作用下的非谐粘弹性响应至关重要。现有的获取这些性质的方法通常依赖于显式的非平衡变形协议，这可能计算成本高昂或难以控制。本研究旨在推导这些非线性动态模量所缺失的涨落表达式。

方法论
作者从相互作用粒子系统的微观运动方程出发，推导了所需的表达式。

1. 运动方程： 推导利用了 DOLLS/SLLOD 运动方程。作者将分析限制在无旋变形（即速度梯度张量 $L$ 为对称张量）的条件下，在此条件下 DOLLS 和 SLLOD 表述是一致的。这使得可以采用哈密顿框架，并辅以粒子动量与流动速度场之间的耦合项。
2. 拉回变量： 为了处理应变的时间依赖性，作者采用了相空间变量（坐标 $\vec{R}_a$ 和动量 $\vec{P}_a$）的“拉回”表示。该变换将当前的变形状态映射回参考构型，确保即使在突然施加应变时，变量仍保持连续。
3. 非平衡分布： 作者通过演化初始正则平衡分布来构建非平衡概率分布函数。他们将应力张量的非平衡平均值表示为对平衡系综的平均，并以相点演化引起的能量变化（$\Delta H_\Gamma$）的指数作为权重。
4. 微扰展开： 应力张量和能量变化针对拉格朗日有限应变张量 $E$ 进行幂级数展开。作者将贡献分离为线性和非线性项，并识别出“玻恩 - 动能”张量（$C^{BK}$ 和 $N^{BK}$），它们代表源自势能和动能项的瞬时弹性常数。
5. 关联函数推导： 通过将上述展开式代入非平衡应力平均表达式，并保留至应变的二阶项，作者推导出了动态模量的显式公式，这些公式表现为涉及应力张量和玻恩 - 动能项的时间相关函数。

主要贡献与结果
本文的主要成果是推导出了线性和非线性动态弹性模量的闭式涨落表达式：

• 线性动态模量（$C_{ijkl}$）： 作者恢复了已知的线性粘弹性涨落公式（公式 45），该公式将模量表示为玻恩 - 动能张量的平衡态平均值与应力张量时间关联函数的组合。
• 非线性动态模量（$N_{ijklmn}$）： 核心贡献是推导了六阶非线性动态模量张量的涨落公式（公式 46）。该表达式涉及：
  • 非线性玻恩 - 动能张量（$N^{BK}$）的平衡态平均值。
  • 应力张量与线性玻恩 - 动能张量之间的交叉关联。
  • 涉及应力张量乘积以及不同时刻玻恩 - 动能张量之差的更高阶关联。
  • 与 $\beta^2$（温度平方的倒数）成比例的项，涉及三重应力关联。

推导出的张量满足特定的对称性质：玻恩 - 动能张量同时具有次要对称性和主要对称性，而完整的时间相关动态模量仅满足次要对称性（非线性张量还具有特定的置换对称性）。

意义与主张
本文声称，这些推导出的公式构成了先前工作的自然延伸，将准静态线性、准静态非线性和线性动态模量的独立结果统一到一个连贯的框架中。

作者强调了这些结果的实用价值：

• 仅基于平衡态的计算： 这些公式允许完全从平衡态 MD 模拟中计算材料的全时间相关非谐响应。这消除了对显式非平衡变形运行的需求。
• 适用性： 这种方法对于非平衡协议成本高昂、难以控制或在微小应变下信噪比差的系统尤为有利。
• 可及性： 所需的关联函数仅涉及应力、玻恩 - 动能项及其时间关联的平衡态平均值，所有这些在标准平衡态模拟中均易于获取。

作者指出，虽然当前工作局限于无旋流动，但该形式体系为未来扩展到旋转流动和高阶模量，以及应用于聚合物和生物材料等非晶固体奠定了基础。他们还指出，数值实施并与直接非平衡模拟进行验证是必要的下一步。