24小时热门版块排行榜    

查看: 1150  |  回复: 3

dgwfly

新虫 (正式写手)

[求助] 势函数的编写问题 已有1人参与

老板让用fortran做分子动力学的模拟。之前有Lennard-Jones势能函数简单的NVE和NVT系宗的程序。但现在不知道怎么用S-W势函数或者其它描述固体的势函数来编写程序?主要是气固两相的模拟。求助相关Fortran代码的例子
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)

本帖仅楼主可见
2楼2016-03-04 12:11:05
已阅   申请模拟EPI   回复此楼   编辑   查看我的主页

lsloneil

专家顾问 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
Stillinger-Webber是三体势,意味着相比起两体势,要多一重循环,比方说你计算LJ势的module里,是
do i=1, N
    do j=i, N
        计算i与j之间相互作用
    end do
end do
对于SW,循环就变成
do i=1, N
    do j=i, N
         do k=j,N
          计算i,j,k的三体势能
         end do
     end do
end do

SW因为是三体势,计算起来比二体势会慢不少。如果你用的程序不是优化得非常好,为了追求速度的话建议用现成的开源软件包,像LAMMPS, DL_POLY里面都有现成的SW势可以用。如果你想参考代码,可以下载LAMMPS(C++)或者DL_POLY (Fortran),看它们的代码是怎么做的。
3楼2016-03-04 14:19:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lvxinyuabc

木虫 (正式写手)

材料大数据

感觉很高深,希望与你交流下。
专心科研,分享成长心路历程。关注个人博客www.yffslxy.com,共同鉴证成长之路.
4楼2016-03-05 19:50:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dgwfly 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +3 gouxfjh 2026-08-28 3/150 2026-08-29 01:01 by jurkat.1640
[基金申请] 国自然评审意见 +11 wangmingqi 2026-08-28 15/750 2026-08-29 00:17 by fangyl2005
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +13 晴天加油 2026-08-26 14/700 2026-08-28 18:36 by huagongfeihu
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 有没有仍没收到信息的 +3 德尚中行 2026-08-27 3/150 2026-08-28 15:24 by rongshuxia
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
请填处理意见