We demonstrate in this work that the use of metasurfaces provides a viable strategy to largely tune and enhance near-field radiative heat transfer between extended structures. In particular, using a rigorous coupled wave analysis, we predict that Si-based metasurfaces featuring two-dimensional periodic arrays of holes can exhibit a room-temperature near-field radiative heat conductance much larger than any unstructured material to date. We show that this enhancement, which takes place in a broad range of separations, relies on the possibility to largely tune the properties of the surface plasmon polaritons that dominate the radiative heat transfer in the near-field regime. DOI: 10.1103/PhysRevLett.118.203901 Thermal radiation is one of the most ubiquitous physical phenomena. In recent years, there has been a renewed interest in this topic due to the confirmation of the prediction that radiative heat transfer can be drastically enhanced for bodies separated by small gaps [1,2]. This occurs when the gap is smaller than the thermal wavelength (9.6 μm at room temperature), and it is due to the contribution of evanescent waves that dominate the near-field regime. The fact that this near-field radiative heat transfer (NFRHT) can overcome the far-field limit set by the Stefan-Boltzmann law has now been verified in a variety of experiments exploring different materials, geometrical shapes, and gaps ranging from micrometers to a few nanometers [3][4][5][6][7][8][9][10][11][12][13][14][15][16][17]. These experiments have also triggered off the hope that NFRHT could have an impact in different thermal technologies [18] such as thermophotovoltaics [19], heat-assisted magnetic recording [20,21], scanning thermal microscopy [22][23][24], nanolithography [25], thermal management [26,27], or coherent thermal sources [28,29].In this context, the question on the fundamental limits of NFRHT is attracting a lot of attention [30]. So far, the largest NFRHT enhancements in extended structures have been reported for polar dielectrics (SiC, SiO 2 , SiN, etc.), in which the NFRHT is dominated by surface phonon polaritons (SPhPs) [31,32]. There has not been any report of an extended structure that has a heat transfer coefficient exceeding that between two planar polar dielectric surfaces, and that includes metamaterials like hyperbolic ones [33,34]. In an attempt to tune NFRHT, several calculations of NFRHT between periodic metallic nanostructures in both 1D [35][36][37][38] and 2D [39] have been reported. These calculations have shown some degree of tunability and a NFRHT enhancement over the corresponding material without nanostructuration. However, the reported NFRHT in these structures is still smaller than in the case of parallel plates made of polar dielectrics. There have also been theoretical studies of the NFRHT between photonic crystals and periodic metamaterials made of dielectrics [40][41][42] that show how the radiative properties can be enhanced with respect to the bulk counterpart. However, the resulting NFRHTs are agai...