[论文解读] Single-cone real-space finite difference schemes for the Dirac von Neumann equation

Magdalena Schreilechner, W. Pötz|arXiv (Cornell University)|Mar 6, 2015

Topological Materials and Phenomena参考文献 10被引用 3

一句话总结

本文提出两种用于求解(2+1)D狄拉克哈密顿量冯·诺依曼方程的单锥实空间有限差分格式，采用交错网格以消除费米子加倍问题。bra-ket格式通过分别对密度算符的ket和bra分量应用哈密顿量来保持正定性，而直接格式则一步计算密度矩阵的时间导数；两种格式均确保单锥色散关系且稳定，适用于拓扑绝缘体表面态，并可扩展至林德布拉德方程和格林函数形式。

ABSTRACT

Two finite difference schemes for the numerical treatment of the von Neumann equation for the (2+1)D Dirac Hamiltonian are presented. Both utilize a single-cone staggered space-time grid which ensures a single-cone energy dispersion to formulate a numerical treatment of the mixed-state dynamics within the von Neumann equation. The first scheme executes the time-derivative according to the product rule for "bra" and "ket" indices of the density operator. It therefore directly inherits all the favorable properties of the difference scheme for the pure-state Dirac equation and conserves positivity. The second scheme proposed here performs the time-derivative in one sweep. This direct scheme is investigated regarding stability and convergence. Both schemes are tested numerically for elementary simulations using parameters which pertain to topological insulator surface states. Application of the schemes to a Dirac Lindblad equation and real-space-time Green's function formulations are discussed.

研究动机与目标

• 解决在(2+1)D狄拉克系统中模拟混合态动力学的数值挑战，特别是针对拓扑绝缘体表面态。
• 克服标准有限差分格式在处理一阶时空导数时固有的费米子加倍问题。
• 开发适用于冯·诺依曼方程的高效、稳定且保持正定性的数值格式，且计算复杂度随网格尺寸线性增长。
• 通过Peierls替换在交错格点上实现电磁势的规范不变处理。
• 通过林德布拉德方程将框架扩展至开放量子系统，通过实空间时间格林函数形式实现非平衡输运建模。

提出的方法

• 采用单锥交错时空网格，确保能量-动量色散关系中仅存在一个狄拉克锥。
• 通过分别对密度算符的左（bra）和右（ket）指标应用哈密顿量，并利用对易子的乘积法则，实现bra-ket格式。
• 开发一种直接格式，通过单次扫描计算密度矩阵的时间导数，避免分别传播bra和ket分量。
• 在交错网格上使用单个晶格间距的中心差分格式，以保持精度并避免费米子加倍。
• 应用Peierls替换，以规范不变方式引入电磁矢势。
• 通过与拓扑绝缘体表面态相关的参数进行数值验证，并分析稳定性与收敛性。

实验结果

研究问题

• RQ1 针对(2+1)D狄拉克系统中的冯·诺依曼方程，有限差分格式能否在保持正定性的同时避免费米子加倍？
• RQ2 bra-ket格式的结构如何从纯态狄拉克方程中继承稳定性和正定性？
• RQ3 直接有限差分格式在密度矩阵演化中的稳定性和收敛性特性如何？
• RQ4 所提出的格式能否扩展至通过林德布拉德方程描述的开放量子系统动力学？
• RQ5 该框架如何适配实空间时间格林函数形式，以实现非平衡输运的建模？

主要发现

• bra-ket格式由于其乘积法则的实现方式，从纯态狄拉克方程中继承了正定性和稳定性。
• 直接格式在每时间步中矩阵运算更少，计算效率更高，但其收敛半径相比bra-ket格式有所减小。
• 两种格式均利用交错网格消除了费米子加倍，并确保能量色散中仅存在一个狄拉克锥。
• 在与拓扑绝缘体表面态相关的参数下，数值测试表明两种格式均稳定且具有收敛性。
• Peierls替换可规范不变地应用于交错网格，且在电磁势作用下仍保持格式的稳定性。
• 对林德布拉德方程和格林函数形式的扩展是可行的，从而能够模拟耗散效应和非平衡输运。