[论文解读] Partition sum of thermal, under-constrained systems

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

Elasticity and Material Modeling被引用 1

一句话总结

该论文为接近无热刚性转变的热力学、欠约束系统构建了首个基于原理的配分函数分析框架。通过将弹性性质（如各向同性张力 $ t $ 和剪切模量 $ G $）以应变 $ \varepsilon $、剪切 $ \gamma $ 和温度 $ T $ 表示，提出了三个通用参数：熵弹性刚度 $ \kappa_S $、能量弹性刚度 $ \kappa_E $ 和耦合参数 $ b_\varepsilon $，统一了聚合物网络、膜结构和顶点模型的力学行为。

ABSTRACT

Athermal (i.e. zero-temperature) under-constrained systems are typically floppy, but they can be rigidified by the application of external strain. Following our recently developed analytical theory for the athermal limit, here and in the companion paper, we extend this theory to under-constrained systems at finite temperatures. Close to the athermal transition point, we derive from first principles the partition sum for a broad class of under-constrained systems, from which we obtain analytic expressions for elastic material properties such as isotropic tension $t$ and shear modulus $G$ in terms of isotropic strain $\varepsilon$, shear strain $\gamma$, and temperature $T$. 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. We provide analytical expressions for these parameters based on the microscopic structure of the system. Our work unifies the physics of systems as diverse as polymer fibers & networks, membranes, and vertex models for biological tissues.

研究动机与目标

• 将欠约束系统的无热理论扩展至有限温度。
• 从第一原理出发，推导热力学、欠约束系统在刚性转变附近的配分函数。
• 统一聚合物网络、膜结构和顶点模型等不同系统在热涨落下的力学行为。
• 识别并解析表达控制弹性响应的三个通用参数（$ \kappa_S $、$ \kappa_E $、$ b_\varepsilon $）。
• 提供各向同性张力 $ t $ 和剪切模量 $ G $ 关于 $ \varepsilon $、$ \gamma $ 和 $ T $ 的解析表达式。

提出的方法

• 利用统计力学对具有 $ N_{\text{dof}} > N_s $ 的广义弹簧网络在周期性边界条件下进行配分函数的形式推导。
• 应用刚性弹簧极限与有限刚度处理方法，以模拟热涨落和应变诱导的刚性。
• 通过对弹簧长度函数 $ L_i(\vec{R}, V, \gamma) $ 引入解析性与齐次性条件，确保在转变点附近的数学可处理性。
• 通过线性代数识别零模态与自应力状态，高阶刚性（K阶模态）采用微扰方法处理。
• 在无热转变点附近推导配分函数，此时系统处于临界刚性状态。
• 利用配分函数及其导数，解析计算弹性模量的热力学平均值。

实验结果

研究问题

• RQ1 如何从第一原理出发，推导热力学、欠约束系统在无热刚性转变附近的配分函数？
• RQ2 各向同性张力 $ t $ 和剪切模量 $ G $ 关于应变 $ \varepsilon $、剪切 $ \gamma $ 和温度 $ T $ 的解析表达式是什么？
• RQ3 三个通用参数——熵弹性刚度 $ \kappa_S $、能量弹性刚度 $ \kappa_E $ 和应变耦合 $ b_\varepsilon $——如何从微观网络结构中涌现？
• RQ4 在何种程度上，聚合物网络、膜结构和顶点模型等不同系统的力学行为可以被统一于单一热力学框架之下？
• RQ5 与无热极限相比，热涨落的引入如何改变刚性转变的行为？

主要发现

• 通过第一原理，成功推导出热力学、欠约束系统在无热刚性转变点附近的配分函数。
• 各向同性张力 $ t $ 和剪切模量 $ G $ 被解析表达为各向同性应变 $ \varepsilon $、剪切应变 $ \gamma $ 和温度 $ T $ 的函数，除 $ \kappa_S $、$ \kappa_E $ 和 $ b_\varepsilon $ 外无其他自由参数。
• 三个关键参数 $ \kappa_S $、$ \kappa_E $ 和 $ b_\varepsilon $ 均源自微观网络结构，从而可实现弹性响应的通用预测。
• 该理论统一了聚合物纤维/网络、膜结构以及生物组织的顶点模型在热涨落下的力学行为，提供统一框架。
• 刚性弹簧极限与有限刚度极限所得配分函数并不相同，证实有限弹簧刚度对物理一致性至关重要。
• 该框架能够捕捉由外加应变诱导的预应力与自应力状态所引发的刚性出现，即使在有限温度下亦成立。