| 查看: 1142 | 回复: 3 | ||||
[求助]
势函数的编写问题 已有1人参与
|
| 老板让用fortran做分子动力学的模拟。之前有Lennard-Jones势能函数简单的NVE和NVT系宗的程序。但现在不知道怎么用S-W势函数或者其它描述固体的势函数来编写程序？主要是气固两相的模拟。求助相关Fortran代码的例子 |
回复此楼» 收录本帖的淘帖专辑推荐
 编程 |
» 猜你喜欢
93BebMhtakh前后11位开头都是大写
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有86人回复
-
感觉是下周放榜了
已经有11人回复
-
filecode=后面第一个是大写字母
已经有10人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 有关编程语言的选择
已经有4人回复
- 势函数问题
已经有8人回复
- 使用Multiwfn绘制原子轨道图形、研究原子壳层结构及相对论效应的影响
已经有17人回复
- 【资源】著名物理学家吴大猷教授的理论物理丛书下载【珍藏版】
已经有475人回复
- 回答一些关于Multiwfn的疑问以及未来Multiwfn的发展打算
已经有37人回复
- 【sobereva个人文集】制作动画分析电子结构特征
已经有29人回复
- 【分享】量化软件 - 输运相关 (zz)
已经有13人回复
匿名
用户注销 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1383.8
- 散金: 210
- 红花: 4
- 帖子: 283
- 在线: 67小时
- 虫号: 0
- 注册: 2013-05-18
- 性别: GG
- 专业: 社会心理学
 本帖仅楼主可见 |
2楼2016-03-04 12:11:05
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)，看它们的代码是怎么做的。 |
3楼2016-03-04 14:19:06
lvxinyuabc
木虫 (正式写手)
材料大数据
- 应助: 3 (幼儿园)
- 金币: 3083.9
- 散金: 5
- 红花: 2
- 帖子: 520
- 在线: 211.4小时
- 虫号: 1043016
- 注册: 2010-06-17
- 专业: 金属基复合材料
4楼2016-03-05 19:50:23
