For systems with periodic boundary condition , plane wave basis is usually adopted. While for molecular systems, gaussian basis set is normally adopted. For periodic systems, the gaussian basis is transformed from real lattice space to momentum space (or k space) through Bloch function. My question is that what is the advantage of the plane wave basis? And a relevant question is that how to determine the number of k points required to converge the energy and density matrix for a periodic system?
In my opinion a solid, brief overview of DFT is given here.
The pro/cons of plane waves and other basis sets are discussed and I will list them here in case the link goes dead.
- Fourier coefficients stored in regular grid.
- Efficient FFT algorithms between r- and G-space representation.
- O(N^2) scaling on CPU
- Complete and orthonormal basis set.
- Not atom-centered -> unbiased.
- Systematically improvable by increasing the cut-off of the Fourier coefficients.
- Large set of basis coefficients. Hamiltonian cannot be stored.
- Sharp nodes of wave functions of core electrons are very expensive. Need pseudo-potential.
- Vacuum as expensive as atoms.
@Susi Lehtola Thanks for sharing the link. I agree that the linked question is more or less the same as the question here. However, I think most of the answers is limited to DFT. I know there are electronic packages that could do calculation for periodic systems with coupled cluster method.
For simplified model Hamiltonian, conventional solid state physicist also developed second quantization techniques in momentum space.
Like classical spin wave approach for ferromagnetism
Also more recently, quantum monte carlo approach is developed to study Hubbard model
All these method adopt a plane-wave basis. I think in the real lattice space, effect of delocalization might not be well captured. But these effect might only consequential for same special system (correspond to long wave length limit in k space). For normal systems, the interaction might be decay within first few nearest neighbours (especially for ground state).
My question is that how is the performance of the plane-wave basis for these post-Hartree-Fock approaches?