| 4 | 1/1 | 返回列表 |
| 查看: 993 | 回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +12 | zjl的號 2026-03-16 | 17/850 |
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +3 | vv迷 2026-03-22 | 3/150 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
[基金申請] 山東省面上項目限額評審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
[考研] 求調(diào)劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研]
|
然11 2026-03-19 | 4/200 |
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|