多个势函数使用方法
lammps吧
全部回复
仅看楼主
level 3
鞋钉ww 楼主
各位大神请帮忙指点一下,感激不尽!!!
我的系统有两种金属A和B,我的目的是A-A之间使用纯金属势函数A.eam.fs,B-B之间使用纯金属势函数B.eam.fs,A-B之间使用合金势函数A-B.eam.alloy,请问以下命令可以办到么?
pair_style hybrid eam/alloy eam/fs eam/fs
pair_coeff * * eam/alloy A-B.eam.alloy A B
pair_coeff * * eam/fs 1 A.eam.fs A NULL
pair_coeff * * eam/fs 2 B.eam.fs NULL B
使用pair_style hybrid命令:eam/alloy A-B.eam.alloy A B中定义的A-A势会被eam/fs 1 A.eam.fs A NULL中定义的A-A势覆盖掉,这样A-A势就只来自于势函数A.eam.fs,同理B-B也是一样。
使用pair_style hybrid/overlay:eam/alloy A-B.eam.alloy A B中定义的A-A势会和eam/fs 1 A.eam.fs A NULL中定义的A-A势叠加在一起,导致A-A势的重复,同理B-B也是一样。
我的理解正确么?请大家讨论!谢谢!
2019年06月17日 07点06分 1
level 9
悬,你这样的话三个eam里面的势能的配位数的选取可能会有问题
2019年06月19日 11点06分 2
level 1
楼主问题解决了么?
2019年10月23日 01点10分 3
level 3
pair_coeff * * eam/fs 1 A.eam.fs A NULL
pair_coeff * * eam/fs 2 B.eam.fs NULL B
这两行命令 NULL位置不同,分别表示什么意思
2020年03月27日 08点03分 4
pair_style hybrid eam/cd eam/alloy pair_coeff * * FeCr.cdeam Fe Cr * * NiAl.eam.alloy Ni Al 这个要怎么改
2020年03月27日 09点03分
level 3
鞋钉ww 楼主
来挖一下自己的坟,写的可能有点啰嗦,希望能对之后的同志们起到一定点作用。
如果仅仅想要实现我上面列出目的,我给出的方法是正确的,也就是说应该使用pair_style hybrid命令,而不是pair_style hybrid/overlay。原因也和我说的一样。我简单解释一下上述用法的含义:
使用pair_style hybrid命令进行多种势函数混合有一个重要特点:在描述同种相互作用时,靠后的势函数会覆盖靠前的势函数。比如,在我的代码中,第二句表示了使用A-B.eam.alloy势函数来描述A-A/A-B/B-B三种相互作用,第三句表示了使用A.eam.fs势函数来描述A-A相互作用,所以第三句和第二句的作用有部分重复,实际上第二句中的A-A相互作用被第三局覆盖了,计算时起作用的是第三句。同理,第四句的作用是覆盖第二句中的B-B相互作用。按照这种用法,就可以实现【A-A之间使用纯金属势函数A.eam.fs,B-B之间使用纯金属势函数B.eam.fs,A-B之间使用合金势函数A-B.eam.alloy】。从上面可以看出,这三句代码的先后顺序是非常重要的,不同顺序会导致不同的效果,需要根据实际应用加以区分。
2020年05月30日 05点05分 5
大佬,我准备做的al-graphite的合金,如果用MS建模,建模的时候需要考虑力场吗?还是在lammps中实现力场势函数设置
2020年06月03日 21点06分
回复 伯努利方程丶🐶 :lammps即可,ms获得坐标就行
2020年11月19日 11点11分
level 1
如果用comb3势函数来表示两种元素见的作用,该怎么使用呢?目前能看到的例子都是eam或者lj,耽搁使用comb3势函数倒是可以,用来组合一直报错不合法。
2021年01月11日 02点01分 6
level 2
请问多元合金的势函数可以用混合势吗,eam和meam可以混用吗,就是假如是有1.A-B-C-D-E以及2.B-C-D-E-M的多主元势函数,但是我要需要的是A-B-C-D-M的的势函数,可以分别用1和2中的一部分相互作用来表示我需要的势函数吗?
2021年03月19日 01点03分 7
楼主,您解决了吗?
2021年05月12日 02点05分
我理解是可以混用的,因为我用的三种势函数也是不同的[哈哈]
2021年05月12日 15点05分
但是后者问题我可能答不了,你看看有没有肯能把所有相互作用关系覆盖全,如果覆盖全了还有可能可以用
2021年05月12日 15点05分
我之前那样组合实在是因为没有合适的势函数可以用了,自己也不会构建。我不知道你是不是在做高熵合金,你可以问问一些论文作者,我觉得也不是每种合金都有现成的势函数可以用的,看看别人是怎么做的[吐舌]
2021年05月12日 15点05分
level 1
厉害,大佬,弄了一天,没想到在这里找到了,十分感谢
2021年12月14日 11点12分 9
level 1
就是请问楼主一下,系统含有A、B、C三种原子,只用一个ABC.eam.fs势函数可不可以?能不能通过pair_coeff设定从势函数ABC.eam.fs调用(A-A)的相互作用势?
2021年12月14日 11点12分 10
1