24小时热门版块排行榜    

查看: 2613  |  回复: 4

hyliao

铜虫 (小有名气)

[求助] lammps 吸附 cnt与ni

本人金币不多,还请大家原谅。

请大家帮下忙,这段时间在做cnt对ni的吸附。但是计算结果总是有问题,也不知道出在哪里?希望大家多多指教。
计算吸附:我是根据参考文件的方法。先分别计算cnt与ni的能量,然后计算ni吸附在cnt上的能量(只有一个ni原子)。(建模是ni在cnt六元环上方1.5A)
其中cnt,CNT+NI是从ms中导出来的。计算后E(CNT)+E(NI)-E(CNT+NI)居然小于零,不易吸附。这与参考文献的结果完全相反啊!
这个问题不解决,后面的工作就很难进行。希望大家帮忙看看问题出在哪里?
#lammps for cnt_ni

units  metal
boundary s s p
atom_style atomic
read_data SWNTcenter.data
pair_style hybrid airebo 3.4 1 1 lj/cut 9.0
pair_coeff * * airebo CH.airebo C NULL
pair_coeff 1 2 lj/cut 0.0688 3.57
pair_coeff 2 2 none
run  0
timestep 0.001
velocity all create 0.00000001 100000000 rot yes dist gaussian
#unfix  1
min_style cg
minimize 1.0e-4 1.0e-6 10000 10000
min_style sd
minimize 0 0 100000 1000000
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

allanjames

新虫 (初入文坛)

【答案】应助回帖


hyliao(金币+3): in文件就这些了,我只用了静力学,能量最小化。cutoff 9 是我从文献中选的。 2011-06-07 16:43:38
ghcacj(金币+1): 谢谢 2011-06-07 20:49:42
in 文件不全阿,cutoff 9 是不是大了些
思考才能进步
2楼2011-06-07 16:13:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hyliao

铜虫 (小有名气)

我这几天考虑了一下,我的data文件中,只有c和ni 的坐标和质量。而lammps的in文件中,考虑单个ni 原子的吸附时。由于cnt取得较长。所以可以不考虑ni与ni 的相互作用。最后只有cnt的碳原子之间的airobo作用 和ni 与cnt间的lj 作用。如果这样的话。似乎只是吸附原子的质量 是唯一的变量。吸附也只与质量有关了。这个过程中,没有考虑原子的电子特性。这当然不合理。
请问大家一下,这个问题该如何解决。
3楼2011-06-07 16:53:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

J_destiny___

铁虫 (初入文坛)

请问楼主 这个问题解决了没?已经解决的话,最后是怎么解决的哈
4楼2012-10-29 09:18:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hkuieagle

新虫 (初入文坛)

楼主解决了吗?
5楼2012-12-05 10:17:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hyliao 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 7lpolszZVXgi 2026-08-14 6/300 2026-08-15 04:21 by 4wMiSEwB6436
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 4/200 2026-08-15 04:17 by 4wMiSEwB6436
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 02:21 by 4wMiSEwB6436
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +11 Tide man 2026-08-10 12/600 2026-08-15 02:12 by home3163
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 01:28 by 4wMiSEwB6436
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 7lpolszZVXgi 2026-08-14 4/200 2026-08-15 01:08 by 4wMiSEwB6436
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +5 Ldrop2023 2026-08-13 5/250 2026-08-14 18:48 by ssxclkj
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
信息提示
请填处理意见