24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

dgwfly

新虫 (正式写手)

应助: 0 (幼儿园)
金币: 683.9
散金: 258
红花: 9
帖子: 543
在线: 202小时
虫号: 3454332
注册: 2014-10-04
性别: MM
专业: 数论

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

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

回复此楼

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

编程

» 猜你喜欢

2025年遐想已经有4人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复
自然科学基金委宣布启动申请书“瘦身提质”行动已经有4人回复
求个博导看看已经有18人回复

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

» 本主题相关价值贴推荐，对您同样有帮助:

有关编程语言的选择已经有4人回复
势函数问题已经有8人回复
使用Multiwfn绘制原子轨道图形、研究原子壳层结构及相对论效应的影响已经有17人回复
【资源】著名物理学家吴大猷教授的理论物理丛书下载【珍藏版】已经有475人回复
回答一些关于Multiwfn的疑问以及未来Multiwfn的发展打算已经有37人回复
【sobereva个人文集】制作动画分析电子结构特征已经有29人回复
【分享】量化软件 - 输运相关 (zz) 已经有13人回复

1楼 2016-03-04 09:49:56

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

lsloneil

专家顾问 (正式写手)

专家经验: +49
模拟EPI: 2
应助: 148 (高中生)
金币: 8640.9
散金: 507
红花: 30
帖子: 505
在线: 598.4小时
虫号: 556471
注册: 2008-05-11
性别: GG
专业: 理论和计算化学
管辖: 分子模拟

【答案】应助回帖

感谢参与，应助指数 +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)，看它们的代码是怎么做的。

赞一下