| 3 | 1/1 | 返回列表 |
| 查看: 766 | 回復(fù): 2 | ||
上官若軒新蟲 (初入文壇)
|
[求助]
求fortran編程小程序 已有1人參與
|
| 我是做分子動(dòng)力學(xué)模擬的,剛剛開始學(xué)習(xí)fortran還不怎么會(huì)編程,老師讓編寫靜電相互作用和范德華相互作用的小程序,求各位大俠幫幫忙!感激不盡! |
木蟲 (小有名氣)
|
如果只是簡單做一下,不去深究如何才能更準(zhǔn)確的處理長程庫侖力的話,你可以這樣來做 (1) 推導(dǎo)出靜電相互作用和范德華相互作用的力的表達(dá)式 這一部分可以參考我的課程講義: http://kczx.suda.edu.cn/mdmc.html lectures>Lecture_slides>3.pdf (2) 如果你有一定的Fortran編程基礎(chǔ),可以直接到第三步,不然可以參考我在課程網(wǎng)站 Outreach 版塊提供的Fortran教程。 (3) 我在教學(xué)網(wǎng)站的Hands-on版塊提供了一個(gè)很基本的分子動(dòng)力學(xué)模擬Lennard-JOnes粒子的程序(lj_nve.f90),可以供你參考,在此基礎(chǔ)上,你需要加上靜電力,另外如果體系較大的話,你可能還需要添加Neighbour list等加快計(jì)算力的技巧。 (4) 如果靜電力的處理只是簡單的用一個(gè)大的cutoff distance, 那編程的難度和處理LJ勢(shì)差不多,如果要需要采用更合理的算法,如Ewald Summation, PPPM等,那你需要去學(xué)習(xí)相關(guān)知識(shí) 。 (5) 如果你需要添加neighbour list, 你可以參考Frenkel 和Smit的 "understanding Molecular Simulation"的Appendix F, 這個(gè)時(shí)候計(jì)算力的模塊自然也需要相應(yīng)改變。 另外, Furio Ercolessi在他的網(wǎng)站上提供了一些MD的 F90程序,也可以供你參考 : http://www.fisica.uniud.it/~ercolessi/md/f90/ |

木蟲 (著名寫手)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 274求調(diào)劑 +5 | S.H1 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 材料專業(yè)求調(diào)劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8 | Ncdx123456 2026-03-13 | 9/450 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +7 | rare12345 2026-03-18 | 7/350 |
|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +6 | 薛云鵬 2026-03-13 | 6/300 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號(hào) 2026-03-16 | 16/800 |
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +16 | 不知道叫什么! 2026-03-15 | 18/900 |
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|