| 3 | 1/1 | 返回列表 |
| 查看: 770 | 回復: 2 | |||
[求助]
求fortran編程小程序 已有1人參與
|
| 我是做分子動力學模擬的,剛剛開始學習fortran還不怎么會編程,老師讓編寫靜電相互作用和范德華相互作用的小程序,求各位大俠幫幫忙!感激不盡! |
木蟲 (小有名氣)
|
如果只是簡單做一下,不去深究如何才能更準確的處理長程庫侖力的話,你可以這樣來做 (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/ |

木蟲 (著名寫手)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 294求調劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|---|---|---|---|---|
|
[考研] 323求調劑 +3 | 洼小桶 2026-03-18 | 3/150 |
|
|
[考研] 288求調劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研] 一志愿北京化工大學0703化學318分,有科研經歷,求調劑 +4 | 一瓶苯甲酸 2026-03-14 | 4/200 |
|
|
[考研] 279分求調劑 一志愿211 +10 | chaojifeixia 2026-03-19 | 11/550 |
|
|
[考研] 一志愿吉林大學材料學碩321求調劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
[考研] 求調劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
[考研] 086500 325 求調劑 +3 | 領帶小熊 2026-03-19 | 3/150 |
|
|
[考研] 環(huán)境工程調劑 +9 | 大可digkids 2026-03-16 | 9/450 |
|
|
[考研] 265求調劑 +8 | 梁梁校校 2026-03-17 | 8/400 |
|
|
[考研] 0856調劑,是學校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 286分人工智能專業(yè)請求調劑愿意跨考! +3 | lemonzzn 2026-03-17 | 4/200 |
|
|
[考研] 0703化學 305求調劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
[考研] 0854可跨調劑,一作一項核心論文五項專利,省、國級證書40+數一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 301求調劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[考研] 308求調劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 一志愿蘇州大學材料工程(085601)專碩有科研經歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研] 070303 總分349求調劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 0856專碩279求調劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 289求調劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|