24小时热门版块排行榜    

查看: 1943  |  回复: 15

shuting1994

新虫 (小有名气)

[求助] 建立氮化硼模型但是dump出来的结构却是碳原子石墨烯结构怎么回事? 已有2人参与

这是我的结构建模:
atom_style atomic
lattice custom 4.35 a1 1 0 0 a2 0 0.57735 0 a3 0 0 10 basis 0.167 0.25 0.5 basis 0.33 0.75 0.5 basis 0.67 0.75 0.5 basis 0.83 0.25 0.5
region box block 0 30 0 400 0 45 units box
create_box 2 box
mass 1 10
mass 2 14
create_atoms 2 region box basis 1 2 basis 2 1 basis 3 1 basis 4 2
pair_style tersoff
pair_coeff * *BNC B N
以上是建模,按道理氮化硼就是把石墨烯的碳原子换成交替的氮和硼,所以我在石墨烯的lattice custom基础上定义了2种原子,并且给出氮和硼质量,pair_style 取BNC.tersoff中的B和N,可是为什么dump原子信息一直就只有一种原子,并且显示是碳呢?哪位高手知道?非常感谢!

发自小木虫Android客户端
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gredream

禁虫 (正式写手)

本帖内容被屏蔽

» 本帖已获得的红花(最新10朵)

2楼2016-12-08 12:15:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abinitio

版主 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
pair_style tersoff  

这个tersoff是按C来写的吧

» 本帖已获得的红花(最新10朵)

3楼2016-12-08 13:50:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shuting1994

新虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by gredream at 2016-12-08 12:15:34
dump怎么写的?

我写dump 1 all custom 5000 bn.lammpstrj id element x y z的时候,输出id element x y x分别是  1(2 3 ……) C(一直显示C)  x的坐标 y的坐标 z的坐标,而当我写dump 1 all atom  1 bn.lammpstrj输出,id type xs ys zs 分别是1(2 3 ……) 1(此处1 2交换) x分子坐标 y 分子坐标 z分子坐标 虽然dump atom好像有2个type但是再看质量,都是输出的12,就算他有2个type ,也都是碳原子的意思,很苦恼,改了好长时间了,没有找出问题在哪

发自小木虫Android客户端
4楼2016-12-08 14:12:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shuting1994

新虫 (小有名气)

送红花一朵
引用回帖:
3楼: Originally posted by abinitio at 2016-12-08 13:50:01
pair_style tersoff  

这个tersoff是按C来写的吧

啊!原来是这样吗,我也看了手册里pair_coeff命令那块,实在没有找到问题,请问这种BCN三原子势函数,我只想取B和N应该怎么定义呢?

发自小木虫Android客户端
5楼2016-12-08 14:15:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abinitio

版主 (著名写手)

【答案】应助回帖

引用回帖:
5楼: Originally posted by shuting1994 at 2016-12-08 14:15:24
啊!原来是这样吗,我也看了手册里pair_coeff命令那块,实在没有找到问题,请问这种BCN三原子势函数,我只想取B和N应该怎么定义呢?
...

pair_coeff * *BNC.tersoff B N C

这样写试试,我也想不出其它原因了

» 本帖已获得的红花(最新10朵)

6楼2016-12-08 14:28:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shuting1994

新虫 (小有名气)

送红花一朵
引用回帖:
6楼: Originally posted by abinitio at 2016-12-08 14:28:26
pair_coeff * *BNC.tersoff B N C

这样写试试,我也想不出其它原因了...

谢谢!我试了改成你说的,直接报错了,我又想,如果是tersoff问题,我直接把tersoff文件里含碳的部分全删了,结果输出的还是碳原子,觉得可能是前面结构部分出问题了,但是又找不到,打算学ms去了,借用其他软件建模生成data试试

发自小木虫Android客户端
7楼2016-12-08 14:59:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abinitio

版主 (著名写手)

引用回帖:
7楼: Originally posted by shuting1994 at 2016-12-08 14:59:23
谢谢!我试了改成你说的,直接报错了,我又想,如果是tersoff问题,我直接把tersoff文件里含碳的部分全删了,结果输出的还是碳原子,觉得可能是前面结构部分出问题了,但是又找不到,打算学ms去了,借用其他软件建 ...

pair_coeff * *BNC B N   这种没有..tersoff 的不错?你将BNC.tersoff文件改名字了?

那在pair_coeff * *BNC B N C呢?就是把C也加到后面

» 本帖已获得的红花(最新10朵)

8楼2016-12-08 15:18:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shuting1994

新虫 (小有名气)

送红花一朵
引用回帖:
8楼: Originally posted by abinitio at 2016-12-08 15:18:38
pair_coeff * *BNC B N   这种没有..tersoff 的不错?你将BNC.tersoff文件改名字了?

那在pair_coeff * *BNC B N C呢?就是把C也加到后面...

没有改名字,我先试了下pair_coeff * * BNC.tersoff B N C 直接报错了,然后我以文本格式打开BNC.tersoff,把含有碳的势函数全部删除了,文件的名字没有改变,此时运行出来的结构依然是石墨烯的

发自小木虫Android客户端
9楼2016-12-08 15:26:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shuting1994

新虫 (小有名气)

问题解决了,是dump的原因,lammps的dump默认原子都为碳,需要设置一下dump_modify

发自小木虫Android客户端
10楼2016-12-10 20:24:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 shuting1994 的主题更新
信息提示
请填处理意见