Power Laws Used to Extrapolate the Coupled Cluster Correlation Energy to the Thermodynamic Limit.

Abstract

Recent calculations using coupled cluster on solids have raised the discussion of using a N-1/3 power law to fit the correlation energy when extrapolating to the thermodynamic limit, an approach which differs from the more commonly used N-1 power law, which is, for example, often used by quantum Monte Carlo methods. In this paper, we present one way to reconcile these viewpoints. Coupled cluster doubles calculations were performed on uniform electron gases reaching system sizes of 922 electrons for an extremely wide range of densities (0.1 < rs < 100.0) to study how the correlation energy approaches the thermodynamic limit. The data were corrected for the basis set incompleteness error and use a selected twist angle approach to mitigate the finite size error from shell filling effects. Analyzing these data, we initially find that a power law of N-1/3 appears to fit the data better than a N-1 power law in the large system size limit. However, we provide an analysis of the transition structure factor showing that N-1 still applies to large system sizes and that the apparent N-1/3 power law occurs only at low N.