title: Solving the parquet equations for the Hubbard model beyond weak coupling.
abstract: We find that imposing crossing symmetry in the iteration process considerably extends the range of convergence for solutions of the parquet equations for the Hubbard model. When crossing symmetry is not imposed, the convergence of both simple iteration and more complicated continuous loading (homotopy) methods is limited to high temperatures and weak interactions. We modify the algorithm to impose the crossing symmetry without increasing the computational complexity. We also imposed time reversal and a subset of the point group symmetries, but they did not further improve the convergence. We elaborate the details of the latency hiding scheme which can significantly improve the performance in the computational implementation. With these modifications, stable solutions for the parquet equations can be obtained by iteration more quickly even for values of the interaction that are a significant fraction of the bandwidth and for temperatures that are much smaller than the bandwidth. This may represent a crucial step towards the solution of two-particle field theories for correlated electron models.
authors: Department of Physics and Astronomy, Louisiana State University, Baton Rouge, Louisiana 70803, USA.
authors: Ka-Ming
authors: KM
authors: Tam
authors: H
authors: H
authors: Fotso
authors: S-X
authors: SX
authors: Yang
authors: Tae-Woo
authors: TW
authors: Lee
authors: J
authors: J
authors: Moreno
authors: J
authors: J
authors: Ramanujam
authors: M
authors: M
authors: Jarrell
mesh_terms: D000465:Algorithms; D003198:Computer Simulation; D004553:Electric Conductivity; D008962:Models, Theoretical; D009716:Numerical Analysis, Computer-Assisted