Correlation-driven multigap d-wave superconductivity in Anderson lattice model

We present the full Gutzwiller-wave-function solution of the Anderson lattice model in two dimensions that leads to the correlation-driven multigap superconducting (SC) ground state with the dominant $d_{x^2-y^2}$-wave symmetry. The results are consistent with the principal properties of the heavy-fermion superconductor CeCoIn$_5$. We regard the pairing mechanism as universal and thus applicable to other Ce-based heavy-fermion compounds. Additionally, a gain in kinetic energy in the SC state takes place, as is also the case for high-temperature superconductors.