计算化学公社

标题: 计算油水两相界面张力的数据处理求助 [打印本页]

作者
Author: perfectasdzz 时间: 2025-9-10 11:11
标题: 计算油水两相界面张力的数据处理求助
各位大佬们，我想计算油水两相界面张力，构建了一个油-水-油的模型，界面的面积为40埃*40埃，计算界面张力的相关代码为：
compute stress all stress/atom NULLcompute cc1 all chunk/atom bin/1d z lower 0.3 units box
fix ave all ave/chunk 1000 200 200000 cc1 c_stress[1] c_stress[2] c_stress[3] norm sample file stress.profile

目前得到的数据主要有以下几项Chunk Coord1 Ncount c_stress[1] c_stress[2] c_stress[3]，请问我该如何按照文献中的公式进行计算呢，我是把数据中xx、yy、zz方向单原子的应力张量乘上原子数量Ncount再除以区块体积（40*40*0.3），得到Pxx、Pyy、Pzz，然后再带入文献中的公式进行计算，但结果偏差较大，请问正确的数据计算处理方法是什么呢？

作者
Author: baoluyao 时间: 2025-9-12 21:51
你用的公式和你的命令不对应。你的公式是用体系的总体压强进行计算，你的in文件中的命令使用压强分布进行计算。你得到的结果需要用压强分布的积分进行计算，你取找一下公式就懂了。

作者
Author: perfectasdzz 时间: 2025-9-15 09:41

baoluyao 发表于 2025-9-12 21:51
你用的公式和你的命令不对应。你的公式是用体系的总体压强进行计算，你的in文件中的命令使用压强分布进行计 ...

谢谢老师，我是计算了每个区块的应力张量然后积分计算的，具体的计算过程和公式确实有些差别，主要的问题是在lammps生成的原始数据处理上，因为compute stress/atom 计算的是单个原子的应力，需要变化到整个区块上的应力分量，就是这一步的计算，我没弄清楚。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)