[论文解读] Pyqcm: An open-source Python library for quantum cluster methods

Théo N. Dionne, Alexandre Foley+2 more|arXiv (Cornell University)|May 29, 2023

Physics of Superconductivity and Magnetism参考文献 41被引用 7

一句话总结

tldr: Pyqcm 是一个实现 CPT、VCA 和 CDMFT 的 Python/C++ 库，使用精确对角化求解器处理哈伯德型模型，提供 Python 接口且采用 GPL 许可。

ABSTRACT

Pyqcm is a Python/C++ library that implements a few quantum cluster methods with an exact diagonalization impurity solver. Quantum cluster methods are used in the study of strongly correlated electrons to provide an approximate solution to Hubbard-like models. The methods covered by this library are Cluster Perturbation Theory (CPT), the Variational Cluster Approach (VCA) and Cellular (or Cluster) Dynamical Mean Field Theory (CDMFT). The impurity solver (the technique used to compute the cluster's interacting Green function) is exact diagonalization from sparse matrices, using the Lanczos algorithm and variants thereof. The core library is written in C++ for performance, but the interface is in Python, for ease of use and inter-operability with the numerical Python ecosystem. The library is distributed under the GPL license.

研究动机与目标

• Provide an open-source implementation of quantum cluster methods (CPT, VCA, CDMFT).
• Deliver an exact diagonalization impurity solver based on sparse matrices and Lanczos algorithms.
• Maintain a Python interface for ease of use and interoperability with numerical Python ecosystems.
• Offer a GPL-licensed framework suitable for studying Hubbard-like models.
• Describe architecture, features, and practical usage of PyQCM.

提出的方法

• Use exact diagonalization from sparse matrices with Lanczos variants to compute cluster Green functions.
• Represent the lattice problem via cluster tiling and inter-cluster coupling, formulating CPT as a matrix equation G^{-1}(ω)=G_c^{-1}(ω)-V.
• Implement partial Fourier transforms and multiple representations to handle cluster and super-lattice structures.
• Provide CPT, VCA, and CDMFT workflows with corresponding self-energy and Green function manipulations.
• Incorporate periodization schemes to obtain fully k-dependent Green functions from cluster quantities.
• Present a Python/C++ architecture that integrates with the numerical Python ecosystem.

实验结果

研究问题

• RQ1 How can CPT, VCA, and CDMFT be implemented efficiently with an exact diagonalization solver?
• RQ2 How do partial Fourier transforms and cluster-periodization affect the translation properties and computed Green functions?
• RQ3 What are the practical considerations and limitations of applying CPT, VCA, and CDMFT to Hubbard-like models using an ED solver?
• RQ4 How does the PyQCM library balance performance (C++) with usability (Python) for quantum cluster computations?

主要发现

• CPT provides a Green function rooted in cluster quantities via G^{-1}(ω)=G_c^{-1}(ω)-V(ω); it is exact in the non-interacting limit and in the strong-coupling limit.
• The library implements CPT, VCA, and CDMFT within a unified framework, enabling approximate lattice Green functions with momentum resolution suitable for ARPES comparisons.
• Periodization schemes (G-scheme, M-scheme) are discussed as methods to recover fully translation-invariant Green functions from cluster data, with different physical implications.
• The exact diagonalization impurity solver relies on Lanczos-type algorithms for ground-state and Green function computations, operating on sparse matrices for performance.
• PyQCM provides a Python interface atop a C++ core to combine computational efficiency with ease of use and interoperability with SciPy/NumPy.
• The architecture supports extended interactions and Hartree approximations within the CPT/CDMFT/VCA framework.