CP2K implements the Gaussian and Augmented Planewaves (GAPW) approach for all-electron calculations. My understanding is that the GAPW method involves using atom-centered Gaussian type orbitals to represent rapidly varying electron density around nuclei, and planewaves for the electron density in the slowly varying interstitial regions.
My question is: does use of the GAPW method mostly eliminate basis set superposition error effects, for example bond length contraction during optimization? I would think so, since the core regions probably will not overlap in most systems, but I may be misunderstanding how GAPW works.