| 1 | 1/1 | 返回列表 |
| 查看: 375 | 回復(fù): 0 | ||
解子林新蟲 (小有名氣)
|
[求助]
結(jié)構(gòu)因子 分子模擬
|
|
哪位能幫忙分析分析我計算出來的結(jié)構(gòu)因子為什么在大q時不趨近于1?我調(diào)整過qx、qy、dq的值,可是沒什么作用,大伙給點意見,謝謝! do i = 1 , conf_2 mm=1 do j = 1 , 9 read(1,*) enddo do j = 1 , all_n read(1,*) num_o , type_o , X_o , Y_o , Z_o if (type_o==1)then !!!!!!!!水分子的坐標(biāo) xx(mm) = X_o yy(mm) = Y_o zz(mm) = Z_o mm=mm+1 !!!!!!!水分子的數(shù)目 endif end do do nx=1,600 do ny=1,600 do nz=1,600 kr=sqrt((nx*qx)**2+(ny*qy)**2+(nz*qz)**2)!!!!!!!q的模 ! if ((kr>9).and.(kr<62))then bin=1+int(kr/0.5) nhist(bin)=nhist(bin)+1 cossum=0.0 sinsum=0.0 mm1=0 do l = 1 , mm-1 rx=0.1*xx(l) ry=0.1*yy(l) rz=0.1*zz(l) cossum=cossum+cos(nx*qx*rx+ny*qy*ry+nz*qz*rz) sinsum=sinsum+sin(nx*qx*rx+ny*qy*ry+nz*qz*rz) mm1=mm1+1 end do hist(bin) =hist(bin) + (cossum**2+sinsum**2)/real(mm-1) write(*,*) nhist(bin),hist(bin),mm1,mm,i !endif enddo enddo enddo end do !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! do k = 1 ,600 if(hist(k)/=0) then g(k) = real(hist(k))/real(nhist(k)) write(2,*) k*0.5,g(k) else g(k)=0 write(2,*) k*0.5,g(k) endif end do |
找到一些相關(guān)的精華帖子,希望有用哦~
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工 +4 | 安全上岸! 2026-03-16 | 4/200 |
|
|---|---|---|---|---|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
[考研] 293求調(diào)劑 +5 | 世界首富 2026-03-11 | 5/250 |
|
|
[考研] 274求調(diào)劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 286求調(diào)劑 +3 | lemonzzn 2026-03-16 | 5/250 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +3 | 張vvvv 2026-03-15 | 5/250 |
|
|
[考研] 機械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 一志愿華中師范071000,325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +7 | Ncdx123456 2026-03-13 | 8/400 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy! 2026-03-13 | 6/300 |
|
|
[考研] 材料工程,326分,求調(diào)劑 +6 | KRSLSR 2026-03-10 | 6/300 |
|
|
[考研] 336求調(diào)劑 +6 | Iuruoh 2026-03-11 | 6/300 |
|
|
[考研] (081700)化學(xué)工程與技術(shù)-298分求調(diào)劑 +12 | 11啦啦啦 2026-03-11 | 35/1750 |
|
|
[考研] 315求調(diào)劑 +9 | 小羊小羊_ 2026-03-11 | 10/500 |
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|