- 1、本文档共84页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
《分子动力学》PPT课件
蛙跳格式 计算机分子模拟方法 第三章、分子动力学方法 蛙跳格式: 计算机分子模拟方法 第三章、分子动力学方法 蛙跳格式: 计算机分子模拟方法 第三章、分子动力学方法 蛙跳格式: 热科学中的计算机分子模拟方法 第三章、分子动力学方法 初始位型:面心立方FCC NC=(FLOAT(N)/4.)**(1./3.) XL=1/(FLOAT(NC) Y=.5*XL R(1)=0,0,0 R(2)=0,Y,Y R(3)=Y,0,Y R(4)=Y,Y,0 M=0 DO 10 I=1,NC DO 10 J=1,NC DO 10 K=1,NC DO 11 IJ=1,4 RX(IJ+M)=RX(IJ)+XL*(K-1) RY(IJ+M)=RY(IJ)+XL*(J-1) RZ(IJ+M)=RZ(IJ)+XL*(I-1) 11 CONTINUE M=M+4 10 CONTINUE 计算机分子模拟方法 第三章、分子动力学方法 XL NC=4 Y Z X 程序F3(N,E,V): (2)量纲转换 SIGMA=(DENS/REAL(N))**(1.0/3.0) RCUT=RCUT*SIGMA DT=DT*SIGMA DENS=DENS/(SIGMA**3) 计算机分子模拟方法 第三章、分子动力学方法 程序F3(N,E,V): (3)长程校正: SR3=(SIGMA/RCUT)**3 SR9=SR3**3 SIGCUB=SIGMA**3 VLRC=(8.0/9.0)*PI*DENS*SIGCUB*REAL(N) *(SR9-3.0*SR3) WLRC=(16.0/9.0)*PI*DENS*SIGCUB*REAL(N) *(2.0*SR9-3.0*SR3) 计算机分子模拟方法 第三章、分子动力学方法 程序F3(N,E,V): (4)算法 算法启动: CALL FORCE(-DT,SIGMA,RCUT,NEWV,NEWVC,NEWW) CALL MOVE(-DT) CALL FORCE(-DT,SIGMA,RCUT,V,VC,W) CALL FORCE(DT,SIGMA,RCUT,V,VC,W) CALL KINET(OLDK) CALL MOVE(DT) CALL FORCE(DT,SIGMA,RCUT,NEWV,NEWVC,NEWW) CALL KINET(NEWK) 计算机分子模拟方法 第三章、分子动力学方法 程序F3(N,E,V): (4)算法 差分格式: SR2=SIGSQ/RIJSQ VIJ=4.0*(SR12-SR6) WIJ=48.0*(SR12-0.5*SR6) VELIJ=WIJ*DT/RIJSQ DVX=VELIJ*RXIJ … … DXI=VXI+DVX … … VX(J)=VX(J)-DVX … … V=V+VIJ W=W+WIJ CALL MOVE(DT) 计算机分子模拟方法 第三章、分子动力学方法 程序F3(N,E,V): (4)算法 最小映像: DO 1000 I=1,N-1 DO 900 J=I+1,N RXIJ=RXI-RX(J) RXIJ=RXIJ-ANINT(RXIJ) IF(ABS(RXIJ).LT.RCUT)THEN RYIJ=RYI-RY(J) RYIJ=RYIJ-ANINT(RYIJ) RIJSQ=RXIJ**2+RYIJ**2 IF(RIJSQ.LT.RCUTSQ)THEN RZIJ=RZI-RZ(J) RZIJ=RZIJ-ANINT(RZIJ) RIJSQ=RIJSQ+RZIJ**2 计算机分子模拟方法 第三章、分子动力学方法 0.6 0.4 RXIJ=RXIJ-ANINT(RXIJ) Example: -0.4=0.6-1 0.4=0.4-0 程序F3(N,E,V): (5)物理量的计算: 动能和势能: K=(OLDK+NEWK)/2+(OLDV-2*V+NEWV)/8 WHER
您可能关注的文档
- 《内科科室会议》ppt课件.pptx
- 《内科护理学课件》PPT课件.ppt
- 《内经》气病学说与糖尿病的分期辨治ppt课件.ppt
- 《内经》“十二刺”ppt课件.pptx
- 《内经选读痹论》ppt课件.ppt
- 《内经选读笔记》ppt课件.ppt
- 《内经选读绪论》PPT课件.ppt
- 《内衣销售培训》PPT课件.ppt
- 《内衣销售技巧培训》PPT课件.ppt
- 《内踝骨折》ppt课件.ppt
- 国际标准 IEC 61158-5-17:2007 EN Industrial communication networks - Fieldbus specifications - Part 5-17: Application layer service definition - Type 17 elements 工业通信网络 - 现场总线规范 - 第5-17部分:应用层服务定义 - 类型17元素.pdf
- 国际标准 IEC 60939-2-2:2004 EN_D 完整抑制无线电干扰滤波器单元-第2-2部分:空白详细规范-用于抑制电磁干扰的被动滤波器单元-仅需进行安全测试的滤波器(安全测试) Complete filter units for radio interference suppression - Part 2-2: Blank detail specification - Passive filter uits for electromagnetic interference .pdf
- 国际标准 IEC 60939-2-2:2004 EN_D Complete filter units for radio interference suppression - Part 2-2: Blank detail specification - Passive filter uits for electromagnetic interference suppression - Filters for which safety tests are required (safety.pdf
- 国际标准 IEC 60670-23:2006 EN-FR Boxes and enclosures for electrical accessories for household and similar fixed electrical installations - Part 23: Particular requirements for floor boxes and enclosures 家用和类似固定电气安装用电器配件的盒子与封闭装置——第23部分:地板箱和封闭装置的要求.pdf
- 国际标准 IEC 60670-23:2006 EN-FR 家用和类似固定电气安装用电器配件的盒子与封闭装置——第23部分:地板箱和封闭装置的要求 Boxes and enclosures for electrical accessories for household and similar fixed electrical installations - Part 23: Particular requirements for floor boxes and enclosures.pdf
- 2024版完整的货物运输合同书.doc
- 2024版无财产的离婚协议书书模板.doc
- 2024版委托融资租赁合同书书.doc
- 2024版汪怡的离婚协议书书.doc
- 2024版委托活动代理服务协议书书.doc
文档评论(0)