Efficient Quasiparticle Determination beyond the Diagonal Approximation via Random Compression

Annabelle Canestraight,Xiaohe Lei,Khaled Ibrahim,Vojtech Vlcek
2023-09-27
Abstract:Calculations of excited states in Green's function formalism often invoke the diagonal approximation, in which the quasiparticle states are taken from a mean-field calculation. Here, we extend the stochastic approaches applied in the many-body perturbation theory and overcome this limitation for large systems in which we are interested in a small subset of states. We separate the problem into a core subspace, whose coupling to the remainder of the system environment is stochastically sampled. This method is exemplified on computing hole injection energies into CO$_2$ on an extended gold surface with nearly 3000 electrons. We find that in the extended system, the size of the problem can be compressed up to $95\%$ using stochastic sampling. This result provides a way forward for self-consistent stochastic methods and determining Dyson orbitals in large systems.
Computational Physics,Materials Science
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to overcome the limitations of the diagonal approximation when dealing with the quasiparticle state calculation in large - scale systems, especially when only a small part of the states in the system are of interest. Specifically, the paper proposes a new stochastic compression technique. Through this method, the quasiparticle states in extended systems can be efficiently determined, especially for large - scale systems containing thousands of electrons. The key to this method lies in dividing the system into a core subspace and the remaining space. The states in the core subspace are processed by the standard stochastic many - body perturbation theory (MBPT), while the states in the remaining space are processed by introducing additional random sampling. This can not only significantly reduce the computational cost but also maintain the computational accuracy, especially when dealing with the interaction between molecules and surfaces, such as the charge injection energy calculation of CO₂ on the gold surface. Through this method, researchers can obtain more accurate quasiparticle energies and Dyson orbitals while maintaining computational efficiency.