| 查看: 728 | 回复: 2 | ||
[求助]
求fortran编程小程序已有1人参与
|
| 我是做分子动力学模拟的，刚刚开始学习fortran还不怎么会编程，老师让编写静电相互作用和范德华相互作用的小程序，求各位大侠帮帮忙！感激不尽！ |
回复此楼» 猜你喜欢
导师想让我从独立一作变成了共一第一
已经有9人回复
-
博士读完未来一定会好吗
已经有23人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
申请2026年博士
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
关于FORTRAN使用DAT文件
已经有5人回复- 急：求用Fortran或者其他语言编写一个程序，处理下面的数据！！！！！！！！！！
已经有37人回复
- fortran编译中$符号的使用
已经有4人回复
- 求一个简单的读写文件的小程序。。。
已经有4人回复
- 求直接模拟蒙特卡罗方法的程序，包括二维和三维的用fortran编写的，谢谢
已经有9人回复
- fortran程序在SSH计算中心上运行的问题！
已经有13人回复
- 编一个fortran程序
已经有17人回复
- 关于fortran各种版本的运算速度的问题
已经有3人回复
- 求助：用Fortran写一个关于晶粒个数统计和编号的程序。卡了我好久了。
已经有4人回复
- fortran程序open文件时的简单问题
已经有3人回复
- 请教一个fortran小程序编译出错的问题，谢谢
已经有9人回复
- 如何用fortran编写高斯白噪声程序
已经有6人回复
- 【求助】Fortran语言赋值问题？
已经有3人回复
- 写了一个fortran90的小程序，编译通不过，请大侠帮忙
已经有59人回复
- 【求助】如何写FORTRAN程序实现求平均最近邻距离
已经有34人回复
- 【求助】用fortran怎么描述下面这样的情况，我是新手，请大家帮忙
已经有8人回复
- 【求助】求Fortran90编译器？
已经有11人回复
- 【求助】这样的公式，在fortran里该怎么写了，谢谢【求助成功】
已经有4人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
上官若轩: 金币+20, ★★★★★最佳答案, 很详细，谢谢帮助！ 2014-11-04 11:21:20
感谢参与，应助指数 +1
上官若轩: 金币+20, ★★★★★最佳答案, 很详细，谢谢帮助！ 2014-11-04 11:21:20
|
如果只是简单做一下，不去深究如何才能更准确的处理长程库仑力的话，你可以这样来做 (1) 推导出静电相互作用和范德华相互作用的力的表达式 这一部分可以参考我的课程讲义： http://kczx.suda.edu.cn/mdmc.html lectures>Lecture_slides>3.pdf (2) 如果你有一定的Fortran编程基础，可以直接到第三步，不然可以参考我在课程网站 Outreach 版块提供的Fortran教程。 (3) 我在教学网站的Hands-on版块提供了一个很基本的分子动力学模拟Lennard-JOnes粒子的程序（lj_nve.f90），可以供你参考，在此基础上，你需要加上静电力，另外如果体系较大的话，你可能还需要添加Neighbour list等加快计算力的技巧。 (4) 如果静电力的处理只是简单的用一个大的cutoff distance, 那编程的难度和处理LJ势差不多，如果要需要采用更合理的算法，如Ewald Summation， PPPM等，那你需要去学习相关知识 。 (5) 如果你需要添加neighbour list, 你可以参考Frenkel 和Smit的 "understanding Molecular Simulation"的Appendix F， 这个时候计算力的模块自然也需要相应改变。 另外， Furio Ercolessi在他的网站上提供了一些MD的 F90程序，也可以供你参考 ： http://www.fisica.uniud.it/~ercolessi/md/f90/ |
2楼2014-10-31 09:33:19
戴世杰
木虫 (著名写手)
- 应助: 103 (高中生)
- 金币: 4018.4
- 散金: 1099
- 红花: 18
- 帖子: 1218
- 在线: 745.7小时
- 虫号: 1275484
- 注册: 2011-04-23
- 专业: 食品加工技术
3楼2014-11-02 09:48:27