Towhee-经验分享

以下是zhang_jaj 的问题:
你好，我想用towhee做NVT下argon气在surface面上nucleation过程。首先我想把potential选为LJ，surface与particles也是LJ作用，

1）问题是这个surface是需要用particles来make吗，如果是的话那么构成它的particles该怎么固定

2 ）或者，直接用towhee里面的LJ9-3 wall吗，想知道这种wall是与particles怎么作用的，它不是用小粒子构成的吗？

3）模拟完成后，液体与固体的结构判断该怎么做，有什么程序代码吗？

4）这个nucleation的rate该怎么算（是单位时间单位体积内产生临街晶核的个数），我想知道程序中怎么实现，我没发现有那个参数相关的，或者是需要自己来加东西吗？

5）还有这个nucleation barrier该怎么来实现？

6）towhee不能并行，那么其实一个任务就是在一个process上处理的；如果我的个人电脑比服务器的主频高，那么是不是用台式机比服务器更快呢（因为台式机是4核的，它是不是能自己把towhee的任务来用4核一起处理？）？

7）希望您能给点关于towhee代码学习的经验，就是对于我要做的这个方向的。

这是我的回答:
实在不敢当, 我现在也是学习阶段, 因为没做过类似的体系, 对于你的问题只能有一些建议.

1）问题是这个surface是需要用particles来make吗，如果是的话那么构成它的particles该怎么固定

可以使用particle模型来模拟, 固定比较简单, 在move设置中不让这些粒子移动就行了.

2 ）或者，直接用towhee里面的LJ9-3 wall吗，想知道这种wall是与particles怎么作用的，它不是用小粒子构成的吗？

这个我以前没有用过, 不过文献当中确实有人这样做. 这种势能的使用情况可以参考例子Fris_Walls, 另外上网搜了一下, 这几个帖子和文章可能对你有用.
http://jcp.aip.org/resource/1/jcpsa6/v67/i5/p2384_s1
http://jcp.aip.org/resource/1/jcpsa6/v73/i8/p4050_s1
http://lammps.sandia.gov/threads/msg24356.html
http://www.sklogwiki.org/SklogWi ... ard-Jones_potential
这种势是粒子势能转化而来的, 第四个参考文献有推导过程, 主要的优势是可以大大提高速度. lammps 的帖子跟你的情况有点像, 使用MD做的, 可以借鉴.

3）模拟完成后，液体与固体的结构判断该怎么做，有什么程序代码吗？

这个主要是分析RDF吧, 代码在 utils 中有, 可以根据你的情况进行修改.

4）这个nucleation的rate该怎么算（是单位时间单位体积内产生临街晶核的个数），我想知道程序中怎么实现，我没发现有那个参数相关的，或者是需要自己来加东西吗？

这个是需要后处理才能得到的吧. towhee本身不是计算成核过程的. 应该需要编写代码实现. 不过MC 如何计算时间呢? 如果用MC步数是不是会有问题? 是不是用MD更合适一点?

5）还有这个nucleation barrier该怎么来实现？

这个应该是进行不同温度点的模拟, 作速率-温度关系图, 然后通过拟合热力学公式得到吧.

6）towhee不能并行，那么其实一个任务就是在一个process上处理的；如果我的个人电脑比服务器的主频高，那么是不是用台式机比服务器更快呢（因为台式机是4核的，它是不是能自己把towhee的任务来用4核一起处理？）？

towhee的不能用多个核并行计算一个任务, 而事实上也不需要. 一般来说, 模拟多个任务可以等价于模拟一个大任务, MC的本质决定它是这样的性质, MD可不行. 所以尽量还是使用towhee提供的并行模式. 注意选择不同的随机种子.

7）希望您能给点关于towhee代码学习的经验，就是对于我要做的这个方向的。

关于代码, 其实我也没有太多经验. 我主要是对结果进行后处理. towhee 本身是一个general Monte Carlo code, 不是专门为某些体系设计的, 但是对于你的体系应该是可以胜任的. 主要针对体系的代码, 差别也主要是建模和分析结果两个过程, 核心代码都是类似的. 我个人认为如果时间有限, 不要陷到源代码中去.当然能够多了解一些代码对于软件的使用也是非常有好处的.
towhee 的代码由于需要我简单的读过一些, 还是比较清楚的. 可以尝试从总体上去学习下.
对于结果的后处理的代码, 建议从修改towhee 提供的代码开始. 使用语言以自己最熟悉的为准, 如果都凑合, 推荐 fortran90.
多学习 linux bash 脚本也是很有好处的, 我现在简单的数据处理基本都用它.
最后强调, 不要陷到代码中去, 能够完成自己的用途就可以了, 越简单越直接越好. 多多思考跟你体系相关的问题, 多查文献.

///希望能多多向你学习towhee知识，导师不管，自己一个人摸索。担心毕业十大问题啊，谢谢了

大家基本都差不多, 坚定信心, 可以做的很好. 其实模拟计算这个方向比起做实验还是比较透明的, <