[论文解读] Nonlinear Elasticity from Atomistic Mechanics

Christoph Ortner, Florian Theil|arXiv (Cornell University)|Feb 17, 2012

Advanced Mathematical Modeling in Engineering参考文献 21被引用 7

一句话总结

本文通过分析静态与动态问题的原子-连续耦合，严格建立了非线性弹性中柯西-伯恩近似的一致收敛界。在晶格稳定性和快速衰减相互作用假设下，证明了在小数据条件下，原子解以 $\varepsilon^2$ 阶误差收敛于柯西-伯恩解，涵盖位移梯度与能量，方法基于应力的逐点一致性与一种局域化技术。

ABSTRACT

We present sharp convergence results for the Cauchy--Born approximation of general classical atomistic interactions, for both static and dynamic problems, for small data.

研究动机与目标

• 严格证明一般多体原子相互作用在非线性弹性中的柯西-伯恩近似。
• 通过将收敛结果从静态问题扩展到动态问题（演化方程），弥合理论理解的空白。
• 在最小物理假设下，为弹性静力学与弹性动力学建立精确的误差估计。
• 通过反例证明晶格稳定性是收敛的必要条件，即不稳定性会导致解无界增长。

提出的方法

• 使用在布拉维晶格 $\varepsilon\Lambda$ 上缩放的多体势 $V$ 形式化原子能 $E^\varepsilon$，其中 $\Lambda = \mathbb{Z}^d$。
• 定义连续柯西-伯恩能量密度 $W(\mathbf{F}) = V(\{\mathbf{F}\rho\}_{\rho \in \Lambda \setminus \{0\}})$，建立离散与连续运动学之间的联系。
• 采用局域化技术将离散对象连续表示，实现弱形式化并进行应力比较。
• 建立原子应力与柯西-伯恩模型中皮奥拉-基尔霍夫应力之间的逐点二阶一致性。
• 应用能量估计与谱分析，以 $L^2$ 与 $H^1$ 范数界定原子解与柯西-伯恩解之间的差异。
• 利用近似特征函数证明晶格不稳定性会导致指数增长，从而证明稳定性假设的必要性。

实验结果

研究问题

• RQ1 在何种条件下，柯西-伯恩近似会收敛于静态非线性弹性问题的原子解？
• RQ2 柯西-伯恩模型能否准确近似由牛顿运动方程控制的原子系统的动力学行为？
• RQ3 在小晶格间距极限 $\varepsilon \to 0$ 下，原子解与柯西-伯恩解之间的收敛速率精确为何？
• RQ4 晶格稳定性是否为收敛的必要条件，还是可在更弱假设下实现收敛？
• RQ5 多体势的结构（如成对、三体、四体项）如何影响柯西-伯恩近似的相容性与收敛性？

主要发现

• 对于静态问题，在小外力与初始数据条件下，原子解与柯西-伯恩极小化器之间梯度的差异被有界于 $C\varepsilon^2$。
• 对于动态问题，在与 $\varepsilon$ 无关的时间区间 $[0,T]$ 内，梯度的 $L^2$-范数解差异被有界于 $C\varepsilon^2$。
• 收敛速率 $\varepsilon^2$ 是最优的，其来源为皮奥拉-基尔霍夫应力张量的逐点二阶一致性。
• 晶格稳定性不仅充分而且必要；若无稳定性，原子系统将存在指数增长的解。
• 通过不稳定晶格（如 $a_1 = -1, a_2 = 1/2$）的反例表明，原子能可能无下界，导致局部极小化器失效，并在动力学中引发不稳定性。
• 该分析涵盖一般多体势，包括具有二面角与扭转项的势能，只要其满足置换不变性以及稳定性和衰减条件。