lammps应力应变曲线与实验的差值巨大问题

saberve

各位老师好
我在使用lammps跑合金的MD单轴拉伸，绘制应力应变曲线。将lammps中pressure的单位换算为GPa后, 屈服强度在0.9GPa左右。但是一些实验上的文章的结果普遍在100~200MPa。查阅了不同体系，如Fe,Cu等文章，发现分子动力学的文章一般结果都是在几个GPa,但是实验上的数值都只有几百MPa，差距非常大。
如果计算晶格常数，聚合能，弹性常数，弹性模量等，与实验上吻合的很好。

想请问下各位老师，为什么会有如此大的差别，是lammps的pressure定义的问题么?

chenzhe

各种原因都有，如力场精度问题、计算参数问题。更可能是宏观实验材料不纯、存在缺陷，这对宏观材料是几乎不可能避免的，但是MD模拟几乎是无缺陷的晶体，实际上就得拉断金属键，这可比沿着缺陷拉，外力要大太多了。