| 4 | 1/1 | 返回列表 |
| 查看: 990 | 回復(fù): 3 | ||||
[求助]
勢函數(shù)的編寫問題 已有1人參與
|
| 老板讓用fortran做分子動力學(xué)的模擬。之前有Lennard-Jones勢能函數(shù)簡單的NVE和NVT系宗的程序。但現(xiàn)在不知道怎么用S-W勢函數(shù)或者其它描述固體的勢函數(shù)來編寫程序?主要是氣固兩相的模擬。求助相關(guān)Fortran代碼的例子 |
編程 |
用戶注銷 (小有名氣)
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗(yàn): +49 |
|
Stillinger-Webber是三體勢,意味著相比起兩體勢,要多一重循環(huán),比方說你計算LJ勢的module里,是 do i=1, N do j=i, N 計算i與j之間相互作用 end do end do 對于SW,循環(huán)就變成 do i=1, N do j=i, N do k=j,N 計算i,j,k的三體勢能 end do end do end do SW因?yàn)槭侨w勢,計算起來比二體勢會慢不少。如果你用的程序不是優(yōu)化得非常好,為了追求速度的話建議用現(xiàn)成的開源軟件包,像LAMMPS, DL_POLY里面都有現(xiàn)成的SW勢可以用。如果你想?yún)⒖即a,可以下載LAMMPS(C++)或者DL_POLY (Fortran),看它們的代碼是怎么做的。 |
木蟲 (正式寫手)
材料大數(shù)據(jù)

| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +3 | niko- 2026-03-22 | 3/150 |
|
|
[考研] 269專碩求調(diào)劑 +6 | 金恩貝 2026-03-21 | 6/300 |
|
|
[考研] 291 求調(diào)劑 +3 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 3/150 |
|
|
[考研] 考研調(diào)劑 +4 | 來好運(yùn)來來來 2026-03-21 | 4/200 |
|
|
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考! +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 310求調(diào)劑 +3 | baibai1314 2026-03-16 | 3/150 |
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 6/300 |
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
[考研]
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 環(huán)境工程調(diào)劑 +9 | 大可digkids 2026-03-16 | 9/450 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|