Liouvillian 插值能修正团簇动力学平均场中非物理的自能吗?

Liouvillian 插值在构造上保持因果性,可对团簇 DMFT 自能进行后处理。

直接答案

按论文检验的范围,答案是肯定的:Liouvillian 插值不再插值依赖频率的函数,而是插值从团簇自能的连分式展开中得到的、不依赖频率的 Liouvillian 矩阵元,因此插值后的自能在构造上保持因果性。作者在一维 Hubbard 模型的 CDMFT 计算,以及二维 Hubbard 模型的四格点 CDMFT 和四 patch DCA 计算上演示了该方法;在后者中,它给出了看似连成闭合曲面的 Fermi 弧与 Luttinger 弧。

1篇文献引用

本文由 WisPaper 驱动的搜索和论文分析生成。

插值要消除的结构性缺陷

动力学平均场理论的团簇扩展能捕捉单格点 DMFT 遗漏的短程非局域关联,代价是自能带有特定类型的非物理性质。团簇动力学平均场理论(CDMFT)把一个实空间格点团簇与自洽浴耦合,这会破坏平移对称性,因此其实空间自能不满足平移不变性。动力学团簇近似(DCA)则把布里渊区粗粒化成若干 patch,这会在自能的动量依赖中引入人为的不连续。

这两个缺陷都不是物理结论,而是团簇嵌入方式造成的假象。标准补救办法是后处理:对 CDMFT 做周期化,对 DCA 做插值,论文把两者统称为插值。目标是从只定义在少数格点或 patch 上的团簇数据出发,重建出光滑且平移不变的格点物理量。

插值 G、Sigma 或 M 为何可能破坏因果性

已有方案分别插值格林函数 G(omega)、自能 Sigma(omega) 或累积量 M(omega)。这三者都是依赖频率的函数,论文指出它们都带有问题性质,有时包括违反因果性。原因在于,插值是一种跨越团簇或 patch 边界施加的光滑化操作,而这一操作本身并不保证物理响应函数必须具备的解析结构。

实际失效在空穴掺杂的方格子 Hubbard 模型中可见:传统方法难以可靠捕捉空穴口袋这类精细结构。因此问题不只是形式上的:本应修复团簇结果的方案,可能恰好扭曲读者最想分辨的动量空间特征。

改为插值不依赖频率的 Liouvillian 矩阵元

论文的提议是插值另一个对象。Liouvillian L 是 Fock 空间中所有线性算符构成的空间里时间平移的生成元,其矩阵元可由推迟自能的连分式展开(CFE)得到。作者插值的是 L 的单粒子不可约部分中不依赖频率的矩阵元,再把插值后的矩阵元代回自能的 CFE。

有两个性质使这些矩阵元比 G、Sigma 或 M 更适合插值。第一,它们不依赖频率,插值过程因此不会触及承载解析结构的频率轴。第二,论文证明 CFE 的组成成分比传统方案所用的函数具有更局域的傅里叶展开,也就是在实空间中通常具有更短的作用范围。由于因果结构由 CFE 承载,而不是由被插值的量承载,所得自能在构造上保持因果性;论文在第 11 页讨论了一个小例外。作者还指出,该方法在技术上并不比其它方案更费力,唯一额外的技术成分就是 CFE。

一维检验中随团簇尺寸更快收敛

第一项检验使用一维 Hubbard 模型的 CDMFT 计算,这是一个可以改变团簇尺寸的受控场景。论文报告,L 插值随团簇尺寸的收敛速度快于 Sigma、G 或 M 插值。这是对同一批团簇数据上不同方案收敛速度的比较,而不是关于一维模型物理的结论;它表明该方法比标准替代方案更高效地逼近大团簇极限。

二维赝能隙区中的 Fermi 弧与 Luttinger 弧

物理应用是对二维 Hubbard 模型赝能隙区做四 patch DCA 计算,用数值重整化群求解,同时用精确对角化做四格点 CDMFT 计算。论文的预览图把带有布里渊区不连续性的原始 DCA 结果,与基于 Sigma、M 和 L 的插值结果并列展示。论文指出,L 插值是最能令人信服地模拟原始不连续团簇结果的方案。

主要物理结果由此重建而来:取决于掺杂,L 插值可以给出 Fermi 弧与 Luttinger 弧,二者合起来看似构成一个闭合曲面,并且两者之间似乎存在光滑过渡。“看似”在这里有实质含义。闭合曲面和光滑过渡都是插值所得量的特征,论文把它们作为该方案的产物呈现,而不是作为模型被独立测量的性质。

该方法没有解决的问题

论文明确说明,其主要目标是方法层面的:详细描述 L 插值的概念和计算步骤,二维应用只是用来展示其能力。论文有意不在 CDMFT 与 DCA 之间作出裁决,并指出具体问题可能有理由偏好其中一种。作者建议,在两者计算成本相当且都可行的情况下,例如四格点 CDMFT 对四 patch DCA,系统比较所得 Fermi 面性质会很有价值。

演示范围也有限。示例覆盖一维 Hubbard 模型的 CDMFT,以及二维 Hubbard 模型的四格点 CDMFT 和四 patch DCA;要处理 d 波超导与反铁磁所需的对称破缺自能,需要推广该表述,论文将其列为未来工作。这里的插值是 DMFT 自洽完成后施加的后处理步骤,而不是在团簇 DMFT 循环内部进行的操作。

Liouvillian 插值的适用判断

论文支持的判断范围窄而具体。当目标是从团簇 DMFT 数据中获得谱特征的良好动量分辨率,并且在插值 G、Sigma、M 还是 Liouvillian 之间做选择时,证据偏向 Liouvillian:它在构造上保持因果性,在一维 CDMFT 检验中随团簇尺寸收敛更快,并在四 patch DCA 的赝能隙计算中给出其它方案未能令人信服地重现的 Fermi 弧与 Luttinger 弧。这一选择的代价是多一个技术成分,即连分式展开。

因此,重点不是对 CDMFT 与 DCA 孰优孰劣下结论,作者对此保持开放;重点是团簇数据已经存在之后,应当插值哪个量。对读者最有分量的主张,是因果性保证与赝能隙区重建出的弧结构这两者的结合,因为旧方案正是在这里明显吃力。

本文引用的文献

1

Liouvillian interpolation of the self-energy of cluster dynamical mean-field theories

插值由团簇自能的连分式展开得到的、单粒子不可约 Liouvillian 中不依赖频率的矩阵元,可得到在构造上保持因果性的插值自能;在一维 CDMFT 中,它随团簇尺寸的收敛快于 Sigma、G 或 M 插值,并在二维 Hubbard 模型的四 patch DCA 计算中给出 Fermi 弧与 Luttinger 弧。