Accurate prediction of the fundamental gap in insulators is a challenge for electronic structure theories. Electronic correlation and finite-size effects can significantly change the fundamental gap. Using quantum Monte Carlo, the authors show here that the fundamental gap can be calculated as the derivative discontinuity of the energy with respect to the electron density in the grand-canonical ensemble. Further, the finite-size error can be understood and corrected using information encoded in the electronic structure factor. The presented methodology applies for any explicitly correlated method.
Accurate prediction of the fundamental gap in insulators is a challenge for electronic structure theories. Electronic correlation and finite-size effects can significantly change the fundamental gap. Using quantum Monte Carlo, the authors show here that the fundamental gap can be calculated as the derivative discontinuity of the energy with respect to the electron density in the grand-canonical ensemble. Further, the finite-size error can be understood and corrected using information encoded in the electronic structure factor. The presented methodology applies for any explicitly correlated method.
A systematically converging series of ab initio, postdensity functional, size-consistent, electron-correlated approximations is desired for predictive computing of electronic band structures of insulating, semiconducting, and metallic solids. A series that meets all of these desiderata (except the applicability to metals) is ab initio many-body Green’s function theory based on Gaussian-type orbital (GTO) basis sets. Here, its leading-order approximation, the second-order Green’s function (GF2) method in the diagonal and frequency-independent approximations with the aug-cc-pVDZ basis set, is applied to the fundamental band gaps of three semiconductors (diamond, silicon, and silicon carbide in the zincblende structure) using cluster models. Corrections are made to the basis set incompleteness errors by the explicit correlation (F12) ansatz (GF2-F12) for the valence band edges. The crystals are modeled as surface-passivated clusters of increasing sizes, whose wave functions are expanded by up to 2709 GTO basis functions. Immense computational costs of these calculations are overcome by the highly scalable stochastic algorithm of the Monte Carlo GF2-F12 method, whose operation cost per state increases only as a cubic power of system size, which has a tiny memory footprint and easily achieves near-perfect parallel efficiency on thousands of CPUs or on hundreds of GPUs. The correlated, F12-corrected highest-occupied and lowest-unoccupied molecular orbital energy (HOMO–LUMO) gap is 5.78 ± 0.07 eV for C 87 H 76 as compared with the experimental value of the fundamental (indirect) band gap of bulk diamond at 5.48 eV. The correlated, F12-corrected HOMO–LUMO gaps for Si 75 H 76 and Si 32 C 43 H 76 are 2.56 ± 0.15 and 3.50 ± 0.12 eV, respectively, which are expected to decrease further with increasing cluster sizes. The experimental fundamental (indirect) band gaps of bulk silicon and silicon carbide are 1.17 and 2.42 eV, respectively.
We review the use of continuum quantum Monte Carlo (QMC) methods for the calculation of energy gaps from first principles, and present a broad set of excited-state calculations carried out with the variational and fixed-node diffusion QMC methods on atoms, molecules, and solids. We propose a finite-size-error correction scheme for bulk energy gaps calculated in finite cells subject to periodic boundary conditions. We show that finite-size effects are qualitatively different in two-dimensional materials, demonstrating the effect in a QMC calculation of the band gap and exciton binding energy of monolayer phosphorene. We investigate the fixed-node errors in diffusion Monte Carlo gaps evaluated with Slater-Jastrow trial wave functions by examining the effects of backflow transformations, and also by considering the formation of restricted multideterminant expansions for excited-state wave functions. For several molecules, we examine the importance of structural relaxation in the excited state in determining excited-state energies. We study the feasibility of using variational Monte Carlo with backflow correlations to obtain accurate excited-state energies at reduced computational cost, finding that this approach can be valid. We find that diffusion Monte Carlo gap calculations can be performed with much larger time steps than are typically required to converge the total energy, at significantly diminished computational expense, but that in order to alleviate fixed-node errors in calculations on solids the inclusion of backflow correlations is sometimes necessary.
Abstract Using the discontinuity of the chemical potential as a function of excess charge, the fundamental gaps for finite systems and the band gaps of extended solids are determined within reduced density matrix functional theory. We also present the necessary and sufficient conditions for the one-body reduced density matrix of a system with fractional charge to be ensemble N-representable. The performance of most modern day reduced density matrix functionals is assessed for the gaps and the correlation energy of finite systems. Our results show that for finite systems the PNOF, BBC3, and power functionals yield very accurate correlation energies while for a correct description of the fundamental gap the removal of self-interaction terms is essential. For extended solids we find that the power functional captures the correct band gap behavior for conventional semiconductors as well as strongly correlated Mott insulators, where a gap is obtained in absence of any magnetic ordering.
Discontinuities / Chemical Potential / Reduced Density Matrix Functional Theory Using the discontinuity of the chemical potential as a function of excess charge, the fundamental gaps for finite systems and the band gaps of extended solids are determined within reduced density matrix functional theory. We also present the necessary and sufficient conditions for the one-body reduced density matrix of a system with fractional charge to be ensemble N -representable. The performance of most modern day reduced density matrix functionals is assessed for the gaps and the correlation energy of finite systems. Our results show that for finite systems the PNOF, BBC3, and power functionals yield very accurate correlation energies while for a correct description of the fundamental gap the removal of self-interaction terms is essential. For extended solids we find that the power functional captures the correct band gap behavior for conventional semiconductors as well as strongly correlated Mott insulators, where a gap is obtained in absence of any magnetic ordering.
Ab initio quantum Monte Carlo methods, in principle, allow for the calculation of exact properties of correlated many-electron systems but are, in general, limited to the simulation of a finite number of electrons N under periodic boundary conditions. Therefore, an accurate theory of finite-size effects is indispensable to bridge the gap to realistic applications in the thermodynamic limit. In this work, we revisit the uniform electron gas at finite temperature, as it is relevant to contemporary research, e.g., in the field of warm dense matter. In particular, we present a new scheme to eliminate finite-size effects both in the static structure factor S(q) and in the interaction energy v, which is based on the density response formalism. We demonstrate that this method often allows us to obtain v in the thermodynamic limit within a relative accuracy of ∼0.2% from as few as N = 4 electrons without any empirical choices or knowledge of results for other values of N. Finally, we evaluate the applicability of our method upon increasing the density parameter rs and decreasing the temperature T.
Quantum Monte Carlo (QMC) is one of the most promising methods for solving quantum many-body problems. QMC combines known analytical properties of wave functions, results from other methods, and stochastic techniques, into a powerful tool for investigation of the electronic structure of real materials. Furthermore, the method has a favorable scaling ( O[N3] in the number of electrons), perfect scalability on parallel architectures, and can be applied to a wide range of systems. Our work has focused both on expanding the capabilities of QMC for calculation of quantities beyond energies, such as interatomic forces, and on application of the method to challenging problems in nanostructure materials research. For the calculation of interatomic forces we have implemented a finite difference correlated sampling method. The correlated sampling enables us to avoid the problem of statistical noise in evaluating the energy differences and allows us to obtain forces with an average 1% error. We have analyzed the effect of the method's main approximations on the estimate of the forces, both by theoretical arguments and by performing a variety of tests on small systems. Furthermore, we have used this newly developed method for estimating equilibrium geometries in a problem where other methods have failed. QMC enabled us to treat the electron correlation effects with high accuracy and to provide accurate predictions in a number of interesting applications. In particular, for carbon rings of C4N+2 stoichiometry we have resolved a long-standing problem regarding the competition between the aromaticity and second order Jahn-Teller effects. We have found that the dominant mechanism at small and intermediate sizes, which stabilizes the bond angle and bond length alternated geometries, is the secondorder Jahn-Teller effect while the aromatic isomer is always found to be a transition state. These high accuracy calculations which have involved multi-reference wave functions, extensive variational optimizations, and the use of massively parallel platforms demonstrate clearly the unique capabilities of the QMC approach.