24小时热门版块排行榜    

查看: 2810  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[基金申请] 为什么国自然不能直接公布 +5 bjdxyxy 2026-08-26 5/250 2026-09-01 17:54 by 小伟大博士
[文学芳草园] 梦想 +5 myrtle 2026-08-26 7/350 2026-09-01 15:18 by myrtle
[基金申请] 学科评审组评审是指会评吗? +4 瞬息宇宙 2026-08-31 4/200 2026-09-01 14:58 by jiaoxg
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见