blends模块计算相互作用参数（F-H参数）求助！

学！！！

各位老师好，我在计算聚合物间的相容性。考虑采用MS中blends模块计算相互作用参数，遇到了几个问题向大家求助。1.如何考虑A组分的质量分数对相互作用参数的影响？
我现在是根据质量分数，分别建立了B组分不同分子数的模型作为Screen，以A组分一条长链作为base，不知道这样计算是否合理。
2.作为screen的分子模型是否应该先进行结构优化、退火、NPT、NVT等一系列操作，以获得接近真实的性能参数？
还请各位老师帮我解惑，谢谢！

学！！！

“We have verified that alternative estimations of x based in the direct evaluation of the intermolecular energy between a pair of units averaged over a statistical sample generated through Monte Carlo simulations ([10]; as incorporated in the MS module ‘blends’) seem to be considerably less precise in the present mixture, depending largely on small details in the forcefield definitions”
————Amirhossein Ahmadi & Juan J. Freire (2008): Molecular dynamics simulation study of compatibility for the polyvinylmethylether/polystyrene mixture, Molecular Simulation, 34:10-15, 1253-1258

学！！！

naoki 发表于 2021-3-11 12:09
我记得blends模块就是简单的在A旁边放置一定数量的B分子，放的分子数是自己指定的，它怎么放是程序内部完 ...

谢谢老师！
我看了一下文献，也是准备放弃使用blends模块进行计算了
后面怎么计算我再看看文献琢磨琢磨

naoki

学！！！ 发表于 2021-3-11 09:40
谢谢老师！
1.我是看文献里有这么计算的，寻思模仿一下，那就将就着看一下效果吧~
2.因为要考虑一定数 ...

我记得blends模块就是简单的在A旁边放置一定数量的B分子，放的分子数是自己指定的，它怎么放是程序内部完成的，所以MD步骤好像没啥用。可以看看M$的help怎么描述blends的。
如果你要计算聚合物和小分子的相互作用，可以考虑用GROMACS直接做动力学模拟然后分组提取相互作用能。如果你想获得溶解度参数，LAMMPS跑动力学可以获取A或B的内聚能，可由此计算溶解度参数和Flory-Huggins参数。

学！！！

naoki 发表于 2021-3-11 09:24
1.我感觉这个模块不好考虑质量分数的影响，也许可以用不同数量B凑合下。。。
2.结构优化就够了

谢谢老师！
1.我是看文献里有这么计算的，寻思模仿一下，那就将就着看一下效果吧~
2.因为要考虑一定数量的B，所以我考虑可能需要将晶胞NPT压实后再做计算。初始密度建立晶胞时会有比较大的自由空间。但我也不确定blends在计算过程中是不是这个原理。。。

naoki

1.我感觉这个模块不好考虑质量分数的影响，也许可以用不同数量B凑合下。。。
2.结构优化就够了