The correlated random-phase-approximation (CRPA) method is applied to study the ground-state and excited-state properties of an anyon gas. We perform a systematic 1/N expansion to the calculation of elementary excitations, where N is the statistical parameter. We find that the leading-(1/N)-order contributions are the same as those obtained by calculations of the RPA response functions. We then consider corrections to the RPA by a perturbation of the exchange interactions, and show that these higher-order corrections are important at large wave vectors. Our calculation gives an excitation spectrum in good agreement with that obtained by the exact solution of a small system.
The correlated random-phase-approximation (CRPA) method is applied to study the ground-state and excited-state properties of an anyon gas. We perform a systematic 1/N expansion to the calculation of elementary excitations, where N is the statistical parameter. We find that the leading-(1/N)-order contributions are the same as those obtained by calculations of the RPA response functions. We then consider corrections to the RPA by a perturbation of the exchange interactions, and show that these higher-order corrections are important at large wave vectors. Our calculation gives an excitation spectrum in good agreement with that obtained by the exact solution of a small system.
Based on the random phase approximation (RPA), RPA renormalization [J. E. Bates and F. Furche, J. Chem. Phys. 139, 171103 (2013)] is a robust many-body perturbation theory that works for molecules and materials because it does not diverge as the Kohn-Sham gap approaches zero. Additionally, RPA renormalization enables the simultaneous calculation of RPA and beyond-RPA correlation energies since the total correlation energy is the sum of a series of independent contributions. The first-order approximation (RPAr1) yields the dominant beyond-RPA contribution to the correlation energy for a given exchange-correlation kernel, but systematically underestimates the total beyond-RPA correction. For both the homogeneous electron gas model and real systems, we demonstrate numerically that RPA renormalization beyond first order converges monotonically to the infinite-order beyond-RPA correlation energy for several model exchange-correlation kernels and that the rate of convergence is principally determined by the choice of the kernel and spin polarization of the ground state. The monotonic convergence is rationalized from an analysis of the RPA renormalized correlation energy corrections, assuming the exchange-correlation kernel and response functions satisfy some reasonable conditions. For spin-unpolarized atoms, molecules, and bulk solids, we find that RPA renormalization is typically converged to 1 meV error or less by fourth order regardless of the band gap or dimensionality. Most spin-polarized systems converge at a slightly slower rate, with errors on the order of 10 meV at fourth order and typically requiring up to sixth order to reach 1 meV error or less. Slowest to converge, however, open-shell atoms present the most challenging case and require many higher orders to converge.
The exactly solvable Luttinger model for a system of one-dimensional, massless, two-component fermions interacting via general two-body forces is used as a test case for a systematically generated sequence of decoupling approximations commonly referred to as random-phase approximations (RPA) of which the simple RPA is the lowest. It is found, irrespective of both the nature and strength of the interaction, that the exact (collective) spectrum, as well as the exact correlation energy, can already be obtained in the simple RPA, and furthermore, that the next higher RPA gives a vanishing correction to these quantities. It is also shown, quite in contradistinction to the exact result, that the free-gas momentum distribution nk persists in the second RPA as well as in the first, thereby eliminating the higher RPA as a possible approach to the calculation of nk. Thus, in the context of this particular model, although the first few orders of RPA are highly effective for calculating the excitation spectrum and correlation energy, they are sensibly powerless in their ability to predict any variation in structure beyond that of the free-gas distribution, a result which most probably holds as well in the realistic three-dimensional case.
The random phase approximation (RPA), the renormalized RPA (RRPA), and the self-consistent RPA (SCRPA) methods are applied to calculate the wave functions of the ground and excited states of an exactly solvable model. The approximated wave functions are expanded in the basis of the exact solutions. It is found that, when the RPA collapses, the RPA wave functions are orthogonal to the exact solutions while the RRPA and SCRPA ones have small but finite overlaps with the exact results. In spite of the apparently good agreement between the results of the RRPA, the SCRPA, and the exact solution, for the energy of the first excited state beyond the point of collapse, it is found that these approximations do not correctly describe the exact wave functions.
The second random-phase approximation (SRPA) is the simplest and most natural extension of the RPA. It enlarges the space of the elementary modes introduced to describe the collective states by adding 2 particle -2 hole excitations to the 1 particle -1 hole ones of the RPA. In deriving the SRPA equations, use is made, as in the RPA, of the so-called quasi-boson approximation (QBA) where expectation values in the ground state of the system are approximated by their values in the uncorrelated reference state. This, however, has been shown to imply a degree of approximation worse than that in the RPA. It is, therefore, necessary to improve the QBA by considering a reference state which contains some correlations. Having in mind to perform such calculations for realistic systems, we consider a simple extension of the SRPA in which the reference state contains 2 particle - 2 hole correlations. The quality of such an extension is tested by applying it to a solvable three-level model and found to be good.
It is well known that the second-order ground-state correlation energy calculated in the usual random-phase approximation (RPA) is off by a factor of 2. It has also been noted that there is a factor-of-2 discrepancy in the correlation corrections to the Hartree-Fock occupation probabilities. The reason for these phenomena is traced to the treatment of the exchange terms in the interaction matrix element. An illustration based on a simple single multipole interaction shows that the factor of 2 arises from a coherent sum of small higher-order terms. Results based on the consideration of one or a few eigenmodes at a time are seen to be free of this difficulty. The correlation energy, occupation-probability corrections, and equations of motion are reexamined in the Beliaev-Zelevinsky boson expansion, as an illustration of clearly defined higher-order treatments.
The correlation energy of an inhomogeneous electron gas is ruled by the density correlation function L. It obeys a Bethe-Salpeter equation (BSE), whose kernel is the functional derivative of the particle self-energy. It is solved on equal footing for two-particle and single-particle Green's functions within Hedin's GW approximation. According to the Hartree and exchange-correlation contributions to the kernel, the correlation expression decays into random-phase-approximation (RPA) -like and screened-exchange-like parts. The actual treatment of the screened potential W in the BSE for L leads to different second-order screened exchange corrections (SOSEX). Including correlation in the kernel even yields approximations beyond SOSEX. The progress beyond RPA is illustrated for spinless and spin-polarized homogeneous electron gases, indicating that the screening treatment is substantial.