版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

秋之愚

新蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 8145.6
紅花: 1
帖子: 719
在線: 72.3小時(shí)
蟲(chóng)號(hào): 2538185
注冊(cè): 2013-07-09
性別: GG
專業(yè): 工程地質(zhì)學(xué)

[求助] 徑向分布函數(shù)怎么畫(huà) 已有1人參與

求助各位大神，我在做Lennard-Jones勢(shì)流體的分子動(dòng)力學(xué)模擬時(shí)，能量圖基本上已經(jīng)畫(huà)出來(lái)了，可徑向分布函數(shù)的計(jì)算怎么弄都弄不出來(lái)，能有人幫幫我嗎？在給出來(lái)的這個(gè)程序中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

回復(fù)此樓

» 猜你喜歡

一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑已經(jīng)有6人回復(fù)
268求調(diào)劑已經(jīng)有7人回復(fù)
302求調(diào)劑已經(jīng)有11人回復(fù)
306求調(diào)劑已經(jīng)有4人回復(fù)
材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) 已經(jīng)有4人回復(fù)
求調(diào)劑已經(jīng)有6人回復(fù)
313求調(diào)劑已經(jīng)有3人回復(fù)
南昌大學(xué)材料專碩311分求調(diào)劑已經(jīng)有6人回復(fù)
316求調(diào)劑已經(jīng)有6人回復(fù)
346求調(diào)劑[0856] 已經(jīng)有7人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

怎樣用forcite畫(huà)徑向分布函數(shù) 已經(jīng)有14人回復(fù)
關(guān)于徑向分布函數(shù)和配位數(shù)公式中的密度已經(jīng)有3人回復(fù)
徑向分布函數(shù)與團(tuán)簇結(jié)構(gòu)關(guān)系已經(jīng)有5人回復(fù)
如何解釋徑向分布函數(shù)與r之間的關(guān)系？已經(jīng)有8人回復(fù)
徑向分布函數(shù)的問(wèn)題已經(jīng)有13人回復(fù)
vmd-做徑向分布函數(shù) 已經(jīng)有6人回復(fù)
forcite徑向分布函數(shù)g(r)物理意義已經(jīng)有8人回復(fù)
請(qǐng)教　徑向分布函數(shù)（ＲＤＦ）　問(wèn)題已經(jīng)有12人回復(fù)
如何做徑向分布函數(shù)和長(zhǎng)程序LRO？已經(jīng)有4人回復(fù)
請(qǐng)問(wèn)大家做徑向分布函數(shù)的時(shí)候結(jié)果有沒(méi)有趨近于1��？已經(jīng)有7人回復(fù)
【求助】怎么從結(jié)果中提取徑向分布函數(shù)、鍵長(zhǎng)、鍵角已經(jīng)有17人回復(fù)
【討論】陳正隆書(shū)上徑向分布函數(shù)的公式是否有誤？已經(jīng)有10人回復(fù)
【求助】應(yīng)用MS計(jì)算徑向分布函數(shù) 已經(jīng)有15人回復(fù)
【求助成功】徑向分布函數(shù)如何分析？已經(jīng)有14人回復(fù)

快樂(lè)學(xué)習(xí)，學(xué)習(xí)快樂(lè)

1樓 2014-08-13 21:45:27

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

秋之愚

新蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 8145.6
紅花: 1
帖子: 719
在線: 72.3小時(shí)
蟲(chóng)號(hào): 2538185
注冊(cè): 2013-07-09
性別: GG
專業(yè): 工程地質(zhì)學(xué)

怎么沒(méi)人回我啊

回復(fù)此樓

快樂(lè)學(xué)習(xí)，學(xué)習(xí)快樂(lè)

2樓2014-08-13 22:42:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dummy1

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 13 (小學(xué)生)
金幣: 5841.7
散金: 32
紅花: 1
帖子: 312
在線: 295.8小時(shí)
蟲(chóng)號(hào): 270440
注冊(cè): 2006-08-06
專業(yè): 理論和計(jì)算化學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
秋之愚: 金幣+5, ★★★很有幫助 2014-08-14 19:01:36

It is easy to get the RDF with VMD if you have the trajectory of your system in XYZ format.

贊一下(1人)

回復(fù)此樓

3樓2014-08-14 02:45:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

秋之愚

新蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 8145.6
紅花: 1
帖子: 719
在線: 72.3小時(shí)
蟲(chóng)號(hào): 2538185
注冊(cè): 2013-07-09
性別: GG
專業(yè): 工程地質(zhì)學(xué)

