Search

The exploitation of space group symmetries in numerical calculations of periodic crystalline solids accelerates calculations and provides physical insight. We present results for a space-group symmetry adaptation of electronic structure calculations within the finite-temperature self-consistent GW method along with an efficient parallelization scheme on accelerators. Our implementation employs the simultaneous diagonalization of the Dirac characters of the orbital representation. Results show that symmetry adaptation in self-consistent many-body codes results in substantial improvements of the runtime, and that block diagonalization on top of a restriction to the irreducible wedge results in additional speedup.

The many-body Green's function provides access to electronic properties beyond density functional theory level in ab inito calculations. In this manuscript, we propose a deep learning framework for predicting the finite-temperature Green's function in atomic orbital space, aiming to achieve a balance between accuracy and efficiency. By predicting the self-energy matrices in Lehmann representation using an equivariant message passing neural network, our method respects its analytical property and the $E(3)$ equivariance. The Green's function is obtained from the predicted self-energy through Dyson equation with target total number of electrons. We present proof-of-concept benchmark results for both molecules and simple periodic systems, showing that our method is able to provide accurate estimate of physical observables such as energy and density of states based on the predicted Green's function.

The real-time contour formalism for Green's functions provides time-dependent information of quantum many-body systems. In practice, the long-time simulation of systems with a wide range of energy scales is challenging due to both the storage requirements of the discretized Green's function and the computational cost of solving the Dyson equation. In this manuscript, we apply a real-time discretization based on a piece-wise high-order orthogonal-polynomial expansion to address these issues. We present a superconvergent algorithm for solving the real-time equilibrium Dyson equation using the Legendre spectral method and the recursive algorithm for Legendre convolution. We show that the compact high order discretization in combination with our Dyson solver enables long-time simulations using far fewer discretization points than needed in conventional multistep methods. As a proof of concept, we compute the molecular spectral functions of H$_2$, LiH, He$_2$ and C$_6$H$_4$O$_2$ using self-consistent second-order perturbation theory and compare the results with standard quantum chemistry methods as well as the auxiliary second-order Green's function perturbation theory method.

We study the fluctuations responsible for pairing in the $d$-wave superconducting state of the two-dimensional Hubbard model at intermediate coupling within a cluster dynamical mean-field theory with a numerically exact quantum impurity solver. By analyzing how momentum and frequency dependent fluctuations generate the $d-$wave superconducting state in different representations, we identify antiferromagnetic fluctuations as the pairing glue of superconductivity both in the underdoped and the overdoped regime. Nevertheless, in the intermediate coupling regime, the predominant magnetic fluctuations may differ significantly from those described by conventional spin-fluctuation theory.

We study the contribution of the electron-spin fluctuation coupling to the superconducting state of the two dimensional Hubbard model within dynamical cluster approximation (DCA) using a numerical exact continuous time Monte Carlo solver. By analyzing the frequency dependence of the self energy, we show that only about half of the superconductivity can be attributed to a "pairing glue" arising from treating spin fluctuations as a pairing boson in the standard one-loop theory.

We investigate charge fluctuations in the two-dimensional Hubbard model as a function of doping, interaction strength, next-nearest-neighbor hopping, and temperature within the eight-site dynamical cluster approximation. In the regime of intermediate interaction strengths, we find that $d$ density wave fluctuations, which had previously been postulated using analytical arguments, are present and strong but cannot be interpreted as the cause of the pseudogap in the model, due to their evolution with doping and interaction strength. For all parameters away from half filling, the charge fluctuations investigated, including $d$ density wave fluctuations, are weaker than $d$-wave superconducting fluctuations., Comment: 6 figures

Quantum many-body systems in thermal equilibrium can be described by the imaginary time Green's function formalism. However, the treatment of large molecular or solid ab inito problems with a fully realistic Hamiltonian in large basis sets is hampered by the storage of the Green's function and the precision of the solution of the Dyson equation. We present a Legendre-spectral algorithm for solving the Dyson equation that addresses both of these issues. By formulating the algorithm in Legendre coefficient space, our method inherits the known faster-than-exponential convergence of the Green's function's Legendre series expansion. In this basis, the fast recursive method for Legendre polynomial convolution, enables us to develop a Dyson equation solver with quadratic scaling. We present benchmarks of the algorithm by computing the dissociation energy of the helium dimer He$_2$ within dressed second-order perturbation theory. For this system, the application of the Legendre spectral algorithm allows us to achieve an energy accuracy of $10^{-9} E_h$ with only a few hundred expansion coefficients.

We compute the dynamical charge susceptibility in the two-dimensional Hubbard model within the dynamical cluster approximation. In order to understand the connection between charge susceptibility and pseudogap, we investigate the momentum, doping, and temperature dependence. We find that as a function of frequency, the dynamical charge susceptibility is well represented by a single peak at a characteristic frequency. It shows little momentum or temperature dependence, while the doping dependence is more evident, and no clear signature of the pseudogap is observed. Data for the doping evolution of the static susceptibility and for fluctuation diagnostics are presented. Our susceptibilities should be directly measurable in future Momentum-resolved electron energy-loss spectroscopy experiments., Comment: 7 pages, 6 figures