| 4 | 1/1 | 返回列表 |
| 查看: 1000 | 回復(fù): 3 | ||||
[求助]
勢函數(shù)的編寫問題 已有1人參與
|
| 老板讓用fortran做分子動力學(xué)的模擬。之前有Lennard-Jones勢能函數(shù)簡單的NVE和NVT系宗的程序。但現(xiàn)在不知道怎么用S-W勢函數(shù)或者其它描述固體的勢函數(shù)來編寫程序?主要是氣固兩相的模擬。求助相關(guān)Fortran代碼的例子 |
編程 |
用戶注銷 (小有名氣)
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +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ōu)化得非常好,為了追求速度的話建議用現(xiàn)成的開源軟件包,像LAMMPS, DL_POLY里面都有現(xiàn)成的SW勢可以用。如果你想?yún)⒖即a,可以下載LAMMPS(C++)或者DL_POLY (Fortran),看它們的代碼是怎么做的。 |
木蟲 (正式寫手)
材料大數(shù)據(jù)

| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 328求調(diào)劑 +7 | 嗯滴的基本都 2026-03-27 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
[考研] 286求調(diào)劑 +3 | 丟掉懶惰 2026-03-27 | 6/300 |
|
|
[考研] 266分,求材料冶金能源化工等調(diào)劑 +5 | 哇呼哼呼哼 2026-03-27 | 7/350 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個名咋這么難 2026-03-27 | 4/200 |
|
|
[考研] 考研調(diào)劑 +4 | Sanmu-124 2026-03-26 | 4/200 |
|
|
[考博] 26申博 +3 | 加油沖啊! 2026-03-26 | 3/150 |
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
[考研] 316求調(diào)劑 +5 | 江辭666 2026-03-26 | 5/250 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 一志愿華理,數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 10/500 |
|
|
[考研] 271求調(diào)劑 +6 | 生如夏花… 2026-03-22 | 6/300 |
|
|
[考研] 機械學(xué)碩310分,數(shù)一英一,一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
[考研] 297求調(diào)劑 +6 | 田洪有 2026-03-26 | 6/300 |
|
|
[考研] 打過很多競賽,085406控制工程300分,求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
[考研] 考研調(diào)劑 +6 | 來好運來來來 2026-03-21 | 7/350 |
|
|
[考研] 284求調(diào)劑 +15 | Zhao anqi 2026-03-22 | 15/750 |
|
|
[考研] 一志愿南航材料專317分求調(diào)劑 +5 | 炸呀炸呀炸薯條 2026-03-23 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|