| 查看: 1132 | 回复: 3 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
势函数的编写问题 已有1人参与
|
||||
| 老板让用fortran做分子动力学的模拟。之前有Lennard-Jones势能函数简单的NVE和NVT系宗的程序。但现在不知道怎么用S-W势函数或者其它描述固体的势函数来编写程序？主要是气固两相的模拟。求助相关Fortran代码的例子 |
回复此楼» 收录本帖的淘帖专辑推荐
 编程 |
» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有41人回复
-
filecode
已经有9人回复
-
请教兼职经验
已经有3人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有5人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有13人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有11人回复
-
读博的好处
已经有4人回复
-
filecode
已经有17人回复
-
FileCode能看出啥？
已经有32人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有15人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 有关编程语言的选择
已经有4人回复
- 势函数问题
已经有8人回复
- 使用Multiwfn绘制原子轨道图形、研究原子壳层结构及相对论效应的影响
已经有17人回复
- 【资源】著名物理学家吴大猷教授的理论物理丛书下载【珍藏版】
已经有475人回复
- 回答一些关于Multiwfn的疑问以及未来Multiwfn的发展打算
已经有37人回复
- 【sobereva个人文集】制作动画分析电子结构特征
已经有29人回复
- 【分享】量化软件 - 输运相关 (zz)
已经有13人回复
lvxinyuabc
木虫 (正式写手)
材料大数据
- 应助: 3 (幼儿园)
- 金币: 3083.9
- 散金: 5
- 红花: 2
- 帖子: 520
- 在线: 211.4小时
- 虫号: 1043016
- 注册: 2010-06-17
- 专业: 金属基复合材料
4楼2016-03-05 19:50:23
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
