 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
小木蟲論壇-學術科研互動平臺 » 計算模擬區(qū) » 分子模擬 » MS » 徑向分布函數(shù)怎么畫
| 5 | 1/1 | 返回列表 |
| 查看: 2541 | 回復: 5 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
秋之愚新蟲 (正式寫手)
|
[求助]
徑向分布函數(shù)怎么畫 已有1人參與
|
|
|
求助各位大神，我在做Lennard-Jones勢流體的分子動力學模擬時，能量圖基本上已經(jīng)畫出來了，可徑向分布函數(shù)的計算怎么弄都弄不出來，能有人幫幫我嗎？在給出來的這個程序中switch在主程序中該怎么調(diào)用��？？？ subroutine gr(switch,box,density,rx,ry,rz,g) dimension rx(108),ry(108),rz(108),g(125) double precision rx,ry,rz,g double precision delg,pi,box,r22,rr,density,nid integer ngr,ig,nhis,switch nhis=125 delg = box/(2*nhis) pi=3.1415926 ! if (switch.eq.0)then ngr=0 !delg=box/(2*nhis) do i=1,nhis g(i)=0 enddo ! elseif (switch.eq.1)then ngr=ngr+1 do i=1,107 do j=i+1,108 xr=rx(i)-rx(j) xr=xr-box*nint(xr/box) yr=ry(i)-ry(j) yr=yr-box*nint(yr/box) zr=rz(i)-rz(j) zr=zr-box*nint(zr/box) r22=xr*xr+yr*yr+zr*zr rr=sqrt(r22) if(rr.lt.box/2) then ig=int(rr/delg) g(ig)=g(ig)+2 endif enddo enddo ! elseif(switch.eq.2) then do i=1,nhis r=delg*(i+0.5) vb=((i+1)**3-i**3)*delg**3 nid=(4/3)*pi*vb*density g(i)=g(i)/(ngr*108.0*nid) write(*,*) nid,r,g(i) enddo ! endif return end |
回復此樓» 猜你喜歡
材料科學與工程求調(diào)劑
已經(jīng)有8人回復
-
一志愿大連理工大學材料求調(diào)劑
已經(jīng)有3人回復
-
一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑
已經(jīng)有9人回復
-
材料求調(diào)劑 一志愿哈工大總分298分，前三科223分
已經(jīng)有8人回復
-
274求調(diào)劑
已經(jīng)有3人回復
-
合肥區(qū)域性重點一本招收調(diào)劑
已經(jīng)有4人回復
-
085600材料與化工調(diào)劑
已經(jīng)有10人回復
-
22408 359分調(diào)劑
已經(jīng)有7人回復
-
抱歉
已經(jīng)有4人回復
-
085602化工求調(diào)劑（331分）
已經(jīng)有8人回復
» 本主題相關商家推薦: (我也要在這里推廣)
» 本主題相關價值貼推薦，對您同樣有幫助:
- 怎樣用forcite畫徑向分布函數(shù)
已經(jīng)有14人回復
- 關于徑向分布函數(shù)和配位數(shù)公式中的密度
已經(jīng)有3人回復
- 徑向分布函數(shù)與團簇結構關系
已經(jīng)有5人回復
- 如何解釋徑向分布函數(shù)與r之間的關系？
已經(jīng)有8人回復
- 徑向分布函數(shù)的問題
已經(jīng)有13人回復
- vmd-做徑向分布函數(shù)
已經(jīng)有6人回復
- forcite徑向分布函數(shù)g(r)物理意義
已經(jīng)有8人回復
- 請教 徑向分布函數(shù)（ＲＤＦ） 問題
已經(jīng)有12人回復
- 如何做徑向分布函數(shù)和長程序LRO？
已經(jīng)有4人回復
- 請問大家做徑向分布函數(shù)的時候結果有沒有趨近于1��？
已經(jīng)有7人回復
- 【求助】怎么從結果中提取徑向分布函數(shù)、鍵長、鍵角
已經(jīng)有17人回復
- 【討論】陳正隆書上徑向分布函數(shù)的公式是否有誤？
已經(jīng)有10人回復
- 【求助】應用MS計算徑向分布函數(shù)
已經(jīng)有15人回復
- 【求助成功】徑向分布函數(shù)如何分析？
已經(jīng)有14人回復
dummy1
木蟲 (正式寫手)
- 應助: 13 (小學生)
- 金幣: 5841.7
- 散金: 32
- 紅花: 1
- 帖子: 312
- 在線: 295.8小時
- 蟲號: 270440
- 注冊: 2006-08-06
- 專業(yè): 理論和計算化學
5樓2014-08-16 02:33:13
dummy1
木蟲 (正式寫手)
- 應助: 13 (小學生)
- 金幣: 5841.7
- 散金: 32
- 紅花: 1
- 帖子: 312
- 在線: 295.8小時
- 蟲號: 270440
- 注冊: 2006-08-06
- 專業(yè): 理論和計算化學
3樓2014-08-14 02:45:40
秋之愚
新蟲 (正式寫手)
- 應助: 0 (幼兒園)
- 金幣: 8145.6
- 紅花: 1
- 帖子: 719
- 在線: 72.3小時
- 蟲號: 2538185
- 注冊: 2013-07-09
- 性別: GG
- 專業(yè): 工程地質(zhì)學
4樓2014-08-14 19:01:25
秋之愚
新蟲 (正式寫手)
- 應助: 0 (幼兒園)
- 金幣: 8145.6
- 紅花: 1
- 帖子: 719
- 在線: 72.3小時
- 蟲號: 2538185
- 注冊: 2013-07-09
- 性別: GG
- 專業(yè): 工程地質(zhì)學
6樓2014-08-16 16:09:07
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 307求調(diào)劑 +17 | 超級伊昂大王 2026-03-24 | 18/900 |
|
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學材料與化工（085600）296求調(diào)劑 +25 | 稻妻小編 2026-03-26 | 25/1250 |
|
|
|
[考研] 083000環(huán)境科學與工程調(diào)劑，總分281 +3 | 橙子（勝意） 2026-03-30 | 3/150 |
|
|
|
[考研] 346求調(diào)劑 一志愿070303有機化學 +7 | 蘿卜燉青菜 2026-03-28 | 8/400 |
|
|
|
[考研] 一志愿鄭州大學，080500學碩，總分317分求調(diào)劑 +10 | 舉個栗子oi 2026-03-24 | 11/550 |
|
|
|
[考研] 311求調(diào)劑 +7 | 藍月亮亮 2026-03-30 | 7/350 |
|
|
|
[考研] 070300化學354求調(diào)劑 +7 | 101次希望 2026-03-28 | 7/350 |
|
|
|
[考博] 26申博自薦 +6 | whh869393 2026-03-24 | 6/300 |
|
|
|
[考研] 0856材料化工調(diào)劑 總分330 +14 | zhubinhao 2026-03-27 | 14/700 |
|
|
|
[考研] 調(diào)劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
|
[考研] 071000生物學求調(diào)劑，初試成績343 +7 | 小小甜面團 2026-03-25 | 7/350 |
|
|
|
[考研] 312，生物學求調(diào)劑 +3 | 小譯同學abc 2026-03-28 | 3/150 |
|
|
|
[考研] 322求調(diào)劑 +5 | 舊吢 2026-03-24 | 5/250 |
|
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
|
[考研] 安徽大學專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級，六級521，求調(diào)劑 +4 | 美味可樂雞翅 2026-03-26 | 4/200 |
|
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
|
[考研] 一志愿 南京郵電大學 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
|
[考研] 打過很多競賽，085406控制工程300分，求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|