应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » 其他 » 势函数的编写问题
41/1返回列表
查看: 1166  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

dgwfly

新虫 (正式写手)

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

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

» 收录本帖的淘帖专辑推荐

编程

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2016-03-04 09:49:56
已阅   回复此楼   关注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 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