24小时热门版块排行榜

返回列表

解子林

新虫 (小有名气)

应助: 2 (幼儿园)
金币: 810.1
散金: 21
帖子: 67
在线: 31.8小时
虫号: 1697306
注册: 2012-03-17
专业: 有机合成

[求助] structure factor

哪位能帮忙分析分析我计算出来的结构因子为什么在大q时不趋近于1？我调整过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       !!!!!!!!水分子的坐标
            xx(mm) = X_o
            yy(mm) = Y_o
            zz(mm) = Z_o
                              mm=mm+1                      !!!!!!!水分子的数目
                              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

回复此楼