24小时热门版块排行榜    

查看: 2053  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 070300化学求调剂 +3 苑豆豆 2026-03-20 3/150 2026-03-22 14:48 by ColorlessPI
[考研] 求调剂院校信息 +5 CX 330 2026-03-21 5/250 2026-03-22 11:32 by 杨杨杨紫
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5 脱颖而出 2026-03-16 16/800 2026-03-22 11:20 by 脱颖而出
[考研] 化学工程321分求调剂 +18 大米饭! 2026-03-15 22/1100 2026-03-21 20:20 by HH领袖
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 279分求调剂 一志愿211 +14 chaojifeixia 2026-03-19 15/750 2026-03-21 13:24 by zhukairuo
[考研] 265求调剂 +3 Jack?k?y 2026-03-17 3/150 2026-03-21 03:17 by JourneyLucky
[考研] 化学求调剂 +4 临泽境llllll 2026-03-17 5/250 2026-03-21 02:23 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 289求调剂 +6 怀瑾握瑜l 2026-03-20 6/300 2026-03-20 20:30 by 学员8dgXkO
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
信息提示
请填处理意见