模拟计算 || 界面张力计算脚本（配合Forcite使用）

本期推送我在MS官方售后论坛里面淘来的界面张力计算脚本

本文的脚本，来自于https://bioviacommunity.force.com/在不登陆账号的情况下，可以自由下载查看的内容

计算界面张力的脚本，脚本作者是Jason DeJoannis，就是写了xlink脚本的那位大佬。对xlink交联脚本感兴趣的同学可以点击Forcite交联脚本

该脚本需要用户自定义的数据在脚本中有着明显的提示，这里就不多赘述了。针对脚本使用前的体系该如何模拟，Jason在论坛中给出了自己的回答，这里很有必要跟大家说一下。

这里可以计算液-液界面和液-真空界面的界面张力，使用的公式在前面的注释部分给出了，如果需要引用文献可以引用J Chem Phys B, 2007, 11, 7812; Molecular Simulation, 2007, 33, 27

*如果计算液-真空界面，就像下图一样。一定要使用NVT系综进行计算。为啥不能用NPT？这个大家建议建模如下图然后尝试一样NPT系综模拟，这样你得到的宝贵经验会指导你以后的模拟计算。不想亲自模拟的，可以在本文的角落处找到我给出的问题答案。

*如果是研究液-液界面，那么需要先进行NPT作为平衡过程然后NVT生产过程收集轨迹文件。

无论哪种界面，都要注意，要保持层层分明，一不小心搞成乳液了就没办法算了。

阅读原文：界面张力计算脚本（配合Forcite使用）