| 查看: 418 | 回复: 0 | ||
[求助]
结构因子 分子模拟
|
|
哪位能帮忙分析分析我计算出来的结构因子为什么在大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 |
回复此楼» 猜你喜欢
什么时候开奖？
已经有4人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有27人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有3人回复
-
2027广东省杰青
已经有9人回复
-
朋友圈看到的
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有110人回复
-
明天放榜？
已经有5人回复
-
时间戳变了，能看出什么问题？
已经有11人回复
-
重要消息，中午系统在维护
已经有10人回复
-
投稿咨询
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
找到一些相关的精华帖子，希望有用哦~
转录因子结合位点
已经有3人回复
- 怎样用forcite画径向分布函数
已经有14人回复
- The_Art_of_Molecular_Dynamics_Simulation
已经有32人回复
- 配体pdb文件中的温度因子
已经有3人回复
- 关于DNA连接的纳米粒子自组装的分子动力学模拟，求交流和分享
已经有15人回复
- 太阳能电池中影响迁移率和短路电流密度的分子结构因素有哪些？希望回帖者附带相关文献
已经有4人回复
- 影响因子2以上生物杂志都要原始数据吗？
已经有10人回复
- Arrhenius公式求频率因子与活化能
已经有3人回复
- 2013年最新公布2012年度SCI杂志IF影响因子.xls
已经有55人回复
- 2013最新发布SCI影响因子（完整版）
已经有22人回复
- 高分子期刊影响因子1.0左右的期刊
已经有6人回复
- 径向分布函数与静态结构因子是什么关系啊？？
已经有4人回复
- 两道生化试题求解答。求答案
已经有3人回复
- 什么软件可以比较快速的模拟反应速率的指前因子和活化能
已经有6人回复
- 分子模拟的统计力学理论基础
已经有6人回复
- 【讨论】用VMD作的RDF图片不知有没有问题
已经有16人回复
- 【求助】fullprof 精修，指标化-欢迎参与讨论
已经有20人回复
- 【请教】X射线衍射确定化合物结构
已经有10人回复
- 【求助】S(Q)－－G(r)转换问题
已经有7人回复
- 【求助】structure factor S(q)的问题
已经有5人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
