[论文解读] Physical insights from imaginary-time density--density correlation functions

Tobias Dornheim, Zhandos A. Moldabekov+3 more|arXiv (Cornell University)|Sep 6, 2022

Advanced Chemical Physics Studies被引用 5

一句话总结

本文证明了虚时密度-密度关联函数 $ F(\boldsymbol{q},\tau) $ 包含与动态结构因子 $ S(\boldsymbol{q},\tilde{\tau}) $ 相同的物理信息，使得无需解析延拓即可直接提取温度和准粒子激发能等关键性质。基于均匀电子气的从头算路径积分蒙特卡罗数据，作者表明非平凡特征如旋子模式在 $ F(\boldsymbol{q},\tau) $ 中清晰显现，并建立了 $ S(\boldsymbol{q},\tilde{\tau}) $ 的频率矩与 $ F(\boldsymbol{q},\tau) $ 在 $ \tau = 0 $ 附近的曲率之间的直接联系，为实频率分析提供了一种稳健的替代方法。

ABSTRACT

The accurate theoretical description of the dynamic properties of correlated quantum many-body systems such as the dynamic structure factor $S(\mathbf{q},ω)$ constitutes an important task in many fields. Unfortunately, highly accurate quantum Monte Carlo methods are usually restricted to the imaginary time domain, and the analytic continuation of the imaginary time density--density correlation function $F(\mathbf{q},τ)$ to real frequencies is a notoriously hard problem. In this work, we argue that no such analytic continuation is required as $F(\mathbf{q},τ)$ contains, by definition, the same physical information as $S(\mathbf{q},ω)$, only in an unfamiliar representation. Specifically, we show how we can directly extract key information such as the temperature or quasi-particle excitation energies from the $τ$-domain, which is highly relevant for equation-of-state measurements of matter under extreme conditions. As a practical example, we consider \emph{ab initio} path integral Monte Carlo results for the uniform electron gas (UEG), and demonstrate that even nontrivial processes such as the \emph{roton feature} of the UEG at low density straightforwardly manifest in $F(\mathbf{q},τ)$. In fact, directly working in the $τ$-domain is advantageous for many reasons and holds the enticing promise for unprecedented agreement between theory and experiment.

研究动机与目标

• 证明虚时密度-密度关联函数 $ F(\boldsymbol{q},\tau) $ 所含物理信息与动态结构因子 $ S(\boldsymbol{q},\tilde{\tau}) $ 相同，但以更易处理的形式呈现。
• 表明温度和准粒子激发能等关键物理量可直接从 $ F(\boldsymbol{q},\tau) $ 提取，而无需转换至实频率。
• 建立 $ S(\boldsymbol{q},\tilde{\tau}) $ 的频率矩与 $ F(\boldsymbol{q},\tau) $ 在 $ \tau = 0 $ 处曲率之间的定量关系，从而提供新的诊断工具。
• 通过均匀电子气的从头算路径积分蒙特卡罗模拟验证该方法，尤其在低密度下准确捕捉旋子特征。
• 倡导将 $ F(\boldsymbol{q},\tau) $ 作为状态方程和热密物质诊断中实频率响应函数的实用且精确的替代方案。

提出的方法

• 作者推导了虚时涨落-耗散定理，表明静态密度响应函数 $ \chi(\boldsymbol{q},0) $ 与 $ F(\boldsymbol{q},\tau) $ 在虚时区间 $ \tau \in [0,\beta] $ 上的积分直接相关，关系式为 $ \chi(\boldsymbol{q},0) = -n \int_0^\beta F(\boldsymbol{q},\tau) d\tau $。
• 他们建立了 $ S(\boldsymbol{q},\omega) $ 的频率矩与 $ F(\boldsymbol{q},\tau) $ 在 $ \tau = 0 $ 附近的泰勒展开系数之间的数学联系，特别是二阶导数（曲率），其编码了第二阶频率矩。
• 该方法依赖于从头算路径积分蒙特卡罗（PIMC）模拟，用于生成有限温度下不同密度的均匀电子气的 $ F(\boldsymbol{q},\tau) $。
• 作者分析了 $ F(\boldsymbol{q},\tau) $ 的函数形式，特别是其曲率和衰减行为，以推断如动态结构因子中的旋子最小值等物理特征。
• 他们比较了不同波矢和温度下 $ F(\boldsymbol{q},\tau) $ 的行为，表明非平凡特征如旋子在 $ \tau $-域中可直接观测。
• 该方法通过直接在 $ \tau $-表示中工作，避免了病态的解析延拓问题，充分利用了量子蒙特卡罗在虚时中的稳定性和准确性。

实验结果

研究问题

• RQ1 是否可直接从虚时密度-密度关联函数 $ F(\boldsymbol{q},\tau) $ 提取温度和准粒子激发能等物理信息，而无需解析延拓？
• RQ2 $ S(\boldsymbol{q},\omega) $ 的频率矩与 $ F(\boldsymbol{q},\tau) $ 在 $ \tau = 0 $ 处的曲率之间有何关系？
• RQ3 在不转换至实频率的情况下，均匀电子气中的非平凡特征（如旋子模式）是否可在 $ F(\boldsymbol{q},\tau) $ 中清晰识别？
• RQ4 $ F(\boldsymbol{q},\tau) $ 在热密物质和状态方程测量的诊断中，能在多大程度上作为 $ S(\boldsymbol{q},\omega) $ 的可靠且实用的替代方案？
• RQ5 高阶虚时关联函数在捕捉非线性响应特性及三体、四体关联方面起到什么作用？

主要发现

• 虚时密度-密度关联函数 $ F(\boldsymbol{q},\tau) $ 包含与动态结构因子 $ S(\boldsymbol{q},\omega) $ 相同的物理信息，使得许多物理洞察无需解析延拓即可获得。
• $ F(\boldsymbol{q},\tau) $ 在 $ \tau = 0 $ 处的曲率与 $ S(\boldsymbol{q},\omega) $ 的第二阶频率矩直接相关，建立了 $ \tau $-域行为与实频率谱矩之间的定量联系。
• 在低密度下，均匀电子气中的旋子特征在 $ F(\boldsymbol{q},\tau) $ 中清晰可见，表明非平凡谱特征可直接在 $ \tau $-表示中显现。
• 静态密度响应函数 $ \chi(\boldsymbol{q},0) $ 可直接通过 $ F(\boldsymbol{q},\tau) $ 在 $ \tau \in [0,\beta] $ 上的积分计算得出，验证了虚时涨落-耗散定理的有效性。
• 直接在 $ \tau $-域工作避免了解析延拓的数值不稳定性，实现了对电子-电子静态结构因子等频率积分性质的稳健提取。
• 通过PIMC模拟可获得的高阶虚时关联函数可用于探测非线性响应函数和多体关联，为量子多体系统诊断开辟了新途径。