24小时热门版块排行榜    

查看: 1320  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

李杜醉酒

新虫 (小有名气)


已领完
lammps 计算原子邻表
领取红包 (小木虫手机app专属红包)

扫一扫,下载小木虫客户端

lammps中可以通过 cna/atom 命令计算每个原子的邻近数,但是我要计算每个原子邻近原子对应的 id,该如何计算。假设只计算完美晶体每个原子第一邻近原子 id,所有原子的邻近数一样,考虑周期性结构

发自小木虫Android客户端
回复此楼

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

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

李杜醉酒

新虫 (小有名气)


送红花一朵
引用回帖:
5楼: Originally posted by brucefan at 2018-04-12 03:31:45
自己写个程序算。如果粒子数小于一万,用简单的二重循环即可。如果粒子数大于一万,可考虑用cell list方法。前者几十行代码,后者也只有一两百行。

谢谢老师,目前粒子数有60w,目前我自己编的是si完美晶体,先找出某个晶胞内所有原子邻近原子id,再按照编号对比找出其它原子邻近原子编号,最后再对边界邻近原子进行修改。但是如果一开始晶体中有缺陷,比如位错,其原子序号规律就无法对比了。因此需要每个原子都需要进行一次邻近原子编号搜索,其时间复杂度O(n?)。计算量很大。请问cell list方法是?

发自小木虫Android客户端
6楼2018-04-12 10:49:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
abinitio: 金币+5, 应助指数+1, 推荐前者,简单粗暴,可靠有效 2018-04-12 08:11:54
自己写个程序算。如果粒子数小于一万,用简单的二重循环即可。如果粒子数大于一万,可考虑用cell list方法。前者几十行代码,后者也只有一两百行。

[ 发自手机版 http://muchong.com/3g ]

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

5楼2018-04-12 03:31:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
6楼: Originally posted by 李杜醉酒 at 2018-04-12 10:49:38
谢谢老师,目前粒子数有60w,目前我自己编的是si完美晶体,先找出某个晶胞内所有原子邻近原子id,再按照编号对比找出其它原子邻近原子编号,最后再对边界邻近原子进行修改。但是如果一开始晶体中有缺陷,比如 ...

https://en.m.wikipedia.org/wiki/Cell_lists

[ 发自手机版 http://muchong.com/3g ]
7楼2018-04-12 13:27:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
pandongzi2楼
2018-04-11 19:25   回复  
李杜醉酒(金币+1): 谢谢参与
已获得1个金币 发自小木虫IOS客户端
51727963楼
2018-04-11 19:26   回复  
李杜醉酒(金币+2): 谢谢参与
已获得2个金币 发自小木虫IOS客户端
对长4楼
2018-04-11 19:26   回复  
李杜醉酒(金币+2): 谢谢参与
已获得2个金币 发自小木虫IOS客户端
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +7 瞬息宇宙 2026-09-01 14/700 2026-09-02 17:33 by ma0526
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
[基金申请] 学科评审组评审是指会评吗? +5 瞬息宇宙 2026-08-31 5/250 2026-09-02 10:13 by 雨冰共舞
[基金申请] 面上合作单位盖章 +6 ssyjh 2026-08-27 9/450 2026-09-01 20:01 by huagongfeihu
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +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 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 面上意见出来了 +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 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知? +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 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +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
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
信息提示
请填处理意见