k空间布点与积分.doc

Outline I k 空间布点与 积分 (一) 北京市计算中 心 云平台 姜骏 ?k 空间布点与 物理量 1 ?k 空间布点与物理量 特殊点布点与 积分方法 2 特殊点布点与积分方法 北京市计算中心 云平台 姜骏 k 空间布点与积分 (一) Outline k 空间布点与 积分 (一) 北京市计算中 心 云平台 姜骏 ?k 空间布点与 物理量 1 ?k 空间布点与物理量 特殊点布点与 积分方法 2 特殊点布点与积分方法 北京市计算中心 云平台 姜骏 k 空间布点与积分 (一) ?k 空间布点与 Fermi 面的确定 k 空间布点与 积分 (一) 北京市计算中 心 云平台 姜骏 ?k 空间布点与 物理量 特殊点布点与 积分方法 (a) Brillouin Zone of Cubic lattice (b) Band Structure of SrSnO 3 在固体能带理论中，能量色散关系 ε(?k) 表示能量在倒空间中分布， 其中量子数 ?k (晶体动量) 描述平移对称性 北京市计算中心 云平台 姜骏 k 空间布点与积分 (一) ?k 空间布点与 Fermi 面的确定 k 空间布点与 积分 (一) 北京市计算中 心 云平台 姜骏 ?k 空间布点与 物理量 特殊点布点与 积分方法 Figure: The relation between unfolded-Band and the Brillouin-zone. 周期体系的 Fermi 能级和 Fermi 面的确定: { 导体: 价电子在 Brillouin-zone 部分填充 半导体-绝缘体: 价电子在 Brillouin-zone 完全填充 北京市计算中心 云平台 姜骏 k 空间布点与积分 (一) ?k 空间布点与 Fermi 面的确定 k 空间布点与 积分 (一) 北京市计算中 心 云平台 姜骏 ?k 空间布点与 物理量 特殊点布点与 积分方法 Figure: The Fermi-surface of Li, Mg, Cu in the ?rst Brillouin-zone. Fermi 面的形状: 最高占据能带折叠到第一 Brillouin-zone 围成的区域 要确定 Fermi 面的精细结构，特别是对于金属和导体体系，必须在 整个 Brillouin-zone 取足够多的采样点 北京市计算中心 云平台 姜骏 k 空间布点与积分 (一) ?k 空间积分与物理量 k 空间布点与 积分 (一) 北京市计算中 心 云平台 姜骏 ?k 空间布点与 物理量 特殊点布点与 积分方法 与 Fermi 面的确定类似，周期体系中所有单粒子期望值可表示为整 个 Brillouin-zone 内占据态的矩阵元的积分 一般地，如果已知 Brillouin-zone 某点 ?k 的能带指标为 n 的波函数 本征态 Ψn(?k) 和本征值 ?n(?k)，算符 X 的期望值 ? X? 是矩阵元 Xn(?k) = ?Ψn(?k)|X|Ψn(?k)? 在倒空间全部占据能带的求和 n ∫ G ∑ 1 ?X? = √V d3kXn(?k)f(εn(?k)) VG 其中 VG 是第一 Brillouin-zone 体积，f(ε) 是占据分布函数 实际计算中，Brillouin-zone 的 ?k 点数是有限的 ?X? = ∑ Xn(?kj)w ?kj n j,n ?k 点数目决定了电子结构和物理量的的精度与计算量 北京市计算中心 云平台 姜骏 k 空间布点与积分 (一) ?k 空间布点方案 k 空间布点与 积分 (一) 北京市计算中 心 云平台 姜骏 1 简单分布函数 Fermi-Dirac 分布函数 ?k 空间布点与 物理量 特殊点布点与 积分方法 1 f(ε) = e(ε? μ)/kT + 1 μ)/kT + 1 其中 μ 是化学势，k 是 Boltzmann 常数，T 是 温度参数 Figure: The Fermi-Dirac Distribution. Gaussian 分布函数 f(ε) = 1 2 σ√2π e? (ε?μ) 2σ2 其中 μ 是化学势，σ 是 展宽参数 Figure: The Gaussian Distribution. 北京市计算中心 云平台 姜骏 k 空间布点与积分 (一) ?k 空间布点方案 k 空间布点与 积分 (一) 北京市计算中 心 云平台 姜骏 ?k 空间布点与 物理量 2 特殊点方法 (Special-point scheme) 这是一种相对高效的积分方法，通过选取少量有代表性的 ?k 点，即可获得较高的计算精度，这些 ?k 点被称为“平均值 点”或“特殊点” 特殊点布点与 积分方法 特殊点方法对导体的收敛性较差 3 四面