求解Hubbard线性系统快速稳定算法.pdfVIP

2014年 2月 计 算 数 学 第 36卷第 1期 Feb．，2014 MATHEMATICA NUMERICA SINICA Vo1．36，No．1 求解 Hubbard线性系统的快速稳定算法木1) 闵文超 黎 稳 柯日焕 (华南师范大学数学科学学院，广州510631) 摘 要 求解Hubbard线性系统是材料物理中DQMC／HQMC模型的核心问题之一，本文讨论了 Hubbard矩阵的结构，然后给出了当能量参数 U=0的情况下快速稳定求解 Hubbard线性系统 的算法．数值实验说明了方法的有效性． 关键词：DQMC模型；HQMC模型；Hubbard矩阵；Hubbard线性系统；线性系统数值解 的快速稳定算法 MR (2000)主题分类：65F10 1．引 言 作为量子蒙特卡罗数值仿真的重要应用，DQMC／HQMC模型在材料科学里扮演着重要 的角色，能够探究出固体材料中电子通过怎样的相互作用产生诸如磁性、超导性、导电性等物 理性质 [14J．但介于模型求解的难度大，一定程度上限制了该模型的实际应用价值，只能分析 拥有几百个电子的小型固体系统．然而即使在这样的限制下，该模型仍能准确地预测出固体 系统的大量物理性质．由此可见，如果突破模型求解的难关，DQMC／HQMC模型将大有所为． 为了简化这两个模型，我们首先引入Hubbard矩阵M： B1 ∈CⅣL 肌 ， (1．1) 其中 ·I是一个N ×N阶的单位矩阵； ·Bl是一个N ×N 阶的如下形式的矩阵 岛 =BDI， B ： et△rK D e ” (h¨ ． ， (1．2) t代表电子跃迁的参数； 2013年 1月4日收到． )基金项 目：国家自然科学基金项 目11271144)、广东省自然科学基金项 目(s2012010009985， s201301o0112530)、教育部博士点基金博导类研究项 目(20104407110001)和华南师范大学 2013年研究生创 新基金项 目(2013kyjj047)资助． 2 计 算 数 学 2014年 △丁= ／， 代表逆温度， 代表将逆温度分段的数 目； 盯=+或者 一； 代表与能量参数U有关的参数，在DQMC模型中，=cosh～e ，在HQMC模 型中，u=(2)△丁； (c)是一个N XN 阶的对角矩阵， V~(h1)=diag(hll，hi2 ．，hlN) 其中 “是随机变量，通常被称为Hubbard—Stratonovish域， K=(kij)用来描述处于第 i个和第J个电子层上电子跃迁情况的邻接矩阵 ]gij-~- {。1 时