Enhancing robustness and efficiency of density matrix embedding theory via semidefinite programming and local correlation potential fitting

Density matrix embedding theory (DMET) is a powerful quantum embedding method for solving strongly correlated quantum systems. Theoretically, the performance of a quantum embedding method should be limited by the computational cost of the impurity solver. However, the practical performance of DMET i...

Full description

Saved in:
Bibliographic Details
Published in:Physical review. B 2020-08, Vol.102 (8), p.1, Article 085123
Main Authors: Wu, Xiaojie, Lindsey, Michael, Zhou, Tiangang, Tong, Yu, Lin, Lin
Format: Article
Language:English
Subjects:
Computational efficiency
Computing costs
Computing time
Correlation
Density
Embedding
Iterative methods
Matrix methods
Numerical stability
Optimization
Robustness (mathematics)
Self consistent fields
Semidefinite programming
Two dimensional models
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Density matrix embedding theory (DMET) is a powerful quantum embedding method for solving strongly correlated quantum systems. Theoretically, the performance of a quantum embedding method should be limited by the computational cost of the impurity solver. However, the practical performance of DMET is often hindered by the numerical stability and the computational time of the correlation potential fitting procedure, which is defined on a single-particle level. Of particular difficulty are cases in which the effective single-particle system is gapless or nearly gapless. To alleviate these issues, we develop a semidefinite programming (SDP) based approach that can significantly enhance the robustness of the correlation potential fitting procedure compared to the traditional least squares fitting approach. We also develop a local correlation potential fitting approach, which allows one to identify the correlation potential from each fragment independently in each self-consistent field iteration, avoiding any optimization at the global level. We prove that the self-consistent solutions of DMET using this local correlation potential fitting procedure are equivalent to those of the original DMET with global fitting. We find that our combined approach, called L-DMET, in which we solve local fitting problems via semidefinite programming, can significantly improve both the robustness and the efficiency of DMET calculations. We demonstrate the performance of L-DMET on the 2D Hubbard model and the hydrogen chain. We also demonstrate with theoretical and numerical evidence that the use of a large fragment size can be a fundamental source of numerical instability in the DMET procedure.
ISSN:2469-9950
2469-9969
DOI:10.1103/PhysRevB.102.085123