| 查看: 5631 | 回复: 44 | ||||
| 【奖励】 本帖被评价4次，作者043114076增加金币 3.2 个 | ||||
[资源]
关于towhee 回答zhang_jaj的问题
|
||||
|
以下是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知识，导师不管，自己一个人摸索。担心毕业十大问题啊，谢谢了 大家基本都差不多, 坚定信心, 可以做的很好. 其实模拟计算这个方向比起做实验还是比较透明的, 重复文献一般没有问题, 技术问题一般论坛上都会回答. 我也是在学习的过程, 大家多多讨论, 相互学习. [ Last edited by 043114076 on 2012-4-26 at 09:55 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 towhee精华 | towhee的运用讨论！！ |
» 本帖已获得的红花（最新10朵）
» 猜你喜欢
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
2025冷门绝学什么时候出结果
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
关于分子模块求助的一些建议
已经有29人回复- towhee软件的安装问题
已经有4人回复
- Towhee能模拟气体吸附等温线吗？
已经有10人回复
- towhee编译失败
已经有9人回复
- 运行 TOWHEE 时出现的问题 ，希望有人帮我解决一下，非常感谢
已经有8人回复
- 使用amber力场，一些参数无法获得
已经有5人回复
- 想做粒子和固体（金属）相互作用，不知道lammps能否实现？
已经有6人回复
- 求问分子模拟断检问题。。。
已经有5人回复
- LAMMPS并行安装问题
已经有4人回复
- MS 进行溶液分子动力学计算时在哪里设置力场参数
已经有5人回复
- 关于气体液化的问题~~~~~
已经有13人回复
- 求用到AMBER force field的开源软件
已经有10人回复
- 水分子在外场下键角、键长的变化？
已经有4人回复
- ms模型导入lammps
已经有8人回复
- lammps里的随机种子是如何确定的？
已经有12人回复
- lammps里log文件里一些数据的单位？
已经有10人回复
- 使用Towhee做NVT模拟时，气液两盒子的hmatrix 大小该怎么设置？ 有什么原则或公式么？
已经有7人回复
- 关于气液平衡的简单问题
已经有24人回复
- 石墨烯 lammps
已经有11人回复
- 石墨烯和下层的SiO2会成键吗？
已经有7人回复
- 使用towhee进行GEMC模拟相平衡的MC运动定义问题..
已经有15人回复
- 【讨论】关于爆炸盐/浴盐/沐浴球的配方
已经有4人回复
- 【求助】弱问有关AMBER远程调用的问题？
已经有3人回复
- 【求助】分子的LJ势能参数的拟合
已经有8人回复
2楼2012-04-26 19:15:47
|
你好，我是按照你之前一个帖子直接编译的并行towhee-7版的，但是这个里面没有那个towheebase的目录（应该是6版以前才有的把）。在做练习的时候，我把towhee_input文件的forcefield改成了自己的绝对路径 /public/software/towhee/Forcefield/() 然后towhee_parallel如下 1 towhee.log /public/software/towhee/Examples/Canonical/Amber_Villin 之后执行了mpirun -np 4 towhee 显示出来的为 Number ofjobs:0 Simulations complete Total time:0.000059 seconds 而且没有towhee.log文件 还有官网说安装需要的autoconf与automake文件是不是可以不用？ example中的makefile文件用来编译什么的？ 非常感谢！（因为一些原因实验室不能上网，只能晚上回来求助） [ Last edited by zhang_jaj on 2012-4-26 at 20:25 ] |
3楼2012-04-26 20:24:20
4楼2012-04-26 21:53:51
|
你好，您知道towhee用intel的编译器好呢还是PGI的好，因为我看官方给的是PGI的。还有我在运行例子的时候，提交大于3个任务并行时（如16），从倒数第三个就会显示出 Job # 14 Insufficient entres in towhee_parallel file 后面就没了，但是前13个job都顺利。不知道是怎么回事。 当小于三时也会有最后一个job显示Insufficient entres in towhee_parallel file 还有想问下，mpich2的并行环境好像不能编译是吗？ 还有想问下，如果在./configrue时，不给CC，F77等赋值，编译好的不会是用GNU的编译器吧？ 麻烦了，不甚感激！ |
5楼2012-05-04 22:00:07
送鲜花一朵|
你好，想再请教下您。关于towhee参数设置问题： 1）.pm*的值是不是按照从小到大进行排列的呢。example中有个是 pmtraat 1.0 pmtracm 0.67 pmrotate 1.0 这种情况，当要判断是进行哪种move时，是不是只进行pmtraat呢？ 2）.给体系加上一个LJ-9-3 wall时，里面有一个参数是ljfrho（the number density of atoms in the integrateed wall potential）是wall上与体系内分子作用的点密度吗？还有ljfsig与ljfeps都是人为的指定大小吗？ 非常感谢！ |
6楼2012-05-24 20:14:55
7楼2012-05-24 22:56:08
|
你好，还想向您请教个问题关于hmatrix参数的。mannul上说这个坐标矩阵可以不正交，但是必须是right hand rule的。我在想如果用其他非正交的坐标矩阵时，程序在运行过程中是不是速度就低了很多？？？（因为个人认为正交坐标里各个原子间进行距离计算时要方便的多（如我在一个面上固定一个单层石墨烯）） 还有个问题是，对于NVT系宗，如果我先让N=1000，T=50K后，该怎么选取box的大小呢（气体），这时不同的box大小应该就对应不同的P吧。是利用PV=NRT吗？这个box的大小是依据我想要的P用那个公式来调节吗？，还是normal atmosphere下算个大约的box大小？ 非常感谢！！！  |
9楼2012-06-02 21:41:54
10楼2012-06-03 13:33:26
11楼2012-06-03 22:48:22
12楼2012-06-05 22:28:48
|
这个看了还是不太明白。按照mannul应该是把wall当做一种一般性的surface，而这个ljfrho就是这个面上的点分布密度。还有下面的那个sigma和epislon都是这个wall与体系中粒子的作用参数。（这个参数该怎么定义呢），还有如果把这个wall当做Si之类的来计算ljfrho，那么这些Si原子的排列该怎么确定(或者说这是一个无序的surface吗)？你之前给的那个文章链接我一直下不了，不知道你有吗，可以给我一份吗（zhjjaj@163.com）？关于LJ9-3的文章非常感谢！http://jcp.aip.org/resource/1/jc ... _s1?isAuthorized=no |
13楼2012-06-06 23:01:36
|
我从头讲好了。 我们一般模拟时使用的是全粒子势能比如UFF，amber等等，这种势能每个原子作为独立的粒子，有自己的力场参数。这样在模拟时需要计算每一对粒子的相互作用能。如果我们的体系是气体或者液体分子和一个表面的相互作用，当然也可以使用全粒子势能参数去模拟，但是我们也可以把表面所有粒子的力场参数积分使其变成一个整体。9-3 wall 势能是从 12-6 粒子势能进行空间积分得到的。积分推导的过程见：http://www.sklogwiki.org/SklogWi ... ard-Jones_potential , 建议自己可以简单推导一下。这样做的好处是大大减少了粒子的数目（表面当成了一个具有特定形状的粒子），从而减少了计算量。推导过程中，原来12-6中的参数sigma ，epsilon意义没有变化，直接继承下来，另外新加入了一个参数ljfrhho，和体系的组成相关，即体相单位体积内的粒子数目。 我只是以Si和SiO2 为例，你要根据你的体系去得到sigma epsilon和ljfrhho。请问你的体系中wall的组成是什么物质？ 费了这么大劲，不知讲清楚了没有。 那个文献我也没有，你可以求助一下 |
14楼2012-06-07 09:13:50
|
感觉你讲了后，确实对这个势有了很多了解（之前特别迷糊）。我想用单层graphene和双层graphene试着模拟Ar的nucleation过程。这个2D的如单层的a为2.46（Angstrom），那我算的时候是不是应该就是2/（2.46*[sqrt（3）/2]*2.46）=0.38.... 再比如一个fcc的（111）面，加入他的a=2，那么计算出来的是不是应该就是4/（2^3）=0.5？这样的话就跟选fcc的哪个切面无关了呀（或者我还是理解错了  ）？还有我以为加上一个LJ9-3wall（如graphene）就是一个虚拟的wall；但如果我再在其他位置加一个同样组成的如graphene-wall时，并且给出这个graphene中各原子的坐标，而且这个wall不是LJ9-3的势。或者说在towhee中模拟的只要是加了一个用来吸附的表面，那么它就应该是当做LJ9-3wall来处理吗。这样他的原子排列还是不好确定（如模拟外延生长时要选不同的面）？ 所以加入一个自己引入的面是不是应该和加一个LJ9-3是完全两个概念吧？ 还有，以你的经验，你觉得像graphene有近1000作用的碳原子，这个得当做一个独立的大分子来建立初始文件，你觉得应该用哪种方式建呢，这样就得在input_style下面有1000个作用的vibration手动参数输入，感觉太恐怖了。 最近在这几个问题上实在搞不定，就只能想到麻烦你了 |
15楼2012-06-07 12:44:42
16楼2012-06-07 13:42:50
17楼2012-06-07 17:01:28
18楼2012-06-07 17:07:53
|
非常感谢啊，我今天下午去把这个大分子上的每个原子当做独立的分子，程序成功运行了。我想了下，其实每个c原子在力场都是一个独立体对待的，所以固定不动时把它当做很多歌分子也是可行的（下午给我乐的，呵呵）。 其实这个LJ-9-3我发现它就是一个弥散的surface吧，我那段话说的意思主要是为了确定这个surface表面的原子排布不好确定，现在我想通了，它应该没有这个原子排列的概念吧。 希望你能给我举个例子，比如我分别用bcc与fcc的铁的来做这个LJ-9-3 wall，加入他们的晶格常数都是2（Angtrosm），那么这个ljfrho分别是怎么求的，值时多少？ 而ljfrho下面那个ljgsig与ljfeps的值，是不是需要重新建个模型用Fe的势与体系的势进行一个cycle，然后从towhee-out中看Fe与体系中各原子的sigma和epsilon值？ 还有towhee mannul中关于lj-9-3介绍下面那个公式里的sigma与epsilon应该就是ljfntypes下面的一系列ljfsig与ljfeps值吧，这样的话，如果体系有n种与surface作用的原子，那么在积分过程中就应该有n中lj-9-3的积分式（参数不同）吧 关于那个pdb2towhee的我也搞定了，用":1, $s/ATOM /HETATM /g"，然后就可以了。 |
19楼2012-06-07 18:39:00
20楼2012-06-07 22:15:20
21楼2012-06-07 22:30:14
|
关于涉及到石墨烯的建模，因为石墨烯是固定的，所以我按照你说的那种方法把每个原子当做一个独立的分子来建模。这样存在着问题： 1.由于突然分子数多了所以对于inix， iniy， iniz 那三个值就需要增加，但是不知道这三个值到底是什么意思，把盒子划分开有意义吗（因为原子坐标都是从towhee_coords中读取的）？ 2.当我用c60试验的时候，这种方法一切正常，但是当我用涉及到1000多碳原子的石墨烯用这方法时，运行什么的都正常，但是在输出文件中有问题：信息只输出到关于各原子间势那块就没了，也就是显示了LJ公式，然后是sigma epsilon这一列，下面没东西了。。。不知道怎么回事，但是程序确实还在运行中，而且还不停有box_*.pdb文件产生，好生纳闷啊 希望你能帮我排除下故障，谢谢啊 |
22楼2012-06-08 12:19:11
23楼2012-06-08 18:46:49
★★★★★ 五星级,优秀推荐
|
楼主你好，我也是刚刚接触到towhee软件的，我是做三元萃取体系气液平衡的，现在需要模拟纯组分噻吩(thiophene)的气液平衡情况，在模拟之前看了一下要用的力场文件，发现里面缺少一部分我所需要的键长键角文件，请问那个力场文件是不是没有显示全部可以处理的键长键角数据呢？（因为之前读过一篇文献说towhee钟的Trappe-EH力场可以处理N原子的问题，但是我进入到towhee 7.0.2版本的那个力场文件后发现里面根本没有N原子的数据。所以才有这种疑问##）。 如果真的没有我该自己添加力场的数据么？需要的话是简单地从别的联合原子力场中找到N原子数据，然后粘贴我现在用的这个力场中吗？ 小弟初学，麻烦你了！ |
24楼2012-06-11 17:58:17
25楼2012-06-12 09:09:15
送鲜花一朵|
楼主师兄你好，我又有几个问题想问一下你。 1.两个力场的混合规则不同，可不可以用什么方法给其中一个导入一种混合规则呢？如果可以怎么做？ 2.对于手动输入电荷时，如果关于力场的描述文件中的Coulombic interactions项下有电荷值，我是不是就该输入此电荷值，如果Coulombic interactions没有值，我是不是就不输入电荷值呢，即为0.0d0？入Dubb2004势，用来建立SiO2分子时？ 3,在建立一个分子时，如果有多个units的atom组成，那么对于towhee_coords文件里面是不是应该按照input_style下的unit值对应的把一种atom输入完再输入另一种呢？ 4，您知道这个程序的结果文件pdb或者movie文件里怎么判断材料是fcc还是bcc之类的吗，或者您知道那个软件有这个功能吗？ 5.（补充的）对于不同的力场文件，如argon的epsilon值都不相同，是不是意味着argon在固体情况下，对于不同的力场他的晶格常数不同呢？ 6，对于涉及同位素的，就拿H和D来说吧，如果我又H的力场文件，可不可以把H的mass修改一下当做D的立场文件来用呢？（对于LJ势） 非常感谢！ [ Last edited by zhang_jaj on 2012-7-3 at 22:01 ] |
26楼2012-07-03 21:41:04
27楼2012-07-04 10:45:34
|
你好： 1.我不知道混合规则不同的能不能用，所以就没用，下午可以去试一试 2.也就是说，只要我的体系不想考了库伦作用，我就可以都不写电荷值吧 3.问题解决，因为当时我导出来的sio2的pdb文件里si和o是间隔的写的，对于towhee需要把si的坐标写完，再写o的坐标。 4.因为我做的是气体低温结晶过程，所以需要判断到底是冻结成了什么结构，这个很郁闷啊 6.关于这个我知道对于vasp是可以这么修改的，但是对于LJ势的力场文件，因为不知道参数和原子的质量有没有关系，所以不知道这么修改得到的同位素力场可不可以用。 问题：还想问一下，不知道您对形核过程了解吗，我不知道该怎么在模拟过程中计算体系的形核势垒。如果说是用成核时的体系总能量与初始的总能量相减的话，那么这个势垒是不是特别依赖与初始模型的原子排布（因为我是看您的教程用packmol生成的文件，所以不知道这个生成的排列对于计算势垒时是否具有一般性） towhee里面有umbrella sampling和flux forwar sampling两种抽样吗？如果没有想要加的话是在编译后的程序里加呢，还是加了之后重新编译呢？ 非常感谢！！！！ |
28楼2012-07-04 12:41:52
29楼2012-07-04 19:36:01
30楼2012-07-04 20:36:21
31楼2012-07-05 10:52:47
33楼2012-07-12 12:24:16
34楼2012-07-12 18:50:44
35楼2012-07-12 23:07:02
36楼2012-07-13 09:28:27
37楼2012-07-13 13:55:05
38楼2012-07-13 15:34:09
39楼2012-07-13 15:38:04
8楼2012-05-25 06:55:15
32楼2012-07-05 12:01:35
40楼2012-09-16 10:35:10
41楼2012-09-17 15:22:48
42楼2012-09-18 12:41:22
43楼2012-11-29 20:18:22
44楼2013-10-19 21:45:54
45楼2014-11-10 15:04:53