24小时热门版块排行榜    

查看: 2519  |  回复: 8

HitYingph

新虫 (初入文坛)

[交流] He-C力场选择,单壁碳纳米管坍塌 已有4人参与

我想做一个单壁碳纳米管里面有氦原子存在,在外界施加载荷使得单壁碳纳米管坍塌的分子动力学模拟,基于LAMMPS,关于He与C的力场选择什么了?大家有什么建议么,我是刚入门的,在这里求问,谢谢各位。
回复此楼

» 猜你喜欢

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

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

fcwu

木虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
你的势函数是怎么设置的,一般来说,气体和碳的作用比较弱,不会有这种情况出现的。

发自小木虫Android客户端
2楼2017-07-13 00:28:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuntao118

新虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
一般来说,He和H都是LJ势函数作用吧?若是的话?直接用LJ共混率就行了!

发自小木虫Android客户端
3楼2017-07-13 07:04:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

HitYingph

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by chuntao118 at 2017-07-13 07:04:07
一般来说,He和H都是LJ势函数作用吧?若是的话?直接用LJ共混率就行了!

你好,谢谢你的回答!
但是我预备C-C之间采用airebo势,C-He 和He-He之间采用LJ/Class2势,这样该如何写相关的程序呢,实在是新手,不好意思。
4楼2017-07-13 15:46:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuntao118

新虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
有势函数,程序不会写吗?对不同粒子分类,根据分类这就行了!若是软件,我不懂!我不用软件!

发自小木虫Android客户端

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

5楼2017-07-13 17:33:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

HitYingph

新虫 (初入文坛)

送红花一朵
引用回帖:
5楼: Originally posted by chuntao118 at 2017-07-13 17:33:03
有势函数,程序不会写吗?对不同粒子分类,根据分类这就行了!若是软件,我不懂!我不用软件!

boundary s s s
lattice fcc 1.0
timestep 0.0002
atom_modify map array
atom_style        full
read_data        ${simname}.lammps05
pair_style hybrid airebo 3.0 lj/cut/coul/long 10
pair_coeff 1 1 CH.airebo C
pair_coeff 1 1 lj/cut/coul/long 0.147999 3.617048
pair_coeff 2 2 lj/cut/coul/long 0.020307 2.555988
pair_modify shift yes mix geometric  

你好,这是我力场定义的方式,其中1指的是C原子,2指的是He原子,这样写是不是对的呢?主要关于这里的Hybrid的用法我不是很确定!!!
6楼2017-07-14 09:55:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kkwzsm

银虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
6楼: Originally posted by HitYingph at 2017-07-14 09:55:28
boundary s s s
lattice fcc 1.0
timestep 0.0002
atom_modify map array
atom_style        full
read_data        ${simname}.lammps05
pair_style hybrid airebo 3.0 lj/cut/coul/long 10
pair_coeff 1 1 CH.airebo C ...

pair_style 后半部分不能用lj/cut/coul/long吧,因为C和He都不带电量,还有下面的pair_coeff 1 1 感觉出现两次?
吃的苦中苦,方为人上人。。。。
7楼2017-07-14 11:40:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

HitYingph

新虫 (初入文坛)

引用回帖:
7楼: Originally posted by kkwzsm at 2017-07-14 11:40:07
pair_style 后半部分不能用lj/cut/coul/long吧,因为C和He都不带电量,还有下面的pair_coeff 1 1 感觉出现两次?...

是 所以我只想用lj/class 2,lj/cut/coul/long是因为用MS导出生成data文件的时候里面默认是有参数的。我想用lj/class 2,但是我找不到C和He的69势不带电的参数,也是之前我问你这样的参数该怎么获取。下面的确是写错了,应该是分别定义12 和22,谢谢你的回复!!!
8楼2017-07-14 14:59:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gredream

禁虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
本帖内容被屏蔽

9楼2017-07-17 23:57:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 HitYingph 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见