| 5 | 1/1 | 返回列表 |
| 查看: 768 | 回復(fù): 2 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
上官若軒新蟲 (初入文壇)
|
[求助]
求fortran編程小程序 已有1人參與
|
|
| 我是做分子動(dòng)力學(xué)模擬的,剛剛開(kāi)始學(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/ |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩英一數(shù)二306 +5 | z1z2z3879 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑 +8 | 妮妮ninicgb 2026-03-15 | 12/600 |
|
|
[考研] 0703化學(xué)調(diào)劑,求各位老師收留 +10 | 秋有木北 2026-03-14 | 10/500 |
|
|
[考研] 330求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
[考研] 321求調(diào)劑 +3 | 何潤(rùn)采123 2026-03-18 | 3/150 |
|
|
[考研] 328求調(diào)劑,英語(yǔ)六級(jí)551,有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-16 | 10/500 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +7 | rare12345 2026-03-18 | 7/350 |
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 290求調(diào)劑 +6 | 孔志浩 2026-03-12 | 11/550 |
|
|
[基金申請(qǐng)] 國(guó)自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 0703一志愿211 285分求調(diào)劑 +5 | ly3471z 2026-03-13 | 5/250 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼?~+123456 2026-03-14 | 3/150 |
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|