引用回帖:

3樓: Originally posted by dummy1 at 2014-08-14 02:45:40
It is easy to get the RDF with VMD if you have the trajectory of your system in XYZ format.

恩，謝謝你，這個(gè)我已經(jīng)看懂了，順便問(wèn)一下LJ粒子系統(tǒng)的速度相關(guān)函數(shù)該怎么理解。。。

贊一下

回復(fù)此樓

快樂(lè)學(xué)習(xí)，學(xué)習(xí)快樂(lè)

4樓2014-08-14 19:01:25

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dummy1

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 13 (小學(xué)生)
金幣: 5841.7
散金: 32
紅花: 1
帖子: 312
在線: 295.8小時(shí)
蟲(chóng)號(hào): 270440
注冊(cè): 2006-08-06
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

4樓: Originally posted by 秋之愚 at 2014-08-14 19:01:25
恩，謝謝你，這個(gè)我已經(jīng)看懂了，順便問(wèn)一下LJ粒子系統(tǒng)的速度相關(guān)函數(shù)該怎么理解。。。...

速度相關(guān)函數(shù)其實(shí)就是自相關(guān)函數(shù)，他是描述t時(shí)刻的和零時(shí)刻的速度的相關(guān)性，一般長(zhǎng)時(shí)間以后都是趨于零。通常用來(lái)計(jì)算體系的光譜或者擴(kuò)散系數(shù)。

贊一下

回復(fù)此樓

5樓2014-08-16 02:33:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

秋之愚

新蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 8145.6
紅花: 1
帖子: 719
在線: 72.3小時(shí)
蟲(chóng)號(hào): 2538185
注冊(cè): 2013-07-09
性別: GG
專業(yè): 工程地質(zhì)學(xué)

引用回帖:

5樓: Originally posted by dummy1 at 2014-08-16 02:33:13
速度相關(guān)函數(shù)其實(shí)就是自相關(guān)函數(shù)，他是描述t時(shí)刻的和零時(shí)刻的速度的相關(guān)性，一般長(zhǎng)時(shí)間以后都是趨于零。通常用來(lái)計(jì)算體系的光譜或者擴(kuò)散系數(shù)。...

恩，謝謝啦，現(xiàn)在對(duì)這些概念有些基本的了解了，能留個(gè)聯(lián)系方式嗎，我這才剛剛開(kāi)始，以后肯定還有很多問(wèn)題會(huì)請(qǐng)教你的.......

贊一下

回復(fù)此樓

快樂(lè)學(xué)習(xí)，學(xué)習(xí)快樂(lè)

6樓2014-08-16 16:09:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主秋之愚的主題更新

返回列表

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 346求調(diào)劑[0856] +4	WayneLim327 2026-03-16	7/350	2026-03-21 04:02 by JourneyLucky
[考研] 初始318分求調(diào)劑（有工作經(jīng)驗(yàn)） +3	1911236844 2026-03-17	3/150	2026-03-21 02:33 by JourneyLucky
[考研] 求調(diào)劑 +3	Ma_xt 2026-03-17	3/150	2026-03-21 02:05 by JourneyLucky
[考研] 307求調(diào)劑 +10	冷笙123 2026-03-17	10/500	2026-03-21 01:54 by JourneyLucky
[考研] 271材料工程求調(diào)劑 +8	.6lL 2026-03-18	8/400	2026-03-21 00:58 by JourneyLucky
[考研] 317求調(diào)劑 +8	申子申申 2026-03-19	13/650	2026-03-21 00:09 by 劉國(guó)森
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑 +3	我想要讀研究生 2026-03-20	3/150	2026-03-20 20:13 by JourneyLucky
[考研] 求調(diào)劑 +3	eation27 2026-03-20	3/150	2026-03-20 19:32 by JourneyLucky
[考研] 材料與化工專碩調(diào)劑 +7	heming3743 2026-03-16	7/350	2026-03-20 19:31 by zhukairuo
[考研] 環(huán)境工程調(diào)劑 +9	大可digkids 2026-03-16	9/450	2026-03-20 17:38 by 醉在風(fēng)里
[考研] 招收調(diào)劑碩士 +4	lidianxing 2026-03-19	12/600	2026-03-20 12:25 by lidianxing
[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 085600材料與化工求調(diào)劑 +6	緒幸與子 2026-03-17	6/300	2026-03-19 13:27 by houyaoxu
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無(wú)懈可擊111
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper