24小时热门版块排行榜    

查看: 1029  |  回复: 2

dry6211

新虫 (初入文坛)

[求助] 碳纳米管设成rigid

研究CNT在Ar剪切流体中的运动情况,希望时间步长能够取大一点儿,这样计算的快一些,因此考虑将CNT设成刚体,用fix rigid,但是Ar原子出现了lost atoms,程序关键部分如下:
fix               6 cnt rigid single
compute           mobile ar temp/partial 1 0 1                         (计算ar的温度,流动方向y的速度不算)
velocity          ar create 300.0 341233 temp mobile
fix           1 nrigid nve                                                           (nrigid代表不是CNT的流体部分,用NVE)
fix           2 ar temp/rescale 100 300.0 300.0 5.0 1.0
fix_modify        2 temp mobile

另外:
neighbor           2.0 bin
neigh_modify       delay 5 exclude group cnt cnt
timestep          0.01                                                        (10fs,这个时间对于纯Ar流体的计算是合理的)

原本考虑C-C键相互作用时timestep是0.001(1fs),没问题,现在CNT是刚体时设成0.001也行,但是增大就不行了,另外这个neighbor怎么设我也不太清楚
小弟求各位大牛帮忙!!!!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
chaizhm: 金币+5, 谢谢~ 2012-09-18 09:06:27
lost atoms一般来说有如下原因:
1)时间步长太大
2)初始构型不合理
3)力场不对
4)另外一个原因就是体系不稳定,而调控温度的间隔过大,也就是在fix 1 nvt temp t1 t2 t3中的t3太大。


另外要确认你的单位。
2楼2012-09-17 23:57:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dry6211

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by jiaoyixiong at 2012-09-17 23:57:04
lost atoms一般来说有如下原因:
1)时间步长太大
2)初始构型不合理
3)力场不对
4)另外一个原因就是体系不稳定,而调控温度的间隔过大,也就是在fix 1 nvt temp t1 t2 t3中的t3太大。


另外要确认你的单 ...

1.时间步长有意取大,如果rigid建立合理就应该没问题
2.初始构型肯定对,在别的问题中算过
3.力场很有可能有问题,可是现在看不出来,还请指教
4.没用到nvt
还请版主再指点一下
3楼2012-09-18 00:17:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dry6211 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 XlxB7q11Bc3x 2026-09-01 3/150 2026-09-01 21:52 by Y3OMSsg8u9zj
[基金申请] 面上合作单位盖章 +6 ssyjh 2026-08-27 9/450 2026-09-01 20:01 by huagongfeihu
[基金申请] 为什么国自然不能直接公布 +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
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +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
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见