Mc04最新版-绪论及蒙特卡罗方法题库.ppt

* PROGRAM MC_gc DO icycl = 1, ncycl ran= RANF(Iseed)* (NPART+NEXC) + 1 IF (ran.LE.NPART)THEN CALL MCMOVE ELSE CALL MCEXC ENDIF IF (MOD(icycl,nsamp).EQ.0) CALL SAMPLE END DO END * SUBROUTINE MCEXCH IF (RANF(Iseed).LE.0.5D0) THEN o = int(NPART* RANF) + 1 CALL ENER arg = NPART *EXP(-BETA*eno)/( ZZ*vol) zz= EXP(-BETA*eno)/c IF (RANF(Iseed).LT.arg) THEN X(o) = X(NPART) NPART = NPART - 1 ENDIF ELSE XN = RANF*BOX CALL ENERI arg = (ZZ*vol )X*EXP(-BETA*enn)/((NPART+1) IF (RANF(Iseed).LT.arg) THEN X(NPART+1) =XN NPART = NPART +1 END IF END IF RETURN END * Variational quantum MC simulations Green’s function MC simulations Path-integral MC simulations Quantum lattice models Topics not Covered * Critical Slowing Down Cluster Algorithms - + + + + + + + + + + + + + + + + - - - - - - - - - - - - - - - - - - - * Other Applications Application in image processing Light transport in computer graphics Neutron transport Growth models (dynamic MC) Computational finance Computational biology/Bioinformatics Strongly correlated systems and quantum field theory 人有了知识，就会具备各种分析能力， 明辨是非的能力。 所以我们要勤恳读书，广泛阅读， 古人说“书中自有黄金屋。 ”通过阅读科技书籍，我们能丰富知识， 培养逻辑思维能力； 通过阅读文学作品，我们能提高文学鉴赏水平， 培养文学情趣； 通过阅读报刊，我们能增长见识，扩大自己的知识面。 有许多书籍还能培养我们的道德情操， 给我们巨大的精神力量， 鼓舞我们前进。 * * * * 格子模型的基本做法是将空间离散化，即链单元只能取空间中某些人为规定的格点。显然格子链在细节上与真实链有较大的区别，但高分子的许多统计性质并不依赖于链模型的细节。 1. Lattice models of Polymers (1) random walk 三维立方格子上的random walk链也常称为“六选择模型” 1 2 3 * （1）将第一个链节固定在坐标原点上，并设格子的边长为1。 以二维方格子为例：来说明格子模型的具体步骤 （2）产生（0，3）整数序列的随机数。 （3）由随机数的值按预先规定的规则，决定下一链节所达的格点。 （4）重复（2）、（3）直至所需的链长n，并记录Rn。 （5）重复（1）~（4）直至达到所需的分子数M。 * (2) self-avoiding walk （1）将第一个链节固定在坐标原点上，并设格子的边长为1 （2）产生（0，3）整数序列的随机数。 （3）由随机数的值按预先规定的规则，找到下一格点。如果此格点未被占据，则为下一链节所达的格点。否则转（2）。 （4）重复（2）、（3）直至所需的链长n，并记录Rn。 （5）重复（1）~（4）直至达到所需的分子数M。 * 抽样的成功率低，如何改进：偏倚抽样；跨大步法等。 其它格子链，多键模型 the atomistic models 2. ‘Continuous’ polymer models The bead models 1 2 3 * 统计力学的一个重要任务就是利用模型来计算体系的宏观平均性质，对polymer来说，均方末端距(mean square end-to-end distance)和均方回转半