马尔可夫链理论和Monte Carlo 取样及实现.pptVIP

马尔可夫链理论和Monte Carlo 取样的实现 Monte Carlo 取样不直接使用 P(X) ,而是以某种方式取样, 大量的样本最终符合所需的分布. Monte Carlo利用转移矩阵 , 从当前 x 生成一个状态 y . 马尔可夫链和Monte Carlo 马尔可夫链是一个简单的随机过程, 给定转移矩阵 W,可以得到平衡分布 P. Monte Carlo 是马尔可夫链的计算机实现. Monte Carlo中, P 给定, 我们需要寻找 W 使得 P = P W. 马尔可夫链的收敛条件(充要) 各态历经 对所有 n nmax, 所有 x 和 x’ 细致平衡 Metropolis 算法 (1953) Metropolis 算法取 T 为一个对称的转移矩阵 The Paper (7500 citations from 1988 to 2003) The Calculation Number of particles N = 224 Monte Carlo sweep ≈ 60 Each sweep took 3 minutes on MANIAC Each data point took 5 hours MANIAC the Computer and the Man The Ising Model 感兴趣的物理量 平均能量 E(s) 比热，利用下面的公式: 3. 磁化 M = |∑isi| 磁化率 5. Binder 4 阶累积量 自旋相关函数 时间相关函数 2D Ising 模型的比热 有限尺寸标度与相变 有限尺寸标度与相变 有限尺寸标度与相变 有限尺寸标度与相变 有限尺寸标度与相变 有限尺寸标度与相变 有限尺寸标度与相变 有限尺寸标度与相变 有限尺寸标度与相变 有限尺寸标度与相变 有限尺寸标度与相变 有限尺寸标度与相变 有限尺寸标度与相变 有限尺寸标度与相变 有限尺寸标度与相变 有限尺寸标度与相变 Finite-Size Scaling Shift of Tc Accurate Exponent Ratio 有限尺寸标度与相变 有限尺寸标度与相变 有限尺寸标度与相变 有限尺寸标度与相变 有限尺寸标度与相变 有限尺寸标度与相变 有限尺寸标度与相变 磁化随MCS的变化，有限系统的各态历经演示。2D Ising 模型，L＝64 相变点会发生什么? 在二级相变点(现一般称为临界点, 二级相变也称为连续相变), 系统的比热, 磁化率等物理量发散, 发散的根源在于相关长度发散. 这样, 前面的分析将不再成立. 在临界点, 热力学量具有很有趣的标度行为: 以单轴各向异性磁体为例(理论上近似以Ising模型描述). 在T! Tc 时, 引进约化温度: 比热: 磁化率: 自发磁化: 关联长度: 当t=0时, 磁化与外场之间有: ?, ?, ?, ?, ? 称为临界指数, 是反映临界点本质的量. 这些指数的来源与关联长度的发散有关. 还有一个指数?, 这里不做介绍. 当系统有限时, 不会有发散! Monte Carlo只能计算有限体系, 如何从Monte Carlo计算获得相变和临界现象的信息? 如何找临界点? 如何算临界指数指数? 为此, 我们考虑有限体系的标度理论. 基本出发点: 设有限系统的线度为L, 则关联长度最多为L, ? L; 发散的物理量不再发散, 而是宽度为?T 的峰. TTC 时的序参量分布 序参量的有限尺寸效应 磁化率的有限尺寸效应 T? TC 时的序参量分布 TTC 时的序参量分布 考虑磁化率: 在临界点附近有一个峰, 最高处的温度定义为有限系统的临界温度Tc(L), 峰宽为 ?T. 当 L ! 1 时, Tc(L)! Tc(1), ?T! 0 假定: 磁化率的极大值: 磁化在Tc(L)不为0, 而是: Fisher 假定: 在临界点只有一个重要的特征长度, 即关联长度. 由于 , 磁化率峰的宽度对应于 , 于是 由此得到: From A M Ferrenberg and D P Landau, Phys Rev B 44 (1991) 5081 From J S Wang, R H Swendsen, and R Kotecky, Phys Rev. B, 42 (1990) 2465. Binder 累积量: 性质1: T Tc, L?, 远离临界区域 性质2: T Tc,