| 3 | 1/1 | 返回列表 |
| 查看: 767 | 回復(fù): 2 | |||
上官若軒新蟲 (初入文壇)
|
[求助]
求fortran編程小程序 已有1人參與
|
| 我是做分子動(dòng)力學(xué)模擬的,剛剛開始學(xué)習(xí)fortran還不怎么會(huì)編程,老師讓編寫靜電相互作用和范德華相互作用的小程序,求各位大俠幫幫忙!感激不盡! |
木蟲 (小有名氣)
|
如果只是簡(jiǎn)單做一下,不去深究如何才能更準(zhǔn)確的處理長(zhǎng)程庫(kù)侖力的話,你可以這樣來(lái)做 (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) 如果靜電力的處理只是簡(jiǎn)單的用一個(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ā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +4 | Ymlll 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 材料專碩英一數(shù)二306 +4 | z1z2z3879 2026-03-18 | 4/200 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +6 | 慕寒mio 2026-03-16 | 6/300 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +4 | ioodiiij 2026-03-17 | 4/200 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家,不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|