| 查看: 8709 | 回复: 20 | ||||
neweroica木虫 (著名写手)
|
[交流]
【求助】如何固定某几个原子不动？ 已有3人参与
|
|||
|
刚接触LAMMPS，有一个问题想请教大家： 我想在整个MD模拟过程中，固定某几个原子始终不动，请问该如何实现？ 具体用什么命令和关键词？ 谢谢！！！ |
回复此楼» 收录本帖的淘帖专辑推荐
 分子模拟 |
» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请帮忙看一下lammps中是否有这个势函数？
已经有3人回复
- lammps中原子固定
已经有4人回复
- 求助 为什么我安装的MS中画不上原子或者是分子环呢
已经有3人回复
- 求助如何固定原子坐标进行优化，快搞疯了
已经有16人回复
- 求助修改说明的写法
已经有4人回复
- 计算过程中出现警告
已经有5人回复
- 如何 用高斯 固定部分原子 进行优化
已经有4人回复
- 请问castep计算的时候如何来固定部分原子
已经有3人回复
- vasp如何优化指定的几个原子
已经有9人回复
- Help: dn/dc value for poly(2,3-dimethyl-1,3-butadiene)
已经有5人回复
- 【求助】怎么强行将4个原子固定在一个面上？
已经有6人回复
- 【求助】请高人指点，如何人为固定原子占有率
已经有8人回复
- 【讨论】局部优化，固定不了原子
已经有16人回复
- 【求助】结构优化时如何将原子位置固定在原胞内
已经有9人回复
- 【求助】表面吸附原子的位置问题！
已经有9人回复
- 【求助】怎样固定原子和距离
已经有14人回复
- 【求助】计算表面能要不要固定中间的原子
已经有12人回复
- 【求助】请问这两个A问题怎么样解决？
已经有11人回复
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+5,VIP+0):谢谢 12-15 09:19
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+5,VIP+0):谢谢 12-15 09:19
|
这个有多种方法可以实现。 首先你把这几个原子放在一个group里面（比如group1），把其它原子放在另外的group（比如group2)里面。然后，你至少可以选择如下几种方法： 1. 运动方程不对这几个原子积分, 只对group2积分（但我不能确定这个时候group1里面的原子对其它原子的力的作用是否还存在（你可以作一个测试，对于一个构型手工计算一下这几个原子旁边的某一原子的势能，和程序输出的该原子的势能比较一下——用computer pe/atom命令） 2. 运动方程对所有原子积分，但是对group1里面的原子，每一步都标定它的速度，让它们的速度（或温度）为零。 3. 运动方程对所有原子积分，令group1里面的原子的初速度为零，然后用 fix setforce命令，把它们的受力搞成零。 [ Last edited by 老虎大王 on 2009-12-15 at 09:03 ] |
2楼2009-12-15 09:01:58
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
3楼2009-12-15 09:12:48
neweroica
木虫 (著名写手)
- 模拟EPI: 1
- 应助: 53 (初中生)
- 金币: 5027.5
- 散金: 301
- 红花: 31
- 帖子: 2489
- 在线: 217.1小时
- 虫号: 36775
- 注册: 2004-02-11
- 专业: 理论和计算化学
|
非常感谢老虎大王！！！ 1. setforce好像不管用。在run之后，原子位置还是随着timestep发生变化。 2. run 100000 every 1 set group s vx 0 vy 0 vz 0管用。 但我还有个问题： 因为我将每个分子都是为一个刚体，但同时又想对系统施加NVT，所以我用了两个fix来实现： fix 1 all nvt 300.0 300.0 100.0 fix 2 all rigid group 2 mol1 mol2 但是运行后却出现警告： WARNING: One or more atoms are time integrated more than once 请问这个问题该怎么解决？ |
4楼2009-12-15 19:31:25
neweroica
木虫 (著名写手)
- 模拟EPI: 1
- 应助: 53 (初中生)
- 金币: 5027.5
- 散金: 301
- 红花: 31
- 帖子: 2489
- 在线: 217.1小时
- 虫号: 36775
- 注册: 2004-02-11
- 专业: 理论和计算化学
5楼2009-12-15 19:32:19
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+6,VIP+0):谢谢 12-16 14:48
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+6,VIP+0):谢谢 12-16 14:48
|
我没有研究过分子系统，只研究熔盐一类的离子体系和金属一类的原子系统，呵呵。现在我看你的问题是两个： A，所有的分子都是rigid body。我其实不大明白你为什么要把分子都搞成刚体，其实分子内可以有相互作用的啊。这在Lammps中也可以实现。 B，有一些分子始终不动。 呵呵。这两个结合起来，有一点点麻烦。但在Lammps中应该也可以实现。 1. 用setforce, 一定要把初始速度先搞成零。这个命令我也没有用过，我现在不能确定它是只管当前步，还是管整个Run, 你测试一下看。刚才我又看了一下手册，手册上说这个命令可以用来固定原子。 2， Fix rigid 这个命令隐含了对所涉及的刚体中原子进行运动方程的恒能积分，不能再对它用其它的积分命令了，你仔细看说明书就知道了。这个命令的下面有两个“Important note", 你一读就知道。 3， 你如果所有的原子都是rigid body，如上所说，你设置fix rigid命令的时候已经同时积分了方程，但你又想用Nvt，那么这个时候你就只能舍弃HOOVER热浴，而改用Berendsen方法（加个 fix temp/berendsen命令，这个命令只标定速度，不积分方程），或者加 fix temp/rescale命令。 这两个命令不但可以控温，而且可以像我前面所说的一样，把你不想让它动的那些原子的温度（速度）设成零。 4，题外话： 你要用XMD，也有很多办法，比如你把你的双体势做成EAM势函数表，但只有双体势部分有数值，而把Dens和Embed部分都搞成零不就行了，主要是Embed搞成零。这样做，只是没有用Ewald求和，但你是列表形式，算起来也很快的，而且只要你把表中数据点搞得密一些，计算的精度也可以保证。你可以这样搞一下，作个常规计算，与Lammps比较一下。比如用同样的构型，只计算一步，看看能量和力是否相同（主要是能量要相同，力的计算可能由于MD算法不同而略有不同）。 [ Last edited by 老虎大王 on 2009-12-16 at 13:53 ] |
6楼2009-12-16 13:50:29
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
7楼2009-12-16 17:55:57
neweroica
木虫 (著名写手)
- 模拟EPI: 1
- 应助: 53 (初中生)
- 金币: 5027.5
- 散金: 301
- 红花: 31
- 帖子: 2489
- 在线: 217.1小时
- 虫号: 36775
- 注册: 2004-02-11
- 专业: 理论和计算化学
8楼2009-12-18 03:29:07
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
9楼2009-12-18 15:48:26
neweroica
木虫 (著名写手)
- 模拟EPI: 1
- 应助: 53 (初中生)
- 金币: 5027.5
- 散金: 301
- 红花: 31
- 帖子: 2489
- 在线: 217.1小时
- 虫号: 36775
- 注册: 2004-02-11
- 专业: 理论和计算化学
10楼2009-12-18 19:35:29
