PIMC方法的应用.pptVIP

* * ——PIMC方法的应用 刘松芬 胡北来 南开大学物理科学院 高温、高压状态下的稠密等离子体热力学特性的研究 Ⅰ 背景介绍 近些年来，费米系统的热动力学性质在很多领域受到人们越来越大的关注。包括等离子体物理、天体物理、凝聚态和核物质。其中人们最感兴趣的是在库仑和量子效应关联条件下的费米液体、金属氢、等离子体相变和束缚态等等。从理论上讲, PIMC方法可以很好的描述这一领域。而且在费米系统中应用这种方法已经取得了明显的进步。但是，它存在着一个本质上的困难，即由于费米子置换的反对称性引起的费米符号问题。 1 N个全同粒子的体系的薛定谔方程 的解是： 即体系的本征能量等于各粒子本征能量之和，体系的本征函数等于各单粒子哈密顿算符的本征函数之积。这样，多粒子体系可以归结为解单粒子薛定谔方程。 对费米系统，体系的波函数可以写成如下的行列式形式： 展开行列式后，每项都具有 (1) 式的形式，因而它是体系薛定谔方程的解。交换任何两粒子，也就是在行列式中两列互相调换，这使行列式改变符号。所以费米体系的波函数是反对称的。即存在N!/2项正的和N!/2项负的，这就是费米符号问题。 我们所感兴趣的是力学量算符 在体系中的平均值，即该力学量算符在系统中的期待值，采用密度（或统计）矩阵的方法，是最有效、最简单的方法。 2 假设系统由一组归一化的纯粹态 表示，相应于 的归一化权为 。再设 可按某一正交归一化的完备基矢集 展开，即 则力学量算符 在纯粹态 中的期待值为： 力学量 的统计平均值为 3 引入密度矩阵 则 的统计平均值可以简单的表示为 因此 的总平均值可由密度矩阵计算出来，而密度矩阵又可通过组成混和态的纯粹态 及相应的统计权重表达出来。 的对角元具有简单而直接的物理意义，即是基矢 在混和态中出现的几率。 同时，如果所考虑的算符 具有非负的本征值，则其平均值 必定也是非负的。因此， 必定为正值。 4 ★ B.Militzer等人提出的 Restrict Path integral Monte Carlo （RPIMC）[1,2] ★ V.S.Filinov 、M.Bonitz、V.Ebeling、V.E.Fortov 等人提出的 direct fermionic PIMC （DPIMC）[3-5] 为了解决这一问题，人们采用了不同的方法，其中以下面两个组最具代表： 费米系统中，由于置换的反对称性，可能会造成密度矩阵为负的情况，这也即为费米符号问题。 5 ★ B.Militzer提出的RPIMC方法，其中引入了节面的概念，通过限制积分的路径使密度矩阵的贡献始终为正，但它只适用于某些特殊的情况。而我们也只知道少数几种密度矩阵的节面，并且这些节面很近似。 ★ V.S.Filinov 、M.Bonitz、V.Ebeling、V.E.Fortov 等人提出的DPIMC，是把N粒子系统的密度算符用新的路径积分表示，而粒子置换的反对称效应则完全包括在一个矩阵里，从而克服了费米符号问题，但V.S.Filinov等人提出的模型只适用于弱耦合及弱简并区域。 在此基础上，我们以 DPIMC 为起点，通过改变其物理模型, 有效的计算了强简并耦合区域的氢等离子体的一些热力学性质。 6 Ⅱ :数学模型 以电子 和质子 组成的稠密氢等离子体的密度矩阵在量子系统中，一般不存在确切的密度矩阵。但密度矩阵遵循简单的组合规律。 重复上面的步骤插入n组完备集，当n—∞时，即为高温近似下的密度矩阵，且这时的密度矩阵等同于没有近似条件下的密度矩阵。 7 对每一个小密度矩阵而言,在高温近似下,n→∞可以近似为 这里的U(q)为粒子间的相互作用能。 8 因此，一般情况下，量子系统的密度矩阵可由高温近似下的路径积分来表示。考虑费米子的自旋效应，则密度矩阵可以写成： 其中 是系统的哈密顿算符 且 为粒子间的相互作用 这意味着粒子由 坐标的费米循环所表示。 其中q和r表示质子和电子的坐标。自旋作用由自旋矩阵S来表示，而置换作用由置换算符 来表示。它作用在电子坐标和自旋投影上,由 的奇偶性求和。费米系统中,求和包括 项正的和 项负的。这就是费米符号问题。由于质子的质量远大于电子，所以它的置换作用可以不考虑。 9 10 对每一个小密度矩阵,由于高温近似下, 可以近似为 is the free-particle d