24小时热门版块排行榜    

查看: 773  |  回复: 2

上官若轩

新虫 (初入文坛)

[求助] 求fortran编程小程序 已有1人参与

我是做分子动力学模拟的,刚刚开始学习fortran还不怎么会编程,老师让编写静电相互作用和范德华相互作用的小程序,求各位大侠帮帮忙!感激不尽!
回复此楼

» 猜你喜欢

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

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

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

ktwyw

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

戴世杰

木虫 (著名写手)

引用回帖:
2楼: Originally posted by ktwyw at 2014-10-31 09:33:19
如果只是简单做一下,不去深究如何才能更准确的处理长程库仑力的话,你可以这样来做
(1) 推导出静电相互作用和范德华相互作用的力的表达式
这一部分可以参考我的课程讲义:
http://kczx.suda.edu.cn/mdmc.html
...

赞一个

[ 发自小木虫客户端 ]
3楼2014-11-02 09:48:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 上官若轩 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 298求调剂 +7 上岸6666@ 2026-03-20 7/350 2026-03-23 10:24 by 呆呆师姐
[考研] 293求调剂 +3 涛涛Wjt 2026-03-22 5/250 2026-03-22 22:21 by jiangpengfei
[考研] 318求调剂 +4 plum李子 2026-03-21 7/350 2026-03-22 14:17 by ColorlessPI
[考研] 286分人工智能专业请求调剂愿意跨考! +4 lemonzzn 2026-03-17 8/400 2026-03-21 22:49 by lemonzzn
[考研] 考研调剂 +3 呼呼?~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 一志愿南大,0703化学,分数336,求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +11 吃吃吃才有意义 2026-03-19 11/550 2026-03-21 18:23 by 学员8dgXkO
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-18 3/150 2026-03-21 18:17 by 学员8dgXkO
[考研] 材料与化工(0856)304求 B区 调剂 +3 邱gl 2026-03-21 3/150 2026-03-21 13:47 by lature00
[考研] 279分求调剂 一志愿211 +14 chaojifeixia 2026-03-19 15/750 2026-03-21 13:24 by zhukairuo
[考研] 一志愿中国石油大学(华东) 本科齐鲁工业大学 +3 石能伟 2026-03-17 3/150 2026-03-21 02:22 by JourneyLucky
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 广西大学家禽遗传育种课题组2026年硕士招生(接收计算机专业调剂) +3 123阿标 2026-03-17 3/150 2026-03-20 15:58 by 飞行琦
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
信息提示
请填处理意见