玻尔百科现代材料科学的核心在于一个基本问题：无数单个原子的集体相互作用如何产生我们观察到和工程应用的宏观性质，例如钢的刚度或流体的粘度？在连续介质力学中，应力被简单地定义为单位面积上的力，但当在原子尺度上观察时，这个概念变得异常复杂，因为在原子尺度上物质大部分是空的空间。维里应力张量是一个优雅的理论框架，它解决了这个悖论，为连接离散的、微观的原子世界与连续的、宏观的材料世界提供了必要的桥梁。
本文将阐释维里应力张量的强大功能和广泛适用性。首先，在“原理与机制”部分，我们将探讨其理论基础，从基本原理出发推导该张量，并剖析其动能和势能分量。我们还将考察其在计算机模拟中计算的实际精妙之处。然后，在“应用与跨学科联系”部分，我们将遍览其多样化的应用，从预测材料失效和工程模拟，到探索生物学和天体物理学的前沿，揭示其作为贯穿现代科学的统一概念。
科学理解常常涉及连接不同尺度的描述。我们可以观察一根钢梁，并用连续介质力学的概念（如密度和应力）来描述它。但我们知道，钢梁并非连续介质；它是一个由原子组成的熙攘城市，一个由微小、振动的粒子稀疏占据的广阔空旷空间，这些粒子被无形的力线维系在一起。一个深刻的问题是：我们如何将这两种图像联系起来？原子的集体舞蹈如何产生我们称之为应力的宏观性质？
想象一下你正站在一座桥上。你感觉到脚下坚实的地面，但你知道它大部分是空的。这种坚实感来自于原子间巨大的力，抵抗着你体重带来的压缩。在宏观世界中，应力被定义为单位面积上的力。如果我们能用一个数学平面将桥切开，应力就是切面一侧对另一侧施加的力，除以切面的单位面积。
但是，在我们的原子城市中如何定义这个概念呢？一个平面切过它，大部分会切过真空。这时，一个更基本的概念——动量通量——为我们提供了解决方案。应力与动量的流动密切相关。再次想象我们的数学平面。动量可以通过两种方式穿过它：
动能贡献： 原子可以实际飞越该平面，随身携带其动量 ()。这就像冰雹敲打屋顶。冰雹的持续撞击会产生压力。类似地，原子的热运动产生了应力张量的动能部分。
势能贡献： 力可以跨越该平面作用。想象两个原子位于我们平面的两侧，被一个力连接。这个力在没有任何物质实际穿过平面的情况下，“传递”了动量。这是应力中更微妙、且在凝聚态物质中通常占主导地位的部分。
维里应力张量正是量化这两种贡献的绝妙数学工具，它连接了原子世界和连续介质世界。
有几种方法可以推导出应力公式，但也许最优雅且最具启发性的路径始于一个简单的问题：如果我们轻轻挤压我们的原子盒子，需要做多少功？
这种方法通过虚功原理直接与连续介质力学中的应力定义相联系。在连续介质力学中，使材料发生一个微小的对称应变张量 所需的单位体积功为 ，其中 是柯西应力张量。因此，一个体积为 的系统的总功为 。
现在，让我们从原子角度来看这个问题。一个均匀应变 使模拟盒子变形，并将每个位于位置 的原子移动一个微小位移 。系统势能的变化量 是克服作用在每个粒子上的力所做的功的总和：
其中 是作用在粒子 上的总力。通过一些张量代数， 这一项可以被重写为一个张量缩并，得到：
这里， 表示张量积或并矢积。现在我们得到了同一个功 的两个表达式。将它们相等，我们便获得了一个纯粹的洞见时刻：
由于这对于任何微小应变 都必须成立，并且知道柯西应力 必须是对称的，我们得到了一个惊人普适的应力张量构型（势能）部分的表达式：
这一项与维里有关。这个推导之所以优美，是因为我们没有对力 的性质做任何假设。它们可以是简单的对偶力，也可以是反应体系或机器学习势中出现的极其复杂的多体力。这个公式无论如何都成立。这证明了基本原理的统一力量。
让我们来整合完整的图像。总应力张量包括动能和势能两部分。在力学中，拉伸为正是一种常见约定，据此我们得到：
等等！ 你可能会在统计力学教科书中看到一个不同的公式，通常称为压力张量 ，其符号是不同的。这是一个源于不同约定的臭名昭著的混淆点。在热力学中，压力对于压缩为正，而在固体力学中，应力对于拉伸为正。对于一个处于平衡状态的各向同性系统，两者通过 根本地联系在一起。我们这里将坚持使用力学约定，但必须意识到这种符号差异。
第一项是动能应力。负号表示原子的热运动总是产生压力（一种压应力），这在拉伸为正的约定中对应一个负值。这里的 是原子 相对于局部原子群平均速度 的特异速度。认为这一项只对气体重要是错误的。在固体中，原子在不断振动，这种振动——即热运动——携带动量并对应力做出贡献。
第二项是构型应力，也称为维里应力。我们在这里为常见的对偶中心力情况重写了它，其中 是原子 对原子 的作用力，而 。这一项是所有相互作用原子对的总和。其贡献的符号取决于力的性质：排斥力（）导致压缩（负应力），而吸引力（）导致拉伸（正应力）。并矢积 优雅地捕捉了此应力贡献的方向和大小。对于中心力（沿着连接粒子的直线作用），力矢量 平行于位置矢量 。这带来一个奇妙的结果：由此产生的张量 是对称的，这意味着总应力张量也是对称的，正如连续介质力学中角动量守恒所要求的那样。
一个源于势能平移不变性的重要性质是，系统中所有内力的总和为零：。一个巧妙的推论是，维里 的值不依赖于我们坐标系原点的选择，这使其成为一个物理上稳健的量。
维里表达式不仅仅是一个理论上的奇珍；它是在几乎每一次分子动力学模拟中计算压力和应力的主力。但在实践中应用它，会揭示出更多的精妙之处和更深层的美。
如果我们想知道的不是整个盒子里的应力，而只是一个微小区域的应力，比如说在一个纳米级裂纹的尖端，该怎么办？我们不能直接使用上面的公式。我们需要一个局域版本。Irving-Kirkwood-Hardy 公式提供了答案。我们定义一个小的数学控制体。然后，通过对控制体内部所有原子的贡献（对于动能部分）以及至关重要的、所有穿过控制体边界的原子间键的贡献求和，来计算局域应力。这个方法正确地划分了应力，并允许我们创建材料在纳米尺度上的应力图，这是纳米力学的一个强大工具。
维里应力的准确性完全取决于力 的准确性。我们对力所做的任何近似都将反映在应力中。
短程截断： 对于许多势，力随距离迅速衰减。为了节省计算时间，我们通常会简单地忽略超过某个截断半径 的相互作用。但是，如果我们只是突然截断力，就会产生一个非物理的不连续性。这个看似微小的“罪过”会带来严重后果，给计算出的性质（如弹性模量）引入系统误差。一种更精巧的方法是使用一个平滑的开关函数，将力和势平缓地降至零，这能得到更准确的结果。
长程力： 那么像 这样缓慢衰减的静电相互作用呢？我们不能简单地截断它们。物理学家使用一个巧妙的数学技巧，称为Ewald 求和，它将计算分为在实空间中计算的短程部分和在傅里叶波的“倒易”空间中计算的长程部分。维里原理要求我们保持一致：如果力有实空间和倒易空间部分，那么应力张量也必须有来自这两部分的贡献。忘记倒易空间项是一个常见且严重的错误，会导致不正确的压力和模拟中的假象。
复杂势： 现代模拟经常使用非常复杂的势，其中相互作用不是固定的。在反应力场中，两个原子之间的“键级”或每个原子上的电荷可以根据局部环境动态变化。人们可能会担心这需要在维里应力中添加复杂的修正项。但在这里，大自然是仁慈的。一个与Hellmann-Feynman 定理相关的深刻结果表明，因为这些内部变量（如电荷）被持续优化以保持能量最小化，它们的导数不会显式地出现在最终的力表达式中。标准的维里公式，当与正确的总力一起使用时，仍然有效！
约束力： 在模拟聚合物或水时，我们经常将化学键建模为固定长度的刚性棒。强制实现这种刚性的力称为约束力。它们对应力有贡献吗？绝对有。它们是系统内力的一部分。维里框架可以优雅地扩展以包含它们，其中对应力的贡献可以用施加约束的拉格朗日乘子来表示。
因此，维里应力张量是一个深刻而强大的概念。它提供了连接原子和力的微观世界与材料和力学的宏观世界的基本纽带。它展示了无数粒子在物理定律支配下的集体推拉，如何产生了我们日常体验到的世界的稳健性质。它从第一性原理的推导以及对最复杂和现代模拟技术的适应性，是物理学统一性与优雅之美的美丽例证。
既然我们已经熟悉了维里应力张量的机制——它是什么以及它从何而来——我们就可以开始一段更激动人心的旅程。我们将不再问“是什么？”，而是问“有什么用？”这个优雅的数学对象有什么好处？答案，正如物理学中常有的情况一样，其广度令人惊叹。维里应力是一把万能钥匙，解开了那些乍一看似乎毫无关联的领域中的秘密。它提供了一种统一的语言来描述物质内部的推与拉，无论这些物质是一块钢，是活细胞内的液体，还是一颗垂死恒星的奇异核心。这是物理定律统一性的完美例证。
让我们从一些坚实的东西开始，一些你可以敲击的东西。我们如何描述一种材料的性质，比如它的刚度或强度？我们当然可以把一个真实的样品带到实验室，对它进行物理上的拉伸、压缩或扭转。但是，如果我们想从原子层面预测这些性质呢？这时，维里应力张量就成了计算物理学家的主要工具。
想象一下，我们建立了一个小型晶体固体的计算机模型，一个由原子组成的完美小晶格，通过它们之间的相互作用力相互作用。我们不能把这个虚拟晶体放进真实的机器里，但我们可以做一些类似的事情：我们可以使它的模拟盒子变形。如果我们对盒子施加一个微小的剪切应变——把它从一个正方形压成一个轻微的菱形——里面的原子就会被 jostle。它们之间的距离和角度发生变化，作为回应，它们之间的力也会调整。维里应力张量正是测量系统对这种变形的集体内部响应的工具。应力张量的非对角分量会改变，而剪切应力的变化量与所施加的剪切应变之比，直接给出了材料的剪切模量——它抵抗扭曲的能力。 我们实质上进行了一次虚拟力学测试，并测量了一个宏观的、真实世界的弹性性质，这个性质直接源于微观力。
但如果我们拉得更用力呢？每种材料都有其断裂点。我们也可以模拟这个过程。通过施加一个逐渐增大的拉伸应变——沿着一个轴拉伸我们的虚拟晶体——我们可以观察到拉伸应力（通常称为“负压”）的累积。原子被拉得离它们愉快的、低能量的位置越来越远。在某个点上，恢复力不再能增加。应力达到一个最大值，然后随着进一步的应变开始减小。应力-应变曲线上的这个峰值标志着材料的理想拉伸强度——它在灾难性失效前可以被拉伸的理论极限。 这是原子键集体放弃的时刻。再一次，维里应力是我们的向导，描绘了材料从温和的弹性一直到剧烈断裂的响应过程。
到目前为止，我们一直将维里应力用作一种被动的测量设备。但它在现代科学中的作用往往要主动得多。在分子动力学模拟的世界里，它是引擎本身的一个关键组成部分。
许多模拟旨在模仿压力恒定而非体积恒定的实验条件。想象一下在一个向大气开放的烧杯中发生的化学反应。为了模拟这一点，我们需要一个“恒压器”，这是一种计算算法，可以动态调整模拟盒子的大小和形状以维持目标压力。恒压器如何知道是应该扩大还是缩小盒子？它持续监测内部压力，这只是维里应力张量的迹。如果内部压力太高，恒压器就扩大盒子；如果太低，就收缩它。
对于更复杂的情况，比如研究固体中的相变，我们可能需要控制整个应力张量，而不仅仅是压力。各向异性恒压器，如著名的 Parrinello-Rahman 方法，使用完整的维里应力张量来控制模拟盒子的长度和角度，使其能够响应内力而改变形状。 维里应力张量就像一个复杂的、六分量的传感器，向模拟的控制系统提供信息。
这使得维里应力成为一个必不可少的诊断工具。我们的模拟是否适当平衡了？它是否进入了稳定状态？我们可以通过观察维里应力来找出答案。如果它的时间平均值与我们的目标值相符，并且不再显示任何系统性漂移，我们就可以确信我们的系统已经达到了力学平衡。 它是模拟的“生命体征”，告诉我们我们的虚拟实验是否健康。
科学中的一大挑战是连接不同尺度。我们知道物质是由离散的原子组成的，但对于许多工程问题，将其视为连续介质更为有用。维里应力在两个世界之间提供了一座强大而优雅的桥梁。
虽然全局应力张量为整个模拟盒子提供一个值，但在一个非均匀系统中——比如一个包含裂纹、晶界或纳米颗粒的材料——应力并非均匀。它在某些地方更高，在另一些地方更低。我们如何捕捉这一点？答案在于将维里局域化。我们可以根据作用在每个原子上的力，将总维里的一部分分配给该原子。这种“单原子”应力随后可以利用数学核函数在空间中“抹开”，就像用海绵蘸取颜料一样。其结果是一个连续的、空间分辨的应力场。 这个场以优美的细节向我们展示了应力如何围绕缺陷流动或在裂纹尖端集中。这个概念是旨在将原子模拟的准确性与连续介质力学的效率相结合的多尺度建模方法的基石。
这座桥梁甚至可以一直延伸到量子世界。当力是用量子力学（通过 Hellmann-Feynman 定理）计算时，同样的公式允许我们构建一个反映底层电子结构的局域应力场，将量子和力学现实统一在同一幅画面中。
更值得注意的是，这座桥梁是双向的。在机器学习势这一前沿领域，科学家训练神经网络来近似材料复杂的势能面。一个好的势不仅要正确地得到能量和力，还要正确地得到应力。因此，从参考量子力学模拟中计算出的维里应力被用作训练过程中的一个目标。神经网络被明确地教导要重现正确的应力，确保最终得到的模型不仅准确，而且在力学上是稳健的。[@problem_-id:102291] 在这里，维里应力不是模型的输出，而是构建模型本身的一个关键输入。
我们的讨论主要集中在处于或接近力学平衡的系统。但统计力学的真正力量往往在于其对涨落和远离平衡系统的描述。在这里，维里应力同样扮演着主角。
考虑一个处于热平衡状态的液体。原子在不停地运动，瞬时维里应力在其平均值（对于非对角分量，平均值为零）附近剧烈涨落。这些涨落不仅仅是随机噪声。它们包含了关于系统性质的深刻信息。Green-Kubo 关系，非平衡统计力学的伟大成就之一，告诉我们像粘度这样的输运系数与相应涨落通量的时间自相关函数的积分有关。对于剪切粘度——衡量流体抗流能力的指标——相关的通量就是剪切应力。通过测量维里剪切应力的涨落“记住”其过去值的时间长度，我们可以直接计算流体的粘度。 这是一个深刻的联系：微观的、短暂的应力涨落之舞决定了宏观的、稳态的粘度性质。
维里概念甚至可以被推向活性物质领域——由自驱动个体组成的系统，如游泳的细菌或合成的微型机器人。这些系统本质上是远离平衡的，因为每个粒子都是一个小引擎，不断消耗能量并将其转化为运动。恒定的自驱动力对系统的总应力有贡献。这种贡献被称为“游泳应力”，是活性的一种独特的力学特征。它产生了一种在平衡系统中没有对应物的“游泳压力”，并可导致有趣的集体行为。通过扩展维里形式体系，我们可以剖析生命和活性系统的力学状态，为统计物理学和生物学之间架起一座桥梁。
从无限小，让我们做最后一次飞跃，到不可想象的巨大。维里定理起源于天体物理学，它将像恒星或星系这样的自引力系统的动能和引力势能联系起来，决定其稳定性。但是，当恒星内部的物质不是简单的气体时，会发生什么？
考虑一个旋转的中子星，宇宙中最极端的物体之一。其内部可能是一种中子超流体。要使这种流体旋转，它必须形成一个密集的、规则的量子化涡旋线阵列，这些微小的量子力学漩涡都与恒星的旋转轴对齐。每条涡旋线都处于张力之下，就像一根被拉伸的橡皮筋。这个充满张力的涡旋线阵列在恒星内部产生了一种内部的、各向异性的应力——它在“赤道”上向外推，但在“两极”沿线向内拉。这种可以被表述为一种维里应力的应力，对恒星的整体能量平衡做出了贡献。要理解中子星的结构和稳定性，必须考虑到这种量子力学应力。 在这里，我们看到了一个令人惊叹的尺度汇合：超流性的量子力学产生了微观应力，而这种应力通过维里定理，对质量是地球数百万倍的天体宏观结构产生了影响。
从晶体的刚度，到液体的粘度，再到恒星的稳定性，维里应力张量已被证明是一个具有非凡力量和广泛适用性的概念。它证明了在物理学中，一个精心阐述的思想可以照亮惊人多样的现象，揭示自然世界深刻而美丽的统一性。