Hubbard模型中格林函数的数值计算.pdfVIP

第34卷第1期 湖北大学学报( 自然科学版) V61．34 N0．1 20 12年3月 J our nal of Hube i Uni ver s i t y( Nat ur al Sci enc e) Ma r ．。2 012 文章编号 ：l o oO一2 3 5( 2 012) 0 1—0 048一04 Hubbar d模 型 中 格 林 函 数 的 数 值 计 算 朱家 昆，肖海波，舒崧 ( 湖北大学物理学与电子技术学院，湖北武汉430062) 摘要从Hubbar d模型出发，由动态平均场论方法结合量子蒙特卡洛模拟数值求解系统的格林函数，重 点讨论在顺磁序和反磁序情形下格林函数的数值结果，给出在不 同温度和不同Hubb ar d 能下格林函数的变 化情形，从格林函数的变化町以分析得到在温度降低或Hubbar d 能升高时，系统会发生从金属到绝缘体的 Mot t 转变 ． 关键词 动态平均场论；格林 函数 ；Hubba r d模型；量子蒙特卡洛方法 中图分类号0469 文献标志码A 高温氧化物超导体理论和实验的研究是当今材料科学研究领域的一个热 门方向，其所涉及到的物 理 内涵异常丰富．随着掺杂程度的变化，高温氧化物将展现出反铁磁性、半导导电行为、超导电性、强关 联电子系统以及新型金属行为等⋯ ．传统的BCS理论 已经无法解释高温氧化物 的超导电性 ，于是在高 温氧化物超导体被发现之后的20余年问，出现了众多的理论流派 【2 。3j ，尤以Ander s on 的非费米液体高 温超导理论最为著名 ，它认为高温超导体的特殊性质是由于组成材料的微观粒子之 间的强关联作用 引 起的，而高温超导体则是一个强关联体系L4j ．对于该强关联体系，可 以用Hubbar d模型来描述， Hubba r d模型在处理具体的强关联系统 当中计算比较复杂，要用到一些数值的方法 ．本文中主要采用动 态平均场论方法结合量子蒙特卡洛模拟( QMC) 研究Hubba r d模型的数值求解问题． 1 Hubbar d模型和动态平均场论 由于高温铜氧化物超导体的超导电性发生于绝缘体一金属相变 附近，Ander son认为可用一个近半 填充的单带Hubbar d模 型描述该系统L5j ： H一一∑￡。( ％+ ％+ ％+ ％) +u ∑疗i 十ni +一∥∑f 。+ ％ ( 1) ( 打) ．d l 白 其中行i —ci +f i 为粒子数算符，％+和。为原子i 中电子的产生与湮没算符，盯表示 自旋．￡一幻( 若 i √为最近邻) 对应于能带半宽度，U则是格位i 上 的库仑排斥能( 即Hu bbar d能) ，卢为化学势． Hubb ar d早期T作提出了一个非常接近于Mot t 初始思想的相变构想，尝试着给出了一个关联系 统的能带描述，以氢原子为例，设想有 ．N个氢原子逐渐靠近而排列成 晶格。当每个原子只有一个电子 时，相当于中性氢 原子态Ho ，当每个原子有两个电子 时，电子之间的库仑排斥作用，使它们之间有正的 Hubbar d能U，这时相当于氢的负离子态H一．若f 。表示第一个电子的能量，f 。+ 【，表示第二个电子的 能量 ．但当氢原子 间相互靠近时，能级拓展为能带，分别称为下Hubb ar d带和上Hubb ar d带．当相邻原 子 电子波函数重叠很小时，能带宽度很窄，上、下Hubbar d 带是分离的，下Hubbar d带是满带，上 Hubbar d带是空带，呈现出绝缘体性质．当原子逐渐靠近，上、下Hubbar d带发生交叠，而且都变成部分 填充的能带，呈现出金属电导的性质，这种由上 、下Hubbar d带引起的金属一绝缘体转变，称为Mot t 转 变L 61． 收稿 日期 ：20 1l —01一 14 基金项 目：国家 自然科学基金( 109050 18) 资助 作者简介：朱家昆( 19 9 一) ．男。硕士，讲师 第l 期 朱家 昆等：Hubb ar