[论文解读] Generic elasticity of thermal, under-constrained systems

Cheng-Tai Lee, Matthias Merkel|arXiv (Cornell University)|Apr 14, 2023

Elasticity and Material Modeling被引用 2

一句话总结

本文通过统一熵弹性与能量弹性，为热力学、欠约束系统（如聚合物网络和顶点模型）的弹性性质建立了一套通用的解析理论。该理论仅使用三个参数——熵刚度（κS）、能量刚度（κE）和应变耦合参数（bε）——推导出各向同性张力与剪切模量的温度、应变和微观结构依赖表达式，揭示了这些贡献项的行为如同串联弹簧，从而解释了在零应变下弹性性质呈现 T^{1/2} 依赖关系的观测现象。

ABSTRACT

Athermal (i.e. zero-temperature) under-constrained systems are typically floppy, but they can be rigidified by the application of external strain, which is theoretically well understood. Here and in the companion paper, we extend this theory to finite temperatures for a very broad class of under-constrained systems. In the vicinity of the athermal transition point, we derive from first principles expressions for elastic properties such as isotropic tension $t$ and shear modulus $G$ on temperature $T$, isotropic strain $\varepsilon$, and shear strain $\gamma$, which we confirm numerically. These expressions contain only three parameters, entropic rigidity $\kappa_S$, energetic rigidity $\kappa_E$, and a parameter $b_\varepsilon$ describing the interaction between isotropic and shear strain, which can be determined from the microstructure of the system. Our results imply that in under-constrained systems, entropic and energetic rigidity interact like two springs in series. This also allows for a simple explanation of the previously numerically observed scaling relation $t\sim G\sim T^{1/2}$ at $\varepsilon=\gamma=0$. Our work unifies the physics of systems as diverse as polymer fibers & networks, membranes, and vertex models for biological tissues.

研究动机与目标

• 开发一种适用于热力学、欠约束系统（超越无温极限）的通用解析框架，以描述其弹性性质。
• 统一系统在无温转变点附近熵与能量对刚度的贡献。
• 解释在零应变下通过数值模拟观测到的各向同性张力与剪切模量呈现 T^{1/2} 依赖关系的成因。
• 证明弹性性质仅依赖于三个由微观结构导出的参数：κS、κE 和 bε。
• 通过在不同应变与温度条件下对随机切割的三角形网络进行数值模拟，验证该理论。

提出的方法

• 在弹簧无限刚硬的极限下，通过相空间体积分析推导出仅由熵弹性主导的弹性性质。
• 通过熵弹性与能量弹性贡献的类串联叠加方式，将上述结果与先前的无温结果（仅能量弹性起作用）相结合。
• 引入三个关键参数：熵刚度（κS）、能量刚度（κE）和应变耦合参数（bε），这些参数均可从网络的微观结构中提取。
• 利用统计力学从第一性原理出发，推导出各向同性张力 t(ε, γ, T) 与剪切模量 G(ε, γ, T) 的解析表达式。
• 通过使用 Metropolis–Hastings 采样和应力张量平均的二维三角形弹簧网络蒙特卡洛模拟进行验证。
• 将参数拟合至模拟数据，并通过相容性矩阵的奇异值分解与特征值分析，独立地从网络结构中预测这些参数。

实验结果

研究问题

• RQ1 在热力学、欠约束系统中，各向同性与剪切弹性模量如何随温度 T、各向同性应变 ε 和剪切应变 γ 变化？
• RQ2 在 ε = γ = 0 时，数值模拟中观测到的弹性性质呈现 T^{1/2} 依赖关系的物理起源是什么？
• RQ3 如何在统一框架中系统地结合熵与能量对刚度的贡献？
• RQ4 是否可以使用单一通用理论描述多种系统（如聚合物网络、薄膜、顶点模型）的弹性响应？
• RQ5 这三个关键参数（κS、κE、bε）在多大程度上可直接从系统的微观结构中预测？

主要发现

• 在 ε = γ = 0 时，各向同性张力 t 与剪切模量 G 均呈现 t ∼ G ∼ T^{1/2} 的依赖关系，其成因在于熵刚度与能量刚度的协同作用。
• 从模拟数据中提取的三个参数：κE ≈ 0.161，κS/kB ≈ 0.0629，bε ≈ 0.835，可对 ε、γ 和 T 跨越多个数量级的条件下，定量预测 t 与 G。
• 通过网络结构（基于相容性矩阵的奇异值分解与特征值计数）推导出的理论 κE 与 κS/kB 值，与拟合值的偏差分别仅为 3% 和 0.1%。
• 描述各向同性与剪切应变之间耦合的参数 bε，其预测值与网络几何结构的吻合度达 2%，证实了其物理相关性。
• 该理论解释了熵刚度与能量刚度的行为如同串联弹簧，为多种系统提供了统一的物理解释。
• 该模型在广泛条件范围内成功复现了张力与剪切模量的数值数据，验证了其普适性与预测能力。